Whole-genome sequencing analysis of volvocine green algae reveals the molecular genetic basis for the diversity and evolution of sex.

This review describes the development of evolutionary studies of sex based on the volvocine lineage of green algae, which was facilitated by whole-genome analyses of both model and non-model species. Volvocine algae, which include Chlamydomonas and Volvox species, have long been considered a model group for experimental studies investigating the evolution of sex. Thus, whole-genomic information on the sex-determining regions of volvocine algal sex chromosomes has been sought to elucidate the molecular genetic basis of sex evolution. By 2010, whole genomes were published for two model species in this group, Chlamydomonas reinhardtii and Volvox carteri. Recent improvements in sequencing technology, particularly next-generation sequencing, allowed our studies to obtain complete genomes for non-model, but evolutionary important, volvocine algal species. These genomes have provided critical details about sex-determining regions that will contribute to our understanding of the diversity and evolution of sex.

Multicellularity and increasing Reynolds number impact on the evolutionary shift in flash-induced ciliary response in Volvocales.

BACKGROUND: Volvocales in green algae have evolved by multicellularity of Chlamydomonas-like unicellular ancestor. Those with various cell numbers exist, such as unicellular Chlamydomonas, four-celled Tetrabaena, and Volvox species with different cell numbers (~1,000, ~5,000, and ~10,000). Each cell of these organisms shares two cilia and an eyespot, which are used for swimming and photosensing. They are all freshwater microalgae but inhabit different fluid environments: unicellular species live in low Reynolds-number (Re) environments where viscous forces dominate, whereas multicellular species live in relatively higher Re where inertial forces become non-negligible. Despite significant changes in the physical environment, during the evolution of multicellularity, they maintained photobehaviors (i.e., photoshock and phototactic responses), which allows them to survive under changing light conditions. RESULTS: In this study, we utilized high-speed imaging to observe flash-induced changes in the ciliary beating manner of 27 Volvocales strains. We classified flash-induced ciliary responses in Volvocales into four patterns: "1: temporal waveform conversion", "2: no obvious response", "3: pause in ciliary beating", and "4: temporal changes in ciliary beating directions". We found that which species exhibit which pattern depends on Re, which is associated with the individual size of each species rather than phylogenetic relationships. CONCLUSIONS: These results suggest that only organisms that acquired different patterns of ciliary responses survived the evolutionary transition to multicellularity with a greater number of cells while maintaining photobehaviors. This study highlights the significance of the Re as a selection pressure in evolution and offers insights for designing propulsion systems in biomimetic micromachines.

Two species of the green algae Volvox sect. Volvox from the Japanese ancient lake, Lake Biwa.

Volvox sect. Volvox is a group of green algae with unique morphological features (thick cytoplasmic bridges between somatic cells and spiny zygote walls) and a worldwide distribution. Despite research interest in the diversity of organisms in ancient lakes, Volvox sect. Volvox from ancient lakes worldwide has not been identified to the species level. Here, we established clonal cultures of two species of this group originating from Lake Biwa, an ancient lake in Japan, and performed identification based on morphological and molecular data. One was identified as Volvox kirkiorum based on the nuclear ribosomal DNA internal spacer region (ITS) sequence, bisexual (monoicous or monoecious) spheroids, and zygote morphology. The other showed genetic separation from related species based on the secondary structure of the ITS and results of phylogenetic analysis of a combined data set from the nuclear actin gene, ITS, and two plastid genes (large subunit of RuBisCO and photosystem II CP43 apoprotein gene); it represented a new phylogenetic lineage within Volvox sect. Volvox, suggesting possible endemism in Lake Biwa. This species produced bisexual spheroids with different zygote morphology and zygote number from other species with bisexual spheroids in Volvox sect. Volvox. Therefore, Volvox biwakoensis Nozaki et H. Yamaguchi sp. nov. is described herein. This is the first endemic species of the genus Volvox described from an ancient lake.

Fossil-calibrated molecular clock data enable reconstruction of steps leading to differentiated multicellularity and anisogamy in the Volvocine algae.

