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.

ABSTRACT

The unicellular green alga Chlamydomonas reinhardtii is a choice reference system for the study of photosynthesis, 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 from manually-curated and community-generated gene lists. We extracted 400 high-confidence cilia-related genes at the intersection of multiple co-expressed lists, illustrating the power of our simple method. Surprisingly, Chlamydomonas experiments did not cluster according to an obvious pattern, suggesting an underappreciated variable during sample collection. One possible source of variation may stem from the strong clustering of nuclear genes as a function of their diurnal phase, even in samples collected in constant conditions, indicating substantial residual synchronization in batch cultures. We provide a step-by-step guide into the analysis of co-expression across Chlamydomonas transcriptome datasets to help foster gene function discovery.

One-sentence summary

we reveal co-expression potential between Chlamydomonas genes and describe strong synchronization of cells grown in batch cultures as a possible source of underappreciated variation.

The photosynthetic unicellular alga Chlamydomonas (Chlamydomonas reinhardtii) is a versatile reference for algal biology because of its ease of culture in the laboratory. Genomic and systems biology approaches have previously described transcriptome responses to environmental changes using bulk data, thus representing the average behavior from pools of cells. Here, we apply single-cell RNA sequencing (scRNA-seq) to probe the heterogeneity of Chlamydomonas cell populations under three environments and in two genotypes differing by the presence of a cell wall. First, we determined that RNA can be extracted from single algal cells with or without a cell wall, offering the possibility to sample natural algal communities. Second, scRNA-seq successfully separated single cells into nonoverlapping cell clusters according to their growth conditions. Cells exposed to iron or nitrogen deficiency were easily distinguished despite a shared tendency to arrest photosynthesis and cell division to economize resources. Notably, these groups of cells not only recapitulated known patterns observed with bulk RNA-seq but also revealed their inherent heterogeneity. A substantial source of variation between cells originated from their endogenous diurnal phase, although cultures were grown in constant light. We exploited this result to show that circadian iron responses may be conserved from algae to land plants. We document experimentally that bulk RNA-seq data represent an average of typically hidden heterogeneity in the population.

Nuclear pore complexes (NPCs) span the nuclear envelope and mediate nucleocytoplasmic exchange. They are a hallmark of eukaryotes and deeply rooted in the evolutionary origin of cellular compartmentalization. NPCs have an elaborate architecture that has been well studied in vertebrates. Whether this architecture is unique or varies significantly in other eukaryotic kingdoms remains unknown, predominantly due to missing in situ structural data. Here, we report the architecture of the algal NPC from the early branching eukaryote Chlamydomonas reinhardtii and compare it to the human NPC. We find that the inner ring of the Chlamydomonas NPC has an unexpectedly large diameter, and the outer rings exhibit an asymmetric oligomeric state that has not been observed or predicted previously. Our study provides evidence that the NPC is subject to substantial structural variation between species. The divergent and conserved features of NPC architecture provide insights into the evolution of the nucleocytoplasmic transport machinery.

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.

In photosynthetic eukaryotes, thousands of proteins are translated in the cytosol and imported into the chloroplast through the concerted action of two translocons-termed TOC and TIC-located in the outer and inner membranes of the chloroplast envelope, respectively. The degree to which the molecular composition of the TOC and TIC complexes is conserved over phylogenetic distances has remained controversial. Here, we combine transcriptomic, biochemical, and genetic tools in the green alga Chlamydomonas (Chlamydomonas reinhardtii) to demonstrate that, despite a lack of evident sequence conservation for some of its components, the algal TIC complex mirrors the molecular composition of a TIC complex from Arabidopsis thaliana. The Chlamydomonas TIC complex contains three nuclear-encoded subunits, Tic20, Tic56, and Tic100, and one chloroplast-encoded subunit, Tic214, and interacts with the TOC complex, as well as with several uncharacterized proteins to form a stable supercomplex (TIC-TOC), indicating that protein import across both envelope membranes is mechanistically coupled. Expression of the nuclear and chloroplast genes encoding both known and uncharacterized TIC-TOC components is highly coordinated, suggesting that a mechanism for regulating its biogenesis across compartmental boundaries must exist. Conditional repression of Tic214, the only chloroplast-encoded subunit in the TIC-TOC complex, impairs the import of chloroplast proteins with essential roles in chloroplast ribosome biogenesis and protein folding and induces a pleiotropic stress response, including several proteins involved in the chloroplast unfolded protein response. These findings underscore the functional importance of the TIC-TOC supercomplex in maintaining chloroplast proteostasis.

