Eukaryotic Sexual Reproduction Evoked "with a Little Help from My Friends".

Bacteria and eukaryotes interact in many ways-from the microbiome that educates the mammalian immune system and enhances nutrition to relationships that are commensal, symbiotic, or parasitic. Now in an unexpected twist, King and colleagues have expanded the repertoire of prokaryotic influence over eukaryotic physiology to include mating.

Green Algal Models for Multicellularity.

The repeated evolution of multicellularity across the tree of life has profoundly affected the ecology and evolution of nearly all life on Earth. Many of these origins were in different groups of photosynthetic eukaryotes, or algae. Here, we review the evolution and genetics of multicellularity in several groups of green algae, which include the closest relatives of land plants. These include millimeter-scale, motile spheroids of up to 50,000 cells in the volvocine algae; decimeter-scale seaweeds in the genus Ulva (sea lettuce); and very plantlike, meter-scale freshwater algae in the genus Chara (stoneworts). We also describe algae in the genus Caulerpa, which are giant, multinucleate, morphologically complex single cells. In each case, we review the life cycle, phylogeny, and genetics of traits relevant to the evolution of multicellularity, and genetic and genomic resources available for the group in question. Finally, we suggest routes toward developing these groups as model organisms for the evolution of multicellularity.

Subscaling of a cytosolic RNA binding protein governs cell size homeostasis in the multiple fission alga Chlamydomonas.

Coordination of growth and division in eukaryotic cells is essential for populations of proliferating cells to maintain size homeostasis, but the underlying mechanisms that govern cell size have only been investigated in a few taxa. The green alga Chlamydomonas reinhardtii (Chlamydomonas) proliferates using a multiple fission cell cycle that involves a long G1 phase followed by a rapid series of successive S and M phases (S/M) that produces 2n daughter cells. Two control points show cell-size dependence: the Commitment control point in mid-G1 phase requires the attainment of a minimum size to enable at least one mitotic division during S/M, and the S/M control point where mother cell size governs cell division number (n), ensuring that daughter distributions are uniform. tny1 mutants pass Commitment at a smaller size than wild type and undergo extra divisions during S/M phase to produce small daughters, indicating that TNY1 functions to inhibit size-dependent cell cycle progression. TNY1 encodes a cytosolic hnRNP A-related RNA binding protein and is produced once per cell cycle during S/M phase where it is apportioned to daughter cells, and then remains at constant absolute abundance as cells grow, a property known as subscaling. Altering the dosage of TNY1 in heterozygous diploids or through mis-expression increased Commitment cell size and daughter cell size, indicating that TNY1 is a limiting factor for both size control points. Epistasis placed TNY1 function upstream of the retinoblastoma tumor suppressor complex (RBC) and one of its regulators, Cyclin-Dependent Kinase G1 (CDKG1). Moreover, CDKG1 protein and mRNA were found to over-accumulate in tny1 cells suggesting that CDKG1 may be a direct target of repression by TNY1. Our data expand the potential roles of subscaling proteins outside the nucleus and imply a control mechanism that ties TNY1 accumulation to pre-division mother cell size.

A conserved RWP-RK transcription factor VSR1 controls gametic differentiation in volvocine algae.

Volvocine green algae are a model for understanding the evolution of mating types and sexes. They are facultatively sexual, with gametic differentiation occurring in response to nitrogen starvation (-N) in most genera and to sex inducer hormone in Volvox. The conserved RWP-RK family transcription factor (TF) MID is encoded by the minus mating-type locus or male sex-determining region of heterothallic volvocine species and dominantly determines minus or male gametic differentiation. However, the factor(s) responsible for establishing the default plus or female differentiation programs have remained elusive. We performed a phylo-transcriptomic screen for autosomal RWP-RK TFs induced during gametogenesis in unicellular isogamous Chlamydomonas reinhardtii (Chlamydomonas) and in multicellular oogamous Volvox carteri (Volvox) and identified a single conserved ortho-group we named Volvocine Sex Regulator 1 (VSR1). Chlamydomonas vsr1 mutants of either mating type failed to mate and could not induce expression of key mating-type-specific genes. Similarly, Volvox vsr1 mutants in either sex could initiate sexual embryogenesis, but the presumptive eggs or androgonidia (sperm packet precursors) were infertile and unable to express key sex-specific genes. Yeast two-hybrid assays identified a conserved domain in VSR1 capable of self-interaction or interaction with the conserved N terminal domain of MID. In vivo coimmunoprecipitation experiments demonstrated association of VSR1 and MID in both Chlamydomonas and Volvox. These data support a new model for volvocine sexual differentiation where VSR1 homodimers activate expression of plus/female gamete-specific-genes, but when MID is present, MID-VSR1 heterodimers are preferentially formed and activate minus/male gamete-specific-genes.

