Mutation of negative regulatory gene CEHC1 encoding an FBXO3 protein results in normoxic expression of HYDA genes in Chlamydomonas reinhardtii.

Oxygen is known to prevent hydrogen production in Chlamydomonas, both by inhibiting the hydrogenase enzyme and by preventing the accumulation of HYDA-encoding transcripts. We developed a screen for mutants showing constitutive accumulation of HYDA1 transcripts in the presence of oxygen. A reporter gene required for ciliary motility, placed under the control of the HYDA1 promoter, conferred motility only in hypoxic conditions. By selecting for mutants able to swim even in the presence of oxygen we obtained strains that express the reporter gene constitutively. One mutant identified a gene encoding an F-box only protein 3 (FBXO3), known to participate in ubiquitylation and proteasomal degradation pathways in other eukaryotes. Transcriptome profiles revealed that the mutation, termed cehc1-1 , leads to constitutive expression of HYDA1 and other genes regulated by hypoxia, and of many genes known to be targets of CRR1, a transcription factor in the nutritional copper signaling pathway. CRR1 was required for the constitutive expression of the HYDA1 reporter gene in cehc1-1 mutants. The CRR1 protein, which is normally degraded in Cu-supplemented cells, was stabilized in cehc1-1 cells, supporting the conclusion that CEHC1 acts to facilitate the degradation of CRR1. Our results reveal a novel negative regulator in the CRR1 pathway and possibly other pathways leading to complex metabolic changes associated with response to hypoxia.

Distinct function of Chlamydomonas CTRA-CTR transporters in Cu assimilation and intracellular mobilization.

Successful acclimation to copper (Cu) deficiency involves a fine balance between Cu import and export. In the green alga Chlamydomonas reinhardtii, Cu import is dependent on a transcription factor, Copper Response Regulator 1 (CRR1), responsible for activating genes in Cu-deficient cells. Among CRR1 target genes are two Cu transporters belonging to the CTR/COPT gene family (CTR1 and CTR2) and a related soluble protein (CTR3). The ancestor of these green algal proteins was likely acquired from an ancient chytrid and contained conserved cysteine-rich domains (named the CTR-associated domains, CTRA) that are predicted to be involved in Cu acquisition. We show by reverse genetics that Chlamydomonas CTR1 and CTR2 are canonical Cu importers albeit with distinct affinities, while loss of CTR3 did not result in an observable phenotype under the conditions tested. Mutation of CTR1, but not CTR2, recapitulates the poor growth of crr1 in Cu-deficient medium, consistent with a dominant role for CTR1 in high-affinity Cu(I) uptake. On the other hand, the overaccumulation of Cu(I) (20 times the quota) in zinc (Zn) deficiency depends on CRR1 and both CTR1 and CTR2. CRR1-dependent activation of CTR gene expression needed for Cu over-accumulation can be bypassed by the provision of excess Cu in the growth medium. Over-accumulated Cu is sequestered into the acidocalcisome but can become remobilized by restoring Zn nutrition. This mobilization is also CRR1-dependent, and requires activation of CTR2 expression, again distinguishing CTR2 from CTR1 and consistent with the lower substrate affinity of CTR2. ONE SENTENCE SUMMARY: Regulation of Cu uptake and sequestration by members of the CTR family of proteins in Chlamydomonas.

Scarcity of fixed carbon transfer in a model microbial phototroph-heterotroph interaction.

While the green alga Chlamydomonas reinhardtii has long served as a reference organism, few studies have interrogated its role as a primary producer in microbial interactions. Here, we quantitatively investigated C. reinhardtii's capacity to support a heterotrophic microbe using the established coculture system with Mesorhizobium japonicum, a vitamin B12-producing α-proteobacterium. Using stable isotope probing and nanoscale secondary ion mass spectrometry (nanoSIMS), we tracked the flow of photosynthetic fixed carbon and consequent bacterial biomass synthesis under continuous and diel light with single-cell resolution. We found that more 13C fixed by the alga was taken up by bacterial cells under continuous light, invalidating the hypothesis that the alga's fermentative degradation of starch reserves during the night would boost M. japonicum heterotrophy. 15NH4 assimilation rates and changes in cell size revealed that the carbon transferred was insufficient for balanced growth of M. japonicum cells, which instead underwent reductive division. However, despite this sign of starvation, M. japonicum still supported a B12-dependent C. reinhardtii mutant. Finally, we showed that bacterial proliferation could be supported solely by the algal lysis that occurred in coculture, highlighting the role of necromass in carbon cycling. Collectively, these results reveal the scarcity of fixed carbon in this microbial trophic relationship, demonstrate B12 exchange even during bacterial starvation, and underscore the importance of quantitative approaches for assessing metabolic coupling in algal-bacterial interactions.