The unicellular green alga Chlamydomonas reinhardtii displays metabolic flexibility in response to a changing environment. We analyzed expression patterns of its three genomes in cells grown under light-dark cycles. Nearly 85% of transcribed genes show differential expression, with different sets of transcripts being up-regulated over the course of the day to coordinate cellular growth before undergoing cell division. Parallel measurements of select metabolites and pigments, physiological parameters, and a subset of proteins allow us to infer metabolic events and to evaluate the impact of the transcriptome on the proteome. Among the findings are the observations that Chlamydomonas exhibits lower respiratory activity at night compared with the day; multiple fermentation pathways, some oxygen-sensitive, are expressed at night in aerated cultures; we propose that the ferredoxin, FDX9, is potentially the electron donor to hydrogenases. The light stress-responsive genes PSBS, LHCSR1, and LHCSR3 show an acute response to lights-on at dawn under abrupt dark-to-light transitions, while LHCSR3 genes also exhibit a later, second burst in expression in the middle of the day dependent on light intensity. Each response to light (acute and sustained) can be selectively activated under specific conditions. Our expression dataset, complemented with coexpression networks and metabolite profiling, should constitute an excellent resource for the algal and plant communities.

Photosynthetic organisms are essential for human life, yet most of their genes remain functionally uncharacterized. Single-celled photosynthetic model systems have the potential to accelerate our ability to connect genes to functions. 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. 78% 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 allow us to place known and previously 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 novel phenotypes and gene functions, including the discovery of three novel components of a defense pathway that counteracts actin cytoskeleton inhibitors released by other organisms. The data also inform phenotype discovery in land plants: mutants in Arabidopsis thaliana genes exhibit similar phenotypes 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.

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 CO₂-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.

Exposing cells to excess metal concentrations well beyond the cellular quota is a powerful tool for understanding the molecular mechanisms of metal homeostasis. Such improved understanding may enable bioengineering of organisms with improved nutrition and bioremediation capacity. We report here that Chlamydomonas reinhardtii can accumulate manganese (Mn) in proportion to extracellular supply, up to 30-fold greater than its typical quota and with remarkable tolerance. As visualized by X-ray fluorescence microscopy and nanoscale secondary ion MS (nanoSIMS), Mn largely co-localizes with phosphorus (P) and calcium (Ca), consistent with the Mn-accumulating site being an acidic vacuole, known as the acidocalcisome. Vacuolar Mn stores are accessible reserves that can be mobilized in Mn-deficient conditions to support algal growth. We noted that Mn accumulation depends on cellular polyphosphate (polyP) content, indicated by 1) a consistent failure of C. reinhardtii vtc1 mutant strains, which are deficient in polyphosphate synthesis, to accumulate Mn and 2) a drastic reduction of the Mn storage capacity in P-deficient cells. Rather surprisingly, X-ray absorption spectroscopy, EPR, and electron nuclear double resonance revealed that only little Mn²⁺ is stably complexed with polyP, indicating that polyP is not the final Mn ligand. We propose that polyPs are a critical component of Mn accumulation in Chlamydomonas by driving Mn relocation from the cytosol to acidocalcisomes. Within these structures, polyP may, in turn, escort vacuolar Mn to a number of storage ligands, including phosphate and phytate, and other, yet unidentified, compounds.