Algal Sex Determination and the Evolution of Anisogamy.

Algae are photosynthetic eukaryotes whose taxonomic breadth covers a range of life histories, degrees of cellular and developmental complexity, and diverse patterns of sexual reproduction. These patterns include haploid- and diploid-phase sex determination, isogamous mating systems, and dimorphic sexes. Despite the ubiquity of sexual reproduction in algae, their mating-type-determination and sex-determination mechanisms have been investigated in only a limited number of representatives. These include volvocine green algae, where sexual cycles and sex-determining mechanisms have shed light on the transition from mating types to sexes, and brown algae, which are a model for UV sex chromosome evolution in the context of a complex haplodiplontic life cycle. Recent advances in genomics have aided progress in understanding sexual cycles in less-studied taxa including ulvophyte, charophyte, and prasinophyte green algae, as well as in diatoms.

Three genomes in the algal genus Volvox reveal the fate of a haploid sex-determining region after a transition to homothallism.

Transitions between separate sexes (dioecy) and other mating systems are common across eukaryotes. Here, we study a change in a haploid dioecious green algal species with male- and female-determining chromosomes (U and V). The genus Volvox is an oogamous (with large, immotile female gametes and small, motile male gametes) and includes both heterothallic species (with distinct male and female genotypes, associated with a mating-type system that prevents fusion of gametes of the same sex) and homothallic species (bisexual, with the ability to self-fertilize). We date the origin of an expanded sex-determining region (SDR) in Volvox to at least 75 Mya, suggesting that homothallism represents a breakdown of dioecy (heterothallism). We investigated the involvement of the SDR of the U and V chromosomes in this transition. Using de novo whole-genome sequences, we identified a heteromorphic SDR of ca 1 Mbp in male and female genotypes of the heterothallic species Volvox reticuliferus and a homologous region (SDLR) in the closely related homothallic species Volvox africanus, which retained several different hallmark features of an SDR. The V. africanus SDLR includes a large region resembling the female SDR of the presumptive heterothallic ancestor, whereas most genes from the male SDR are absent. However, we found a multicopy array of the male-determining gene, MID, in a different genomic location from the SDLR. Thus, in V. africanus, an ancestrally female genotype may have acquired MID and thereby gained male traits.

The landscape of Chlamydomonas histone H3 lysine 4 methylation reveals both constant features and dynamic changes during the diurnal cycle.

Chromatin modifications are epigenetic regulatory features with major roles in various cellular events, yet they remain understudied in algae. We interrogated the genome-wide distribution pattern of mono- and trimethylated histone H3 lysine 4 (H3K4) using chromatin-immunoprecipitation followed by deep-sequencing (ChIP-seq) during key phases of the Chlamydomonas cell cycle: early G1 phase, Zeitgeber Time 1 (ZT1), when cells initiate biomass accumulation, S/M phase (ZT13) when cells are replicating DNA and undergoing mitosis, and late G0 phase (ZT23) when they are quiescent. Tri-methylated H3K4 was predominantly enriched at transcription start sites of the majority of protein coding genes (85%). The likelihood of a gene being marked by H3K4me3 correlated with it being transcribed at some point during the life cycle but not necessarily by continuous active transcription, as exemplified by early zygotic genes, which may remain transcriptionally dormant for thousands of generations between sexual cycles. The exceptions to this rule were around 120 loci, some of which encode non-poly-adenylated transcripts, such as small nuclear RNAs and replication-dependent histones that had H3K4me3 peaks only when they were being transcribed. Mono-methylated H3K4 was the default state for the vast majority of histones that were bound outside of transcription start sites and terminator regions of genes. A small fraction of the genome that was depleted of any H3 lysine 4 methylation was enriched for DNA cytosine methylation and the genes within these DNA methylation islands were poorly expressed. Besides marking protein coding genes, H3K4me3 ChIP-seq data served also as a annotation tool for validation of hundreds of long non-coding RNA genes.