Distinct function of Chlamydomonas CTRA-CTR transporters in Cu assimilation and intracellular mobilization.

Successful acclimation to copper (Cu) deficiency involves a fine balance between Cu import and export. In the unicellular green alga Chlamydomonas reinhardtii , Cu import is dependent on C opper R esponse R egulator 1 (CRR1), the master regulator of Cu homeostasis. Among CRR1 target genes are two Cu transporters belonging to the CTR/COPT gene family ( CTR1 and CTR2 ) and a related soluble cysteine-rich protein (CTR3). The ancestor of these green algal proteins was likely acquired from an ancient chytrid and contained conserved cysteine-rich domains (named the CTR-associated domains, CTRA) that are predicted to be involved in Cu acquisition. We show by reverse genetics that Chlamydomonas CTR1 and CTR2 are canonical Cu importers albeit with distinct affinities, while loss of CTR3 did not result in an observable phenotype under the conditions tested. Mutation of CTR1 , but not CTR2 , recapitulate the poor growth of crr1 in Cu-deficient medium, consistent with a dominant role for CTR1 in high affinity Cu(I) uptake. Notably, the over-accumulation of Cu(I) in Zinc (Zn)-deficiency (20 times the quota) depends on CRR1 and both CTR1 and CTR2. CRR1-dependent activation of CTR gene expression needed for Cu over-accumulation can be bypassed by the provision of excess Cu in the growth medium. Over-accumulated Cu is sequestered into the acidocalcisome but can become remobilized by restoring Zn nutrition. This mobilization is also CRR1-dependent, and requires activation of CTR2 expression, again distinguishing CTR2 from CTR1 and is consistent with the lower substrate affinity of CTR2.

Zn deficiency disrupts Cu and S homeostasis in Chlamydomonas resulting in over accumulation of Cu and Cysteine.

Growth of Chlamydomonas reinhardtii in zinc (Zn) limited medium leads to disruption of copper (Cu) homeostasis, resulting in up to 40-fold Cu over-accumulation relative to its typical Cu quota. We show that Chlamydomonas controls its Cu quota by balancing Cu import and export, which is disrupted in a Zn deficient cell, thus establishing a mechanistic connection between Cu and Zn homeostasis. Transcriptomics, proteomics and elemental profiling revealed that Zn-limited Chlamydomonas cells up-regulate a subset of genes encoding "first responder" proteins involved in sulfur (S) assimilation and consequently accumulate more intracellular S, which is incorporated into L-cysteine, γ-glutamylcysteine, and homocysteine. Most prominently, in the absence of Zn, free L-cysteine is increased ∼80-fold, corresponding to ∼2.8 × 109 molecules/cell. Interestingly, classic S-containing metal binding ligands like glutathione and phytochelatins do not increase. X-ray fluorescence microscopy showed foci of S accumulation in Zn-limited cells that co-localize with Cu, phosphorus and calcium, consistent with Cu-thiol complexes in the acidocalcisome, the site of Cu(I) accumulation. Notably, cells that have been previously starved for Cu do not accumulate S or Cys, causally connecting cysteine synthesis with Cu accumulation. We suggest that cysteine is an in vivo Cu(I) ligand, perhaps ancestral, that buffers cytosolic Cu.

Pumping iron: A multi-omics analysis of two extremophilic algae reveals iron economy management.

Marine algae are responsible for half of the world's primary productivity, but this critical carbon sink is often constrained by insufficient iron. One species of marine algae, Dunaliella tertiolecta, is remarkable for its ability to maintain photosynthesis and thrive in low-iron environments. A related species, Dunaliella salina Bardawil, shares this attribute but is an extremophile found in hypersaline environments. To elucidate how algae manage their iron requirements, we produced high-quality genome assemblies and transcriptomes for both species to serve as a foundation for a comparative multiomics analysis. We identified a host of iron-uptake proteins in both species, including a massive expansion of transferrins and a unique family of siderophore-iron-uptake proteins. Complementing these multiple iron-uptake routes, ferredoxin functions as a large iron reservoir that can be released by induction of flavodoxin. Proteomic analysis revealed reduced investment in the photosynthetic apparatus coupled with remodeling of antenna proteins by dramatic iron-deficiency induction of TIDI1, which is closely related but identifiably distinct from the chlorophyll binding protein, LHCA3. These combinatorial iron scavenging and sparing strategies make Dunaliella unique among photosynthetic organisms.