BACKGROUND: Throughout its nearly four-billion-year history, life has undergone evolutionary transitions in which simpler subunits have become integrated to form a more complex whole. Many of these transitions opened the door to innovations that resulted in increased biodiversity and/or organismal efficiency. The evolution of multicellularity from unicellular forms represents one such transition, one that paved the way for cellular differentiation, including differentiation of male and female gametes. A useful model for studying the evolution of multicellularity and cellular differentiation is the volvocine algae, a clade of freshwater green algae whose members range from unicellular to colonial, from undifferentiated to completely differentiated, and whose gamete types can be isogamous, anisogamous, or oogamous. To better understand how multicellularity, differentiation, and gametes evolved in this group, we used comparative genomics and fossil data to establish a geologically calibrated roadmap of when these innovations occurred. RESULTS: Our ancestral-state reconstructions, show that multicellularity arose independently twice in the volvocine algae. Our chronograms indicate multicellularity evolved during the Carboniferous-Triassic periods in Goniaceae + Volvocaceae, and possibly as early as the Cretaceous in Tetrabaenaceae. Using divergence time estimates we inferred when, and in what order, specific developmental changes occurred that led to differentiated multicellularity and oogamy. We find that in the volvocine algae the temporal sequence of developmental changes leading to differentiated multicellularity is much as proposed by David Kirk, and that multicellularity is correlated with the acquisition of anisogamy and oogamy. Lastly, morphological, molecular, and divergence time data suggest the possibility of cryptic species in Tetrabaenaceae. CONCLUSIONS: Large molecular datasets and robust phylogenetic methods are bringing the evolutionary history of the volvocine algae more sharply into focus. Mounting evidence suggests that extant species in this group are the result of two independent origins of multicellularity and multiple independent origins of cell differentiation. Also, the origin of the Tetrabaenaceae-Goniaceae-Volvocaceae clade may be much older than previously thought. Finally, the possibility of cryptic species in the Tetrabaenaceae provides an exciting opportunity to study the recent divergence of lineages adapted to live in very different thermal environments.

Different transcription of novel, functional long non-coding RNA genes by UV-B in green algae, Volvox carteri.

Long non-coding RNAs (lncRNAs) are identified as important regulatory molecules related to diverse biological processes. In recent years, benefiting from the rapid development of high-throughput sequencing technology, RNA-seq, and analysis methods, more lncRNAs have been identified and discovered in various plant and algal species. However, so far, only limited studies related to algal lncRNAs are available. Volvox carteri f. nagariensis is the best multicellular model organism to study in developmental and evolutionary biology; therefore, studying and increasing information about this species is important. This study identified lncRNAs in the multicellular green algae Volvox carteri and 1457 lncRNAs were reported, using RNA-seq data and with the help of bioinformatics tools and software. This study investigated the effect of low-dose UV-B radiation on changes in the expression profile of lncRNAs in gonidial and somatic cells. The differential expression of lncRNAs was analyzed between the treatment (UV-B) and the control (WL) groups in gonidial and somatic cells. A total of 37 and 26 lncRNAs with significant differential expression in gonidial and somatic cells, respectively, were reported. Co-expression analysis between the lncRNAs and their neighbor protein-coding genes (in the interval of ± 10 Kb) was accomplished. In gonidial cells, 184 genes with a positive correlation and 13 genes with a negative correlation (greater than 0.95), and in somatic cells, 174 genes with a positive correlation, and 18 genes with a negative correlation were detected. Functional analysis of neighboring coding genes was also performed based on gene ontology. The results of the current work may help gain deeper insight into the regulation of gene expression in the studied model organism, Volvox carteri.

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.

Cell Type-Specific Promoters of Volvox carteri for Molecular Cell Biology Studies.

