Prevalence, evolution, and cis-regulation of diel transcription in Chlamydomonas reinhardtii

Dominant constraints on the evolution of rhythmic gene expression

Although the individual transcriptional regulators of the core circadian clock are distinct among different organisms, the autoregulatory feedback loops they form are conserved. This unified design principle explains how daily physiological activities oscillate across species. However, it is unknown whether analogous design principles govern the gene expression output of circadian clocks. In this study, we performed a comparative analysis of rhythmic gene expression in eight diverse species and identified four common distribution patterns of cycling gene expression across these species. We hypothesized that the maintenance of reduced energetic costs constrains the evolution of rhythmic gene expression. Our large-scale computational simulations support this hypothesis by showing that selection against high-energy expenditure completely regenerates all cycling gene patterns. Moreover, we find that the peaks of rhythmic expression have been subjected to this type of selective pressure. The results suggest that selective pressure from circadian regulation efficiently removes unnecessary gene products from the transcriptome, thereby significantly impacting its evolutionary path.

Cis-acting elements involved in the G2/M-phase-specific transcription of the cyclin B gene in the unicellular alga Cyanidioschyzon merolae

M-specific activator (MSA) cis-acting elements have been determined to be involved in the regulation of G2/M-phase-specific transcription in spermatophytes. In this study, the involvement of MSA-core elements in G2/M-phase-specific transcription was examined in the unicellular red alga Cyanidioschyzon merolae. In the C. merolae genome, MSA-core elements do not accumulate specifically in the upstream of mitosis-specific transcriptional genes. Mutations of the four MSA-core elements of the cyclin B gene, which encodes a central factor of the G2-to-M-phase transition, have resulted in the abolishment of transcription or permission of transcription even in the G1 phase. These results suggest that all four MSA-core elements located in the upstream region of cyclin B are involved in G2/M-phase-specific transcription in C. merolae; however, the nature of the involvement of MSA-core elements in G2/M-phase-specific transcription differed among the four elements. Thus, MSA-core-element-mediated G2/M-phase-specific transcription in C. merolae seems to be regulated by a complex mechanism.

Co-expression networks in Chlamydomonas reveal significant rhythmicity in batch cultures and empower gene function discovery

The unicellular green alga Chlamydomonas reinhardtii is a choice reference system for the study of photosynthesis and chloroplast metabolism, cilium assembly and function, lipid and starch metabolism, and metal homeostasis. Despite decades of research, the functions of thousands of genes remain largely unknown, and new approaches are needed to categorically assign genes to cellular pathways. Growing collections of transcriptome and proteome data now allow a systematic approach based on integrative co-expression analysis. We used a dataset comprising 518 deep transcriptome samples derived from 58 independent experiments to identify potential co-expression relationships between genes. We visualized co-expression potential with the R package corrplot, to easily assess co-expression and anti-correlation between genes. We extracted several hundred high-confidence genes at the intersection of multiple curated lists involved in cilia, cell division, and photosynthesis, illustrating the power of our method. Surprisingly, Chlamydomonas experiments retained a significant rhythmic component across the transcriptome, suggesting an underappreciated variable during sample collection, even in samples collected in constant light. Our results therefore document substantial residual synchronization in batch cultures, contrary to assumptions of asynchrony. We provide step-by-step protocols for the analysis of co-expression across transcriptome data sets from Chlamydomonas and other species to help foster gene function discovery.

Regulation of Multiple Fission and Cell-Cycle-Dependent Gene Expression by CDKA1 and the Rb-E2F Pathway in Chlamydomonas

The green alga Chlamydomonas proliferates by "multiple fission": a long G1 with >10-fold cell growth followed by multiple rapid divisions. Cells above a critical size threshold are "committed" and will divide independent of light and further cell growth. The number of divisions carried out depends on the initial size of the committed mother cell. Here, I show that CDKA1, the ortholog of the yeast and animal mitotic inducer CDK1, regulates the critical size for commitment. The Rb/E2F/Dp1 pathway regulates division number as well as commitment size. Epistasis analysis indicated that CDKA1 and Rb/E2F/Dp1 regulate multiple fission by distinct mechanisms. Rb-E2F/Dp1 regulates G1/S gene expression in animals and land plants. Transcriptome analysis showed that mat3 or dp1 disruption altered regulation of a large group of cell-division-associated genes with respect to cell size, but not with respect to synchronization timing. In contrast, cdka1 inactivation disturbed both temporal and cell-size regulation of expression. These defects were enhanced by double inactivation of cdka1 and dp1, suggesting interaction between CDKA1 and the Rb-E2F/Dp1 pathways in regulating cell-cycle-specific gene expression and cell-cycle initiation. In the context of a theoretical model for regulation of Chlamydomonas multiple fission, these results suggest that CDKA1 may promote a switch into a division-competent state, and E2F/Dp1 may promote maintenance of this state.