Optimizing Fe Nutrition for Algal Growth.

Chlamydomonas is an excellent reference system for dissecting the impact of iron (Fe) nutrition on photosynthetic and other metabolisms. The operational definition of four stages of Fe nutrition is described and a guide to the practical use of these stages is offered, specifically the preparation of media and growth of mixotrophic cultures. A key consideration is the impact of carbon metabolism on the expression of Fe-containing enzymes and hence the Fe quota. The absolute concentration of Fe in the medium is less determinative of gene expression than the Fe available on a per-cell basis. In nature, algal cells may transition from Fe-replete to -deficient to -limited during a bloom.

Quantitative elemental imaging in eukaryotic algae.

All organisms, fundamentally, are made from the same raw material, namely the elements of the periodic table. Biochemical diversity is achieved by how these elements are utilized, for what purpose, and in which physical location. Determining elemental distributions, especially those of trace elements that facilitate metabolism as cofactors in the active centers of essential enzymes, can determine the state of metabolism, the nutritional status, or the developmental stage of an organism. Photosynthetic eukaryotes, especially algae, are excellent subjects for quantitative analysis of elemental distribution. These microbes utilize unique metabolic pathways that require various trace nutrients at their core to enable their operation. Photosynthetic microbes also have important environmental roles as primary producers in habitats with limited nutrient supplies or toxin contaminations. Accordingly, photosynthetic eukaryotes are of great interest for biotechnological exploitation, carbon sequestration, and bioremediation, with many of the applications involving various trace elements and consequently affecting their quota and intracellular distribution. A number of diverse applications were developed for elemental imaging, allowing subcellular resolution, with X-ray fluorescence microscopy (XFM, XRF) being at the forefront, enabling quantitative descriptions of intact cells in a non-destructive method. This Tutorial Review summarizes the workflow of a quantitative, single-cell elemental distribution analysis of a eukaryotic alga using XFM.

Cysteine: an ancestral Cu binding ligand in green algae?

Growth of Chlamydomonas reinhardtii in zinc (Zn) limited medium leads to disruption of copper (Cu) homeostasis, resulting in up to 40-fold Cu over-accumulation relative to its typical Cu quota. We show that Chlamydomonas controls its Cu quota by balancing Cu import and export, which is disrupted in a Zn deficient cell, thus establishing a mechanistic connection between Cu and Zn homeostasis. Transcriptomics, proteomics and elemental profiling revealed that Zn-limited Chlamydomonas cells up-regulate a subset of genes encoding "first responder" proteins involved in sulfur (S) assimilation and consequently accumulate more intracellular S, which is incorporated into L-cysteine, γ-glutamylcysteine and homocysteine. Most prominently, in the absence of Zn, free L-cysteine is increased ~80-fold, corresponding to ~ 2.8 × 10 9 molecules/cell. Interestingly, classic S-containing metal binding ligands like glutathione and phytochelatins do not increase. X-ray fluorescence microscopy showed foci of S accumulation in Zn-limited cells that co-localize with Cu, phosphorus and calcium, consistent with Cu-thiol complexes in the acidocalcisome, the site of Cu(I) accumulation. Notably, cells that have been previously starved for Cu do not accumulate S or Cys, causally connecting cysteine synthesis with Cu accumulation. We suggest that cysteine is an in vivo Cu(I) ligand, perhaps ancestral, that buffers cytosolic Cu.

Structural and functional regulation of Chlamydomonas lysosome-related organelles during environmental changes.

Lysosome-related organelles (LROs) are a class of heterogeneous organelles conserved in eukaryotes that primarily play a role in storage and secretion. An important function of LROs is to mediate metal homeostasis. Chlamydomonas reinhardtii is a model organism for studying metal ion metabolism; however, structural and functional analyses of LROs in C. reinhardtii are insufficient. Here, we optimized a method for purifying these organelles from 2 populations of cells: stationary phase or overloaded with iron. The morphology, elemental content, and lysosomal activities differed between the 2 preparations, even though both have phosphorus and metal ion storage functions. LROs in stationary phase cells had multiple non-membrane-bound polyphosphate granules to store phosphorus. Those in iron-overloaded cells were similar to acidocalcisomes (ACs), which have a boundary membrane and contain 1 or 2 large polyphosphate granules to store more phosphorus. We established a method for quantifying the capacity of LROs to sequester individual trace metals. Based on a comparative proteomic analysis of these 2 types of LROs, we present a comprehensive AC proteome and identified 113 putative AC proteins. The methods and protein inventories provide a framework for studying the biogenesis and modification of LROs and the mechanisms by which they participate in regulating metal ion metabolism.

