The DIURNAL project: DIURNAL and circadian expression profiling, model-based pattern matching, and promoter analysis

Widespread inter-individual gene expression variability in Arabidopsis thaliana

A fundamental question in biology is how gene expression is regulated to give rise to a phenotype. However, transcriptional variability is rarely considered although it could influence the relationship between genotype and phenotype. It is known in unicellular organisms that gene expression is often noisy rather than uniform, and this has been proposed to be beneficial when environmental conditions are unpredictable. However, little is known about inter-individual transcriptional variability in multicellular organisms. Using transcriptomic approaches, we analysed gene expression variability between individual Arabidopsis thaliana plants growing in identical conditions over a 24-h time course. We identified hundreds of genes that exhibit high inter-individual variability and found that many are involved in environmental responses, with different classes of genes variable between the day and night. We also identified factors that might facilitate gene expression variability, such as gene length, the number of transcription factors regulating the genes and the chromatin environment. These results shed new light on the impact of transcriptional variability in gene expression regulation in plants.

Robust predictions of specialized metabolism genes through machine learning

Specialized metabolites are critical for plant–environment interactions, e.g., attracting pollinators or defending against herbivores, and are important sources of plant-based pharmaceuticals. However, it is unclear what proportion of enzyme-encoding genes play a role in specialized metabolism (SM) as opposed to general metabolism (GM) in any plant species. This is because of the diversity of specialized metabolites and the considerable number of incompletely characterized pathways responsible for their production. In addition, SM gene ancestors frequently played roles in GM. We evaluate features distinguishing SM and GM genes and build a computational model that accurately predicts SM genes. Our predictions provide candidates for experimental studies, and our modeling approach can be applied to other species that produce medicinally or industrially useful compounds.

Plant specialized metabolism (SM) enzymes produce lineage-specific metabolites with important ecological, evolutionary, and biotechnological implications. Using Arabidopsis thaliana as a model, we identified distinguishing characteristics of SM and GM (general metabolism, traditionally referred to as primary metabolism) genes through a detailed study of features including duplication pattern, sequence conservation, transcription, protein domain content, and gene network properties. Analysis of multiple sets of benchmark genes revealed that SM genes tend to be tandemly duplicated, coexpressed with their paralogs, narrowly expressed at lower levels, less conserved, and less well connected in gene networks relative to GM genes. Although the values of each of these features significantly differed between SM and GM genes, any single feature was ineffective at predicting SM from GM genes. Using machine learning methods to integrate all features, a prediction model was established with a true positive rate of 87% and a true negative rate of 71%. In addition, 86% of known SM genes not used to create the machine learning model were predicted. We also demonstrated that the model could be further improved when we distinguished between SM, GM, and junction genes responsible for reactions shared by SM and GM pathways, indicating that topological considerations may further improve the SM prediction model. Application of the prediction model led to the identification of 1,220 A. thaliana genes with previously unknown functions, each assigned a confidence measure called an SM score, providing a global estimate of SM gene content in a plant genome.

Comparative Genomics Analysis Provides New Insight Into Molecular Basis of Stomatal Movement in Kalanchoë fedtschenkoi

CO₂ uptake and water loss in plants are regulated by microscopic pores on the surface of leaves, called stomata. This enablement of gas exchange by the opening and closing of stomata is one of the most essential processes in plant photosynthesis and transpiration, affecting water-use efficiency (WUE) and thus drought susceptibility. In plant species with crassulacean acid metabolism (CAM) photosynthesis, diel stomatal movement pattern is inverted relative to C₃ and C₄ photosynthesis species, resulting in much higher WUE and drought tolerance. However, little is known about the molecular basis of stomatal movement in CAM species. The goal of this study is to identify candidate genes that could play a role in stomatal movement in an obligate CAM species, Kalanchoë fedtschenkoi. By way of a text-mining approach, proteins were identified in various plant species, spanning C₃, C₄, and CAM photosynthetic types, which are orthologous to proteins known to be involved in stomatal movement. A comparative analysis of diel time-course gene expression data was performed between K. fedtschenkoi and two C₃ species (i.e., Arabidopsis thaliana and Solanum lycopersicum) to identify differential gene expression between the dusk and dawn phases of the 24-h cycle. A rescheduled catalase gene known to be involved in stomatal movement was identified, suggesting a role for H₂O₂ in CAM-like stomatal movement. Overall, these results provide new insights into the molecular regulation of stomatal movement in CAM plants, facilitating genetic improvement of drought resistance in agricultural crops through manipulation of stomata-related genes.

EFFECTOR OF TRANSCRIPTION factors are novel plant-specific regulators associated with genomic DNA methylation in Arabidopsis

Plant-specific EFFECTORS OF TRANSCRIPTION (ET) are characterised by a variable number of highly conserved ET repeats, which are involved in zinc and DNA binding. In addition, ETs share a GIY-YIG domain, involved in DNA nicking activity. It was hypothesised that ETs might act as epigenetic regulators. Here, methylome, transcriptome and phenotypic analyses were performed to investigate the role of ET factors and their involvement in DNA methylation in Arabidopsis thaliana. Comparative DNA methylation and transcriptome analyses in flowers and seedlings of et mutants revealed ET-specific differentially expressed genes and mostly independently characteristic, ET-specific differentially methylated regions. Loss of ET function results in pleiotropic developmental defects. The accumulation of cyclobutane pyrimidine dimers after ultraviolet stress in et mutants suggests an ET function in DNA repair.

Conservation and Diversification of Circadian Rhythmicity Between a Model Crassulacean Acid Metabolism Plant Kalanchoë fedtschenkoi and a Model C3 Photosynthesis Plant Arabidopsis thaliana

Crassulacean acid metabolism (CAM) improves photosynthetic efficiency under limited water availability relative to C3 photosynthesis. It is widely accepted that CAM plants have evolved from C3 plants and it is hypothesized that CAM is under the control of the internal circadian clock. However, the role that the circadian clock plays in the evolution of CAM is not well understood. To identify the molecular basis of circadian control over CAM evolution, rhythmic gene sets were identified in a CAM model plant species (Kalanchoë fedtschenkoi) and a C3 model plant species (Arabidopsis thaliana) through analysis of diel time-course gene expression data using multiple periodicity detection algorithms. Based on protein sequences, ortholog groups were constructed containing genes from each of these two species. The ortholog groups were categorized into five gene sets based on conservation and diversification of rhythmic gene expression. Interestingly, minimal functional overlap was observed when comparing the rhythmic gene sets of each species. Specifcally, metabolic processes were enriched in the gene set under circadian control in K. fedtschenkoi and numerous genes were found to have retained or gained rhythmic expression in K. fedtsechenkoi. Additonally, several rhythmic orthologs, including CAM-related orthologs, displayed phase shifts between species. Results of this analysis point to several mechanisms by which the circadian clock plays a role in the evolution of CAM. These genes provide a set of testable hypotheses for future experiments.

Circadian control of abscisic acid biosynthesis and signalling pathways revealed by genome-wide analysis of LHY binding targets

The LATE ELONGATED HYPOCOTYL (LHY) transcription factor functions as part of the oscillatory mechanism of the Arabidopsis circadian clock. This paper reports the genome-wide analysis of its binding targets and reveals a role in the control of abscisic acid (ABA) biosynthesis and downstream responses. LHY directly repressed expression of 9-cis-epoxycarotenoid dioxygenase enzymes, which catalyse the rate-limiting step of ABA biosynthesis. This suggested a mechanism for the circadian control of ABA accumulation in wild-type plants. Consistent with this hypothesis, ABA accumulated rhythmically in wild-type plants, peaking in the evening. LHY-overexpressing plants had reduced levels of ABA under drought stress, whereas loss-of-function mutants exhibited an altered rhythm of ABA accumulation. LHY also bound the promoter of multiple components of ABA signalling pathways, suggesting that it may also act to regulate responses downstream of the hormone. LHY promoted expression of ABA-responsive genes responsible for increased tolerance to drought and osmotic stress but alleviated the inhibitory effect of ABA on seed germination and plant growth. This study reveals a complex interaction between the circadian clock and ABA pathways, which is likely to make an important contribution to plant performance under drought and osmotic stress conditions.

