A proxitome-RNA-capture approach reveals that processing bodies repress coregulated hub genes.

Cellular condensates are usually ribonucleoprotein assemblies with liquid- or solid-like properties. Because these subcellular structures lack a delineating membrane, determining their compositions is difficult. Here we describe a proximity-biotinylation approach for capturing the RNAs of the condensates known as processing bodies (PBs) in Arabidopsis (Arabidopsis thaliana). By combining this approach with RNA detection, in silico, and high-resolution imaging approaches, we studied PBs under normal conditions and heat stress. PBs showed a much more dynamic RNA composition than the total transcriptome. RNAs involved in cell wall development and regeneration, plant hormonal signaling, secondary metabolism/defense, and RNA metabolism were enriched in PBs. RNA-binding proteins and the liquidity of PBs modulated RNA recruitment, while RNAs were frequently recruited together with their encoded proteins. In PBs, RNAs follow distinct fates: in small liquid-like PBs, RNAs get degraded while in more solid-like larger ones, they are stored. PB properties can be regulated by the actin-polymerizing SCAR (suppressor of the cyclic AMP)-WAVE (WASP family verprolin homologous) complex. SCAR/WAVE modulates the shuttling of RNAs between PBs and the translational machinery, thereby adjusting ethylene signaling. In summary, we provide an approach to identify RNAs in condensates that allowed us to reveal a mechanism for regulating RNA fate.

Binding by the Polycomb complex component BMI1 and H2A monoubiquitination shape local and long-range interactions in the Arabidopsis genome.

Three-dimensional (3D) chromatin organization is highly dynamic during development and seems to play a crucial role in regulating gene expression. Self-interacting domains, commonly called topologically associating domains (TADs) or compartment domains (CDs), have been proposed as the basic structural units of chromatin organization. Surprisingly, although these units have been found in several plant species, they escaped detection in Arabidopsis (Arabidopsis thaliana). Here, we show that the Arabidopsis genome is partitioned into contiguous CDs with different epigenetic features, which are required to maintain appropriate intra-CD and long-range interactions. Consistent with this notion, the histone-modifying Polycomb group machinery is involved in 3D chromatin organization. Yet, while it is clear that Polycomb repressive complex 2 (PRC2)-mediated trimethylation of histone H3 on lysine 27 (H3K27me3) helps establish local and long-range chromatin interactions in plants, the implications of PRC1-mediated histone H2A monoubiquitination on lysine 121 (H2AK121ub) are unclear. We found that PRC1, together with PRC2, maintains intra-CD interactions, but it also hinders the formation of H3K4me3-enriched local chromatin loops when acting independently of PRC2. Moreover, the loss of PRC1 or PRC2 activity differentially affects long-range chromatin interactions, and these 3D changes differentially affect gene expression. Our results suggest that H2AK121ub helps prevent the formation of transposable element/H3K27me1-rich long loops and serves as a docking point for H3K27me3 incorporation.

A gene regulatory network critical for axillary bud dormancy directly controlled by Arabidopsis BRANCHED1.

The Arabidopsis thaliana transcription factor BRANCHED1 (BRC1) plays a pivotal role in the control of shoot branching as it integrates environmental and endogenous signals that influence axillary bud growth. Despite its remarkable activity as a growth inhibitor, the mechanisms by which BRC1 promotes bud dormancy are largely unknown. We determined the genome-wide BRC1 binding sites in vivo and combined these with transcriptomic data and gene co-expression analyses to identify bona fide BRC1 direct targets. Next, we integrated multi-omics data to infer the BRC1 gene regulatory network (GRN) and used graph theory techniques to find network motifs that control the GRN dynamics. We generated an open online tool to interrogate this network. A group of BRC1 target genes encoding transcription factors (BTFs) orchestrate this intricate transcriptional network enriched in abscisic acid-related components. Promoter::ß-GLUCURONIDASE transgenic lines confirmed that BTFs are expressed in axillary buds. Transient co-expression assays and studies in planta using mutant lines validated the role of BTFs in modulating the GRN and promoting bud dormancy. This knowledge provides access to the developmental mechanisms that regulate shoot branching and helps identify candidate genes to use as tools to adapt plant architecture and crop production to ever-changing environmental conditions.