A Pan-plant Protein Complex Map Reveals Deep Conservation and Novel Assemblies

Plants are foundational for global ecological and economic systems, but most plant proteins remain uncharacterized. Protein interaction networks often suggest protein functions and open new avenues to characterize genes and proteins. We therefore systematically determined protein complexes from 13 plant species of scientific and agricultural importance, greatly expanding the known repertoire of stable protein complexes in plants. By using co-fractionation mass spectrometry, we recovered known complexes, confirmed complexes predicted to occur in plants, and identified previously unknown interactions conserved over 1.1 billion years of green plant evolution. Several novel complexes are involved in vernalization and pathogen defense, traits critical for agriculture. We also observed plant analogs of animal complexes with distinct molecular assemblies, including a megadalton-scale tRNA multi-synthetase complex. The resulting map offers a cross-species view of conserved, stable protein assemblies shared across plant cells and provides a mechanistic, biochemical framework for interpreting plant genetics and mutant phenotypes.

Chlamydomonas CHT7 Is Required for an Effective Quiescent State by Regulating Nutrient-Responsive Cell Cycle Gene Expression

COMPROMISED HYDROLYSIS OF TRIACYLGLYCEROLS7 (CHT7) in Chlamydomonas (Chlamydomonas reinhardtii) was previously shown to affect the transcription of a subset of genes during nitrogen (N)-replete growth and following N refeeding. Here, we show that an extensive derepression of genes involved in DNA metabolism and cell cycle-related processes, as well as downregulation of genes encoding oxidoreductases and nutrient transporters, occurs in the cht7 mutant during N deprivation. Cellular mutant phenotypes are consistent with the observed transcriptome misregulation, as cht7 cells fail to properly arrest growth, nuclear replication, and cell division following N deprivation. Reduction in cht7 colony formation following N refeeding is explained by its compromised viability during N deprivation and by the occurrence of abortive divisions during N refeeding. Surprisingly, the largely unstructured C-terminal half of CHT7 with predicted protein binding domains, but not the canonical CXC DNA binding domain, is essential for the ability of CHT7 to form stable complexes and reverse the cellular phenotypes and transcription levels in the cht7 mutant. Hence, although lacking the presumed DNA binding domain, CHT7 modulates the expression of cell cycle genes in response to N availability, which is essential for establishing an effective quiescent state and the coordinated resumption of growth following N refeeding.

Improved recovery of cell-cycle gene expression in Saccharomyces cerevisiae from regulatory interactions in multiple omics data

BACKGROUND

Gene expression is regulated by DNA-binding transcription factors (TFs). Together with their target genes, these factors and their interactions collectively form a gene regulatory network (GRN), which is responsible for producing patterns of transcription, including cyclical processes such as genome replication and cell division. However, identifying how this network regulates the timing of these patterns, including important interactions and regulatory motifs, remains a challenging task.

RESULTS

We employed four in vivo and in vitro regulatory data sets to investigate the regulatory basis of expression timing and phase-specific patterns cell-cycle expression in Saccharomyces cerevisiae. Specifically, we considered interactions based on direct binding between TF and target gene, indirect effects of TF deletion on gene expression, and computational inference. We found that the source of regulatory information significantly impacts the accuracy and completeness of recovering known cell-cycle expressed genes. The best approach involved combining TF-target and TF-TF interactions features from multiple datasets in a single model. In addition, TFs important to multiple phases of cell-cycle expression also have the greatest impact on individual phases. Important TFs regulating a cell-cycle phase also tend to form modules in the GRN, including two sub-modules composed entirely of unannotated cell-cycle regulators (STE12-TEC1 and RAP1-HAP1-MSN4).

CONCLUSION