Early Responses to Severe Drought Stress in the Arabidopsis thaliana Cell Suspension Culture Proteome

Abiotic stresses are considered the most deleterious factor affecting growth and development of plants worldwide. Such stresses are largely unavoidable and trigger adaptive responses affecting different cellular processes and target different compartments. Shotgun proteomic and mass spectrometry-based approaches offer an opportunity to elucidate the response of the proteome to abiotic stresses. In this study, the severe drought or water-deficit response in Arabidopsis thaliana was mimicked by treating cell suspension callus with 40% polyethylene glycol for 10 and 30 min. Resulting data demonstrated that 310 proteins were differentially expressed in response to this treatment with a strict ±2.0-fold change. Over-representation was observed in the gene ontology categories of 'ribosome' and its related functions as well as 'oxidative phosphorylation', indicating both structural and functional drought responses at the cellular level. Proteins in the category 'endocytosis' also show significant enrichment and this is consistent with increased active transport and recycling of membrane proteins in response to abiotic stress. This is supported by the particularly pronounced enrichment in proteins of the endosomal sorting complexes that are required for membrane remodelling. Taken together, the findings point to rapid and complex physiological and structural changes essential for survival in response to sudden severe drought stress.

Genome-wide association studies and expression-based quantitative trait loci analyses reveal roles of HCT2 in caffeoylquinic acid biosynthesis and its regulation by defense-responsive transcription factors in Populus

3-O-caffeoylquinic acid, also known as chlorogenic acid (CGA), functions as an intermediate in lignin biosynthesis in the phenylpropanoid pathway. It is widely distributed among numerous plant species and acts as an antioxidant in both plants and animals. Using GC-MS, we discovered consistent and extreme variation in CGA content across a population of 739 4-yr-old Populus trichocarpa accessions. We performed genome-wide association studies (GWAS) from 917 P. trichocarpa accessions and expression-based quantitative trait loci (eQTL) analyses to identify key regulators. The GWAS and eQTL analyses resolved an overlapped interval encompassing a hydroxycinnamoyl-CoA:shikimate hydroxycinnamoyl transferase 2 (PtHCT2) that was significantly associated with CGA and partially characterized metabolite abundances. PtHCT2 leaf expression was significantly correlated with CGA abundance and it was regulated by cis-eQTLs containing W-box for WRKY binding. Among all nine PtHCT homologs, PtHCT2 is the only one that responds to infection by the fungal pathogen Sphaerulina musiva (a Populus pathogen). Validation using protoplast-based transient expression system suggests that PtHCT2 is regulated by the defense-responsive WRKY. These results are consistent with reports of CGA functioning as an antioxidant in response to biotic stress. This study provides insights into data-driven and omics-based inference of gene function in woody species.

Changes in iron availability in Arabidopsis are rapidly sensed in the leaf vasculature and impaired sensing leads to opposite transcriptional programs in leaves and roots

The OLIGOPEPTIDE TRANSPORTER 3 (OPT3) has recently been identified as a component of the systemic network mediating iron (Fe) deficiency responses in Arabidopsis. Reduced expression of OPT3 induces an over accumulation of Fe in roots and leaves, due in part by an elevated expression of the IRON-REGULATED TRANSPORTER 1. Here we show however, that opt3 leaves display a transcriptional program consistent with an Fe overload, suggesting that Fe excess is properly sensed in opt3 leaves and that the OPT3-mediated shoot-to-root signaling is critical to prevent a systemic Fe overload. We also took advantage of the tissue-specific localization of OPT3, together with other Fe-responsive genes, to determine the timing and location of early transcriptional events during Fe limitation and resupply. Our results show that the leaf vasculature responds more rapidly than roots to both Fe deprivation and resupply, suggesting that the leaf vasculature is within the first tissues that sense and respond to changes in Fe availability. Our data highlight the importance of the leaf vasculature in Fe homeostasis by sensing changes in apoplastic levels of Fe coming through the xylem and relaying this information back to roots via the phloem to regulate Fe uptake at the root level.

Using light to improve commercial value

The plasticity of plant morphology has evolved to maximize reproductive fitness in response to prevailing environmental conditions. Leaf architecture elaborates to maximize light harvesting, while the transition to flowering can either be accelerated or delayed to improve an individual's fitness. One of the most important environmental signals is light, with plants using light for both photosynthesis and as an environmental signal. Plants perceive different wavelengths of light using distinct photoreceptors. Recent advances in LED technology now enable light quality to be manipulated at a commercial scale, and as such opportunities now exist to take advantage of plants' developmental plasticity to enhance crop yield and quality through precise manipulation of a crops' lighting regime. This review will discuss how plants perceive and respond to light, and consider how these specific signaling pathways can be manipulated to improve crop yield and quality.

Diel rewiring and positive selection of ancient plant proteins enabled evolution of CAM photosynthesis in Agave

Background

Crassulacean acid metabolism (CAM) enhances plant water-use efficiency through an inverse day/night pattern of stomatal closure/opening that facilitates nocturnal CO₂ uptake. CAM has evolved independently in over 35 plant lineages, accounting for ~ 6% of all higher plants. Agave species are highly heat- and drought-tolerant, and have been domesticated as model CAM crops for beverage, fiber, and biofuel production in semi-arid and arid regions. However, the genomic basis of evolutionary innovation of CAM in genus Agave is largely unknown.

Results

Using an approach that integrated genomics, gene co-expression networks, comparative genomics and protein structure analyses, we investigated the molecular evolution of CAM as exemplified in Agave. Comparative genomics analyses among C₃, C₄ and CAM species revealed that core metabolic components required for CAM have ancient genomic origins traceable to non-vascular plants while regulatory proteins required for diel re-programming of metabolism have a more recent origin shared among C₃, C₄ and CAM species. We showed that accelerated evolution of key functional domains in proteins responsible for primary metabolism and signaling, together with a diel re-programming of the transcription of genes involved in carbon fixation, carbohydrate processing, redox homeostasis, and circadian control is required for the evolution of CAM in Agave. Furthermore, we highlighted the potential candidates contributing to the adaptation of CAM functional modules.

Conclusions

This work provides evidence of adaptive evolution of CAM related pathways. We showed that the core metabolic components required for CAM are shared by non-vascular plants, but regulatory proteins involved in re-reprogramming of carbon fixation and metabolite transportation appeared more recently. We propose that the accelerated evolution of key proteins together with a diel re-programming of gene expression were required for CAM evolution from C₃ ancestors in Agave.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4964-7) contains supplementary material, which is available to authorized users.

Circadian control of ORE1 by PRR9 positively regulates leaf senescence in Arabidopsis

The circadian clock is involved in aging in animals, where mutations in core clock genes accelerate aging. However, little is known about the relationship between aging and the circadian clock in plants. Using the well-studied process of leaf senescence in Arabidopsis, a higher plant, as a model for aging, we show that the circadian clock has a critical role in regulating the aging process in plants. Specifically, we show that PSEUDO-RESPONSE REGULATOR 9 (PRR9), a core clock component, positively regulates leaf senescence. ORESARA 1 (ORE1), an aging regulator, is controlled by PRR9 via direct transcriptional activation and indirectly by suppressing miR164, a posttranscriptional repressor of ORE1, thus forming a coherent feed-forward regulatory loop.