The Chlamydomonas Genome Project, version 6: Reference assemblies for mating-type plus and minus strains reveal extensive structural mutation in the laboratory.

Five versions of the Chlamydomonas reinhardtii reference genome have been produced over the last two decades. Here we present version 6, bringing significant advances in assembly quality and structural annotations. PacBio-based chromosome-level assemblies for two laboratory strains, CC-503 and CC-4532, provide resources for the plus and minus mating-type alleles. We corrected major misassemblies in previous versions and validated our assemblies via linkage analyses. Contiguity increased over ten-fold and >80% of filled gaps are within genes. We used Iso-Seq and deep RNA-seq datasets to improve structural annotations, and updated gene symbols and textual annotation of functionally characterized genes via extensive manual curation. We discovered that the cell wall-less classical reference strain CC-503 exhibits genomic instability potentially caused by deletion of the helicase RECQ3, with major structural mutations identified that affect >100 genes. We therefore present the CC-4532 assembly as the primary reference, although this strain also carries unique structural mutations and is experiencing rapid proliferation of a Gypsy retrotransposon. We expect all laboratory strains to harbor gene-disrupting mutations, which should be considered when interpreting and comparing experimental results. Collectively, the resources presented here herald a new era of Chlamydomonas genomics and will provide the foundation for continued research in this important reference organism.

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.

Simple steps to enable reproducibility: culture conditions affecting Chlamydomonas growth and elemental composition.

Even subtle modifications in growth conditions elicit acclimation responses affecting the molecular and elemental makeup of organisms, both in the laboratory and in natural habitats. We systematically explored the effect of temperature, pH, nutrient availability, culture density, and access to CO2 and O2 in laboratory-grown algal cultures on growth rate, the ionome, and the ability to accumulate Fe. We found algal cells accumulate Fe in alkaline conditions, even more so when excess Fe is present, coinciding with a reduced growth rate. Using a combination of Fe-specific dyes, X-ray fluorescence microscopy, and NanoSIMS, we show that the alkaline-accumulated Fe was intracellularly sequestered into acidocalcisomes, which are localized towards the periphery of the cells. At high photon flux densities, Zn and Ca specifically over-accumulate, while Zn alone accumulates at low temperatures. The impact of aeration was probed by reducing shaking speeds and changing vessel fill levels; the former increased the Cu quota of cultures, the latter resulted in a reduction in P, Ca, and Mn at low fill levels. Trace element quotas were also affected in the stationary phase, where specifically Fe, Cu, and Zn accumulate. Cu accumulation here depends inversely on the Fe concentration of the medium. Individual laboratory strains accumulate Ca, P, and Cu to different levels. All together, we identified a set of specific changes to growth rate, elemental composition, and the capacity to store Fe in response to subtle differences in culturing conditions of Chlamydomonas, affecting experimental reproducibility. Accordingly, we recommend that these variables be recorded and reported as associated metadata.

Systematic characterization of gene function in the photosynthetic alga Chlamydomonas reinhardtii.

Most genes in photosynthetic organisms remain functionally uncharacterized. Here, using a barcoded mutant library of the model eukaryotic alga Chlamydomonas reinhardtii, we determined the phenotypes of more than 58,000 mutants under more than 121 different environmental growth conditions and chemical treatments. A total of 59% of genes are represented by at least one mutant that showed a phenotype, providing clues to the functions of thousands of genes. Mutant phenotypic profiles place uncharacterized genes into functional pathways such as DNA repair, photosynthesis, the CO2-concentrating mechanism and ciliogenesis. We illustrate the value of this resource by validating phenotypes and gene functions, including three new components of an actin cytoskeleton defense pathway. The data also inform phenotype discovery in land plants; mutants in Arabidopsis thaliana genes exhibit phenotypes similar to those we observed in their Chlamydomonas homologs. We anticipate that this resource will guide the functional characterization of genes across the tree of life.

An optimized ChIP-Seq framework for profiling histone modifications in Chromochloris zofingiensis.

The eukaryotic green alga Chromochloris zofingiensis is a reference organism for studying carbon partitioning and a promising candidate for the production of biofuel precursors. Recent transcriptome profiling transformed our understanding of its biology and generally algal biology, but epigenetic regulation remains understudied and represents a fundamental gap in our understanding of algal gene expression. Chromatin immunoprecipitation followed by deep sequencing (ChIP-Seq) is a powerful tool for the discovery of such mechanisms, by identifying genome-wide histone modification patterns and transcription factor-binding sites alike. Here, we established a ChIP-Seq framework for Chr. zofingiensis yielding over 20 million high-quality reads per sample. The most critical steps in a ChIP experiment were optimized, including DNA shearing to obtain an average DNA fragment size of 250 bp and assessment of the recommended formaldehyde concentration for optimal DNA-protein cross-linking. We used this ChIP-Seq framework to generate a genome-wide map of the H3K4me3 distribution pattern and to integrate these data with matching RNA-Seq data. In line with observations from other organisms, H3K4me3 marks predominantly transcription start sites of genes. Our H3K4me3 ChIP-Seq data will pave the way for improved genome structural annotation in the emerging reference alga Chr. zofingiensis.