The multicellular green alga Volvox carteri has emerged as a valuable model organism for investigating various aspects of multicellularity and cellular differentiation, photoreception and phototaxis, cell division, biogenesis of the extracellular matrix and morphogenetic movements. While a range of molecular tools and bioinformatics resources have been made available for exploring these topics, the establishment of cell type-specific promoters in V. carteri has not been achieved so far. Therefore, here, we conducted a thorough screening of transcriptome data from RNA sequencing analyses of V. carteri in order to identify potential cell type-specific promoters. Eventually, we chose two putative strong and cell type-specific promoters, with one exhibiting specific expression in reproductive cells (gonidia), the PCY1 promoter, and the other in somatic cells, the PFP promoter. After cloning both promoter regions, they were introduced upstream of a luciferase reporter gene. By using particle bombardment, the DNA constructs were stably integrated into the genome of V. carteri. The results of the expression analyses, which were conducted at both the transcript and protein levels, demonstrated that the two promoters drive cell type-specific expression in their respective target cell types. Transformants with considerably diverse expression levels of the chimeric genes were identifiable. In conclusion, the screening and analysis of transcriptome data from RNA sequencing allowed for the identification of potential cell type-specific promoters in V. carteri. Reporter gene constructs demonstrated the actual usability of two promoters. The investigated PCY1 and PFP promoters were proven to be potent molecular tools for genetic engineering in V. carteri.

The Genetics of Fitness Reorganization during the Transition to Multicellularity: The Volvocine regA-like Family as a Model.

The evolutionary transition from single-celled to multicellular individuality requires organismal fitness to shift from the cell level to a cell group. This reorganization of fitness occurs by re-allocating the two components of fitness, survival and reproduction, between two specialized cell types in the multicellular group: soma and germ, respectively. How does the genetic basis for such fitness reorganization evolve? One possible mechanism is the co-option of life history genes present in the unicellular ancestors of a multicellular lineage. For instance, single-celled organisms must regulate their investment in survival and reproduction in response to environmental changes, particularly decreasing reproduction to ensure survival under stress. Such stress response life history genes can provide the genetic basis for the evolution of cellular differentiation in multicellular lineages. The regA-like gene family in the volvocine green algal lineage provides an excellent model system to study how this co-option can occur. We discuss the origin and evolution of the volvocine regA-like gene family, including regA-the gene that controls somatic cell development in the model organism Volvox carteri. We hypothesize that the co-option of life history trade-off genes is a general mechanism involved in the transition to multicellular individuality, making volvocine algae and the regA-like family a useful template for similar investigations in other lineages.

A personal cost of cheating can stabilize reproductive altruism during the early evolution of clonal multicellularity.

Understanding how cooperation evolved and is maintained remains an important and often controversial topic because cheaters that reap the benefits of cooperation without paying the costs can threaten the evolutionary stability of cooperative traits. Cooperation-and especially reproductive altruism-is particularly relevant to the evolution of multicellularity, as somatic cells give up their reproductive potential in order to contribute to the fitness of the newly emerged multicellular individual. Here, we investigated cheating in a simple multicellular species-the green alga Volvox carteri, in the context of the mechanisms that can stabilize reproductive altruism during the early evolution of clonal multicellularity. We found that the benefits cheater mutants can gain in terms of their own reproduction are pre-empted by a cost in survival due to increased sensitivity to stress. This personal cost of cheating reflects the antagonistic pleiotropic effects that the gene coding for reproductive altruism-regA-has at the cell level. Specifically, the expression of regA in somatic cells results in the suppression of their reproduction potential but also confers them with increased resistance to stress. Since regA evolved from a life-history trade-off gene, we suggest that co-opting trade-off genes into cooperative traits can provide a built-in safety system against cheaters in other clonal multicellular lineages.

Cellular organization in lab-evolved and extant multicellular species obeys a maximum entropy law.

The prevalence of multicellular organisms is due in part to their ability to form complex structures. How cells pack in these structures is a fundamental biophysical issue, underlying their functional properties. However, much remains unknown about how cell packing geometries arise, and how they are affected by random noise during growth - especially absent developmental programs. Here, we quantify the statistics of cellular neighborhoods of two different multicellular eukaryotes: lab-evolved 'snowflake' yeast and the green alga Volvox carteri. We find that despite large differences in cellular organization, the free space associated with individual cells in both organisms closely fits a modified gamma distribution, consistent with maximum entropy predictions originally developed for granular materials. This 'entropic' cellular packing ensures a degree of predictability despite noise, facilitating parent-offspring fidelity even in the absence of developmental regulation. Together with simulations of diverse growth morphologies, these results suggest that gamma-distributed cell neighborhood sizes are a general feature of multicellularity, arising from conserved statistics of cellular packing.

Molecular and cellular dynamics of early embryonic cell divisions in Volvox carteri.