Transcriptomic characterization of the response to a microalgae extract in Arabidopsis thaliana and Solanum lycopersicum.

BACKGROUND: The steady world population growth and the current climate emergency crisis demand the development of sustainable methods to increase crop performance and resilience to the abiotic and biotic stresses produced by global warming. Microalgal extracts are being established as sustainable sources to produce compounds that improve agricultural yield, concurrently contributing during their production process to atmospheric CO2 abatement through the photosynthetic activity of microalgae. RESULTS: In the present study, we characterize the transcriptomic response in the model plant Arabidopsis thaliana and the plant of horticultural interest Solanum lycopersicum to the foliar application of a microalgae-based commercial preparation LRM™ (AlgaEnergy, Madrid, Spain). The foliar spray of LRM™ has a substantial effect over both transcriptomes potentially mediated by various compounds within LRM™, including its phytohormone content, activating systemic acquired resistance, possibly mediated by salicylic acid biosynthetic processes, and drought/heat acclimatization, induced by stomatal control and wax accumulation during cuticle development. Specifically, the agronomic improvements observed in treated S. lycopersicum (tomato) plants include an increase in the number of fruits, an acceleration in flowering time and the provision of higher drought resistance. The effect of LRM™ foliar spray in juvenile and adult plants was similar, producing a fast response detectable 2 h from its application that was also maintained 24 h later. CONCLUSION: The present study improves our knowledge on the transcriptomic effect of a novel microalgal extract on crops and provides the first step towards a full understanding of the yield and resistance improvement of crops. © 2024 The Authors. Journal of The Science of Food and Agriculture published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Role of a cryptic tRNA gene operon in survival under translational stress.

As compared to eukaryotes, bacteria have a reduced tRNA gene set encoding between 30 and 220 tRNAs. Although in most bacterial phyla tRNA genes are dispersed in the genome, many species from distinct phyla also show genes forming arrays. Here, we show that two types of arrays with distinct evolutionary origins exist. This work focuses on long tRNA gene arrays (L-arrays) that encompass up to 43 genes, which disseminate by horizontal gene transfer and contribute supernumerary tRNA genes to the host. Although in the few cases previously studied these arrays were reported to be poorly transcribed, here we show that the L-array of the model cyanobacterium Anabaena sp. PCC 7120, encoding 23 functional tRNAs, is largely induced upon impairment of the translation machinery. The cellular response to this challenge involves a global reprogramming of the transcriptome in two phases. tRNAs encoded in the array are induced in the second phase of the response, directly contributing to cell survival. Results presented here show that in some bacteria the tRNA gene set may be partitioned between a housekeeping subset, which constantly sustains translation, and an inducible subset that is generally silent but can provide functionality under particular conditions.

ALGAEFUN with MARACAS, microALGAE FUNctional enrichment tool for MicroAlgae RnA-seq and Chip-seq AnalysiS.