Our findings illustrate the importance of integrating both multiple omics data and regulatory motifs in order to understand the significance regulatory interactions involved in timing gene expression. This integrated approached allowed us to recover both known cell-cycles interactions and the overall pattern of phase-specific expression across the cell-cycle better than any single data set. Likewise, by looking at regulatory motifs in the form of TF-TF interactions, we identified sets of TFs whose co-regulation of target genes was important for cell-cycle expression, even when regulation by individual TFs was not. Overall, this demonstrates the power of integrating multiple data sets and models of interaction in order to understand the regulatory basis of established biological processes and their associated gene regulatory networks.

Molecular Traits of Long Non-protein Coding RNAs from Diverse Plant Species Show Little Evidence of Phylogenetic Relationships

Long non-coding RNAs (lncRNAs) represent a diverse class of regulatory loci with roles in development and stress responses throughout all kingdoms of life. LncRNAs, however, remain under-studied in plants compared to animal systems. To address this deficiency, we applied a machine learning prediction tool, Classifying RNA by Ensemble Machine learning Algorithm (CREMA), to analyze RNAseq data from 11 plant species chosen to represent a wide range of evolutionary histories. Transcript sequences of all expressed and/or annotated loci from plants grown in unstressed (control) conditions were assembled and input into CREMA for comparative analyses. On average, 6.4% of the plant transcripts were identified by CREMA as encoding lncRNAs. Gene annotation associated with the transcripts showed that up to 99% of all predicted lncRNAs for Solanum tuberosum and Amborella trichopoda were missing from their reference annotations whereas the reference annotation for the genetic model plant Arabidopsis thaliana contains 96% of all predicted lncRNAs for this species. Thus a reliance on reference annotations for use in lncRNA research in less well-studied plants can be impeded by the near absence of annotations associated with these regulatory transcripts. Moreover, our work using phylogenetic signal analyses suggests that molecular traits of plant lncRNAs display different evolutionary patterns than all other transcripts in plants and have molecular traits that do not follow a classic evolutionary pattern. Specifically, GC content was the only tested trait of lncRNAs with consistently significant and high phylogenetic signal, contrary to high signal in all tested molecular traits for the other transcripts in our tested plant species.

Kingdom-wide comparison reveals the evolution of diurnal gene expression in Archaeplastida

Plants have adapted to the diurnal light-dark cycle by establishing elaborate transcriptional programs that coordinate many metabolic, physiological, and developmental responses to the external environment. These transcriptional programs have been studied in only a few species, and their function and conservation across algae and plants is currently unknown. We performed a comparative transcriptome analysis of the diurnal cycle of nine members of Archaeplastida, and we observed that, despite large phylogenetic distances and dramatic differences in morphology and lifestyle, diurnal transcriptional programs of these organisms are similar. Expression of genes related to cell division and the majority of biological pathways depends on the time of day in unicellular algae but we did not observe such patterns at the tissue level in multicellular land plants. Hence, our study provides evidence for the universality of diurnal gene expression and elucidates its evolutionary history among different photosynthetic eukaryotes.

Nitrogen-dependent coordination of cell cycle, quiescence and TAG accumulation in Chlamydomonas

Microalgae hold great promises as sustainable cellular factories for the production of alternative fuels, feeds, and biopharmaceuticals for human health. While the biorefinery approach for fuels along with the coproduction of high-value compounds with industrial, therapeutic, or nutraceutical applications have the potential to make algal biofuels more economically viable, a number of challenges continue to hamper algal production systems at all levels. One such hurdle includes the metabolic trade-off often observed between the increased yields of desired products, such as triacylglycerols (TAG), and the growth of an organism. Initial genetic engineering strategies to improve lipid productivity in microalgae, which focused on overproducing the enzymes involved in fatty acid and TAG biosynthesis or inactivating competing carbon (C) metabolism, have seen some successes albeit at the cost of often greatly reduced biomass. Emergent approaches that aim at modifying the dynamics of entire metabolic pathways by engineering of pertinent transcription factors or signaling networks appear to have successfully achieved a balance between growth and neutral lipid accumulation. However, the biological knowledge of key signaling networks and molecular components linking these two processes is still incomplete in photosynthetic eukaryotes, making it difficult to optimize metabolic engineering strategies for microalgae. Here, we focus on nitrogen (N) starvation of the model green microalga, Chlamydomonas reinhardtii, to present the current understanding of the nutrient-dependent switch between proliferation and quiescence, and the drastic reprogramming of metabolism that results in the storage of C compounds following N starvation. We discuss the potential components mediating the transcriptional repression of cell cycle genes and the establishment of quiescence in Chlamydomonas, and highlight the importance of signaling pathways such as those governed by the target of rapamycin (TOR) and sucrose nonfermenting-related (SnRK) kinases in the coordination of metabolic status with cellular growth. A better understanding of how the cell division cycle is regulated in response to nutrient scarcity and of the signaling pathways linking cellular growth to energy and lipid homeostasis, is essential to improve the prospects of biofuels and biomass production in microalgae.