The circadian clock coordinates the daily cyclic rhythm of numerous biological processes by regulating a large portion of the transcriptome. In animals, the circadian clock is involved in aging and senescence, and circadian disruption by mutations in clock genes frequently accelerates aging. Conversely, aging alters circadian rhythmicity, which causes age-associated physiological alterations. However, interactions between the circadian clock and aging have been rarely studied in plants. Here, we investigated potential roles for the circadian clock in the regulation of leaf senescence in plants. Members of the evening complex in Arabidopsis circadian clock, EARLY FLOWERING 3 (ELF3), EARLY FLOWERING 4 (ELF4), and LUX ARRHYTHMO (LUX), as well as the morning component PSEUDO-RESPONSE REGULATOR 9 (PRR9), affect both age-dependent and dark-induced leaf senescence. The circadian clock regulates the expression of several senescence-related transcription factors. In particular, PRR9 binds directly to the promoter of the positive aging regulator ORESARA1 (ORE1) gene to promote its expression. PRR9 also represses miR164, a posttranscriptional repressor of ORE1. Consistently, genetic analysis revealed that delayed leaf senescence of a prr9 mutant was rescued by ORE1 overexpression. Thus, PRR9, a core circadian component, is a key regulator of leaf senescence via positive regulation of ORE1 through a feed-forward pathway involving posttranscriptional regulation by miR164 and direct transcriptional regulation. Our results indicate that, in plants, the circadian clock and leaf senescence are intimately interwoven as are the clock and aging in animals.

Expression of the Intracellular COPT3-Mediated Cu Transport Is Temporally Regulated by the TCP16 Transcription Factor

Copper is an essential element in plants. When scarce, copper is acquired from extracellular environment or remobilized from intracellular sites, through members of the high affinity copper transporters family COPT located at the plasma membrane and internal membrane, respectively. Here, we show that COPT3 is an intracellular copper transporter, located at a compartment of the secretory pathway, that is mainly expressed in pollen grains and vascular bundles. Contrary to the COPT1 plasma membrane member, the expression of the internal COPT3 membrane transporter was higher at 12 h than at 0 h of a neutral photoperiod day under copper deficiency. The screening of a library of conditionally overexpressed transcription factors implicated members of the TCP family in the COPT3 differential temporal expression pattern. Particularly, in vitro, TCP16 was found to bind to the COPT3 promoter and down-regulated its expression. Accordingly, TCP16 was mainly expressed at 0 h under copper deficiency and induced at 12 h by copper excess. Moreover, TCP16 overexpression resulted in increased sensitivity to copper deficiency, whereas the tcp16 mutant was sensitive to copper excess. Both copper content and the expression of particular copper status markers were altered in plants with modified levels of TCP16. Consistent with TCP16 affecting pollen development, the lack of COPT3 function led to altered pollen morphology. Furthermore, analysis of copt3 and COPT3 overexpressing plants revealed that COPT3 function exerted a negative effect on TCP16 expression. Taken together, these results suggest a differential daily regulation of copper uptake depending on the external and internal copper pools, in which TCP16 inhibits copper remobilization at dawn through repression of intracellular transporters.

A 13C isotope labeling method for the measurement of lignin metabolic flux in Arabidopsis stems

BACKGROUND

Metabolic fluxes represent the functional phenotypes of biochemical pathways and are essential to reveal the distribution of precursors among metabolic networks. Although analysis of metabolic fluxes, facilitated by stable isotope labeling and mass spectrometry detection, has been applied in the studies of plant metabolism, we lack experimental measurements for carbon flux towards lignin, one of the most abundant polymers in nature.

RESULTS

We developed a feeding strategy of excised Arabidopsis stems with 13C labeled phenylalanine (Phe) for the analysis of lignin biosynthetic flux. We optimized the feeding methods and found the stems continued to grow and lignify. Consistent with lignification profiles along the stems, higher levels of phenylpropanoids and activities of lignin biosynthetic enzymes were detected in the base of the stem. In the feeding experiments, 13C labeled Phe was quickly accumulated and used for the synthesis of phenylpropanoid intermediates and lignin. The intermediates displayed two different patterns of labeling kinetics during the feeding period. Analysis of lignin showed rapid incorporation of label into all three subunits in the polymers.

CONCLUSIONS

Our feeding results demonstrate the effectiveness of the stem feeding system and suggest a potential application for the investigations of other aspects in plant metabolism. The supply of exogenous Phe leading to a higher lignin deposition rate indicates the availability of Phe is a determining factor for lignification rates.

GEsture: an online hand-drawing tool for gene expression pattern search

Gene expression profiling data provide useful information for the investigation of biological function and process. However, identifying a specific expression pattern from extensive time series gene expression data is not an easy task. Clustering, a popular method, is often used to classify similar expression genes, however, genes with a 'desirable' or 'user-defined' pattern cannot be efficiently detected by clustering methods. To address these limitations, we developed an online tool called GEsture. Users can draw, or graph a curve using a mouse instead of inputting abstract parameters of clustering methods. GEsture explores genes showing similar, opposite and time-delay expression patterns with a gene expression curve as input from time series datasets. We presented three examples that illustrate the capacity of GEsture in gene hunting while following users' requirements. GEsture also provides visualization tools (such as expression pattern figure, heat map and correlation network) to display the searching results. The result outputs may provide useful information for researchers to understand the targets, function and biological processes of the involved genes.

Mutations in EID1 and LNK2 caused light-conditional clock deceleration during tomato domestication

Circadian period and phase of cultivated tomato (Solanum lycopersicum) were changed during domestication, likely adapting the species to its new agricultural environments. Whereas the delayed circadian phase is mainly caused by allelic variation of EID1, the genetic basis of the long circadian period has remained elusive. Here we show that a partial deletion of the clock gene LNK2 is responsible for the period lengthening in cultivated tomatoes. We use resequencing data to phylogenetically classify hundreds of tomato accessions and investigate the evolution of the eid1 and lnk2 mutations along successive domestication steps. We reveal signatures of selection across the genomic region of LNK2 and different patterns of fixation of the mutant alleles. Strikingly, LNK2 and EID1 are both involved in light input to the circadian clock, indicating that domestication specifically targeted this input pathway. In line with this, we show that the clock deceleration in the cultivated tomato is light-dependent and requires the phytochrome B1 photoreceptor. Such conditional variation in circadian rhythms may be key for latitudinal adaptation in a variety of species, including crop plants and livestock.

Temporal transcriptional logic of dynamic regulatory networks underlying nitrogen signaling and use in plants

Our study exploits time—the relatively unexplored fourth dimension of gene regulatory networks (GRNs)—to learn the temporal transcriptional logic underlying dynamic nitrogen (N) signaling in plants. We introduce several conceptual innovations to the analysis of time-series data in the area of predictive GRNs. Our resulting network now provides the “transcriptional logic” for transcription factor perturbations aimed at improving N-use efficiency, an important issue for global food production in marginal soils and for sustainable agriculture. More broadly, the combination of the time-based approaches we develop and deploy can be applied to uncover the temporal “transcriptional logic” for any response system in biology, agriculture, or medicine.

This study exploits time, the relatively unexplored fourth dimension of gene regulatory networks (GRNs), to learn the temporal transcriptional logic underlying dynamic nitrogen (N) signaling in plants. Our “just-in-time” analysis of time-series transcriptome data uncovered a temporal cascade of cis elements underlying dynamic N signaling. To infer transcription factor (TF)-target edges in a GRN, we applied a time-based machine learning method to 2,174 dynamic N-responsive genes. We experimentally determined a network precision cutoff, using TF-regulated genome-wide targets of three TF hubs (CRF4, SNZ, and CDF1), used to “prune” the network to 155 TFs and 608 targets. This network precision was reconfirmed using genome-wide TF-target regulation data for four additional TFs (TGA1, HHO5/6, and PHL1) not used in network pruning. These higher-confidence edges in the GRN were further filtered by independent TF-target binding data, used to calculate a TF “N-specificity” index. This refined GRN identifies the temporal relationship of known/validated regulators of N signaling (NLP7/8, TGA1/4, NAC4, HRS1, and LBD37/38/39) and 146 additional regulators. Six TFs—CRF4, SNZ, CDF1, HHO5/6, and PHL1—validated herein regulate a significant number of genes in the dynamic N response, targeting 54% of N-uptake/assimilation pathway genes. Phenotypically, inducible overexpression of CRF4 in planta regulates genes resulting in altered biomass, root development, and ¹⁵NO₃⁻ uptake, specifically under low-N conditions. This dynamic N-signaling GRN now provides the temporal “transcriptional logic” for 155 candidate TFs to improve nitrogen use efficiency with potential agricultural applications. Broadly, these time-based approaches can uncover the temporal transcriptional logic for any biological response system in biology, agriculture, or medicine.