BACKGROUND: Microalgae are emerging as promising sustainable sources for biofuels, biostimulants in agriculture, soil bioremediation, feed and human nutrients. Nonetheless, the molecular mechanisms underpinning microalgae physiology and the biosynthesis of compounds of biotechnological interest are largely uncharacterized. This hinders the development of microalgae full potential as cell-factories. The recent application of omics technologies into microalgae research aims at unraveling these systems. Nevertheless, the lack of specific tools for analysing omics raw data generated from microalgae to provide biological meaningful information are hampering the impact of these technologies. The purpose of ALGAEFUN with MARACAS consists in providing researchers in microalgae with an enabling tool that will allow them to exploit transcriptomic and cistromic high-throughput sequencing data. RESULTS: ALGAEFUN with MARACAS consists of two different tools. First, MARACAS (MicroAlgae RnA-seq and Chip-seq AnalysiS) implements a fully automatic computational pipeline receiving as input RNA-seq (RNA sequencing) or ChIP-seq (chromatin immunoprecipitation sequencing) raw data from microalgae studies. MARACAS generates sets of differentially expressed genes or lists of genomic loci for RNA-seq and ChIP-seq analysis respectively. Second, ALGAEFUN (microALGAE FUNctional enrichment tool) is a web-based application where gene sets generated from RNA-seq analysis as well as lists of genomic loci from ChIP-seq analysis can be used as input. On the one hand, it can be used to perform Gene Ontology and biological pathways enrichment analysis over gene sets. On the other hand, using the results of ChIP-seq data analysis, it identifies a set of potential target genes and analyses the distribution of the loci over gene features. Graphical representation of the results as well as tables with gene annotations are generated and can be downloaded for further analysis. CONCLUSIONS: ALGAEFUN with MARACAS provides an integrated environment for the microalgae research community that facilitates the process of obtaining relevant biological information from raw RNA-seq and ChIP-seq data. These applications are designed to assist researchers in the interpretation of gene lists and genomic loci based on functional enrichment analysis. ALGAEFUN with MARACAS is publicly available on https://greennetwork.us.es/AlgaeFUN/ .

Comparative transcriptomic analysis unveils interactions between the regulatory CarS protein and light response in Fusarium.

BACKGROUND: The orange pigmentation of the agar cultures of many Fusarium species is due to the production of carotenoids, terpenoid pigments whose synthesis is stimulated by light. The genes of the carotenoid pathway and their regulation have been investigated in detail in Fusarium fujikuroi. In this and other Fusarium species, such as F. oxysporum, deep-pigmented mutants affected in the gene carS, which encodes a protein of the RING-finger family, overproduce carotenoids irrespective of light. The induction of carotenogenesis by light and its deregulation in carS mutants are achieved on the transcription of the structural genes of the pathway. We have carried out global RNA-seq transcriptomics analyses to investigate the relationship between the regulatory role of CarS and the control by light in these fungi. RESULTS: The absence of a functional carS gene or the illumination exert wide effects on the transcriptome of F. fujikuroi, with predominance of genes activated over repressed and a greater functional diversity in the case of genes induced by light. The number of the latter decreases drastically in a carS mutant (1.1% vs. 4.8% in the wild-type), indicating that the deregulation produced by the carS mutation affects the light response of many genes. Moreover, approximately 27% of the genes activated at least 2-fold by light or by the carS mutation are coincident, raising to 40% for an 8-fold activation threshold. As expected, the genes with the highest changes under both regulatory conditions include those involved in carotenoid metabolism. In addition, light and CarS strongly influence the expression of some genes associated with stress responses, including three genes with catalase domains, consistent with roles in the control of oxidative stress. The effects of the CarS mutation or light in the transcriptome of F. oxysporum were partially coincident with those of F. fujikuroi, indicating the conservation of the objectives of their regulatory mechanisms. CONCLUSIONS: The CarS RING finger protein down-regulates many genes whose expression is up-regulated by light in wild strains of the two investigated Fusarium species, indicating a regulatory interplay between the mechanism of action of the CarS protein and the control by light.

The Arabidopsis Polycomb Repressive Complex 1 (PRC1) Components AtBMI1A, B, and C Impact Gene Networks throughout All Stages of Plant Development.

Polycomb Group regulation in Arabidopsis (Arabidopsis thaliana) is required to maintain cell differentiation and allow developmental phase transitions. This is achieved by the activity of three PcG repressive complex 2s (PRC2s) and the participation of a yet poorly defined PRC1. Previous results showed that apparent PRC1 components perform discrete roles during plant development, suggesting the existence of PRC1 variants; however, it is not clear in how many processes these components participate. We show that AtBMI1 proteins are required to promote all developmental phase transitions and to control cell proliferation during organ growth and development, expanding their proposed range of action. While AtBMI1 function during germination is closely linked to B3 domain transcription factors VAL1/2 possibly in combination with GT-box binding factors, other AtBMI1 regulatory networks require participation of different factor combinations. Conversely, EMF1 and LHP1 bind many H3K27me3 positive genes up-regulated in atbmi1a/b/c mutants; however, loss of their function affects expression of a different subset, suggesting that even if EMF1, LHP1, and AtBMI1 exist in a common PRC1 variant, their role in repression depends on the functional context.