Identification of cell-type specific alternative transcripts in the multicellular alga Volvox carteri.

BACKGROUND: Cell type specialization is a hallmark of complex multicellular organisms and is usually established through implementation of cell-type-specific gene expression programs. The multicellular green alga Volvox carteri has just two cell types, germ and soma, that have previously been shown to have very different transcriptome compositions which match their specialized roles. Here we interrogated another potential mechanism for differentiation in V. carteri, cell type specific alternative transcript isoforms (CTSAI). METHODS: We used pre-existing predictions of alternative transcripts and de novo transcript assembly with HISAT2 and Ballgown software to compile a list of loci with two or more transcript isoforms, identified a small subset that were candidates for CTSAI, and manually curated this subset of genes to remove false positives. We experimentally verified three candidates using semi-quantitative RT-PCR to assess relative isoform abundance in each cell type. RESULTS: Of the 1978 loci with two or more predicted transcript isoforms 67 of these also showed cell type isoform expression biases. After curation 15 strong candidates for CTSAI were identified, three of which were experimentally verified, and their predicted gene product functions were evaluated in light of potential cell type specific roles. A comparison of genes with predicted alternative splicing from Chlamydomonas reinhardtii, a unicellular relative of V. carteri, identified little overlap between ortholog pairs with alternative splicing in both species. Finally, we interrogated cell type expression patterns of 126 V. carteri predicted RNA binding protein (RBP) encoding genes and found 40 that showed either somatic or germ cell expression bias. These RBPs are potential mediators of CTSAI in V. carteri and suggest possible pre-adaptation for cell type specific RNA processing and a potential path for generating CTSAI in the early ancestors of metazoans and plants. CONCLUSIONS: We predicted numerous instances of alternative transcript isoforms in Volvox, only a small subset of which showed cell type specific isoform expression bias. However, the validated examples of CTSAI supported existing hypotheses about cell type specialization in V. carteri, and also suggested new hypotheses about mechanisms of functional specialization for their gene products. Our data imply that CTSAI operates as a minor but important component of V. carteri cellular differentiation and could be used as a model for how alternative isoforms emerge and co-evolve with cell type specialization.

Phosphorus Availability Regulates TORC1 Signaling via LST8 in Chlamydomonas.

Target of rapamycin complex 1 (TORC1) is a central regulator of cell growth. It balances anabolic and catabolic processes in response to nutrients, growth factors, and energy availability. Nitrogen- and carbon-containing metabolites have been shown to activate TORC1 in yeast, animals, and plants. Here, we show that phosphorus (P) regulates TORC1 signaling in the model green alga Chlamydomonas (Chlamydomonas reinhardtii) via LST8, a conserved TORC1 subunit that interacts with the kinase domain of TOR. P starvation results in a sharp decrease in LST8 abundance and downregulation of TORC1 activity. A hypomorphic lst8 mutation resulted in decreased LST8 abundance, and it both reduced TORC1 signaling and altered the cellular response to P starvation. Additionally, we found that LST8 levels and TORC1 activity were not properly regulated in a mutant defective in the transcription factor PSR1, which is the major mediator of P deprivation responses in Chlamydomonas. Unlike wild-type cells, the psr1 mutant failed to downregulate LST8 abundance and TORC1 activity when under P limitation. These results identify PSR1 as an upstream regulator of TORC1 and demonstrate that TORC1 is a key component in P signaling in Chlamydomonas.

Cleavage-furrow formation without F-actin in Chlamydomonas.

