- Research article
- Open Access
Does the core circadian clock in the moss Physcomitrella patens (Bryophyta) comprise a single loop?
BMC Plant Biology volume 10, Article number: 109 (2010)
The endogenous circadian clock allows the organism to synchronize processes both to daily and seasonal changes. In plants, many metabolic processes such as photosynthesis, as well as photoperiodic responses, are under the control of a circadian clock. Comparative studies with the moss Physcomitrella patens provide the opportunity to study many aspects of land plant evolution. Here we present a comparative overview of clock-associated components and the circadian network in the moss P. patens.
The moss P. patens has a set of conserved circadian core components that share genetic relationship and gene expression patterns with clock genes of vascular plants. These genes include Myb-like transcription factors PpCCA1a and PpCCA1b, pseudo-response regulators PpPRR1-4, and regulatory elements PpELF3, PpLUX and possibly PpELF4. However, the moss lacks homologs of AtTOC1, AtGI and the AtZTL-family of genes, which can be found in all vascular plants studied here. These three genes constitute essential components of two of the three integrated feed-back loops in the current model of the Arabidopsis circadian clock mechanism. Consequently, our results suggest instead a single loop circadian clock in the moss. Possibly as a result of this, temperature compensation of core clock gene expression appears to be decreased in P. patens.
This study is the first comparative overview of the circadian clock mechanism in a basal land plant, the moss P. patens. Our results indicate that the moss clock mechanism may represent an ancestral state in contrast to the more complex and partly duplicated structure of subsequent land plants. These findings may provide insights into the understanding of the evolution of circadian network topology.
The importance of timing of biological processes to the day-night cycle is reflected in the ubiquity and independent evolution of circadian clock mechanisms in different taxa . The endogenous circadian clock allows the organism to synchronize processes both to daily and seasonal changes. In plants, many processes such as photosynthesis that exhibit a daily rhythm are under the control of a circadian clock. The circadian system is also important to track seasonal changes, e.g. the anticipation of spring or autumn to induce reproduction . Even though circadian clocks in bacteria, fungi, plants and animals have evolved largely independently and thus are composed of different sets of genes, these clocks show a common structure or underlying principle. Recent experimental data and modeling suggests that circadian systems are complex networks of several interlocked feedback loops and interconnected input and output circuits [3, 4]. These loops contain both positive and negative elements on transcriptional and post transcriptional levels [5, 4].
Arabidopsis thaliana (Arabidopsis) has been the main plant model in studies of plant circadian clocks, and the first identified feedback loop contained two closely related and partly redundant Myb-like transcription factors CIRCADIAN CLOCK ASSOCIATED 1 (CCA1) and LATE ELONGATED HYPOCOTYL (LHY) forming a negative feedback loop with TIMING OF CAB EXPRESSION 1 (TOC1 also referred to as PRR1) . Experimental data have since shown that additional loops are needed to explain observed data. Mathematical modeling and empirical studies suggest that CCA1/LHY also form a second loop with PSEUDO RESPONSE REGULATOR 7 (PRR7) and PRR9 which both belong to the same family as TOC1. In addition, an evening phased loop is predicted including TOC1 and a factor that has been suggested to be GIGANTEA (GI) [7–9].
CCA1 and LHY are important regulators of circadian rhythm in Arabidopsis, not only involved in regulating the expression of TOC1, but also in mediating light input to the core oscillator [10–12]. CCA1/LHY regulate several output genes of the circadian clock by binding to the evening element of the promoter of these genes .
TOC1, or PRR1, share domains with the other PRR proteins (PRR9/7/5/3). The pseudo receiver domain is similar to the phospho-accepting domain found in authentic response regulators, and the plant specific CCT motif (CO, COL and TOC1) contains a putative nuclear localization signal and has also been shown to mediate protein-protein interactions [14, 15]. PRR9, PRR7 and PRR5 have, besides regulatory roles close to the Arabidopsis core oscillator, also functions in the response of the clock to both light and temperature [16–18].
Simulations and experimental data suggest that GI is a component of the circadian core oscillator with rhythm defect phenotypes of mutants in both continuous light (LL) and continuous dark (DD) [19, 20]. GI also has an important role in the photoperiodic control of flowering as it induces the circadian expression of CONSTANS (CO) in the late afternoon during long days in Arabidopsis . In addition, GI has a function in temperature compensation of the circadian clock as a functional GI has been shown to extend the range of temperatures at which robust rhythmicity can be maintained .
Further components associated with core clock function in Arabidopsis include LUX ARRYTHMO (LUX), EARLY FLOWERING 4 (ELF4), EARLY FLOWERING 3 (ELF3), and members of the ZTL gene family, ZEITLUPE (ZTL), LOV KELCH PROTEIN 2 (LKP2) and FLAVIN-BINDING, KELCH REPEAT, F-BOX 1 (FKF1). ZTL has been shown to have a functional role in core clock function by controlling the cyclical degradation of TOC1 [22, 23]. This is achieved by a blue light dependent protein interaction with GI, which stabilizes the ZTL protein . LKP2 may have a similar role, although its effect on circadian clock function appears to be smaller [24, 25]. FKF1 seems not to have a role in clock function itself but acts down stream of the core oscillator. The FKF1 protein is also stabilized by protein-protein interaction with GI, and regulates flowering time by targeted degradation of CDF2, a repressor of CO transcription [26–29]. Thus, all members of the ZTL family of proteins have functional roles closely associated to the evening phased loop comprising TOC1 and GI.
The Myb transcription factor LUX ARRHYTHMO, however, has been shown to interact with components of the morning-phased loop. While LUX is required for the expression of the core clock genes LHY and CCA1, the LUX promoter also contains an evening element, indicating that LUX itself is under negative control by LHY and CCA1 . In Arabidopsis, ELF4 has been proposed to have a function in the core clock and recent experimental data suggest that ELF4 likely has two inputs, at both PRR9/PRR7 and GI/LUX, and functions to repress the light-induced expression of these components of both the morning- and evening-phased loops [31, 32]. ELF3 acts in the gating of light input to the core oscillator and is necessary for light induced expression of CCA1/LHY; as such, ELF3 contributes to the resetting of the clock each day [33, 34, 32]. In addition, a recent study has demonstrated the dual role of ELF3 by showing that the protein regulates both circadian rhythm and flowering time in Arabidopsis by interacting with COP1 and thus mediating the cyclical degradation of GI . Although a specific placement of ELF3 in the plant circadian network has not yet been determined, the ELF3 protein seems to interact with components of both the morning- and evening phased circuits of the current three-loop model.
Comparative studies of clock associated components and the circadian systems in photosynthesizing organisms have been carried out in the green algae Chlamydomonas reinhardtii and Ostreococcus tauri [36–38], and in several seed plants, Mesembryanthemum crystallinum , Glycine max , Phaseolus vulgaris , Pisum sativum , Lemna gibba and Lemna paucicostata [43, 44], Oryza sativa [ and references therein], Castanea sativa [46, 47], Populus nigra  and Picea abies [ and unpublished data, our lab]. The consensus result of studies in seed plants is that homologs to core clock genes in Arabidopsis are present and furthermore appear to display a high degree of functional conservation. However, the situation is different in the algae C. reinhardtii and O. tauri. Beside implicated similarities in the phototransduction (CRYPTOCHROMES) and kinase pathways (CKI and CKIIs), that constitute important input and regulatory components of the circadian clock, homology to core clock genes in Arabidopsis is limited to Myb-related and PRR family-like proteins [36–38] (Additional file 1).
As a representative of Bryophytes that separated from the lineage leading to vascular plants more than 400 MYA , Physcomitrella patens provides the opportunity to study many aspects of land plant evolution. Circadian rhythms and photoperiodic responses are known from all categories of plants from algae and mosses to higher plants. In P. patens, as in higher plants, photoperiod affects the switch from vegetative growth to reproduction, as short days (SD) induce sporophyte development in Physcomitrella . Diurnal rhythms in expression have also been observed for some P. patens genes involved in photosynthesis [52, 53]. Furthermore, the presence of CONSTANS-like (COL) genes with a diurnal expression pattern has been reported [54, 4]. These genes are related to the central regulator of flowering time, CO in Arabidopsis, but their involvement in photoperiodic induction of reproduction in P. patens is not supported by available data [[54, 55]; our own unpublished data]. Finally, two homologs to known circadian clock genes in higher plants, PpCCA1a and PpCCA1b, have recently been characterized in P. patens .
To systematically identify putative photoperiod pathway and circadian clock components in P. patens and to study the evolution of the circadian clock, genes associated with the photoperiod pathway in Arabidopsis (Additional file 1) were used in BLAST searches against the genome sequence of selected plant species representing algae, non-vascular plants, non-seed vascular plants as well as seed plants. Identified putative P. patens clock genes were further characterized by phylogenetic analyses and by assessing their temporal expression profiles in photoperiod and free-running conditions.
Database searches in sequenced genomes representing algae, mosses, lycophytes, and angiosperms identified a limited set of homologs to genes implicated in circadian clock in non-land plants, while homologs to a majority of the included Arabidopsis clock associated genes were identified in all land plants (Additional file 1), with a few striking exceptions in moss (see below). These data suggest that a majority of the components in the circadian clock present in higher plants arose with the evolution of land plants. To study the evolution of the circadian clock of land plants in more detail, we focused on P. patens representing the earliest land plant with a sequenced genome, with the aim of identifying the putative core circadian clock genes present early in the evolution land plants. We focused on a set of genes for which there is strong evidence for a circadian clock function in higher plants.
All studies of diurnal and circadian expression of Physcomitrella genes have thus far observed clear cycling under light-dark cycles (LD) and/or in DD, but rapid dampening of the rhythm in LL [[52–56]; our own unpublished data]. This is in contrast to the situation in Arabidopsis, where endogenous expression rhythms can more readily be measured in LL whereas they dampen more quickly in DD [14, 40, 31, 34]. It is perhaps premature to state that the circadian clock in Physcomitrella is dysfunctional in LL, as has been suggested , and more studies of clock behavior in LL are needed. However, in the present study we chose to concentrate on rhythmic expression under LD and DD conditions.
Putative clock genes in P. patens
Single Myb domain genes related to CCA1 and LHY
The transcription factors CCA1 and LHY contain a single Myb-like domain of the well-conserved SHAQKYF motif at the N-terminal [Figure 1A, ]. One ortholog to CCA1/LHY has been identified and characterized in O. sativa  and we identified two predicted proteins in Selaginella moellendorffii that share significant sequence homology to CCA1/LHY (Figure 1B). Two predicted protein sequences containing related single Myb domains were also found in C. reinhardtii. The DNA-binding domain in Chlre4-402780 was identical to the Myb domain in the previously described ROC40 . A putative CCA1 homolog has also recently been identified and characterized in O. tauri . In P. patens, two orthologs have been identified and characterized, PpCCA1a and PpCCA1b [58, 56]. In the present study, PpCCA1a and PpCCA1b show stable circadian expression patterns in LD and for 48 hrs in DD without any apparent dampening (Figure 1C). The rhythmic pattern was statistically significant with COSOPT pMMC-β values at 0.026 for PpCCA1a and 0.030 for PpCCA1b (a listing of all statistical test results can be found in Additional file 2). In addition, phases are well superimposed onto those of AtCCA1 and AtLHY with peaks close to dawn.
Pseudo-response regulators (PRR)
The receiver domain at the N-terminal of the pseudo-response regulators, TOC1/PRR1 and PRR9/7/5/3, differ from the authentic response regulators, the ARR-family, by lacking an aspartate at the phospo-accepting site of the domain . In addition, the PRRs include a CCT motif at the C-terminal of the polypeptide (Figure 2A). Orthologs to TOC1 and the PRR family of proteins have been identified and extensively characterized in O. sativa [60, 45]. One TOC1-like sequence, Selmo1_438647, and its allelic variant, Selmo1_447266 (not shown in tree 2B), are annotated as predicted proteins in the Selaginella genome database, although both gene annotations contain a stop codon before the C-terminal CCT domain. Four additional S. moellendorffii sequences positioned on scaffold_88:201505-204488 and scaffold_62:750390-752881, including their probable allelic variants located on scaffold_42:1313168-13161115 and scaffold_20:1992667-1995113, cluster within the remaining PRR sequences (Figure 2B). These are not annotated as predicted proteins and lack protein ID numbers, because they appear as fused together with other genes. Even if it cannot be determined at present whether these putative proteins are expressed and translated, it appears obvious that the corresponding nucleotide sequences are present in the S. moellendorffii genome. In version 4.0 of the Chlamydomonas genome database, only one predicted amino acid sequence displays resemblance to the domain structure of the TOC1/PRR family of proteins, Chlre4-166515 [Additional file 1]. This protein clusters with a recently identified and described PRR-like protein in another unicellular green alga, Ostreococcus tauri . This protein, Ot-AY740079, also shares domain structure with the TOC1/PRR family of proteins and has a clearly shown function within the O. tauri circadian clock . However, although the O. tauri protein is referred to as a TOC1-homologue, the phylogenetic analysis reveals a distant relationship to the remaining TOC1/PRR family of proteins (Figure 2B). This suggests that the algal proteins constitute sister lineages to the whole TOC1/PRR family, and that it is not possible from presently available data to determine whether these proteins are more closely related to the TOC1-group, or the group consisting of the remaining PRR proteins. We therefore choose to refer to them as PRR family-like (Additional file 1). Four pseudo-response regulator homologs have been identified in P. patens, PpPRR1-4 [Figure 2B, ]. PpPRR1 to PpPRR4 are closely related and form a cluster next to the S. moellendorffii PRRs in the phylogenetic tree (Figure 2B). Thus, no ortholog to TOC1 appears to exist in P. patens. All putative response-regulators in P. patens show stable rhythm under LD conditions, with very slight dampening of amplitude in DD (COSOPT pMMC-β values ranging from 0.019 to 0.049) (Figure 2C, right panel). However, under our sampling scheme, we could not detect any differential expression of phase peaks among PpPRR1 to PpPRR4, similar to the sequential expression of phase displayed by the Arabidopsis quintet of PRR genes [Figure 2C, left panel; ]. Rather, our data suggest that the closely related, and probably recently duplicated PpPRR1 to PpPRR4, have not diverged in terms of their temporal gene expression pattern.
The emergence in vascular plants of a markedly diverged TOC1 sub-family distinct from the other PRR genes, suggests that positive Darwinian selection might have been important in the evolution of TOC1 genes. To test this hypothesis, we estimated dN/dS ratios in a phylogenetic tree including identified PRR and TOC1-like genes from a diverse set of species (Additional file 3). The results suggest that the branch leading to the TOC1 clade has diverged due to adaptive evolution (2Δl = 6.75, p < 0.01) and 21.8% of the sites in the two conserved domains show a significant signal of positive selection (Additional file 3).
In Arabidopsis, GI is a 127-kD nuclear protein consisting of 1173 amino acids. It does not contain any previously described protein domains or motifs. Orthologs to GI have been identified in O. sativa , and in addition, we found two predicted protein sequences in S. moellendorffii, Selmo1-140066 and Selmo1-170553, that displayed significant sequence homology to AtGI (Additional file 1). However, GI appear to be absent in Chlamydomonas reinhardtii and in O. tauri as previously reported [Additional file 1, [36–38]].
Similarly, no searches for GI, on either nucleotide or protein levels, yielded any hits against version 1.1 of the P. patens genome database. We conclude that GI, or any GI-like sequences, are absent in the Physcomitrella genome.
LUX ARRYTHMO (LUX) or PHYTOCLOCK 1 (PCL), contains a Myb-related DNA-binding domain similar to, but distinct from those found in LHY and CCA1 [Figure 3A; [30, 62]]. A LUX ortholog has been identified in O. sativa  and one predicted protein in S. moellendorffii displays sequence homology to AtLUX (Figure 3B). Two proteins in C. reinhardtii, CrROC15 and CrROC75, also include Myb DNA binding motifs that share homology with AtLUX . A search for putative P. patens orthologs revealed four predicted protein sequences that clustered together with OsPCL1 and AtLUX (Figure 3B). Of these putative LUX orthologs, only Phypa-47310 showed a tendency to a circadian expression pattern (Figure 3C). Although the rhythmic expression pattern of Phypa-47310 was not significant (COSOPT pMMC-β = 0.23 and Fischer's exact g p-value = 0.07), a visual inspection of the expression pattern in DD suggests a periodic cycling with a phase peak very similar to that of AtLUX (Figure 3C).
EARLY FLOWERING 4
ELF4 belongs to the DUF1313 family of plant proteins; a group of relatively short proteins with unknown function (Figure 4A). Several gene products in both O. sativa and S. moellendorffii show significant sequence homology to AtELF4 (Figure 4B); however, the expression patterns of the rice sequences were reported not to be under circadian control and thus not considered to be ELF4 orthologs . A predicted protein incorporating the DUF1313 domain is also present in C. reinhardtii [Additional file 1]. In our phylogenetic analysis of the DUF1313 protein sequences, AtELF4 clusters together with AtEFL1, while AtEFL2-4 form a separate and more distantly related lineage (Figure 4B). This grouping is supported by recent detailed analyses of the Arabidopsis DUF1313 proteins, where in addition, it was shown that only AtEFL1 were able to fully complement the elf4 loss-of-function phenotype . One predicted protein in P. patens, Phypa-49622, displays significant sequence homology to the DUF1313 domain and clusters together with the Selaginella sequences (Figure 4B). However, the expression profile for Phypa-49622 displayed very weak amplitude in DD and was out of phase with ELF4 in Arabidopsis (Figure 4C). Present data suggest that Phypa-49622 is not likely to be a functional homolog of AtELF4.
EARLY FLOWERING 3
The ELF3 protein does not contain any known functional domains or motifs. Two gene products in O. sativa, Os01g0566100 and Os06g01426000, show significant sequence homology to AtELF3, however, their genes do not display any circadian expression pattern in rice . We found two sequences in S. moellendorffii with significant homology to ELF3, Selmo1-415241 and Selmo1-411196 including their allelic variants Selmo1-443557 and Selmo1-443909, respectively (not shown in figure 5). No ELF3 like proteins were found in the latest version of the C. reinhardtii database . A search of the P. patens database produced three predicted proteins with significant sequence homology to AtELF3, that all cluster together in the phylogenetic tree (Figure 5A). Of the three putative ELF3-like genes, Phypa-66647 and Phypa-165364 displayed an oscillating pattern in both LD and DD (COSOPT pMMC-β values = 0.025 and 0.16, respectively), with a peak of phase that agreed well with expression data of ELF3 in Arabidopsis (Figure 5B).
The ZEITLUPE family
ZTL, and its homologs LKP2 and FKF1, are unique F-box proteins with a PAS/PAC domain at the N-terminal and a series of Kelch repeats at the C-terminal [26, 25, 63]. The PAS/PAC domain is a blue light receptor also found in PHOTOTROPINS, whereas the F-box mediates ubiquitination of specific target proteins . The Kelch repeats form groups of beta sheets that interact with other proteins . Six separate Kelch-repeats have been identified in AtZTL and AtFKF1 [26, 25, 63]; however, for comparison we display only five repeats in figure 6, which is the result obtained from the SMART database [95, 96]. Three proteins with domain structures similar to the ZTL family in Arabidopsis have been identified in O. sativa, OsZTL1, OsZTL2 and OsFKF1 . Two predicted proteins in S. moellendorffii with high sequence similarity, Selmo1-174189 and Selmo1-185595, also share this typical domain structure. In accordance with previous studies, we could not detect any ZTL family-like proteins or predicted protein sequences in Chlamydomonas reinhardtii or O. tauri [Additional file 1, [36–38]]. Our scan of the Physcomitrella genome revealed several proteins containing the above described functional domains, however, the three domains were never present in the same predicted protein sequence. Instead, the PAS/PAC domain was found in duplicates in several predicted phototropins, e.g. PpPHOTA1, PpPHOTB1 and PpPHOTB2. The duplicate PAS/PAC domains also occurred in isolation, as in Phypa-143199 and Phypa-4514 (Figure 6). F-boxes followed by Kelch-like repeats could be found in two predicted proteins in P. patens, Phypa-131411 and Phypa-14114 but never with the PAS/PAC domain at the N-terminal (Figure 6). Thus, proteins with a domain architecture analogous to the one found in the Arabidopsis ZTL family, in three annotated O. sativa proteins and in two predicted S. moellendorffii protein sequences, appear to be absent in P. patens.
Temperature compensation in P. patens
GI has been implicated in QTL mapping experiments and mutant studies as a strong candidate for mediating temperature compensation in the Arabidopsis circadian clock [66, 21]. It has been shown that GI critically affects clock function under natural conditions by extending the range of temperatures at which circadian rhythmicity can be maintained and period length remains stable. In particular, it would appear that the GI/TOC1 loop participates by regulating the expression levels of CCA1 and LHY at both low (12°C) and high (27°C) temperatures . Since neither GI, nor a GI/TOC1 loop appear to be present in P. patens, it could be hypothesized that the ability to compensate for changes in ambient temperature is either maintained by other components related to the clock mechanism, or is reduced in the moss. To test the temperature response of putative clock genes in P. patens, cultures were sampled in parallel in constant dark (DD) in three different temperatures, 12°C, 17°C and 25°C. Temperature was changed one day before the onset of DD at ZT -24. Gene expression levels were measured for all putative clock genes in P. patens (Figure 7A-J).
A decrease in amplitude and less well defined phase peaks are evident in the expression pattern of all genes at lower and higher temperatures. The broader phase peaks are noticeable as early as the first day in DD, whereas the loss of amplitude is more evident during the second day. These indications of a less robust rhythmic pattern are more pronounced at 12°C where several genes lose any trace of rhythmic expression pattern (Figure 7C-E, H). Although a slight decrease in amplitude can also be seen during the second day at 17°C, most likely due to a lack of photoperiodic input, most genes maintain a well-defined phase peak at this temperature. The observed decay of robust rhythm is reflected in the statistical analysis of the time series data. At 17°C, expression patterns of PpCCA1a-b, PpPRR1-4 and PpELF3L1 all display a significant circadian rhythm as reported by COSOPT (pMMC-β values ≤ 0.05) (Additional file 2). The group of significantly oscillating genes are reduced to PpCCA1a-b and PpPRR1-2 at 25°C, whereas only PpCCA1a-b maintain significant rhythms at 12°C, as reported by either COSOPT or Fisher's exact G test (Additional file 2). Although estimates of period length reported by COSOPT are not to be taken at face value under the present sampling scheme, it is notable that all significant expression patterns cluster together at 17°C with an estimated period length of approximately 22 h (Additional file 4). In summary, time series expression data suggest that the temperature interval of robust circadian rhythm is more narrow in P. patens compared to wild type Arabidopsis. However, the comparatively robust rhythm and stable period lengths of PpCCA1a and PpCCA1b at all temperatures may indicate a degree of temperature compensation of these core components independent of GI or a GI/TOC1 loop.
The moss, P. patens, has a set of conserved clock-associated components, which share genetic relationship and gene expression profiles with clock genes of higher plants. These components include members of the Myb family of transcription factors, PpCCA1a and PpCCA1b, and four pseudo-response regulators, PpPRR1-4 [58, 56]. Further putative clock genes present in the moss include one LUX-like, as well as two ELF3 and possibly one ELF4-like component of the Arabidopsis circadian mechanism (see Additional file 1 for a comparative overview of clock associated genes in Arabidopsis and P. patens). However, we note that P. patens appears to lack any genes coding for proteins orthologous to TOC1, GI or any of the ZTL family of proteins, all of which are present in S. moellendoerffii. Genes orthologous to some or all of these genes have also been reported in a number of seed plants, e.g. the monocotyledons rice [61, 45], two species of Lemna [43, 44], the angiosperms soybean , pea , chestnut , and the gymnosperm Norway spruce .
None of the three core clock genes absent in the moss are present in the green algae C. reinhardtii, as has been previously reported [36, 37] or O. tauri , if as we propose, OtTOC1 constitute a sister lineage to the whole PRR gene-family rather than a TOC1 ortholog. The result of the phylogenetic analysis of the pseudo-response regulator family of proteins implies that a common ancestor of all PRR genes was likely present already in the green algae C. reinhardtii and O. tauri [Figure 2B and Additional file 1, [37, 38]]. This notion seems further supported by the fact that the OtPRR-like sequence is the only pseudo-response regulator to date that display a functional phosphate acceptor aspartyl residue, implying recent divergence from the ARR family of proteins . These data collectively suggests that the three core clock genes have not been lost in the moss, but rather appeared before, or simultaneously, with the emergence of the vascular plants.
The four PRR proteins in P. patens display a high degree of sequence similarity and cluster tightly together among the PRR proteins of higher plants, but are as a group well separated from the branch including the TOC1/PRR1s of these plants (Figure 2B). The relatively recent expansion of the P. patens PRR gene family is reflected in the expression profiles of PpPRR1-4; while the Arabidopsis PRR9/7/3/5 and AtTOC1 display a serial distribution of phase peaks suggesting functional divergence , the PRR orthologs in P. patens have a very similar expression pattern with near identical phase peaks (Figure 2C).
Our data suggest that the TOC1/PRR1 orthologs in S. moellendorffii, O. sativa and Arabidopsis represent a more recent and functionally divergent variant among the two component pseudo-response regulators. A functional divergence of the TOC1-like genes compared the other PRR genes is supported by the elevated dN/dS ratio inferred on the branch separating the PRR and TOC1 genes. It suggests that a PRR gene duplication in the lineage leading to present day vascular plants was followed by positive selection on one of the copies possibly due the recruitment of new genes and features to the circadian clock. It is tempting to speculate that the adaptive evolution of TOC1 is coupled with the occurrence of GI and ZTL with which TOC1 closely interact.
While proteins with a domain architecture analogous to the ZTL family of proteins are also absent in the moss, the functional domains of these proteins do appear in other combinations; these sequences include predicted phototropins and also F-box proteins followed by Kelch repeats, but without the LOV domain at the N-terminal (Figure 6). From sequence homology alone, we have not been able to determine if any of these domains are more closely related to corresponding domains of ZTL like proteins in higher plants. In any case, the emergence of ZTL family-like proteins in S. moellendorffii suggest the formation of this novel protein architecture from the modular assembly of already existing functional domains, possibly through mechanisms such as duplication followed by domain shuffling [67, 68].
No trace of GI or GI-like sequences can be found at the nucleotide or protein levels in available complete genomes other than in S. moellendorffii and in representatives of higher plants (see references above). Although gene duplication is considered one of the most common and important mechanisms for the emergence of novel genes, it is also known that each genome harbors a certain set of genes that cannot be associated with other known genes and that often are of unclear origin . Occasionally, these orphan, or de novo, genes can be traced to ancestral noncoding DNA and intergenic regions [70–72], but there are also instances, e.g. in Drosophila and rice, where the exact origin of functional genes remains obscure [73, 74]. It may be that the additional sequencing of full genomes representing different phylogenetic lineages may help to reveal an origin of GI in intergenic regions of early land plants or algae. However, there is also the rare possibility of horizontal gene transfer, which has been exemplified as occurring between fungi and early plant lineages such as bryophytes and lycophytes . Presently, it would appear that GI is a de novo gene with unclear origin in the land plant lineage.
In view of the current three-loop model of the Arabidopsis circadian clock, the components absent in P. patens, TOC1 and GI, comprise an entire feed-back loop, also referred to as the evening-phased loop, or the evening oscillator [Figure 8A, 7]. In addition, TOC1 acts as the second feed-back component together with CCA1/LHY in the originally described central loop model of the Arabidopsis clock. Thus, while the clock mechanism of higher plants may at least comprise a three-loop design, our findings imply a single-loop for the moss clock (Figure 8B). It is notable that ZTL, LKP2 and FKF1, with known functions closely related to TOC1 and GI in Arabidopsis, also appear not to be present in P. patens. The singular morning-phased loop in the P. patens circadian mechanism would then consist of PpCCA1a/PpCCA1b joined together with one or some of the PpPRRs present in the moss (Figure 8B). Interestingly, such a hypothesized single-loop plant clock is entirely analogous to the "simplified slave version" of the circadian clock described in Arabidopsis roots. In roots, the morning-phased loop is decoupled from the central- and evening-phased loops and the clock appears to run only on the feed-back interaction between CCA1/LHY and PRR9/7 . However, in a recent study it was shown that double disruptants of PpCCA1a and PpCCA1b in the moss still display rhythmic output of the plastid sigma factor PpSIG5 and a member of the PpPRR-quartet . This suggests the presence of further components, and/or possibly additional unknown loops, other than a TOC1/GI loop, that contribute to the maintenance of endogenous circadian rhythm in P. patens. Recent experimental data illustrate a similar situation in the Arabidopsis three-loop model, where triple mutants prr9-10, prr7-11 and toc1-2, still display detectable endogenous rhythm in LHY/CCA1 expression, although no feed-back loop should remain intact . To what extent known clock components ELF3, ELF4 and LUX may contribute to this rhythm in Arabidopsis is still unknown.
In Arabidopsis, the protein GI has a critical effect on robustness of circadian rhythmicity and stability of period length of clock output during prolonged fluctuations in temperature . Our assay of temperature response of putative clock genes in P. patens shows that most genes lose their circadian expression at temperatures 12°C and 25°C in constant dark (Figure 7). The decay of rhythmic expression pattern is most pronounced at 12°C where only PpCCA1a and PpCCA1b maintain significant circadian oscillation (Additional file 2 and 4). One reason for this reduced level of temperature compensation may be the absent GI/TOC1 loop in the moss.
The proposed simpler loop structure of the circadian clock in P. patens may have implications for the evolutionary origins and development of clock network topology. Hallmark features of endogenous circadian oscillators include the ability to maintain period length under constant conditions, to compensate for fluctuations in ambient variables, and yet be able to entrain, or reset, the clock daily according to environmental cues of e.g. light and temperature . It has been pointed out that the complex loop structure of different circadian clock mechanisms may be a result of these inherently conflicting evolutionary aims of simultaneous robustness and sensitivity . Mathematical modeling and simulation studies have shown that the degree of flexibility in an interconnected metabolic system, for example its ability to attain several aims simultaneously, depends less on the number of components in such a system, than on the structure of the network itself. Accordingly, there would be a selective advantage in increasing loop number and complexity, as well as on mechanisms that would enable this, e.g. gene duplication [79, 80]. In this perspective, the circadian clock in P. patens may well represent an antecedent state comprising a single loop structure in contrast to the, in part, duplicated three-loop clock design of subsequent land plants.
Although the Arabidopsis three-loop model has proven to have greater realism than previous models, and provides a valuable framework for comparative studies, it is not yet complete. It was initially pointed out that the model in its proposed version did not account for all known clock associated genes or interactions, and experimental data continues to imply the presence of unknown components and possibly additional loops [7, 81, 77]. We believe that the further study of a seemingly less complex network system not only can provide important insights into the evolution of circadian network topology, but in addition, may help clarifying some of the remaining issues of the circadian clock mechanism of higher model plants.
Plant materials and growth conditions
Physcomitrella patens of the Gransden Wood wild type strain was cultured as protonemal tissue on cellophane covered BCD medium (0.8% agar, 1 mM MgSO4, 1.85 mM KH2PO4, 10 mM KNO3, 45 μM FeSO4, 1 mM CaCl2, 1× Hoagland's Number 2 solution), supplemented with 5 mM ammonium tartrate. All cultures were grown in Sanyo MLR-350 growth chambers supplied with white fluorescent tubes, Toshiba FL40SS W/37 (430-650 nm), at 65-75 μmol m-2 s-1. Light/dark (LD) conditions used were 16 h light/8 h dark, followed by constant dark (DD). Cultures were grown at 20°C, and transferred to temperatures 12°C, 17°C and 25°C one day before the onset of DD as noted in the text for each experiment. Cultures transferred to 17°C were sampled every 4th hour for 72 hours, the first day under LD conditions, followed by two days in DD. Cultures transferred to 12°C and 25°C were sampled every 4th hour for 52 hours, two time points in light followed by 12 time points in DD, as indicated in each figure.
RNA sample preparation and cDNA synthesis
Total RNA was extracted from each time point sample using the RNeasy plant mini kit (Qiagen). For each sample, 0.5 μg of total RNA was reverse transcribed to cDNA using random hexamer primers (Invitrogen) and Superscript III reverse transcriptase (Invitrogen), following the manufacturer's instructions.
Quantitative RT-PCR assay for time-series analyses of gene expression
cDNA samples were diluted 1:50 and amplified with the DyNAmo Flash SYBR Green qPCR Kit (Finnzymes) on a MyiQ Real-Time PCR Detection System (Bio-Rad). The two-step cycling program was as follows: 95°C for 7 min, followed by 40 cycles at 95°C for 15 s and 60°C for 1 min. Melt curve analyses were performed after each amplification to ensure specificity of products. Each cDNA sample was run in duplicates on separate plates. PCR efficiencies were calculated for each amplification with the software LinRegPCR . Any wells showing strongly deviating PCR efficiencies of either target or reference genes were excluded from further analysis. We used transcription level measurements for the P. patens beta-tubulin 1 gene (Phypa_186458), as a reference to normalize target gene transcription levels. The PpTUB1 reference gene displayed consistent amplification over all sample time points and treatments without evidence of circadian rhythmicity (Additional file 5). The CT values of duplicates were averaged, and the difference of the mean CT values for reference and target genes (ΔCT) was calculated for each sample time point and treatment combination. Target and reference primer sequences are listed in Additional file 6. The software COSOPT was used to detect significant circadian expression patterns under constant conditions (DD), under all temperatures, among the genes in the P. patens time series data [83–85]. In addition, a Fisher's exact g test, as implemented in the statistical software R package GeneCycle http://cran.r-project.org/web/packages/GeneCycle/index.html, was used to search for any harmonic pattern with unknown frequency in the time series expression data measured under constant conditions [86, 87]. Comparative expression data for Arabidopsis measured under diurnal and circadian conditions were downloaded from the DIURNAL database [, http://diurnal.cgrb.oregonstate.edu/].
Sequence analyses and reconstruction of phylogenetic relationships
Searches for homologs of Arabidopsis circadian clock components in other plant genomes used the following releases: Physcomitrella patens v1.1 [http://www.cosmoss.org/; ], Oryza sativa http://riceblast.dna.affrc.go.jp, Selaginella moellendorffii v1.0 http://genome.jgi-psf.org/Selmo1/Selmo1.home.html, and Chlamydomonas reinhardtii v4.0 [http://genome.jgi-psf.org/Chlre4/Chlre4.home.html; ]. Search criteria included sequences from the above databases with e-values ≤ 1e-10 and a domain structure similar to the one found in Arabidopsis proteins. Conserved regions of predicted protein sequences were aligned with the software MUSCLE . Protein sequence phylogenies were analyzed with a maximum likelihood model as implemented in PHYML 3.0  and with a Bayesian inference method as implemented in MrBayes 3.1 [93, 94]. The SMART database http://smart.embl-heidelberg.de/ was used for analyses of annotated protein domains and motifs [95, 96]. Test for positive selection using dN/dS (ω) ratios were performed on an alignment of selected PRR genes from angiosperms (Additional file 3). We used the branch-site test of positive selection as implemented in PAML 4.3 . Briefly, the branches of the tree are divided into foreground (the one leading to the TOC1 clade) and background branches (all other branches). Codon sites are separated into four classes with proportions p 0 , p 1 , p 2a and p 2b . Class 0 have ω < 1, and class 1 ω = 1 on both types of branches. Class 2a has ω < 1 on the background branch and ω ≥ 1 on the foreground branch. Finally, class 2b has ω = 1 on the background branch and ω ≥ 1 on the foreground branch. The test for positive selection is then obtained by comparing the likelihood of this model to the same model except that ω is fixed to 1 in site classes 2a and 2b on the foreground branch . Twice the log likelihood difference between the two compared models (2Δl), was compared against χ2with 1 df.
Young MW, Kay SA: Time zones: a comparative genetics of circadian clocks. Nat Rev Genet. 2001, 2: 702-15. 10.1038/35088576.
Yakir E, Hilman D, Harir Y, Green RM: Regulation of output from the plant circadian clock. FEBS J. 2007, 274: 335-45. 10.1111/j.1742-4658.2006.05616.x.
Roenneberg T, Merrow M: The network of time: understanding the molecular circadian system. Curr Biol. 2003, 13: R198-207. 10.1016/S0960-9822(03)00124-6.
Harmer SL: The circadian system in higher plants. Annual review of plant biology. 2009, 60: 357-77. 10.1146/annurev.arplant.043008.092054.
McClung CR: Comes a time. Curr Opin Plant Biol. 2008, 11: 514-20. 10.1016/j.pbi.2008.06.010.
Alabadí D, Oyama T, Yanovsky MJ, Harmon FG, Más P, Kay SA: Reciprocal regulation between TOC1 and LHY/CCA1 within the Arabidopsis circadian clock. Science. 2001, 293: 880-3. 10.1126/science.1061320.
Locke JCW, Kozma-Bognár L, Gould PD, Fehér B, Kevei E, Nagy F, Turner MS, Hall A, Millar AJ: Experimental validation of a predicted feedback loop in the multi-oscillator clock of Arabidopsis thaliana. Mol Syst Biol. 2006, 2: 59-10.1038/msb4100102.
Zeilinger MN, Farr'e EM, Taylor SR, Kay SA, Doyle FJ: A novel computational model of the circadian clock in Arabidopsis that incorporates PRR7 and PRR9. Mol Syst Biol. 2006, 2: 58-10.1038/msb4100101.
Kim WY, Fujiwara S, Suh SS, Kim J, Kim Y, Han L, David K, Putterill J, Nam HG, Somers DE: ZEITLUPE is a circadian photoreceptor stabilized by GIGANTEA in blue light. Nature. 2007, 449: 356-60. 10.1038/nature06132.
Wang ZY, Tobin EM: Constitutive expression of the CIRCADIAN CLOCK ASSOCIATED 1 (CCA1) gene disrupts circadian rhythms and suppresses its own expression. Cell. 1998, 93: 1207-17. 10.1016/S0092-8674(00)81464-6.
Schaffer R, Ramsay N, Samach A, Corden S, Putterill J, Carré IA, Coupland G: The late elongated hypocotyl mutation of Arabidopsis disrupts circadian rhythms and the photoperiodic control of flowering. Cell. 1998, 93: 1219-29. 10.1016/S0092-8674(00)81465-8.
Mizoguchi T, Wheatley K, Hanzawa Y, Wright L, Mizoguchi M, Song HR, Carré IA, Coupland G: LHY and CCA1 are partially redundant genes required to maintain circadian rhythms in Arabidopsis. Dev Cell. 2002, 2: 629-41. 10.1016/S1534-5807(02)00170-3.
Harmer SL, Hogenesch JB, Straume M, Chang HS, Han B, Zhu T, Wang X, Kreps JA, Kay SA: Orchestrated transcription of key pathways in Arabidopsis by the circadian clock. Science. 2000, 290: 2110-3. 10.1126/science.290.5499.2110.
Strayer C, Oyama T, Schultz TF, Raman R, Somers DE, Más P, Panda S, Kreps JA, Kay SA: Cloning of the Arabidopsis clock gene TOC1, an autoregulatory response regulator homolog. Science. 2000, 289: 768-71. 10.1126/science.289.5480.768.
Ben-Naim O, Eshed R, Parnis A, Teper-Bamnolker P, Shalit A, Coupland G, Samach A, Lifschitz E: The CCAAT binding factor can mediate interactions between CONSTANS-like proteins and DNA. Plant J. 2006, 46: 462-76. 10.1111/j.1365-313X.2006.02706.x.
Farré EM, Harmer SL, Harmon FG, Yanovsky MJ, Kay SA: Overlapping and distinct roles of PRR7 and PRR9 in the Arabidopsis circadian clock. Curr Biol. 2005, 15: 47-54. 10.1016/j.cub.2004.12.067.
Nakamichi N, Kita M, Ito S, Yamashino T, Mizuno T: PSEUDO-RESPONSE REGULATORS, PRR9, PRR7 and PRR5, together play essential roles close to the circadian clock of Arabidopsis thaliana. Plant Cell Physiol. 2005, 46: 686-98. 10.1093/pcp/pci086.
Salomé PA, McClung CR: PSEUDO-RESPONSE REGULATOR 7 and 9 are partially redundant genes essential for the temperature responsiveness of the Arabidopsis circadian clock. Plant Cell. 2005, 17: 791-803. 10.1105/tpc.104.029504.
Locke JCW, Southern MM, Kozma-Bognár L, Hibberd V, Brown PE, Turner MS, Millar AJ: Extension of a genetic network model by iterative experimentation and mathematical analysis. Mol Syst Biol. 2005, 1: 2005.0013-10.1038/msb4100018.
Mizoguchi T, Wright L, Fujiwara S, Cremer F, Lee K, Onouchi H, Mouradov A, Fowler S, Kamada H, Putterill J, Coupland G: Distinct roles of GIGANTEA in promoting flowering and regulating circadian rhythms in Arabidopsis. Plant Cell. 2005, 17: 2255-70. 10.1105/tpc.105.033464.
Gould PD, Locke JCW, Larue C, Southern MM, Davis SJ, Hanano S, Moyle R, Milich R, Putterill J, Millar AJ, Hall A: The molecular basis of temperature compensation in the Arabidopsis circadian clock. Plant Cell. 2006, 18: 1177-87. 10.1105/tpc.105.039990.
Más P, Kim WY, Somers DE, Kay SA: Targeted degradation of TOC1 by ZTL modulates circadian function in Arabidopsis thaliana. Nature. 2003, 426: 567-70. 10.1038/nature02163.
Kevei E, Gyula P, Hall A, Kozma-Bognár L, Kim WY, Eriksson ME, Tóth R, Hanano S, Fehér B, Southern MM, Bastow RM, Viczián A, Hibberd V, Davis SJ, Somers DE, Nagy F, Millar AJ: Forward genetic analysis of the circadian clock separates the multiple functions of ZEITLUPE. Plant Physiol. 2006, 140: 933-45. 10.1104/pp.105.074864.
Yasuhara M, Mitsui S, Hirano H, Takanabe R, Tokioka Y, Ihara N, Komatsu A, Seki M, Shinozaki K, Kiyosue T: Identification of ASK and clock-associated proteins as molecular partners of LKP2 (LOV kelch protein 2) in Arabidopsis. J Exp Bot. 2004, 55: 2015-27. 10.1093/jxb/erh226.
Somers DE, Schultz TF, Milnamow M, Kay SA: ZEITLUPE encodes a novel clock-associated PAS protein from Arabidopsis. Cell. 2000, 101: 319-29. 10.1016/S0092-8674(00)80841-7.
Nelson DC, Lasswell J, Rogg LE, Cohen MA, Bartel B: FKF1, a clock-controlled gene that regulates the transition to flowering in Arabidopsis. Cell. 2000, 101: 331-40. 10.1016/S0092-8674(00)80842-9.
Imaizumi T, Tran HG, Swartz TE, Briggs WR, Kay SA: FKF1 is essential for photoperiodic-specific light signalling in Arabidopsis. Nature. 2003, 426: 302-6. 10.1038/nature02090.
Imaizumi T, Schultz TF, Harmon FG, Ho LA, Kay SA: FKF1 F-box protein mediates cyclic degradation of a repressor of CONSTANS in Arabidopsis. Science. 2005, 309: 293-7. 10.1126/science.1110586.
Fornara F, Panigrahi KCS, Gissot L, Sauerbrunn N, Rühl M, Jarillo JA, Coupland G: Arabidopsis DOF transcription factors act redundantly to reduce CONSTANS expression and are essential for a photoperiodic flowering response. Dev Cell. 2009, 17: 75-86. 10.1016/j.devcel.2009.06.015.
Hazen SP, Schultz TF, Pruneda-Paz JL, Borevitz JO, Ecker JR, Kay SA: LUX ARRHYTHMO encodes a Myb domain protein essential for circadian rhythms. Proc Natl Acad Sci USA. 2005, 102: 10387-92. 10.1073/pnas.0503029102.
Doyle MR, Davis SJ, Bastow RM, McWatters HG, Kozma-Bognár L, Nagy F, Millar AJ, Amasino RM: The ELF4 gene controls circadian rhythms and flowering time in Arabidopsis thaliana. Nature. 2002, 419: 74-7. 10.1038/nature00954.
Kolmos E, Nowak M, Werner M, Fischer K, Schwarz G, Mathews S, Schoof H, Nagy F, Bujnicki JM, Davis SJ: Integrating ELF4 into the circadian system through combined structural and functional studies. HSFP Journal. 2009, 3: 350-66.
McWatters HG, Bastow RM, Hall A, Millar AJ: The ELF3 zeitnehmer regulates light signalling to the circadian clock. Nature. 2000, 408: 716-20. 10.1038/35047079.
Covington MF, Panda S, Liu XL, Strayer CA, Wagner DR, Kay SA: ELF3 modulates resetting of the circadian clock in Arabidopsis. Plant Cel l. 2001, 13: 1305-15.
Yu JW, Rubio V, Lee NY, Bai S, Lee SY, Kim SS, Liu L, Zhang Y, Irigoyen ML, Sullivan JA, Zhang Y, Lee I, Xie Q, Paek NC, Deng XW: COP1 and ELF3 control circadian function and photoperiodic flowering by regulating GI stability. Mol Cell. 2008, 32: 617-30. 10.1016/j.molcel.2008.09.026.
Mittag M, Kiaulehn S, Johnson CH: The circadian clock in Chlamydomonas reinhardtii. What is it for? What is it similar to?. Plant Physiology. 2005, 137: 399-409. 10.1104/pp.104.052415.
Matsuo T, Okamoto K, Onai K, Niwa Y, Shimogawara K, Ishiura M: A systematic forward genetic analysis identified components of the Chlamydomonas circadian system. Genes Dev. 2008, 22: 918-30. 10.1101/gad.1650408.
Corellou F, Schwartz C, Motta JP, Djouani-Tahri EB, Sanchez F, Bouget FY: Clocks in the green lineage: comparative functional analysis of the circadian architecture of the picoeukaryote Ostreococcus. Plant Cell. 2009, 21: 3436-49. 10.1105/tpc.109.068825.
Boxall S, Foster J, Bohnert H, Cushman J, Nimmo H, Hartwell J: Conservation and divergence of circadian clock operation in a stress-inducible crassulacean acid metabolism species reveals clock compensation against stress. Plant Physiology. 2005, 137: 969-982. 10.1104/pp.104.054577.
Liu H, Wang H, Gao P, Xu J, Xu T, Wang J, Wang B, Lin C, Fu YF: Analysis of clock gene homologs using unifoliolates as target organs in soybean (Glycine max). J Plant Physiol. 2009, 166: 278-89. 10.1016/j.jplph.2008.06.003.
Kaldis AD, Kousidis P, Kesanopoulos K, Prombona A: Light and circadian regulation in the expression of LHY and Lhcb genes in Phaseolus vulgaris. Plant Mol Biol. 2003, 52: 981-97. 10.1023/A:1025433529082.
Hecht V, Knowles CL, Schoor JKV, Liew LC, Jones SE, Lambert MJM, Weller JL: Pea LATE BLOOMER1 is a GIGANTEA ortholog with roles in photoperiodic flowering, deetiolation, and transcriptional regulation of circadian clock gene homologs. Plant Physiology. 2007, 144: 648-61. 10.1104/pp.107.096818.
Miwa K, Serikawa M, Suzuki S, Kondo T, Oyama T: Conserved expression profiles of circadian clock-related genes in two Lemna species showing long-day and short-day photoperiodic flowering responses. Plant Cell Physiol. 2006, 47: 601-12. 10.1093/pcp/pcj027.
Serikawa M, Miwa K, Kondo T, Oyama T: Functional conservation of clock-related genes in flowering plants: overexpression and RNA interference analyses of the circadian rhythm in the monocotyledon Lemna gibba. Plant Physiology. 2008, 146: 1952-63. 10.1104/pp.107.114611.
Murakami M, Tago Y, Yamashino T, Mizuno T: Comparative overviews of clock-associated genes of Arabidopsis thaliana and Oryza sativa. Plant Cell Physiol. 2007, 48: 110-21. 10.1093/pcp/pcl043.
Ramos A, Pérez-Solís E, Ibañes C, Casado R, Collada C, Gómez L, Aragoncillo C, Allona I: Winter disruption of the circadian clock in chestnut. Proc Natl Acad Sci USA. 2005, 102: 7037-42. 10.1073/pnas.0408549102.
Ibanez C, Ramos A, Acebo P, Contreras A, Casado R, Allona I, Aragoncillo C: Overall alteration of circadian clock gene expression in the chestnut cold response. PLoS ONE. 2008, 3: e3567-10.1371/journal.pone.0003567.
Takata N, Saito S, Saito CT, Nanjo T, Shinohara K, Uemura M: Molecular phylogeny and expression of poplar circadian clock genes, LHY1 and LHY2. New Phytol. 2008
Heuertz M, De Paoli E, Källman T, Larsson H, Jurman I, Morgante M, Lascoux M, Gyllenstrand N: Multilocus patterns of nucleotide diversity, linkage disequilibrium and demographic history of norway spruce [Picea abies (L.) Karst]. Genetics. 2006, 174: 2095-2105. 10.1534/genetics.106.065102.
Kenrick P, Crane PR: The Origin and Early Diversification of Land Plants. A Cladistic Study. 1997, Washington, London: Smithsonian Institution Press
Hohe A, Rensing S, Midner M, Lang D, Reski A: Day length and temperature strongly influence sexual reproduction and expression of a novel MADS-box gene in the moss Physcomitrella patens. Plant Biology. 2002, 4: 595-602. 10.1055/s-2002-35440.
Aoki S, Kato S, Ichikawa K, Shimizu M: Circadian expression of the PpLhcb2 gene encoding a major light-harvesting chlorophyll a/b-binding protein in the moss Physcomitrella patens. Plant Cell Physiol. 2004, 45: 68-76. 10.1093/pcp/pch006.
Ichikawa K, Sugita M, Imaizumi T, Wada M, Aoki S: Differential expression on a daily basis of plastid sigma factor genes from the moss Physcomitrella patens. Regulatory interactions among PpSig5, the circadian clock, and blue light signaling mediated by cryptochromes. Plant Physiology. 2004, 136: 4285-98. 10.1104/pp.104.053033.
Shimizu M, Ichikawa K, Aoki S: Photoperiod-regulated expression of the PpCOL1 gene encoding a homolog of CO/COL proteins in the moss Physcomitrella patens. Biochem Biophys Res Commun. 2004, 324: 1296-301. 10.1016/j.bbrc.2004.09.194.
Zobell O, Coupland G, Reiss B: The family of CONSTANS-like genes in Physcomitrella patens. Plant Biol (Stuttg). 2005, 7: 266-75. 10.1055/s-2005-865621.
Okada R, Kondo S, Satbhai S, Yamaguchi N, Tsukuda M, Aoki S: Functional characterization of CCA1/LHY homolog genes, PpCCA1a and PpCCA1b, in the moss Physcomitrella patens. Plant J. 2009, 60 (3): 551-63. 10.1111/j.1365-313X.2009.03979.x.
Rose A, Meier I, Wienand U: The tomato I-box binding factor LeMYBI is a member of a novel class of myb-like proteins. Plant J. 1999, 20: 641-52. 10.1046/j.1365-313X.1999.00638.x.
Rensing SA, Lang D, Zimmer AD, Terry A, Salamov A, Shapiro H, Nishiyama T, Perroud PF, Lindquist EA, Kamisugi Y, Tanahashi T, Sakakibara K, Fujita T, Oishi K, Shin-I T, Kuroki Y, Toyoda A, Suzuki Y, Hashimoto SI, Yamaguchi K, Sugano S, Kohara Y, Fujiyama A, Anterola A, Aoki S, Ashton N, Barbazuk WB, Barker E, Bennetzen JL, Blankenship R, Cho SH, Dutcher SK, Estelle M, Fawcett JA, Gundlach H, Hanada K, Heyl A, Hicks KA, Hughes J, Lohr M, Mayer K, Melkozernov A, Murata T, Nelson DR, Pils B, Prigge M, Reiss B, Renner T, Rombauts S, Rushton PJ, Sanderfoot A, Schween G, Shiu SH, Stueber K, Theodoulou FL, Tu H, de Peer YV, Verrier PJ, Waters E, Wood A, Yang L, Cove D, Cuming AC, Hasebe M, Lucas S, Mishler BD, Reski R, Grigoriev IV, Quatrano RS, Boore JL: The Physcomitrella genome reveals evolutionary insights into the conquest of land by plants. Science. 2008, 319: 64-9. 10.1126/science.1150646.
Makino S, Kiba T, Imamura A, Hanaki N, Nakamura A, Suzuki T, Taniguchi M, Ueguchi C, Sugiyama T, Mizuno T: Genes encoding pseudo-response regulators: insight into His-to-Asp phosphorelay and circadian rhythm in Arabidopsis thaliana. Plant Cell Physiol. 2000, 41: 791-803.
Murakami M, Ashikari M, Miura K, Yamashino T, Mizuno T: The evolutionarily conserved OsPRR quintet: rice pseudo-response regulators implicated in circadian rhythm. Plant Cell Physiol. 2003, 44: 1229-36. 10.1093/pcp/pcg135.
Hayama R, Izawa T, Shimamoto K: Isolation of rice genes possibly involved in the photoperiodic control of flowering by a fluorescent differential display method. Plant Cell Physiol. 2002, 43: 494-504. 10.1093/pcp/pcf059.
Onai K, Ishiura M: PHYTOCLOCK 1 encoding a novel GARP protein essential for the Arabidopsis circadian clock. Genes Cells. 2005, 10: 963-72. 10.1111/j.1365-2443.2005.00892.x.
Schultz TF, Kiyosue T, Yanovsky M, Wada M, Kay SA: A role for LKP2 in the circadian clock of Arabidopsis. Plant Cell. 2001, 13: 2659-70. 10.1105/tpc.13.12.2659.
Bai C, Sen P, Hofmann K, Ma L, Goebl M, Harper JW, Elledge SJ: SKP1 connects cell cycle regulators to the ubiquitin proteolysis machinery through a novel motif, the F-box. Cell. 1996, 86: 263-74. 10.1016/S0092-8674(00)80098-7.
Adams J, Kelso R, Cooley L: The kelch repeat superfamily of proteins: propellers of cell function. Trends Cell Biol. 2000, 10: 17-24. 10.1016/S0962-8924(99)01673-6.
Edwards KD: Natural Allelic Variation in the Temperature-Compensation Mechanisms of the Arabidopsis thaliana Circadian Clock. Genetics. 2005, 170: 387-400. 10.1534/genetics.104.035238.
Long M, Deutsch M, Wang W, Betrán E, Brunet FG, Zhang J: Origin of new genes: evidence from experimental and computational analyses. Genetica. 2003, 118: 171-82. 10.1023/A:1024153609285.
Babushok DV, Ostertag EM, Kazazian HH: Current topics in genome evolution: molecular mechanisms of new gene formation. Cell Mol Life Sci. 2007, 64: 542-54. 10.1007/s00018-006-6453-4.
Domazet-Loso T, Tautz D: An evolutionary analysis of orphan genes in Drosophila. Genome Res. 2003, 13: 2213-9. 10.1101/gr.1311003.
Levine MT, Jones CD, Kern AD, Lindfors HA, Begun DJ: Novel genes derived from noncoding DNA in Drosophila melanogaster are frequently X-linked and exhibit testis-biased expression. Proc Natl Acad Sci USA. 2006, 103: 9935-9. 10.1073/pnas.0509809103.
Zhou Q, Zhang G, Zhang Y, Xu S, Zhao R, Zhan Z, Li X, Ding Y, Yang S, Wang W: On the origin of new genes in Drosophila. Genome Res. 2008, 18: 1446-55. 10.1101/gr.076588.108.
Heinen TJAJ, Staubach F, Haming D, Tautz D: Emergence of a new gene from an intergenic region. Curr Biol. 2009, 19: 1527-31. 10.1016/j.cub.2009.07.049.
Chen ST, Cheng HC, Barbash DA, Yang HP: Evolution of hydra, a recently evolved testis-expressed gene with nine alternative first exons in Drosophila melanogaster. PLoS Genet. 2007, 3: e107-10.1371/journal.pgen.0030107.
Xiao W, Liu H, Li Y, Li X, Xu C, Long M, Wang S: A rice gene of de novo origin negatively regulates pathogen-induced defense response. PLoS ONE. 2009, 4: e4603-10.1371/journal.pone.0004603.
Richards TA, Soanes DM, Foster PG, Leonard G, Thornton CR, Talbot NJ: Phylogenomic analysis demonstrates a pattern of rare and ancient horizontal gene transfer between plants and fungi. Plant Cell. 2009, 21: 1897-911. 10.1105/tpc.109.065805.
James AB, Monreal JA, Nimmo GA, Kelly CL, Herzyk P, Jenkins GI, Nimmo HG: The circadian clock in Arabidopsis roots is a simplified slave version of the clock in shoots. Science. 2008, 322: 1832-5. 10.1126/science.1161403.
Ito S, Kawamura H, Niwa Y, Nakamichi N, Yamashino T, Mizuno T: A genetic study of the Arabidopsis circadian clock with reference to the TIMING OF CAB EXPRESSION 1 (TOC1) gene. Plant Cell Physiol. 2009, 50: 290-303. 10.1093/pcp/pcn198.
McClung CR: Plant circadian rhythms. Plant Cell. 2006, 18: 792-803. 10.1105/tpc.106.040980.
Rand DA, Shulgin BV, Salazar D, Millar AJ: Design principles underlying circadian clocks. Journal of the Royal Society, Interface/the Royal Society. 2004, 1: 119-30. 10.1098/rsif.2004.0014.
Troein C, Locke JCW, Turner MS, Millar AJ: Weather and seasons together demand complex biological clocks. Current Biology. 2009, 19: 1961-64. 10.1016/j.cub.2009.09.024.
Pruneda-Paz JL, Breton G, Para A, Kay SA: A functional genomics approach reveals CHE as a component of the Arabidopsis circadian clock. Science. 2009, 323: 1481-5. 10.1126/science.1167206.
Ramakers C, Ruijter JM, Lekanne Deprez RH, Moorman AFM: Assumption-free analysis of quantitative real-time polymerase chain reaction (PCR) data. Neuroscience Letters. 2003, 339: 62-66. 10.1016/S0304-3940(02)01423-4.
Straume M: DNA microarray time series analysis: automated statistical assessment of circadian rhythms in gene expression patterning. Meth Enzymol. 2004, 383: 149-66. full_text.
Ceriani MF, Hogenesch JB, Yanovsky M, Panda S, Straume M, Kay SA: Genome-wide expression analysis in Drosophila reveals genes controlling circadian behavior. J Neurosci. 2002, 22: 9305-19.
Panda S, Antoch MP, Miller BH, Su AI, Schook AB, Straume M, Schultz PG, Kay SA, Takahashi JS, Hogenesch JB: Coordinated transcription of key pathways in the mouse by the circadian clock. Cell. 2002, 109: 307-20. 10.1016/S0092-8674(02)00722-5.
Fisher RA: Tests of significance in harmonic analysis. Proc. Soc. A. 1929, 125: 54-59. 10.1098/rspa.1929.0151.
Wichert S, Fokianos K, Strimmer K: Identifying periodically expressed transcripts in microarray time series data. Bioinformatics. 2004, 20: 5-20. 10.1093/bioinformatics/btg364.
Mockler TC, Michael TP, Priest HD, Shen R, Sullivan CM, Givan SA, McEntee C, Kay SA, Chory J: The DIURNAL project: DIURNAL and circadian expression profiling, model-based pattern matching, and promoter analysis. Cold Spring Harbor Symp. Quant. Biol. 2007, 72: 353-363. 10.1101/sqb.2007.72.006.
Lang D, Zimmer AD, Rensing SA, Reski R: Exploring plant biodiversity: the Physcomitrella genome and beyond. Trends Plant Sci. 2008, 13: 542-9. 10.1016/j.tplants.2008.07.002.
Merchant SS, Prochnik SE, Vallon O, Harris EH, Karpowicz SJ, Witman GB, Terry A, Salamov A, Fritz-Laylin LK, Mar'echal-Drouard L, Marshall WF, Qu LH, Nelson DR, Sanderfoot AA, Spalding MH, Kapitonov VV, Ren Q, Ferris P, Lindquist E, Shapiro H, Lucas SM, Grimwood J, Schmutz J, Cardol P, Cerutti H, Chanfreau G, Chen CL, Cognat V, Croft MT, Dent R, Dutcher S, Fernández E, Fukuzawa H, González-Ballester D, González-Halphen D, Hallmann A, Hanikenne M, Hippler M, Inwood W, Jabbari K, Kalanon M, Kuras R, Lefebvre PA, Lemaire SD, Lobanov AV, Lohr M, Manuell A, Meier I, Mets L, Mittag M, Mittelmeier T, Moroney JV, Moseley J, Napoli C, Nedelcu AM, Niyogi K, Novoselov SV, Paulsen IT, Pazour G, Purton S, Ral JP, Riano-Pachón DM, Riekhof W, Rymarquis L, Schroda M, Stern D, Umen J, Willows R, Wilson N, Zimmer SL, Allmer J, Balk J, Bisova K, Chen CJ, Elias M, Gendler K, Hauser C, Lamb MR, Ledford H, Long JC, Minagawa J, Page MD, Pan J, Pootakham W, Ro je S, Rose A, Stahlberg E, Terauchi AM, Yang P, Ball S, Bowler C, Dieckmann CL, Gladyshev VN, Green P, Jorgensen R, Mayfield S, Mueller-Roeber B, Ra jamani S, Sayre RT, Brokstein P, Dubchak I, Goodstein D, Hornick L, Huang YW, Jhaveri J, Luo Y, Martínez D, Ngau WCA, Otillar B, Poliakov A, Porter A, Sza jkowski L, Werner G, Zhou K, Grigoriev IV, Rokhsar DS, Grossman AR: The Chlamydomonas genome reveals the evolution of key animal and plant functions. Science. 2007, 318: 245-50. 10.1126/science.1143609.
Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Research. 2004, 32: 1792-7. 10.1093/nar/gkh340.
Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Systematic Biology. 2003, 52: 696-704. 10.1080/10635150390235520.
Huelsenbeck JP, Ronquist F: MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001, 17: 754-5. 10.1093/bioinformatics/17.8.754.
Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-4. 10.1093/bioinformatics/btg180.
Schultz J, Milpetz F, Bork P, Ponting CP: SMART, a simple modular architecture research tool: identification of signaling domains. Proc Natl Acad Sci USA. 1998, 95: 5857-64. 10.1073/pnas.95.11.5857.
Letunic I, Doerks T, Bork P: SMART 6: recent updates and new developments. Nucleic Acids Research. 2009, D229-32. 10.1093/nar/gkn808. 37 Database
Yang Z: PAML 4: a program package for phylogenetic analysis by maximum likelihood. Molecular Biology and Letters. 2007, 24: 1586-1591. 10.1093/molbev/msm088.
We thank four anonymous reviewers for valuable comments on the manuscript. We thank Kerstin Santesson for running the quantitative RT-PCR assay. UL was funded by the Swedish Research Council and the Swedish Research Council Formas. KH acknowledges financial support from the Nilsson-Ehle foundation.
KH and UL conceived and designed the experiments and wrote the paper. KH carried out bioinformatics searches, phylogenetic reconstructions and the statistical analyses of time series data. TK carried out the PAML analysis. KH, TK, NG and HH carried out the experiments. All authors read, edited and approved the final manuscript.
Electronic supplementary material
Additional file 1: Comparative overview of photoperiodic pathway components associated to clock function in Arabidopsis and P. patens. Comparative overview of photoperiodic pathway components associated to clock function in A. thaliana and putative orthologs in P. patens, followed by the number of putative orthologs in the red alga Cyanidioschyzon merolae, the green algae Ostreococcus tauri and Chlamydomonas reinhardtii and the non-seed vascular plant Selaginella moellendorffii. (DOC 89 KB)
Additional file 2: Statistical analyses of circadian rhythms of gene expression data. Statistical analyses of circadian rhythms in constant darkness (DD) at different temperatures. (DOC 70 KB)
Additional file 3: Test for positive selection. Sequence data and phylogeny used for test of positive selection on PRR genes using PAML 4.3. Likelihood values and parameter estimates. (DOCX 117 KB)
Additional file 4: Period length estimates of gene expression data in different temperatures. Period lengths estimated with COSOPT plotted against temperature for putative clock genes in P. patens. (PDF 260 KB)
Additional file 5: Test of reference gene ( PpTUB1 ) stability. Experimental procedure for measurement of cDNA concentration in all time series and plots of PpTUB1 expression stability over time. (DOCX 385 KB)
Additional file 6: Quantitative RT-PCR primer sequences and time series plots. Names and sequences of primers used in the assay of gene expression and time series plots of gene expression levels including standard deviation from duplicate runs. (PDF 1 MB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Holm, K., Källman, T., Gyllenstrand, N. et al. Does the core circadian clock in the moss Physcomitrella patens (Bryophyta) comprise a single loop?. BMC Plant Biol 10, 109 (2010). https://doi.org/10.1186/1471-2229-10-109
- Circadian Clock
- Clock Gene
- Predict Protein Sequence
- Significant Sequence Homology
- Core Clock Gene