Time is honey: circadian clocks of bees and flowers and how their interactions may influence ecological communities

The interactions between flowering plants and insect pollinators shape ecological communities and provide one of the best examples of coevolution. Although these interactions have received much attention in both ecology and evolution, their temporal aspects are little explored. Here we review studies on the circadian organization of pollination-related traits in bees and flowers. Research, mostly with the honeybee, Apis mellifera, has implicated the circadian clock in key aspects of their foraging for flower rewards. These include anticipation, timing of visits to flowers at specified locations and time-compensated sun-compass orientation. Floral rhythms in traits such as petal opening, scent release and reward availability also show robust daily rhythms. However, in only few studies was it possible to adequately determine whether these oscillations are driven by external time givers such as light and temperature cycles, or endogenous circadian clocks. The interplay between the timing of flower and pollinator rhythms may be ecologically significant. Circadian regulation of pollination-related traits in only few species may influence the entire pollination network and thus affect community structure and local biodiversity. We speculate that these intricate chronobiological interactions may be vulnerable to anthropogenic effects such as the introduction of alien invasive species, pesticides or environmental pollutants.This article is part of the themed issue 'Wild clocks: integrating chronobiology and ecology to understand timekeeping in free-living animals'.

Diurnal down-regulation of ethylene biosynthesis mediates biomass heterosis

Heterosis is widely applied in agriculture; however, the underlying molecular mechanisms for superior performance are not well understood. Ethylene biosynthesis and signaling genes are shown to be down-regulated in Arabidopsis interspecific hybrids. Ethylene is a plant hormone that promotes fruit ripening and maturation but inhibits hypocotyl elongation. Here we report that application of exogenous ethylene could eliminate biomass vigor in Arabidopsis thaliana F1 hybrids, suggesting a negative role of ethylene in heterosis. Ethylene biosynthesis is mediated by the rate-limiting enzyme, 1-aminocyclopropane-1-carboxylate synthase (ACS). Down-regulation of ACS genes led to the decrease of ethylene production, which was associated with the high-vigor F1 hybrids, but not with the low-vigor ones. At the mechanistic level, expression of ACS genes was down-regulated diurnally and indirectly by Circadian Clock Associated 1 (CCA1) during the day and directly by Phyotochrome-Interacting Factor 5 (PIF5) at night. Consistent with the negative role of ethylene in plant growth, biomass vigor was higher in the acs mutants than in wild-type plants, while increasing endogenous ethylene production in the hybridizing parents reduced growth vigor in the hybrids. Thus, integrating circadian rhythms and light signaling into ethylene production is another regulatory module of complex biological networks, leading to biomass heterosis in plants.

ZINC-FINGER interactions mediate transcriptional regulation of hypocotyl growth in Arabidopsis

Light coordinates energy production, growth, and survival throughout plant development. In Arabidopsis, light stimulates transcriptional reprogramming during developmental transitions such as photomorphogenesis and flowering through the action of photoreceptors, transcription factors, and signaling components. Here we assign a function to a member of the zinc-finger homeodomain (ZFHD) transcription factor family in regulating light-induced development. Our findings reveal ZFHD10 to be a missing link in understanding how the recently discovered integrator of light and photoperiodic flowering, TANDEM ZINC-FINGER PLUS3 (TZP), controls the expression of growth-promoting transcriptional regulators via direct association with light-regulated promoter elements. Elucidating how such novel protein complexes coordinate gene expression will allow scientists and breeders to optimize plant growth and development in response to unfavorable environmental conditions.

Integration of environmental signals and interactions among photoreceptors and transcriptional regulators is key in shaping plant development. TANDEM ZINC-FINGER PLUS3 (TZP) is an integrator of light and photoperiodic signaling that promotes flowering in Arabidopsis thaliana. Here we elucidate the molecular role of TZP as a positive regulator of hypocotyl elongation. We identify an interacting partner for TZP, the transcription factor ZINC-FINGER HOMEODOMAIN 10 (ZFHD10), and characterize its function in coregulating the expression of blue-light–dependent transcriptional regulators and growth-promoting genes. By employing a genome-wide approach, we reveal that ZFHD10 and TZP coassociate with promoter targets enriched in light-regulated elements. Furthermore, using a targeted approach, we show that ZFHD10 recruits TZP to the promoters of key coregulated genes. Our findings not only unveil the mechanism of TZP action in promoting hypocotyl elongation at the transcriptional level but also assign a function to an uncharacterized member of the ZFHD transcription factor family in promoting plant growth.

Photoperiodic control of the Arabidopsis proteome reveals a translational coincidence mechanism

Plants respond to seasonal cues such as the photoperiod, to adapt to current conditions and to prepare for environmental changes in the season to come. To assess photoperiodic responses at the protein level, we quantified the proteome of the model plant Arabidopsis thaliana by mass spectrometry across four photoperiods. This revealed coordinated changes of abundance in proteins of photosynthesis, primary and secondary metabolism, including pigment biosynthesis, consistent with higher metabolic activity in long photoperiods. Higher translation rates in the day than the night likely contribute to these changes, via an interaction with rhythmic changes in RNA abundance. Photoperiodic control of protein levels might be greatest only if high translation rates coincide with high transcript levels in some photoperiods. We term this proposed mechanism "translational coincidence", mathematically model its components, and demonstrate its effect on the Arabidopsis proteome. Datasets from a green alga and a cyanobacterium suggest that translational coincidence contributes to seasonal control of the proteome in many phototrophic organisms. This may explain why many transcripts but not their cognate proteins exhibit diurnal rhythms.

Dark period transcriptomic and metabolic profiling of two diverse Eutrema salsugineum accessions

Eutrema salsugineum is a model species for the study of plant adaptation to abiotic stresses. Two accessions of E. salsugineum, Shandong (SH) and Yukon (YK), exhibit contrasting morphology and biotic and abiotic stress tolerance. Transcriptome profiling and metabolic profiling from tissue samples collected during the dark period were used to investigate the molecular and metabolic bases of these contrasting phenotypes. RNA sequencing identified 17,888 expressed genes, of which 157 were not in the published reference genome, and 65 of which were detected for the first time. Differential expression was detected for only 31 genes. The RNA sequencing data contained 14,808 single nucleotide polymorphisms (SNPs) in transcripts, 3,925 of which are newly identified. Among the differentially expressed genes, there were no obvious candidates for the physiological or morphological differences between SH and YK. Metabolic profiling indicated that YK accumulates free fatty acids and long-chain fatty acid derivatives as compared to SH, whereas sugars are more abundant in SH. Metabolite levels suggest that carbohydrate and respiratory metabolism, including starch degradation, is more active during the first half of the dark period in SH. These metabolic differences may explain the greater biomass accumulation in YK over SH. The accumulation of 56% of the identified metabolites was lower in F1 hybrids than the mid-parent averages and the accumulation of 17% of the metabolites in F1 plants transgressed the level in both parents. Concentrations of several metabolites in F1 hybrids agree with previous studies and suggest a role for primary metabolism in heterosis. The improved annotation of the E. salsugineum genome and newly identified high-quality SNPs will permit accelerated studies using the standing variation in this species to elucidate the mechanisms of its diverse adaptations to the environment.

Multiple links between shade avoidance and auxin networks

Auxin has emerged as a key player in the adjustment of plant morphology to the challenge imposed by variable environmental conditions. Shade-avoidance responses, including the promotion of stem and petiole growth, leaf hyponasty, and the inhibition of branching, involve an intimate connection between light and auxin signalling. Low activity of photo-sensory receptors caused by the presence of neighbouring vegetation enhances the activity of PHYTOCHROME INTERACTING FACTORs (PIFs), which directly promote the expression of genes involved in auxin biosynthesis, conjugation, transport, perception, and signalling. In seedlings, neighbour signals increase auxin levels in the foliage, which then moves to the stem, where it reaches epidermal tissues to promote growth. However, this model only partially accounts for shade-avoidance responses (which may also occur in the absence of increased auxin levels), and understanding the whole picture will require further insight into the functional significance of the multiple links between shade and auxin networks.