It is widely believed that cleavage-furrow formation during cytokinesis is driven by the contraction of a ring containing F-actin and type-II myosin. However, even in cells that have such rings, they are not always essential for furrow formation. Moreover, many taxonomically diverse eukaryotic cells divide by furrowing but have no type-II myosin, making it unlikely that an actomyosin ring drives furrowing. To explore this issue further, we have used one such organism, the green alga Chlamydomonas reinhardtii We found that although F-actin is associated with the furrow region, none of the three myosins (of types VIII and XI) is localized there. Moreover, when F-actin was eliminated through a combination of a mutation and a drug, furrows still formed and the cells divided, although somewhat less efficiently than normal. Unexpectedly, division of the large Chlamydomonas chloroplast was delayed in the cells lacking F-actin; as this organelle lies directly in the path of the cleavage furrow, this delay may explain, at least in part, the delay in cytokinesis itself. Earlier studies had shown an association of microtubules with the cleavage furrow, and we used a fluorescently tagged EB1 protein to show that microtubules are still associated with the furrows in the absence of F-actin, consistent with the possibility that the microtubules are important for furrow formation. We suggest that the actomyosin ring evolved as one way to improve the efficiency of a core process for furrow formation that was already present in ancestral eukaryotes.

Evolutionary divergence of the sex-determining gene MID uncoupled from the transition to anisogamy in volvocine algae.

Volvocine algae constitute a unique comparative model for investigating the evolution of oogamy from isogamous mating types. The sex- or mating type-determining gene MID encodes a conserved RWP-RK transcription factor found in either the MT- or male mating locus of dioecious volvocine species. We previously found that MID from the isogamous species Chlamydomonas reinhardtii (CrMID) could not induce ectopic spermatogenesis when expressed heterologously in Volvox carteri females, suggesting coevolution of Mid function with gamete dimorphism. Here we found that ectopic expression of MID from the anisogamous species Pleodorina starrii (PsMID) could efficiently induce spermatogenesis when expressed in V. carteri females and, unexpectedly, that GpMID from the isogamous species Gonium pectorale was also able to induce V. carteri spermatogenesis. Neither VcMID nor GpMID could complement a C. reinhardtii mid mutant, at least partly owing to instability of heterologous Mid proteins. Our data show that Mid divergence was not a major contributor to the transition between isogamy and anisogamy/oogamy in volvocine algae, and instead implicate changes in cis-regulatory interactions and/or trans-acting factors of the Mid network in the evolution of sexual dimorphism.

Inositol polyphosphates and target of rapamycin kinase signalling govern photosystem II protein phosphorylation and photosynthetic function under light stress in Chlamydomonas.

Stress and nutrient availability influence cell proliferation through complex intracellular signalling networks. In a previous study it was found that pyro-inositol polyphosphates (InsP7 and InsP8 ) produced by VIP1 kinase, and target of rapamycin (TOR) kinase signalling interacted synergistically to control cell growth and lipid metabolism in the green alga Chlamydomonas reinhardtii. However, the relationship between InsPs and TOR was not completely elucidated. We used an in vivo assay for TOR activity together with global proteomic and phosphoproteomic analyses to assess differences between wild-type and vip1-1 in the presence and absence of rapamycin. We found that TOR signalling is more severely affected by the inhibitor rapamycin in a vip1-1 mutant compared with wild-type, indicating that InsP7 and InsP8 produced by VIP1 act independently but also coordinately with TOR. Additionally, among hundreds of differentially phosphorylated peptides detected, an enrichment for photosynthesis-related proteins was observed, particularly photosystem II proteins. The significance of these results was underscored by the finding that vip1-1 strains show multiple defects in photosynthetic physiology that were exacerbated under high light conditions. These results suggest a novel role for inositol pyrophosphates and TOR signalling in coordinating photosystem phosphorylation patterns in Chlamydomonas cells in response to light stress and possibly other stresses.

Modulation of CHT7 Complexes during Light/Dark- and Nitrogen-Mediated Life Cycle Transitions of Chlamydomonas.