Photoperiodic control of carbon distribution during the floral transition in Arabidopsis.

Flowering is a crucial process that demands substantial resources. Carbon metabolism must be coordinated with development through a control mechanism that optimizes fitness for any physiological need and growth stage of the plant. However, how sugar allocation is controlled during the floral transition is unknown. Recently, the role of a CONSTANS (CO) ortholog (Cr-CO) in the control of the photoperiod response in the green alga Chlamydomonas reinhardtii and its influence on starch metabolism was demonstrated. In this work, we show that transitory starch accumulation and glycan composition during the floral transition in Arabidopsis thaliana are regulated by photoperiod. Employing a multidisciplinary approach, we demonstrate a role for CO in regulating the level and timing of expression of the GRANULE BOUND STARCH SYNTHASE (GBSS) gene. Furthermore, we provide a detailed characterization of a GBSS mutant involved in transitory starch synthesis and analyze its flowering time phenotype in relation to its altered capacity to synthesize amylose and to modify the plant free sugar content. Photoperiod modification of starch homeostasis by CO may be crucial for increasing the sugar mobilization demanded by the floral transition. This finding contributes to our understanding of the flowering process.

ChlamyNET: a Chlamydomonas gene co-expression network reveals global properties of the transcriptome and the early setup of key co-expression patterns in the green lineage.

BACKGROUND: Chlamydomonas reinhardtii is the model organism that serves as a reference for studies in algal genomics and physiology. It is of special interest in the study of the evolution of regulatory pathways from algae to higher plants. Additionally, it has recently gained attention as a potential source for bio-fuel and bio-hydrogen production. The genome of Chlamydomonas is available, facilitating the analysis of its transcriptome by RNA-seq data. This has produced a massive amount of data that remains fragmented making necessary the application of integrative approaches based on molecular systems biology. RESULTS: We constructed a gene co-expression network based on RNA-seq data and developed a web-based tool, ChlamyNET, for the exploration of the Chlamydomonas transcriptome. ChlamyNET exhibits a scale-free and small world topology. Applying clustering techniques, we identified nine gene clusters that capture the structure of the transcriptome under the analyzed conditions. One of the most central clusters was shown to be involved in carbon/nitrogen metabolism and signalling, whereas one of the most peripheral clusters was involved in DNA replication and cell cycle regulation. The transcription factors and regulators in the Chlamydomonas genome have been identified in ChlamyNET. The biological processes potentially regulated by them as well as their putative transcription factor binding sites were determined. The putative light regulated transcription factors and regulators in the Chlamydomonas genome were analyzed in order to provide a case study on the use of ChlamyNET. Finally, we used an independent data set to cross-validate the predictive power of ChlamyNET. CONCLUSIONS: The topological properties of ChlamyNET suggest that the Chlamydomonas transcriptome posseses important characteristics related to error tolerance, vulnerability and information propagation. The central part of ChlamyNET constitutes the core of the transcriptome where most authoritative hub genes are located interconnecting key biological processes such as light response with carbon and nitrogen metabolism. Our study reveals that key elements in the regulation of carbon and nitrogen metabolism, light response and cell cycle identified in higher plants were already established in Chlamydomonas. These conserved elements are not only limited to transcription factors, regulators and their targets, but also include the cis-regulatory elements recognized by them.

An Evolutionarily Conserved DOF-CONSTANS Module Controls Plant Photoperiodic Signaling.