Computational exploration of cis-regulatory modules in rhythmic expression data using the "Exploration of Distinctive CREs and CRMs" (EDCC) and "CRM Network Generator" (CNG) programs

Understanding the effect of cis-regulatory elements (CRE) and clusters of CREs, which are called cis-regulatory modules (CRM), in eukaryotic gene expression is a challenge of computational biology. We developed two programs that allow simple, fast and reliable analysis of candidate CREs and CRMs that may affect specific gene expression and that determine positional features between individual CREs within a CRM. The first program, "Exploration of Distinctive CREs and CRMs" (EDCC), correlates candidate CREs and CRMs with specific gene expression patterns. For pairs of CREs, EDCC also determines positional preferences of the single CREs in relation to each other and to the transcriptional start site. The second program, "CRM Network Generator" (CNG), prioritizes these positional preferences using a neural network and thus allows unbiased rating of the positional preferences that were determined by EDCC. We tested these programs with data from a microarray study of circadian gene expression in Arabidopsis thaliana. Analyzing more than 1.5 million pairwise CRE combinations, we found 22 candidate combinations, of which several contained known clock promoter elements together with elements that had not been identified as relevant to circadian gene expression before. CNG analysis further identified positional preferences of these CRE pairs, hinting at positional information that may be relevant for circadian gene expression. Future wet lab experiments will have to determine which of these combinations confer daytime specific circadian gene expression.

Photoperiodic Regulation of Shoot Apical Growth in Poplar

Woody perennials adapt their genetic traits to local climate conditions. Day length plays an essential role in the seasonal growth of poplar trees. When photoperiod falls below a given critical day length, poplars undergo growth cessation and bud set. A leaf-localized mechanism of photoperiod measurement triggers the transcriptional modulation of a long distance signaling molecule, FLOWERING LOCUS T (FT). This molecule targets meristem function giving rise to these seasonal responses. Studies over the past decade have identified conserved orthologous genes involved in photoperiodic flowering in Arabidopsis that regulate poplar vegetative growth. However, phenological and molecular examination of key photoperiod signaling molecules reveals functional differences between these two plant model systems suggesting alternative components and/or regulatory mechanisms operating during poplar vegetative growth. Here, we review current knowledge and provide new data regarding the molecular components of the photoperiod measuring mechanism that regulates annual growth in poplar focusing on main achievements and new perspectives.

Environmentally Sensitive Molecular Switches Drive Poplar Phenology

Boreal and temperate woody perennials are highly adapted to their local climate, which delimits the length of the growing period. Moreover, seasonal control of growth-dormancy cycles impacts tree productivity and geographical distribution. Therefore, traits related to phenology are of great interest to tree breeders and particularly relevant in the context of global warming. The recent application of transcriptional profiling and genetic association studies to poplar species has provided a robust molecular framework for investigating molecules with potential links to phenology. The environment dictates phenology by modulating the expression of endogenous molecular switches, the identities of which are currently under investigation. This review outlines the current knowledge of these molecular switches in poplar and covers several perspectives concerning the environmental control of growth-dormancy cycles. In the process, we highlight certain genetic pathways which are affected by short days, low temperatures and cold-induced signaling.

The Kalanchoë genome provides insights into convergent evolution and building blocks of crassulacean acid metabolism

Crassulacean acid metabolism (CAM) is a water-use efficient adaptation of photosynthesis that has evolved independently many times in diverse lineages of flowering plants. We hypothesize that convergent evolution of protein sequence and temporal gene expression underpins the independent emergences of CAM from C3 photosynthesis. To test this hypothesis, we generate a de novo genome assembly and genome-wide transcript expression data for Kalanchoë fedtschenkoi, an obligate CAM species within the core eudicots with a relatively small genome (~260 Mb). Our comparative analyses identify signatures of convergence in protein sequence and re-scheduling of diel transcript expression of genes involved in nocturnal CO2 fixation, stomatal movement, heat tolerance, circadian clock, and carbohydrate metabolism in K. fedtschenkoi and other CAM species in comparison with non-CAM species. These findings provide new insights into molecular convergence and building blocks of CAM and will facilitate CAM-into-C3 photosynthesis engineering to enhance water-use efficiency in crops.

What makes ribosomes tick?

In most organisms, gene expression over the course of the day is under the control of the circadian clock. The canonical clock operates as a gene expression circuit that is controlled at the level of transcription, and transcriptional control is also a major clock output. However, rhythmic transcription cannot explain all the observed rhythms in protein accumulation. Although it is clear that rhythmic gene expression also involves RNA processing and protein turnover, until two years ago little was known in any eukaryote about diel dynamics of mRNA translation into protein. A recent series of studies in animals and plants demonstrated that diel cycles of translation efficiency are widespread across the tree of life and its transcriptomes. There are surprising parallels between the patterns of diel translation in mammals and plants. For example, ribosomal proteins and mitochondrial proteins are under translational control in mouse liver, human tissue culture, and Arabidopsis seedlings. In contrast, the way in which the circadian clock, light-dark changes, and other environmental factors such as nutritional signals interact to drive the cycles of translation may differ between organisms. Further investigation is needed to identify the signaling pathways, biochemical mechanisms, RNA sequence features, and the physiological implications of diel translation.

Guidelines for Genome-Scale Analysis of Biological Rhythms

Genome biology approaches have made enormous contributions to our understanding of biological rhythms, particularly in identifying outputs of the clock, including RNAs, proteins, and metabolites, whose abundance oscillates throughout the day. These methods hold significant promise for future discovery, particularly when combined with computational modeling. However, genome-scale experiments are costly and laborious, yielding “big data” that are conceptually and statistically difficult to analyze. There is no obvious consensus regarding design or analysis. Here we discuss the relevant technical considerations to generate reproducible, statistically sound, and broadly useful genome-scale data. Rather than suggest a set of rigid rules, we aim to codify principles by which investigators, reviewers, and readers of the primary literature can evaluate the suitability of different experimental designs for measuring different aspects of biological rhythms. We introduce CircaInSilico, a web-based application for generating synthetic genome biology data to benchmark statistical methods for studying biological rhythms. Finally, we discuss several unmet analytical needs, including applications to clinical medicine, and suggest productive avenues to address them.

Shining a light on the Arabidopsis circadian clock

The circadian clock provides essential timing information to ensure optimal growth to prevailing external environmental conditions. A major time-setting mechanism (zeitgeber) in clock synchronization is light. Differing light wavelengths, intensities, and photoperiodic duration are processed for the clock-setting mechanism. Many studies on light-input pathways to the clock have focused on Arabidopsis thaliana. Photoreceptors are specific chromic proteins that detect light signals and transmit this information to the central circadian oscillator through a number of different signalling mechanisms. The most well-characterized clock-mediating photoreceptors are cryptochromes and phytochromes, detecting blue, red, and far-red wavelengths of light. Ultraviolet and shaded light are also processed signals to the oscillator. Notably, the clock reciprocally generates rhythms of photoreceptor action leading to so-called gating of light responses. Intermediate proteins, such as Phytochrome interacting factors (PIFs), constitutive photomorphogenic 1 (COP1) and EARLY FLOWERING 3 (ELF3), have been established in signalling pathways downstream of photoreceptor activation. However, the precise details for these signalling mechanisms are not fully established. This review highlights both historical and recent efforts made to understand overall light input to the oscillator, first looking at how each wavelength of light is detected, this is then related to known input mechanisms and their interactions.

FE Controls the Transcription of Downstream Flowering Regulators Through Two Distinct Mechanisms in Leaf Phloem Companion Cells