The Chlamydomonas reinhardtii Compromised Hydrolysis of Triacylglycerols7 (CHT7) protein has been previously implicated in the regulation of DNA metabolism and cell-cycle-related gene expression during nitrogen (N) deprivation, and its predicted protein interaction domains are necessary for function. Here, we examined impacts of the cht7 mutation during the cell division cycle under nutrient deficiency in light-dark synchronized cultures. We explored the potential mechanisms affecting CHT7 complex activities during the cell cycle and N starvation, with a focus on the possible interaction between CHT7 and the C. reinhardtii retinoblastoma tumor suppressor (RB) protein homolog MAT3. Notably, the absence of CHT7 did not negatively impact the synchrony of cell division and cell cycle progression during diel growth. Although the majority of CHT7 and MAT3/RB proteins were observed in separate complexes by blue native-PAGE, the two proteins coimmunoprecipitated both during synchronized growth and following N deprivation, suggesting the presence of low abundance subcomplexes containing CHT7 and MAT3/RB. Furthermore, we observed several phosphorylated isoforms of CHT7 under these conditions. To test the potential role of phosphorylation on the structure and function of CHT7, we performed site-directed mutagenesis of previously identified phosphorylated amino acids within CHT7. These phosphorylated residues were dispensable for CHT7 function, but phosphorylated variants of CHT7 persisted, indicating that yet-unidentified residues within CHT7 are also likely phosphorylated. Based on the interaction of CHT7 and MAT3/RB, we postulate the presence of a low-abundance or transient regulatory complex in C. reinhardtii that may be similar to DREAM-like complexes in other organisms.

Investigating the effect of target of rapamycin kinase inhibition on the Chlamydomonas reinhardtii phosphoproteome: from known homologs to new targets.

Target of rapamycin (TOR) kinase is a conserved regulator of cell growth whose activity is modulated in response to nutrients, energy and stress. Key proteins involved in the pathway are conserved in the model photosynthetic microalga Chlamydomonas reinhardtii, but the substrates of TOR kinase and downstream signaling network have not been elucidated. Our study provides a new resource for investigating the phosphorylation networks governed by the TOR kinase pathway in Chlamydomonas. We used quantitative phosphoproteomics to investigate the effects of inhibiting Chlamydomonas TOR kinase on dynamic protein phosphorylation. Wild-type and AZD-insensitive Chlamydomonas strains were treated with TOR-specific chemical inhibitors (rapamycin, AZD8055 and Torin1), after which differentially affected phosphosites were identified. Our quantitative phosphoproteomic dataset comprised 2547 unique phosphosites from 1432 different proteins. Inhibition of TOR kinase caused significant quantitative changes in phosphorylation at 258 phosphosites, from 219 unique phosphopeptides. Our results include Chlamydomonas homologs of TOR signaling-related proteins, including a site on RPS6 with a decrease in phosphorylation. Additionally, phosphosites on proteins involved in translation and carotenoid biosynthesis were identified. Follow-up experiments guided by these phosphoproteomic findings in lycopene beta/epsilon cyclase showed that carotenoid levels are affected by TORC1 inhibition and carotenoid production is under TOR control in algae.

Synergism between Inositol Polyphosphates and TOR Kinase Signaling in Nutrient Sensing, Growth Control, and Lipid Metabolism in Chlamydomonas.

The networks that govern carbon metabolism and control intracellular carbon partitioning in photosynthetic cells are poorly understood. Target of Rapamycin (TOR) kinase is a conserved growth regulator that integrates nutrient signals and modulates cell growth in eukaryotes, though the TOR signaling pathway in plants and algae has yet to be completely elucidated. We screened the unicellular green alga Chlamydomonas reinhardtii using insertional mutagenesis to find mutants that conferred hypersensitivity to the TOR inhibitor rapamycin. We characterized one mutant, vip1-1, that is predicted to encode a conserved inositol hexakisphosphate kinase from the VIP family that pyrophosphorylates phytic acid (InsP6) to produce the low abundance signaling molecules InsP7 and InsP8 Unexpectedly, the rapamycin hypersensitive growth arrest of vip1-1 cells was dependent on the presence of external acetate, which normally has a growth-stimulatory effect on Chlamydomonas. vip1-1 mutants also constitutively overaccumulated triacylglycerols (TAGs) in a manner that was synergistic with other TAG inducing stimuli such as starvation. vip1-1 cells had reduced InsP7 and InsP8, both of which are dynamically modulated in wild-type cells by TOR kinase activity and the presence of acetate. Our data uncover an interaction between the TOR kinase and inositol polyphosphate signaling systems that we propose governs carbon metabolism and intracellular pathways that lead to storage lipid accumulation.

Identification of Chlamydomonas reinhardtii endogenous genic flanking sequences for improved transgene expression.