Sizing up the cell cycle: systems and quantitative approaches in Chlamydomonas

The unicellular green alga Chlamydomonas provides a simplified model for defining core cell cycle functions conserved in the green lineage and for understanding multiple fission, a common cell cycle variation found in many algae. Systems-level approaches including a recent groundbreaking screen for conditional lethal cell cycle mutants and genome-wide transcriptome analyses are revealing the complex relationships among cell cycle regulators and helping define roles for CDKA/CDK1 and CDKB, the latter of which is unique to the green lineage and plays a central role in mitotic regulation. Genetic screens and quantitative single-cell analyses have provided insight into cell-size control during multiple fission including the identification of a candidate `sizer' protein. Quantitative single-cell tracking and modeling are promising approaches for gaining additional insight into regulation of cellular and subcellular scaling during the Chlamydomonas cell cycle.

The Intracellular Dynamics of Circadian Clocks Reach for the Light of Ecology and Evolution

A major challenge for biology is to extend our understanding of molecular regulation from the simplified conditions of the laboratory to ecologically relevant environments. Tractable examples are essential to make these connections for complex, pleiotropic regulators and, to go further, to link relevant genome sequences to field traits. Here, I review the case for the biological clock in higher plants. The gene network of the circadian clock drives pervasive, 24-hour rhythms in metabolism, behavior, and physiology across the eukaryotes and in some prokaryotes. In plants, the scope of chronobiology is now extending from the most tractable, intracellular readouts to the clock's many effects at the whole-organism level and across the life cycle, including biomass and flowering. I discuss five research areas where recent progress might be integrated in the future, to understand not only circadian functions in natural conditions but also the evolution of the clock's molecular mechanisms.

A G-Box-Like Motif Is Necessary for Transcriptional Regulation by Circadian Pseudo-Response Regulators in Arabidopsis

PSEUDO-RESPONSE REGULATORs (PRRs) play overlapping and distinct roles in maintaining circadian rhythms and regulating diverse biological processes, including the photoperiodic control of flowering, growth, and abiotic stress responses. PRRs act as transcriptional repressors and associate with chromatin via their conserved C-terminal CCT (CONSTANS, CONSTANS-like, and TIMING OF CAB EXPRESSION 1 [TOC1/PRR1]) domains by a still-poorly understood mechanism. Here, we identified genome-wide targets of PRR9 using chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq) and compared them with PRR7, PRR5, and TOC1/PRR1 ChIP-seq data. We found that PRR binding sites are located within genomic regions of low nucleosome occupancy and high DNase I hypersensitivity. Moreover, conserved noncoding regions among Brassicaceae species are enriched around PRR binding sites, indicating that PRRs associate with functionally relevant cis-regulatory regions. The PRRs shared a significant number of binding regions, and our results indicate that they coordinately restrict the expression of target genes to around dawn. A G-box-like motif was overrepresented at PRR binding regions, and we showed that this motif is necessary for mediating transcriptional regulation of CIRCADIAN CLOCK ASSOCIATED 1 and PRR9 by the PRRs. Our results further our understanding of how PRRs target specific promoters and provide an extensive resource for studying circadian regulatory networks in plants.

Dynamic Changes in the Transcriptome and Methylome of Chlamydomonas reinhardtii throughout Its Life Cycle

The green alga Chlamydomonas reinhardtii undergoes gametogenesis and mating upon nitrogen starvation. While the steps involved in its sexual reproductive cycle have been extensively characterized, the genome-wide transcriptional and epigenetic changes underlying different life cycle stages have yet to be fully described. Here, we performed transcriptome and methylome sequencing to quantify expression and DNA methylation from vegetative and gametic cells of each mating type and from zygotes. We identified 361 gametic genes with mating type-specific expression patterns and 627 genes that are specifically induced in zygotes; furthermore, these sex-related gene sets were enriched for secretory pathway and alga-specific genes. We also examined the C. reinhardtii nuclear methylation map with base-level resolution at different life cycle stages. Despite having low global levels of nuclear methylation, we detected 23 hypermethylated loci in gene-poor, repeat-rich regions. We observed mating type-specific differences in chloroplast DNA methylation levels in plus versus minus mating type gametes followed by chloroplast DNA hypermethylation in zygotes. Lastly, we examined the expression of candidate DNA methyltransferases and found three, DMT1a, DMT1b, and DMT4, that are differentially expressed during the life cycle and are candidate DNA methylases. The expression and methylation data we present provide insight into cell type-specific transcriptional and epigenetic programs during key stages of the C. reinhardtii life cycle.

High-Resolution Profiling of a Synchronized Diurnal Transcriptome from Chlamydomonas reinhardtii Reveals Continuous Cell and Metabolic Differentiation

The green alga Chlamydomonas reinhardtii is a useful model organism for investigating diverse biological processes, such as photosynthesis and chloroplast biogenesis, flagella and basal body structure/function, cell growth and division, and many others. We combined a highly synchronous photobioreactor culture system with frequent temporal sampling to characterize genome-wide diurnal gene expression in Chlamydomonas. Over 80% of the measured transcriptome was expressed with strong periodicity, forming 18 major clusters. Genes associated with complex structures and processes, including cell cycle control, flagella and basal bodies, ribosome biogenesis, and energy metabolism, all had distinct signatures of coexpression with strong predictive value for assigning and temporally ordering function. Importantly, the frequent sampling regime allowed us to discern meaningful fine-scale phase differences between and within subgroups of genes and enabled the identification of a transiently expressed cluster of light stress genes. Coexpression was further used both as a data-mining tool to classify and/or validate genes from other data sets related to the cell cycle and to flagella and basal bodies and to assign isoforms of duplicated enzymes to their cognate pathways of central carbon metabolism. Our diurnal coexpression data capture functional relationships established by dozens of prior studies and are a valuable new resource for investigating a variety of biological processes in Chlamydomonas and other eukaryotes.

Cyclin-Dependent Kinase Regulation of Diurnal Transcription in Chlamydomonas

We analyzed global transcriptome changes during synchronized cell division in the green alga Chlamydomonas reinhardtii. The Chlamydomonas cell cycle consists of a long G1 phase, followed by an S/M phase with multiple rapid, alternating rounds of DNA replication and segregation. We found that the S/M period is associated with strong induction of ∼2300 genes, many with conserved roles in DNA replication or cell division. Other genes, including many involved in photosynthesis, are reciprocally downregulated in S/M, suggesting a gene expression split correlating with the temporal separation between G1 and S/M. The Chlamydomonas cell cycle is synchronized by light-dark cycles, so in principle, these transcriptional changes could be directly responsive to light or to metabolic cues. Alternatively, cell-cycle-periodic transcription may be directly regulated by cyclin-dependent kinases. To distinguish between these possibilities, we analyzed transcriptional profiles of mutants in the kinases CDKA and CDKB, as well as other mutants with distinct cell cycle blocks. Initial cell-cycle-periodic expression changes are largely CDK independent, but later regulation (induction and repression) is under differential control by CDKA and CDKB. Deviation from the wild-type transcriptional program in diverse cell cycle mutants will be an informative phenotype for further characterization of the Chlamydomonas cell cycle.

Transcriptional coordination of physiological responses in Nannochloropsis oceanica CCMP1779 under light/dark cycles

Nannochloropsis oceanica CCMP1779 is a marine unicellular stramenopile and an emerging reference species for basic research on oleogenic microalgae with biotechnological relevance. We investigated its physiology and transcriptome under light/dark cycles. We observed oscillations in lipid content and a predominance of cell division in the first half of the dark phase. Globally, more than 60% of the genes cycled in N. oceanica CCMP1779, with gene expression peaking at different times of the day. Interestingly, the phase of expression of genes involved in certain biological processes was conserved across photosynthetic lineages. Furthermore, in agreement with our physiological studies we found the processes of lipid metabolism and cell division enriched in cycling genes. For example, there was tight coordination of genes involved in the lower part of glycolysis, fatty acid synthesis and lipid production at dawn preceding lipid accumulation during the day. Our results suggest that diel lipid storage plays a key role for N. oceanica CCMP1779 growth under natural conditions making this alga a promising model to gain a basic mechanistic understanding of triacylglycerol production in photosynthetic cells. Our data will help the formulation of new hypotheses on the role of cyclic gene expression in cell growth and metabolism in Nannochloropsis.