In the facultative long-day plant Arabidopsis thaliana, FLOWERING LOCUS T (FT), encoding the mobile hormone florigen, plays an essential role in modulating the optimal timing of flowering to ensure reproductive success. Under inductive long-day conditions, the transcription of FT is activated by the CONSTANS (CO)/NUCLEAR FACTOR-Y (NF-Y) protein complex in leaf phloem companion cells. FT is transported to the shoot apical meristem through interaction with florigen transporters, such as SODIUM POTASSIUM ROOT DEFECTIVE 1 (NaKR1). Some regulators involved in photoperiod-dependent FT function have been reported previously; however, the molecular mechanism that coordinates FT protein synthesis and transport efficiently needs to be investigated. The present study examined the role of an Myb-related transcription factor, FE, in the activation of FT gene transcription and FT protein transport. Expression analysis using FE-inducible systems and chromatin immunoprecipitation assays showed that FE directly bound to the FT and NaKR1 promoters and activated the transcription of downstream target genes. FE failed to activate FT expression without CO function, whereas FE-mediated NaKR1 induction was not affected by CO function. Taken together, our data indicate that FE regulates the transcription of FT and florigen transporter genes via different mechanisms.

Cross-species complementation reveals conserved functions for EARLY FLOWERING 3 between monocots and dicots

Plant responses to the environment are shaped by external stimuli and internal signaling pathways. In both the model plant Arabidopsis thaliana (Arabidopsis) and crop species, circadian clock factors are critical for growth, flowering, and circadian rhythms. Outside of Arabidopsis, however, little is known about the molecular function of clock gene products. Therefore, we sought to compare the function of Brachypodium distachyon (Brachypodium) and Setaria viridis (Setaria) orthologs of EARLY FLOWERING 3, a key clock gene in Arabidopsis. To identify both cycling genes and putative ELF3 functional orthologs in Setaria, a circadian RNA-seq dataset and online query tool (Diel Explorer) were generated to explore expression profiles of Setaria genes under circadian conditions. The function of ELF3 orthologs from Arabidopsis, Brachypodium, and Setaria was tested for complementation of an elf3 mutation in Arabidopsis. We find that both monocot orthologs were capable of rescuing hypocotyl elongation, flowering time, and arrhythmic clock phenotypes. Using affinity purification and mass spectrometry, our data indicate that BdELF3 and SvELF3 could be integrated into similar complexes in vivo as AtELF3. Thus, we find that, despite 180 million years of separation, BdELF3 and SvELF3 can functionally complement loss of ELF3 at the molecular and physiological level.

Temporal and spatial transcriptomic and microRNA dynamics of CAM photosynthesis in pineapple

The altered carbon assimilation pathway of crassulacean acid metabolism (CAM) photosynthesis results in an up to 80% higher water-use efficiency than C₃ photosynthesis in plants making it a potentially useful pathway for engineering crop plants with improved drought tolerance. Here we surveyed detailed temporal (diel time course) and spatial (across a leaf gradient) gene and microRNA (miRNA) expression patterns in the obligate CAM plant pineapple [Ananas comosus (L.) Merr.]. The high-resolution transcriptome atlas allowed us to distinguish between CAM-related and non-CAM gene copies. A differential gene co-expression network across green and white leaf diel datasets identified genes with circadian oscillation, CAM-related functions, and source-sink relations. Gene co-expression clusters containing CAM pathway genes are enriched with clock-associated cis-elements, suggesting circadian regulation of CAM. About 20% of pineapple microRNAs have diel expression patterns, with several that target key CAM-related genes. Expression and physiology data provide a model for CAM-specific carbohydrate flux and long-distance hexose transport. Together these resources provide a list of candidate genes for targeted engineering of CAM into C₃ photosynthesis crop species.

Diurnal Cycling Transcription Factors of Pineapple Revealed by Genome-Wide Annotation and Global Transcriptomic Analysis

Circadian clock provides fitness advantage by coordinating internal metabolic and physiological processes to external cyclic environments. Core clock components exhibit daily rhythmic changes in gene expression, and the majority of them are transcription factors (TFs) and transcription coregulators (TCs). We annotated 1,398 TFs from 67 TF families and 80 TCs from 20 TC families in pineapple, and analyzed their tissue-specific and diurnal expression patterns. Approximately 42% of TFs and 45% of TCs displayed diel rhythmic expression, including 177 TF/TCs cycling only in the nonphotosynthetic leaf tissue, 247 cycling only in the photosynthetic leaf tissue, and 201 cycling in both. We identified 68 TF/TCs whose cycling expression was tightly coupled between the photosynthetic and nonphotosynthetic leaf tissues. These TF/TCs likely coordinate key biological processes in pineapple as we demonstrated that this group is enriched in homologous genes that form the core circadian clock in Arabidopsis and includes a STOP1 homolog. Two lines of evidence support the important role of the STOP1 homolog in regulating CAM photosynthesis in pineapple. First, STOP1 responds to acidic pH and regulates a malate channel in multiple plant species. Second, the cycling expression pattern of the pineapple STOP1 and the diurnal pattern of malate accumulation in pineapple leaf are correlated. We further examined duplicate-gene retention and loss in major known circadian genes and refined their evolutionary relationships between pineapple and other plants. Significant variations in duplicate-gene retention and loss were observed for most clock genes in both monocots and dicots.

Plant dual-specificity tyrosine phosphorylation-regulated kinase optimizes light-regulated growth and development in Arabidopsis

Light controls vegetative and reproductive development of plants. For a plant, sensing the light input properly ensures coordination with the ever-changing environment. Previously, we found that LIGHT-REGULATED WD1 (LWD1) and LWD2 regulate the circadian clock and photoperiodic flowering. Here, we identified Arabidopsis YET ANOTHER KINASE1 (AtYAK1), an evolutionarily conserved protein and a member of dual-specificity tyrosine phosphorylation-regulated kinases (DYRKs), as an interacting protein of LWDs. Our study revealed that AtYAK1 is an important regulator for various light responses, including the circadian clock, photomorphogenesis and reproductive development. AtYAK1 could antagonize the function of LWDs in regulating the circadian clock and photoperiodic flowering. By examining phenotypes of atyak1, we found that AtYAK1 regulated light-induced period-length shortening and photomorphogenic development. Moreover, AtYAK1 mediated plant fertility especially under inferior light conditions including low light and short-day length. This study discloses a new regulator connecting environmental light to plant growth.

Crosstalk between diurnal rhythm and water stress reveals an altered primary carbon flux into soluble sugars in drought-treated rice leaves

Plants retain rhythmic physiological responses when adapting to environmental challenges. However, possible integrations between drought conditions and those responses have not received much focus, especially regarding crop plants, and the relationship between abiotic stress and the diurnal cycle is generally not considered. Therefore, we conducted a genome-wide analysis to identify genes showing both diurnal regulation and water-deficiency response in rice (Oryza sativa). Among the 712 drought-responsive genes primary identified, 56.6% are diurnally expressed while 47.6% of the 761 that are down-regulated by drought are also diurnal. Using the β-glucuronidase reporter system and qRT-PCR analyses, we validated expression patterns of two candidate genes, thereby supporting the reliability of our transcriptome data. MapMan analysis indicated that diurnal genes up-regulated by drought are closely associated with the starch-sucrose pathway while those that are down-regulated are involved in photosynthesis. We then confirmed that starch-sucrose contents and chlorophyll fluorescence are altered in a diurnal manner under drought stress, suggesting these metabolic diurnal alterations as a novel indicator to evaluate the drought response in rice leaves. We constructed a functional gene network associated with the starch-sucrose KEGG metabolic pathway for further functional studies, and also developed a regulatory pathway model that includes OsbZIP23 transcription factor.

Heterologous protein-DNA interactions lead to biased allelic expression of circadian clock genes in interspecific hybrids