Chlamydomonas ATX1 is essential for Cu distribution to multiple cupro-enzymes and maintenance of biomass in conditions demanding cupro-enzyme-dependent metabolic pathways.

Copper (Cu) chaperones, of which yeast ATX1 is a prototype, are small proteins with a Cu(I) binding MxCxxC motif and are responsible for directing intracellular Cu toward specific client protein targets that use Cu as a cofactor. The Chlamydomonas reinhardtii ATX1 (CrATX1) was identified by its high sequence similarity with yeast ATX1. Like the yeast homologue, CrATX1 accumulates in iron-deficient cells (but is not impacted by other metal-deficiencies). N- and C-terminally YFP-ATX1 fusion proteins are distributed in the cytoplasm. Reverse genetic analysis using artificial microRNA (amiRNA) to generate lines with reduced CrATX1 abundance and CRISPR/Cpf1 to generate atx1 knockout lines validated a function for ATX1 in iron-poor cells, again reminiscent of yeast ATX1, most likely because of an impact on metalation of the multicopper oxidase FOX1, which is an important component in high-affinity iron uptake. We further identify other candidate ATX1 targets owing to reduced growth of atx1 mutant lines on guanine as a sole nitrogen source, which we attribute to loss of function of UOX1, encoding a urate oxidase, a cupro-enzyme involved in guanine assimilation. An impact of ATX1 on Cu distribution in atx1 mutants is strikingly evident by a reduced amount of intracellular Cu in all conditions probed in this work.

A pathogenic role for histone H3 copper reductase activity in a yeast model of Friedreich's ataxia.

Disruptions to iron-sulfur (Fe-S) clusters, essential cofactors for a broad range of proteins, cause widespread cellular defects resulting in human disease. A source of damage to Fe-S clusters is cuprous (Cu1+) ions. Since histone H3 enzymatically produces Cu1+ for copper-dependent functions, we asked whether this activity could become detrimental to Fe-S clusters. Here, we report that histone H3­mediated Cu1+ toxicity is a major determinant of cellular functional pool of Fe-S clusters. Inadequate Fe-S cluster supply, due to diminished assembly as occurs in Friedreich's ataxia or defective distribution, causes severe metabolic and growth defects in Saccharomyces cerevisiae. Decreasing Cu1+ abundance, through attenuation of histone cupric reductase activity or depletion of total cellular copper, restored Fe-S cluster­dependent metabolism and growth. Our findings reveal an interplay between chromatin and mitochondria in Fe-S cluster homeostasis and a potential pathogenic role for histone enzyme activity and Cu1+ in diseases with Fe-S cluster dysfunction.

Single-cell visualization and quantification of trace metals in Chlamydomonas lysosome-related organelles.

The acidocalcisome is an acidic organelle in the cytosol of eukaryotes, defined by its low pH and high calcium and polyphosphate content. It is visualized as an electron-dense object by transmission electron microscopy (TEM) or described with mass spectrometry (MS)-based imaging techniques or multimodal X-ray fluorescence microscopy (XFM) based on its unique elemental composition. Compared with MS-based imaging techniques, XFM offers the additional advantage of absolute quantification of trace metal content, since sectioning of the cell is not required and metabolic states can be preserved rapidly by either vitrification or chemical fixation. We employed XFM in Chlamydomonas reinhardtii to determine single-cell and organelle trace metal quotas within algal cells in situations of trace metal overaccumulation (Fe and Cu). We found up to 70% of the cellular Cu and 80% of Fe sequestered in acidocalcisomes in these conditions and identified two distinct populations of acidocalcisomes, defined by their unique trace elemental makeup. We utilized the vtc1 mutant, defective in polyphosphate synthesis and failing to accumulate Ca, to show that Fe sequestration is not dependent on either. Finally, quantitation of the Fe and Cu contents of individual cells and compartments via XFM, over a range of cellular metal quotas created by nutritional and genetic perturbations, indicated excellent correlation with bulk data from corresponding cell cultures, establishing a framework to distinguish the nutritional status of single cells.

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.