Chlamydomonas reinhardtii is a unicellular green alga that has attracted interest due to its potential biotechnological applications, and as a model for algal biofuel and energy metabolism. Despite all the advantages that this unicellular alga offers, poor and inconsistent expression of nuclear transgenes remains an obstacle for basic and applied research. We used a data-mining strategy to identify highly expressed genes in Chlamydomonas whose flanking sequences were tested for the ability to drive heterologous nuclear transgene expression. Candidates identified in this search included two ribosomal protein genes, RPL35a and RPL23, and ferredoxin, FDX1, whose flanking regions including promoters, terminators and untranslated sequences could drive stable luciferase transgene expression to significantly higher levels than the commonly used Hsp70A-RBCS2 (AR) hybrid promoter/terminator sequences. The RPL23 flanking sequences were further tested using the zeocin resistance gene sh-ble as a reporter in monocistronic and dicistronic constructs, and consistently yielded higher numbers of zeocin-resistant transformants and higher levels of resistance than AR- or PSAD-based vectors. Chlamydomonas RPL23 sequences also enabled transgene expression in Volvox carteri. Our study provides an additional benchmark for strong constitutive expression of transgenes in Chlamydomonas, and develops a general approach for identifying flanking sequences that can be used to drive transgene expression for any organism where transcriptome data are available.

Probing the global kinome and phosphoproteome in Chlamydomonas reinhardtii via sequential enrichment and quantitative proteomics.

The identification of dynamic protein phosphorylation events is critical for understanding kinase/phosphatase-regulated signaling pathways. To date, protein phosphorylation and kinase expression have been examined independently in photosynthetic organisms. Here we present a method to study the global kinome and phosphoproteome in tandem in a model photosynthetic organism, the alga Chlamydomonas reinhardtii (Chlamydomonas), using mass spectrometry-based label-free proteomics. A dual enrichment strategy targets intact protein kinases via capture on immobilized multiplexed inhibitor beads with subsequent proteolytic digestion of unbound proteins and peptide-based phosphorylation enrichment. To increase depth of coverage, both data-dependent and data-independent (via SWATH, Sequential Windowed Acquisition of All Theoretical Fragment Ion Mass Spectra) mass spectrometric acquisitions were performed to obtain a more than 50% increase in coverage of the enriched Chlamydomonas kinome over coverage found with no enrichment. The quantitative phosphoproteomic dataset yielded 2250 phosphopeptides and 1314 localized phosphosites with excellent reproducibility across biological replicates (90% of quantified sites with coefficient of variation below 11%). This approach enables simultaneous investigation of kinases and phosphorylation events at the global level to facilitate understanding of kinase networks and their influence in cell signaling events.

Autophagic flux is required for the synthesis of triacylglycerols and ribosomal protein turnover in Chlamydomonas.

Autophagy is an intracellular catabolic process that allows cells to recycle unneeded or damaged material to maintain cellular homeostasis. This highly dynamic process is characterized by the formation of double-membrane vesicles called autophagosomes, which engulf and deliver the cargo to the vacuole. Flow of material through the autophagy pathway and its degradation in the vacuole is known as autophagic flux, and reflects the autophagic degradation activity. A number of assays have been developed to determine autophagic flux in yeasts, mammals, and plants, but it has not been examined yet in algae. Here we analyzed autophagic flux in the model green alga Chlamydomonas reinhardtii. By monitoring specific autophagy markers such as ATG8 lipidation and using immunofluorescence and electron microscopy techniques, we show that concanamycin A, a vacuolar ATPase inhibitor, blocks autophagic flux in Chlamydomonas. Our results revealed that vacuolar lytic function is needed for the synthesis of triacylglycerols and the formation of lipid bodies in nitrogen- or phosphate-starved cells. Moreover, we found that concanamycin A treatment prevented the degradation of ribosomal proteins RPS6 and RPL37 under nitrogen or phosphate deprivation. These results indicate that autophagy might play an important role in the regulation of lipid metabolism and the recycling of ribosomal proteins under nutrient limitation in Chlamydomonas.

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.

Non-model model organisms.

Model organisms are widely used in research as accessible and convenient systems to study a particular area or question in biology. Traditionally only a handful of organisms have been widely studied, but modern research tools are enabling researchers to extend the set of model organisms to include less-studied and more unusual systems. This Forum highlights a range of 'non-model model organisms' as emerging systems for tackling questions across the whole spectrum of biology (and beyond), the opportunities and challenges, and the outlook for the future.