Genomic interactions in allopolyploids create expression variation of homoeologous alleles through protein-protein and protein-DNA interactions. However, the molecular basis for this is largely unknown. Here we investigated the protein-protein and protein-DNA interactions among homoeologous transcription factors in the circadian-clock feedback loop, consisting of CCA1 HIKING EXPEDITION (CHE), CIRCADIAN CLOCK ASSOCIATED1 (CCA1), and TIMING OF CAB EXPRESSION1 (TOC1), plus the interaction with a chromatin factor, HISTONE DEACETYLASE1 (HD1). In the allotetraploids formed between A. thaliana (At) and Arabidopsis arenosa (Aa), AtCCA1 is expressed at lower levels than AaCCA1, which could alter clock output traits. The reduced AtCCA1 expressions in the allotetraploids are consistent with the biochemical data that AaCHE showed preferential binding to the AtCCA1 promoter, in which AaCHE interacts with a higher affinity to AtHD1 than AtCHE. AaCHE also showed a higher affinity to TOC1 than AtCHE, consistent with the effect of TOC1 on repressing CCA1. Thus, stronger AaCHE-TOC1 and AaCHE-AtHD1 interactions reduce AtCC1 allelic expression. Our current data suggest a biochemical basis for protein interactions in trans with a preference to the cis-acting elements in heterologous combinations to reduce AtCCA1 expression, while altered CCA1 expression has been shown to affect metabolic and biomass heterosis in interspecific hybrids or allotetraploids.

Comparative Analysis of Vertebrate Diurnal/Circadian Transcriptomes

From photosynthetic bacteria to mammals, the circadian clock evolved to track diurnal rhythms and enable organisms to anticipate daily recurring changes such as temperature and light. It orchestrates a broad spectrum of physiology such as the sleep/wake and eating/fasting cycles. While we have made tremendous advances in our understanding of the molecular details of the circadian clock mechanism and how it is synchronized with the environment, we still have rudimentary knowledge regarding its connection to help regulate diurnal physiology. One potential reason is the sheer size of the output network. Diurnal/circadian transcriptomic studies are reporting that around 10% of the expressed genome is rhythmically controlled. Zebrafish is an important model system for the study of the core circadian mechanism in vertebrate. As Zebrafish share more than 70% of its genes with human, it could also be an additional model in addition to rodent for exploring the diurnal/circadian output with potential for translational relevance. Here we performed comparative diurnal/circadian transcriptome analysis with established mouse liver and other tissue datasets. First, by combining liver tissue sampling in a 48h time series, transcription profiling using oligonucleotide arrays and bioinformatics analysis, we profiled rhythmic transcripts and identified 2609 rhythmic genes. The comparative analysis revealed interesting features of the output network regarding number of rhythmic genes, proportion of tissue specific genes and the extent of transcription factor family expression. Undoubtedly, the Zebrafish model system will help identify new vertebrate outputs and their regulators and provides leads for further characterization of the diurnal cis-regulatory network.

Phosphorylation of Ribosomal Protein RPS6 Integrates Light Signals and Circadian Clock Signals

The translation of mRNA into protein is tightly regulated by the light environment as well as by the circadian clock. Although changes in translational efficiency have been well documented at the level of mRNA-ribosome loading, the underlying mechanisms are unclear. The reversible phosphorylation of RIBOSOMAL PROTEIN OF THE SMALL SUBUNIT 6 (RPS6) has been known for 40 years, but the biochemical significance of this event remains unclear to this day. Here, we confirm using a clock-deficient strain of Arabidopsis thaliana that RPS6 phosphorylation (RPS6-P) is controlled by the diel light-dark cycle with a peak during the day. Strikingly, when wild-type, clock-enabled, seedlings that have been entrained to a light-dark cycle are placed under free-running conditions, the circadian clock drives a cycle of RPS6-P with an opposite phase, peaking during the subjective night. We show that in wild-type seedlings under a light-dark cycle, the incoherent light and clock signals are integrated by the plant to cause an oscillation in RPS6-P with a reduced amplitude with a peak during the day. Sucrose can stimulate RPS6-P, as seen when sucrose in the medium masks the light response of etiolated seedlings. However, the diel cycles of RPS6-P are observed in the presence of 1% sucrose and in its absence. Sucrose at a high concentration of 3% appears to interfere with the robust integration of light and clock signals at the level of RPS6-P. Finally, we addressed whether RPS6-P occurs uniformly in polysomes, non-polysomal ribosomes and their subunits, and non-ribosomal protein. It is the polysomal RPS6 whose phosphorylation is most highly stimulated by light and repressed by darkness. These data exemplify a striking case of contrasting biochemical regulation between clock signals and light signals. Although the physiological significance of RPS6-P remains unknown, our data provide a mechanistic basis for the future understanding of this enigmatic event.

Evolution of Daily Gene Co-expression Patterns from Algae to Plants

Daily rhythms play a key role in transcriptome regulation in plants and microalgae orchestrating responses that, among other processes, anticipate light transitions that are essential for their metabolism and development. The recent accumulation of genome-wide transcriptomic data generated under alternating light:dark periods from plants and microalgae has made possible integrative and comparative analysis that could contribute to shed light on the evolution of daily rhythms in the green lineage. In this work, RNA-seq and microarray data generated over 24 h periods in different light regimes from the eudicot Arabidopsis thaliana and the microalgae Chlamydomonas reinhardtii and Ostreococcus tauri have been integrated and analyzed using gene co-expression networks. This analysis revealed a reduction in the size of the daily rhythmic transcriptome from around 90% in Ostreococcus, being heavily influenced by light transitions, to around 40% in Arabidopsis, where a certain independence from light transitions can be observed. A novel Multiple Bidirectional Best Hit (MBBH) algorithm was applied to associate single genes with a family of potential orthologues from evolutionary distant species. Gene duplication, amplification and divergence of rhythmic expression profiles seems to have played a central role in the evolution of gene families in the green lineage such as Pseudo Response Regulators (PRRs), CONSTANS-Likes (COLs), and DNA-binding with One Finger (DOFs). Gene clustering and functional enrichment have been used to identify groups of genes with similar rhythmic gene expression patterns. The comparison of gene clusters between species based on potential orthologous relationships has unveiled a low to moderate level of conservation of daily rhythmic expression patterns. However, a strikingly high conservation was found for the gene clusters exhibiting their highest and/or lowest expression value during the light transitions.

Order restricted inference for oscillatory systems for detecting rhythmic signals

MOTIVATION

Many biological processes, such as cell cycle, circadian clock, menstrual cycles, are governed by oscillatory systems consisting of numerous components that exhibit rhythmic patterns over time. It is not always easy to identify such rhythmic components. For example, it is a challenging problem to identify circadian genes in a given tissue using time-course gene expression data. There is a great potential for misclassifying non-rhythmic as rhythmic genes and vice versa. This has been a problem of considerable interest in recent years. In this article we develop a constrained inference based methodology called Order Restricted Inference for Oscillatory Systems (ORIOS) to detect rhythmic signals. Instead of using mathematical functions (e.g. sinusoidal) to describe shape of rhythmic signals, ORIOS uses mathematical inequalities. Consequently, it is robust and not limited by the biologist's choice of the mathematical model. We studied the performance of ORIOS using simulated as well as real data obtained from mouse liver, pituitary gland and data from NIH3T3, U2OS cell lines. Our results suggest that, for a broad collection of patterns of gene expression, ORIOS has substantially higher power to detect true rhythmic genes in comparison to some popular methods, while also declaring substantially fewer non-rhythmic genes as rhythmic.

AVAILABILITY AND IMPLEMENTATION

A user friendly code implemented in R language can be downloaded from http://www.niehs.nih.gov/research/atniehs/labs/bb/staff/peddada/index.cfm CONTACT: peddada@niehs.nih.gov.

Utility and Limitations of Using Gene Expression Data to Identify Functional Associations

Gene co-expression has been widely used to hypothesize gene function through guilt-by association. However, it is not clear to what degree co-expression is informative, whether it can be applied to genes involved in different biological processes, and how the type of dataset impacts inferences about gene functions. Here our goal is to assess the utility and limitations of using co-expression as a criterion to recover functional associations between genes. By determining the percentage of gene pairs in a metabolic pathway with significant expression correlation, we found that many genes in the same pathway do not have similar transcript profiles and the choice of dataset, annotation quality, gene function, expression similarity measure, and clustering approach significantly impacts the ability to recover functional associations between genes using Arabidopsis thaliana as an example. Some datasets are more informative in capturing coordinated expression profiles and larger data sets are not always better. In addition, to recover the maximum number of known pathways and identify candidate genes with similar functions, it is important to explore rather exhaustively multiple dataset combinations, similarity measures, clustering algorithms and parameters. Finally, we validated the biological relevance of co-expression cluster memberships with an independent phenomics dataset and found that genes that consistently cluster with leucine degradation genes tend to have similar leucine levels in mutants. This study provides a framework for obtaining gene functional associations by maximizing the information that can be obtained from gene expression datasets.