The response to daylength is a crucial process that evolved very early in plant evolution, entitling the early green eukaryote to predict seasonal variability and attune its physiological responses to the environment. The photoperiod responses evolved into the complex signaling pathways that govern the angiosperm floral transition today. The Chlamydomonas reinhardtii DNA-Binding with One Finger (CrDOF) gene controls transcription in a photoperiod-dependent manner, and its misexpression influences algal growth and viability. In short days, CrDOF enhances CrCO expression, a homolog of plant CONSTANS (CO), by direct binding to its promoter, while it reduces the expression of cell division genes in long days independently of CrCO. In Arabidopsis (Arabidopsis thaliana), transgenic plants overexpressing CrDOF show floral delay and reduced expression of the photoperiodic genes CO and FLOWERING LOCUS T. The conservation of the DOF-CO module during plant evolution could be an important clue to understanding diversification by the inheritance of conserved gene toolkits in key developmental programs.

CONSTANS alters the circadian clock in Arabidopsis thaliana.

Plants are sessile organisms that have acquired highly plastic developmental strategies to adapt to the environment. Among these processes, the floral transition is essential to ensure reproductive success and is finely regulated by several internal and external genetic networks. The photoperiodic pathway, which controls plant response to day length, is one of the most important pathways controlling flowering. In Arabidopsis photoperiodic flowering, CONSTANS (CO) is the central gene activating the expression of the florigen FLOWERING LOCUS T (FT) in the leaves at the end of a long day. The circadian clock strongly regulates CO expression. However, to date, no evidence has been reported regarding a feedback loop from the photoperiod pathway back to the circadian clock. Using transcriptional networks, we have identified relevant network motifs regulating the interplay between the circadian clock and the photoperiod pathway. Gene expression, chromatin immunoprecipitation experiments, and phenotypic analysis allowed us to elucidate the role of CO over the circadian clock. Plants with altered CO expression showed a different internal clock period, measured by daily leaf rhythmic movements. We showed that CO upregulates the expression of key genes related to the circadian clock, such as CCA1, LHY, PRR5, and GI, at the end of a long day by binding to specific sites on their promoters. Moreover, a high number of PRR5-repressed target genes are upregulated by CO, and this could explain the phase transition promoted by CO. The CO-PRR5 complex interacts with the bZIP transcription factor HY5 and helps to localize the complex in the promoters of clock genes. Taken together, our results indicate that there may be a feedback loop in which CO communicates back to the circadian clock, providing seasonal information to the circadian system.

Transcriptomic and Metabolomic Response to High Light in the Charophyte Alga Klebsormidium nitens.

The characterization of the molecular mechanisms, such as high light irradiance resistance, that allowed plant terrestralization is a cornerstone in evolutionary studies since the conquest of land by plants played a pivotal role in life evolution on Earth. Viridiplantae or the green lineage is divided into two clades, Chlorophyta and Streptophyta, that in turn splits into Embryophyta or land plants and Charophyta. Charophyta are used in evolutionary studies on plant terrestralization since they are generally accepted as the extant algal species most closely related to current land plants. In this study, we have chosen the facultative terrestrial early charophyte alga Klebsormidium nitens to perform an integrative transcriptomic and metabolomic analysis under high light in order to unveil key mechanisms involved in the early steps of plants terrestralization. We found a fast chloroplast retrograde signaling possibly mediated by reactive oxygen species and the inositol polyphosphate 1-phosphatase (SAL1) and 3'-phosphoadenosine-5'-phosphate (PAP) pathways inducing gene expression and accumulation of specific metabolites. Systems used by both Chlorophyta and Embryophyta were activated such as the xanthophyll cycle with an accumulation of zeaxanthin and protein folding and repair mechanisms constituted by NADPH-dependent thioredoxin reductases, thioredoxin-disulfide reductases, and peroxiredoxins. Similarly, cyclic electron flow, specifically the pathway dependent on proton gradient regulation 5, was strongly activated under high light. We detected a simultaneous co-activation of the non-photochemical quenching mechanisms based on LHC-like stress related (LHCSR) protein and the photosystem II subunit S that are specific to Chlorophyta and Embryophyta, respectively. Exclusive Embryophyta systems for the synthesis, sensing, and response to the phytohormone auxin were also activated under high light in K. nitens leading to an increase in auxin content with the concomitant accumulation of amino acids such as tryptophan, histidine, and phenylalanine.