Cell division is fundamental to all organisms and the green alga used here exhibits both key animal and plant functions. Specifically, we analyzed the molecular and cellular dynamics of early embryonic divisions of the multicellular green alga Volvox carteri (Chlamydomonadales). Relevant proteins related to mitosis and cytokinesis were identified in silico, the corresponding genes were cloned, fused to yfp, and stably expressed in Volvox, and the tagged proteins were studied by live-cell imaging. We reveal rearrangements of the microtubule cytoskeleton during centrosome separation, spindle formation, establishment of the phycoplast, and generation of previously unknown structures. The centrosomes participate in initiation of spindle formation and determination of spindle orientation. Although the nuclear envelope does not break down during early mitosis, intermixing of cytoplasm and nucleoplasm results in loss of nuclear identity. Finally, we present a model for mitosis in Volvox. Our study reveals enormous dynamics, clarifies spatio-temporal relationships of subcellular structures, and provides insight into the evolution of cell division.

Cell Type-Specific Pherophorins of Volvox carteri Reveal Interplay of Both Cell Types in ECM Biosynthesis.

The spheroidal green algae Volvox carteri serves as a model system to investigate the formation of a complex, multifunctional extracellular matrix (ECM) in a relatively simple, multicellular organism with cell differentiation. The V. carteri ECM is mainly composed of hydroxyproline-rich glycoproteins (HRGPs) and there are diverse region-specific, anatomically distinct structures in the ECM. One large protein family with importance for ECM biosynthesis stands out: the pherophorins. The few pherophorins previously extracted from the ECM and characterized, were specifically expressed by somatic cells. However, the localization and function of most pherophorins is unknown. Here, we provide a phylogenetic analysis of 153 pherophorins of V. carteri and its unicellular relative Chlamydomonas reinhardtii. Our analysis of cell type-specific mRNA expression of pherophorins in V. carteri revealed that, contrary to previous assumptions, only about half (52%) of the 102 investigated pherophorin-related genes show stronger expression in somatic cells, whereas about one-third (34%) of the genes show significant higher expression in reproductive cells (gonidia). We fused two pherophorin genes that are expressed by different cell types to yfp, stably expressed them in Volvox and studied the tagged proteins by live-cell imaging. In contrast to earlier biochemical approaches, this genetic approach also allows the in vivo analysis of non-extractable, covalently cross-linked ECM proteins. We demonstrate that the soma-specific pherophorin SSG185 is localized in the outermost ECM structures of the spheroid, the boundary zone and at the flagellar hillocks. SSG185:YFP is detectable as early as 1.5 h after completion of embryogenesis. It is then present for the rest of the life cycle. The gonidia-specific pherophorin PhG is localized in the gonidial cellular zone 1 ("gonidial vesicle") suggesting its involvement in the protection of gonidia and developing embryos until hatching. Even if somatic cells produce the main portion of the ECM of the spheroids, ECM components produced by gonidia are also required to cooperatively assemble the total ECM. Our results provide insights into the evolution of the pherophorin protein family and convey a more detailed picture of Volvox ECM synthesis.

Genome sequencing of the multicellular alga Astrephomene provides insights into convergent evolution of germ-soma differentiation.

Germ-soma differentiation evolved independently in many eukaryotic lineages and contributed to complex multicellular organizations. However, the molecular genetic bases of such convergent evolution remain unresolved. Two multicellular volvocine green algae, Volvox and Astrephomene, exhibit convergent evolution of germ-soma differentiation. The complete genome sequence is now available for Volvox, while genome information is scarce for Astrephomene. Here, we generated the de novo whole genome sequence of Astrephomene gubernaculifera and conducted RNA-seq analysis of isolated somatic and reproductive cells. In Volvox, tandem duplication and neofunctionalization of the ancestral transcription factor gene (RLS1/rlsD) might have led to the evolution of regA, the master regulator for Volvox germ-soma differentiation. However, our genome data demonstrated that Astrephomene has not undergone tandem duplication of the RLS1/rlsD homolog or acquisition of a regA-like gene. Our RNA-seq analysis revealed the downregulation of photosynthetic and anabolic gene expression in Astrephomene somatic cells, as in Volvox. Among genes with high expression in somatic cells of Astrephomene, we identified three genes encoding putative transcription factors, which may regulate somatic cell differentiation. Thus, the convergent evolution of germ-soma differentiation in the volvocine algae may have occurred by the acquisition of different regulatory circuits that generate a similar division of labor.