There remain genes with no known function even in the most well studied, model species. One common way to hypothesize gene function is based on the assumption that genes with similar expression profiles tend to have similar functions. However, using datasets and biological pathway information from the model plant Arabidopsis thaliana as an example, we discovered that, although genes in the same pathways are functionally related, genes in only a subset of the pathways have highly similar expression patterns. In addition, our ability to hypothesize gene functions based on expression is significantly impacted by how the dataset is processed and combined as well as the methodology used to identify genes with similar expression. Therefore, multiple datasets and methods should be tested to maximize the functional information that we can get based on similarity in gene expression.

Genome-Wide Analysis of Soybean JmjC Domain-Containing Proteins Suggests Evolutionary Conservation Following Whole-Genome Duplication

Histone modifications, such as methylation and demethylation, play an important role in regulating chromatin structure and gene expression. The JmjC domain-containing proteins, an important family of histone lysine demethylases (KDMs), play a key role in maintaining homeostasis of histone methylation in vivo. In this study, we performed a comprehensive analysis of the jumonji C (JmjC) gene family in the soybean genome and identified 48 JmjC genes (GmJMJs) distributed unevenly across 18 chromosomes. Phylogenetic analysis showed that these JmjC domain-containing genes can be divided into eight groups. GmJMJs within the same phylogenetic group share similar exon/intron organization and domain composition. In addition, 16 duplicated gene pairs were formed by a Glycine-specific whole-genome duplication (WGD) event approximately 13 million years ago (Mya). By investigating the expression profiles of these gene pairs in various tissues, we showed that the expression pattern is conserved in the polyploidy-derived JmjC duplicates, demonstrating that the majority of GmJMJs were preferentially retained after the most recent WGD event and suggesting important roles for demethylase duplications in soybean evolution. These results shed light on the evolutionary history of this family in soybean and provide insights into the JmjCs which will be helpful to reveal their functions in controlling soybean development.

ABA-dependent control of GIGANTEA signalling enables drought escape via up-regulation of FLOWERING LOCUS T in Arabidopsis thaliana

One strategy deployed by plants to endure water scarcity is to accelerate the transition to flowering adaptively via the drought escape (DE) response. In Arabidopsis thaliana, activation of the DE response requires the photoperiodic response gene GIGANTEA (GI) and the florigen genes FLOWERING LOCUS T (FT) and TWIN SISTER OF FT (TSF). The phytohormone abscisic acid (ABA) is also required for the DE response, by promoting the transcriptional up-regulation of the florigen genes. The mode of interaction between ABA and the photoperiodic genes remains obscure. In this work we use a genetic approach to demonstrate that ABA modulates GI signalling and consequently its ability to activate the florigen genes. We also reveal that the ABA-dependent activation of FT, but not TSF, requires CONSTANS (CO) and that impairing ABA signalling dramatically reduces the expression of florigen genes with little effect on the CO transcript profile. ABA signalling thus has an impact on the core genes of photoperiodic signalling GI and CO by modulating their downstream function and/or activities rather than their transcript accumulation. In addition, we show that as well as promoting flowering, ABA simultaneously represses flowering, independent of the florigen genes. Genetic analysis indicates that the target of the repressive function of ABA is the flowering-promoting gene SUPPRESSOR OF OVEREXPRESSION OF CONSTANS1 (SOC1), a transcription factor integrating floral cues in the shoot meristem. Our study suggests that variations in ABA signalling provide different developmental information that allows plants to co-ordinate the onset of the reproductive phase according to the available water resources.

Transcript, protein and metabolite temporal dynamics in the CAM plant Agave

Already a proven mechanism for drought resilience, crassulacean acid metabolism (CAM) is a specialized type of photosynthesis that maximizes water-use efficiency by means of an inverse (compared to C₃ and C₄ photosynthesis) day/night pattern of stomatal closure/opening to shift CO₂ uptake to the night, when evapotranspiration rates are low. A systems-level understanding of temporal molecular and metabolic controls is needed to define the cellular behaviour underpinning CAM. Here, we report high-resolution temporal behaviours of transcript, protein and metabolite abundances across a CAM diel cycle and, where applicable, compare the observations to the well-established C₃ model plant Arabidopsis. A mechanistic finding that emerged is that CAM operates with a diel redox poise that is shifted relative to that in Arabidopsis. Moreover, we identify widespread rescheduled expression of genes associated with signal transduction mechanisms that regulate stomatal opening/closing. Controlled production and degradation of transcripts and proteins represents a timing mechanism by which to regulate cellular function, yet knowledge of how this molecular timekeeping regulates CAM is unknown. Here, we provide new insights into complex post-transcriptional and -translational hierarchies that govern CAM in Agave. These data sets provide a resource to inform efforts to engineer more efficient CAM traits into economically valuable C₃ crops.

CGDB: a database of circadian genes in eukaryotes

We report a database of circadian genes in eukaryotes (CGDB, http://cgdb.biocuckoo.org), containing ∼73 000 circadian-related genes in 68 animals, 39 plants and 41 fungi. Circadian rhythm is ∼24 h rhythm in behavioral and physiological processes that exists in almost all organisms on the earth. Defects in the circadian system are highly associated with a number of diseases such as cancers. Although several databases have been established for rhythmically expressed genes, a comprehensive database of cycling genes across phyla is still lacking. From the literature, we collected 1382 genes of which transcript level oscillations were validated using methods such as RT-PCR, northern blot and in situ hybridization. Given that many genes exhibit different oscillatory patterns in different tissues/cells within an organism, we have included information regarding the phase and amplitude of the oscillation, as well as the tissue/cells in which the oscillation was identified. Using these well characterized cycling genes, we have then conducted an orthologous search and identified ∼45 000 potential cycling genes from 148 eukaryotes. Given that significant effort has been devoted to identifying cycling genes by transcriptome profiling, we have also incorporated these results, a total of over 26 000 genes, into our database.

Into the Evening: Complex Interactions in the Arabidopsis Circadian Clock

In Arabidopsis thaliana an assembly of proteins named the evening complex (EC) has been established as an essential component of the circadian clock with conserved functions in regulating plant growth and development. Recent studies identifying EC-regulated genes and EC-interacting proteins have expanded our understanding of EC function. In this review we focus on new progress uncovering how the EC contributes to the circadian network through the integration of environmental inputs and the direct regulation of key clock genes. We also summarize new findings of how the EC directly regulates clock outputs, such as photoperiodic and thermoresponsive growth, and provide new perspectives on future experiments to address unsolved questions related to the EC.

The 5'UTR Intron of Arabidopsis GGT1 Aminotransferase Enhances Promoter Activity by Recruiting RNA Polymerase II

Photorespiration is essential for the detoxification of glycolate and recycling of carbon to the Calvin Benson Bassham cycle. Enzymes participating in the pathway have been identified, and investigations now focus on the regulation of photorespiration by transporters and metabolites. However, regulation of photorespiration on the gene level has not been intensively studied. Here, we show that maximum transcript abundance of Glu:glyoxylate aminotransferase 1 (GGT1) is regulated by intron-mediated enhancement (IME) of the 5' leader intron rather than by regulatory elements in the 5' upstream region. The intron is rich in CT-stretches and contains the motif TGTGATTTG that is highly similar to the IME-related motif TTNGATYTG. The GGT1 intron also confers leaf-specific expression of foreign promoters. Quantitative PCR analysis and GUS activity measurements revealed that IME of the GGT1 5'UTR intron is controlled on the transcriptional level. IME by the GGT1 5'UTR intron was at least 2-fold. Chromatin immunoprecipitation experiments showed that the abundance of RNA polymerase II binding to the intron-less construct is reduced.