Unveiling the underlying molecular basis of astaxanthin accumulation in Haematococcus through integrative metabolomic-transcriptomic analysis.

Astaxanthin is a valuable and highly demanded ketocarotenoid pigment, for which the chlorophycean microalga Haematococcus pluvialis is an outstanding natural source. Although information on astaxanthin accumulation in H. pluvialis has substantially advanced in recent years, its underlying molecular bases remain elusive. An integrative metabolic and transcriptomic analysis has been performed for vegetative Haematococcus cells, grown both under N sufficiency (green palmelloid cells) and under moderate N limitation, allowing concurrent active cell growth and astaxanthin synthesis (reddish palmelloid cells). Transcriptional activation was noticeable in reddish cells of key enzymes participating in glycolysis, pentose phosphate cycle and pyruvate metabolism, determining the adequate provision of glyceraldehyde 3 phosphate and pyruvate, precursors of carotenoids and fatty acids. Moreover, for the first time, transcriptional regulators potentially involved in controlling astaxanthin accumulation have been identified, a knowledge enabling optimization of commercial astaxanthin production by Haematococcus through systems metabolic engineering.

Changes at a Critical Branchpoint in the Anthocyanin Biosynthetic Pathway Underlie the Blue to Orange Flower Color Transition in Lysimachia arvensis.

Anthocyanins are the primary pigments contributing to the variety of flower colors among angiosperms and are considered essential for survival and reproduction. Anthocyanins are members of the flavonoids, a broader class of secondary metabolites, of which there are numerous structural genes and regulators thereof. In western European populations of Lysimachia arvensis, there are blue- and orange-petaled individuals. The proportion of blue-flowered plants increases with temperature and daylength yet decreases with precipitation. Here, we performed a transcriptome analysis to characterize the coding sequences of a large group of flavonoid biosynthetic genes, examine their expression and compare our results to flavonoid biochemical analysis for blue and orange petals. Among a set of 140 structural and regulatory genes broadly representing the flavonoid biosynthetic pathway, we found 39 genes with significant differential expression including some that have previously been reported to be involved in similar flower color transitions. In particular, F3'5'H and DFR, two genes at a critical branchpoint in the ABP for determining flower color, showed differential expression. The expression results were complemented by careful examination of the SNPs that differentiate the two color types for these two critical genes. The decreased expression of F3'5'H in orange petals and differential expression of two distinct copies of DFR, which also exhibit amino acid changes in the color-determining substrate specificity region, strongly correlate with the blue to orange transition. Our biochemical analysis was consistent with the transcriptome data indicating that the shift from blue to orange petals is caused by a change from primarily malvidin to largely pelargonidin forms of anthocyanins. Overall, we have identified several flavonoid biosynthetic pathway loci likely involved in the shift in flower color in L. arvensis and even more loci that may represent the complex network of genetic and physiological consequences of this flower color polymorphism.

H2AK121ub in Arabidopsis associates with a less accessible chromatin state at transcriptional regulation hotspots.

Although it is well established that the Polycomb Group (PcG) complexes maintain gene repression through the incorporation of H2AK121ub and H3K27me3, little is known about the effect of these modifications on chromatin accessibility, which is fundamental to understand PcG function. Here, by integrating chromatin accessibility, histone marks and expression analyses in different Arabidopsis PcG mutants, we show that PcG function regulates chromatin accessibility. We find that H2AK121ub is associated with a less accessible but still permissive chromatin at transcriptional regulation hotspots. Accessibility is further reduced by EMF1 acting in collaboration with PRC2 activity. Consequently, H2AK121ub/H3K27me3 marks are linked to inaccessible although responsive chromatin. In contrast, only-H3K27me3-marked chromatin is less responsive, indicating that H2AK121ub-marked hotspots are required for transcriptional responses. Nevertheless, despite the loss of PcG activities leads to increased chromatin accessibility, this is not necessarily accompanied by transcriptional activation, indicating that accessible chromatin is not always predictive of gene expression.