Characterization and Transformation of reg Cluster Genes in Volvox powersii Enable Investigation of Convergent Evolution of Cellular Differentiation in Volvox.

The evolution of germ-soma cellular differentiation represents a key step in the evolution of multicellular individuality. Volvox carteri and its relatives, the volvocine green algae, provide a model system for studying the evolution of cellular differentiation. In V. carteri, the regA gene controls somatic cell differentiation and is found in a group of paralogs called the reg cluster, along with rlsA, rlsB, and rlsC. However, the developmental program of V. carteri is derived compared to other volvocine algae. Here we examine Volvox powersii which possesses an ancestral developmental program and independent evolution of the Volvox body plan. We sequenced the reg cluster from V. powersii wild-type and a mutant with fewer cells and altered germ-soma ratio. We found that the mutant strain's rlsB gene has a deletion predicted to cause a truncated protein product. We developed a genetic transformation procedure to insert wild-type rlsB into the mutant strain. Transformation did not result in phenotypic rescue, suggesting the rlsB mutation is insufficient for generating the mutant phenotype. The transformation techniques and sequences described here provide essential tools to study V. powersii, a species well suited for studying the evolution of cellular differentiation and convergent evolution of Volvox morphology.

Phylotranscriptomics points to multiple independent origins of multicellularity and cellular differentiation in the volvocine algae.

BACKGROUND: The volvocine algae, which include the single-celled species Chlamydomonas reinhardtii and the colonial species Volvox carteri, serve as a model in which to study the evolution of multicellularity and cellular differentiation. Studies reconstructing the history of this group have by and large relied on datasets of one to a few genes for phylogenetic inference and ancestral character state reconstruction. As a result, volvocine phylogenies lack concordance depending on the number and/or type of genes (i.e., chloroplast vs nuclear) chosen for phylogenetic inference. While multiple studies suggest that multicellularity evolved only once in the volvocine algae, that each of its three colonial families is monophyletic, and that there have been at least three independent origins of cellular differentiation in the group, other studies call into question one or more of these conclusions. An accurate assessment of the evolutionary history of the volvocine algae requires inference of a more robust phylogeny. RESULTS: We performed RNA sequencing (RNA-seq) on 55 strains representing 47 volvocine algal species and obtained similar data from curated databases on 13 additional strains. We then compiled a dataset consisting of transcripts for 40 single-copy, protein-coding, nuclear genes and subjected the predicted amino acid sequences of these genes to maximum likelihood, Bayesian inference, and coalescent-based analyses. These analyses show that multicellularity independently evolved at least twice in the volvocine algae and that the colonial family Goniaceae is not monophyletic. Our data further indicate that cellular differentiation arose independently at least four, and possibly as many as six times, within the volvocine algae. CONCLUSIONS: Altogether, our results demonstrate that multicellularity and cellular differentiation are evolutionarily labile in the volvocine algae, affirming the importance of this group as a model system for the study of major transitions in the history of life.

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.

Phototoxicities Caused by Continuous Light Exposure Were Not Induced in Retinal Ganglion Cells Transduced by an Optogenetic Gene.

The death of photoreceptor cells is induced by continuous light exposure. However, it is unclear whether light damage was induced in retinal ganglion cells with photosensitivity by transduction of optogenetic genes. In this study, we evaluated the phototoxicities of continuous light exposure on retinal ganglion cells after transduction of the optogenetic gene mVChR1 using an adeno-associated virus vector. Rats were exposed to continuous light for a week, and visually evoked potentials (VEPs) were recorded. The intensities of continuous light (500, 1000, 3000, and 5000 lx) increased substantially after VEP recordings. After the final recording of VEPs, retinal ganglion cells (RGCs) were retrogradely labeled with a fluorescein tracer, FluoroGold, and the number of retinal ganglion cells was counted under a fluorescent microscope. There was no significant reduction in the amplitudes of VEPs and the number of RGCs after exposure to any light intensity. These results indicated that RGCs were photosensitive after the transduction of optogenetic genes and did not induce any phototoxicity by continuous light exposure.

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.

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.