H2A monoubiquitination in Arabidopsis thaliana is generally independent of LHP1 and PRC2 activity.

BACKGROUND: Polycomb group complexes PRC1 and PRC2 repress gene expression at the chromatin level in eukaryotes. The classic recruitment model of Polycomb group complexes in which PRC2-mediated H3K27 trimethylation recruits PRC1 for H2A monoubiquitination was recently challenged by data showing that PRC1 activity can also recruit PRC2. However, the prevalence of these two mechanisms is unknown, especially in plants as H2AK121ub marks were examined at only a handful of Polycomb group targets. RESULTS: By using genome-wide analyses, we show that H2AK121ub marks are surprisingly widespread in Arabidopsis thaliana, often co-localizing with H3K27me3 but also occupying a set of transcriptionally active genes devoid of H3K27me3. Furthermore, by profiling H2AK121ub and H3K27me3 marks in atbmi1a/b/c, clf/swn, and lhp1 mutants we found that PRC2 activity is not required for H2AK121ub marking at most genes. In contrast, loss of AtBMI1 function impacts the incorporation of H3K27me3 marks at most Polycomb group targets. CONCLUSIONS: Our findings show the relationship between H2AK121ub and H3K27me3 marks across the A. thaliana genome and unveil that ubiquitination by PRC1 is largely independent of PRC2 activity in plants, while the inverse is true for H3K27 trimethylation.

Evolution of photoperiod sensing in plants and algae.

Measuring day length confers a strong fitness improvement to photosynthetic organisms as it allows them to anticipate light phases and take the best decisions preceding diurnal transitions. In close association with signals from the circadian clock and the photoreceptors, photoperiodic sensing constitutes also a precise way to determine the passing of the seasons and to take annual decisions such as the best time to flower or the beginning of dormancy. Photoperiodic sensing in photosynthetic organisms is ancient and two major stages in its evolution could be identified, the cyanobacterial time sensing and the evolutionary tool kit that arose in green algae and developed into the photoperiodic system of modern plants. The most recent discoveries about the evolution of the perception of light, measurement of day length and relationship with the circadian clock along the evolution of the eukaryotic green lineage will be discussed in this review.

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.

Evolutionary Analysis of DELLA-Associated Transcriptional Networks.

DELLA proteins are transcriptional regulators present in all land plants which have been shown to modulate the activity of over 100 transcription factors in Arabidopsis, involved in multiple physiological and developmental processes. It has been proposed that DELLAs transduce environmental information to pre-wired transcriptional circuits because their stability is regulated by gibberellins (GAs), whose homeostasis largely depends on environmental signals. The ability of GAs to promote DELLA degradation coincides with the origin of vascular plants, but the presence of DELLAs in other land plants poses at least two questions: what regulatory properties have DELLAs provided to the behavior of transcriptional networks in land plants, and how has the recruitment of DELLAs by GA signaling affected this regulation. To address these issues, we have constructed gene co-expression networks of four different organisms within the green lineage with different properties regarding DELLAs: Arabidopsis thaliana and Solanum lycopersicum (both with GA-regulated DELLA proteins), Physcomitrella patens (with GA-independent DELLA proteins) and Chlamydomonas reinhardtii (a green alga without DELLA), and we have examined the relative evolution of the subnetworks containing the potential DELLA-dependent transcriptomes. Network analysis indicates a relative increase in parameters associated with the degree of interconnectivity in the DELLA-associated subnetworks of land plants, with a stronger effect in species with GA-regulated DELLA proteins. These results suggest that DELLAs may have played a role in the coordination of multiple transcriptional programs along evolution, and the function of DELLAs as regulatory 'hubs' became further consolidated after their recruitment by GA signaling in higher plants.