The Effects of De Novo Mutation on Gene Expression and the Consequences for Fitness in Chlamydomonas reinhardtii

Mutation is the ultimate source of genetic variation, the bedrock of evolution. Yet, predicting the consequences of new mutations remains a challenge in biology. Gene expression provides a potential link between a genotype and its phenotype. But the variation in gene expression created by de novo mutation and the fitness consequences of mutational changes to expression remain relatively unexplored. Here, we investigate the effects of >2,600 de novo mutations on gene expression across the transcriptome of 28 mutation accumulation lines derived from 2 independent wild-type genotypes of the green algae Chlamydomonas reinhardtii. We observed that the amount of genetic variance in gene expression created by mutation (Vm) was similar to the variance that mutation generates in typical polygenic phenotypic traits and approximately 15-fold the variance seen in the limited species where Vm in gene expression has been estimated. Despite the clear effect of mutation on expression, we did not observe a simple additive effect of mutation on expression change, with no linear correlation between the total expression change and mutation count of individual MA lines. We therefore inferred the distribution of expression effects of new mutations to connect the number of mutations to the number of differentially expressed genes (DEGs). Our inferred DEE is highly L-shaped with 95% of mutations causing 0-1 DEG while the remaining 5% are spread over a long tail of large effect mutations that cause multiple genes to change expression. The distribution is consistent with many cis-acting mutation targets that affect the expression of only 1 gene and a large target of trans-acting targets that have the potential to affect tens or hundreds of genes. Further evidence for cis-acting mutations can be seen in the overabundance of mutations in or near differentially expressed genes. Supporting evidence for trans-acting mutations comes from a 15:1 ratio of DEGs to mutations and the clusters of DEGs in the co-expression network, indicative of shared regulatory architecture. Lastly, we show that there is a negative correlation with the extent of expression divergence from the ancestor and fitness, providing direct evidence of the deleterious effects of perturbing gene expression.

Moderate high temperature is beneficial or detrimental depending on carbon availability in the green alga Chlamydomonas reinhardtii

High temperatures impair plant growth and reduce agricultural yields, but the underlying mechanisms remain elusive. The unicellular green alga Chlamydomonas reinhardtii is an excellent model to study heat responses in photosynthetic cells due to its fast growth rate, many similarities in cellular processes to land plants, simple and sequenced genome, and ample genetic and genomics resources. Chlamydomonas grows in light by photosynthesis and with externally supplied acetate as an organic carbon source. Understanding how organic carbon sources affect heat responses is important for the algal industry but remains understudied. We cultivated wild-type Chlamydomonas under highly controlled conditions in photobioreactors at 25 °C (control), 35 °C (moderate high temperature), or 40 °C (acute high temperature) with or without constant acetate supply for 1 or 4 day. Treatment at 35 °C increased algal growth with constant acetate supply but reduced algal growth without sufficient acetate. The overlooked and dynamic effects of 35 °C could be explained by induced acetate uptake and metabolism. Heat treatment at 40 °C for more than 2 day was lethal to algal cultures with or without constant acetate supply. Our findings provide insights to understand algal heat responses and help improve thermotolerance in photosynthetic cells.

Two disulfide-reducing pathways are required for the maturation of plastid c-type cytochromes in Chlamydomonas reinhardtii

In plastids, conversion of light energy into ATP relies on cytochrome f, a key electron carrier with a heme covalently attached to a CXXCH motif. Covalent heme attachment requires reduction of the disulfide-bonded CXXCH by CCS5 and CCS4. CCS5 receives electrons from the oxidoreductase CCDA, while CCS4 is a protein of unknown function. In Chlamydomonas reinhardtii, loss of CCS4 or CCS5 yields a partial cytochrome f assembly defect. Here, we report that the ccs4ccs5 double mutant displays a synthetic photosynthetic defect characterized by a complete loss of holocytochrome f assembly. This defect is chemically corrected by reducing agents, confirming the placement of CCS4 and CCS5 in a reducing pathway. CCS4-like proteins occur in the green lineage, and we show that HCF153, a distant ortholog from Arabidopsis thaliana, can substitute for Chlamydomonas CCS4. Dominant suppressor mutations mapping to the CCS4 gene were identified in photosynthetic revertants of the ccs4ccs5 mutants. The suppressor mutations yield changes in the stroma-facing domain of CCS4 that restore holocytochrome f assembly above the residual levels detected in ccs5. Because the CCDA protein accumulation is decreased specifically in the ccs4 mutant, we hypothesize the suppressor mutations enhance the supply of reducing power through CCDA in the absence of CCS5. We discuss the operation of a CCS5-dependent and a CCS5-independent pathway controlling the redox status of the heme-binding cysteines of apocytochrome f.

Natural variation of nutrient homeostasis among laboratory and field strains of Chlamydomonas reinhardtii

Natural variation among individuals and populations exists in all species, playing key roles in response to environmental stress and adaptation. Micro- and macronutrients have a wide range of functions in photosynthetic organisms, and mineral nutrition thus plays a sizable role in biomass production. To maintain nutrient concentrations inside the cell within physiological limits and prevent the detrimental effects of deficiency or excess, complex homeostatic networks have evolved in photosynthetic cells. The microalga Chlamydomonas reinhardtii (Chlamydomonas) is a unicellular eukaryotic model for studying such mechanisms. In this work, 24 Chlamydomonas strains, comprising field isolates and laboratory strains, were examined for intraspecific differences in nutrient homeostasis. Growth and mineral content were quantified in mixotrophy, as full nutrition control, and compared with autotrophy and nine deficiency conditions for macronutrients (-Ca, -Mg, -N, -P, and -S) and micronutrients (-Cu, -Fe, -Mn, and -Zn). Growth differences among strains were relatively limited. However, similar growth was accompanied by highly divergent mineral accumulation among strains. The expression of nutrient status marker genes and photosynthesis were scored in pairs of contrasting field strains, revealing distinct transcriptional regulation and nutrient requirements. Leveraging this natural variation should enable a better understanding of nutrient homeostasis in Chlamydomonas.

CLPB3 is required for the removal of chloroplast protein aggregates and for thermotolerance in Chlamydomonas

In the cytosol of plant cells, heat-induced protein aggregates are resolved by ClpB/Hsp100 family member HSP101, which is essential for thermotolerance. For chloroplast family member CLPB3 this is less clear with controversial reports on its role in conferring thermotolerance. To shed light onto this issue, we have characterized two clpb3 mutants in Chlamydomonas reinhardtii. We show that chloroplast CLPB3 is required for resolving heat-induced protein aggregates containing stromal TIG1 and the small heat shock proteins HSP22E/F in vivo and for conferring thermotolerance under heat stress. Although CLPB3 accumulates to similarly high levels as stromal HSP70B under ambient conditions, we observed no prominent constitutive phenotypes. However, we found decreased accumulation of the ribosomal subunit PRPL1 and increased accumulation of the stromal protease DEG1C in the clpb3 mutants, suggesting that reduction in chloroplast protein synthesis capacity and increase in proteolytic capacity may compensate for loss of CLPB3 function. Under ambient conditions, CLPB3 was distributed throughout the chloroplast but reorganized into stromal foci upon heat stress, which mostly disappeared during recovery. CLPB3 foci were localized next to HSP22E/F, which accumulated largely to the thylakoid membrane occupied area. This suggests a possible role for CLPB3 in disentangling protein aggregates from the thylakoid membrane system.

Assembly Apparatus of Light-Harvesting Complexes: Identification of Alb3.1-cpSRP-LHCP Complexes in the Green Alga Chlamydomonas reinhardtii

The unicellular green alga, Chlamydomonas reinhardtii, contains many light-harvesting complexes (LHCs) associating chlorophylls a/b and carotenoids; the major LHCIIs (types I, II, III and IV) and minor light-harvesting complexes, CP26 and CP29, for photosystem II, as well as nine LHCIs (LHCA1-9), for photosystem I. A pale green mutant BF4 exhibited impaired accumulation of LHCs due to deficiency in the Alb3.1 gene, which encodes the insertase involved in insertion, folding and assembly of LHC proteins in the thylakoid membranes. To elucidate the molecular mechanism by which ALB3.1 assists LHC assembly, we complemented BF4 to express ALB3.1 fused with no, single or triple Human influenza hemagglutinin (HA) tag at its C-terminus (cAlb3.1, cAlb3.1-HA or cAlb3.1-3HA). The resulting complemented strains accumulated most LHC proteins comparable to wild-type (WT) levels. The affinity purification of Alb3.1-HA and Alb3.1-3HA preparations showed that ALB3.1 interacts with cpSRP43 and cpSRP54 proteins of the chloroplast signal recognition particle (cpSRP) and several LHC proteins; two major LHCII proteins (types I and III), two minor LHCII proteins (CP26 and CP29) and eight LHCI proteins (LHCA1, 2, 3, 4, 5, 6, 8 and 9). Pulse-chase labeling experiments revealed that the newly synthesized major LHCII proteins were transiently bound to the Alb3.1 complex. We propose that Alb3.1 interacts with cpSRP43 and cpSRP54 to form an assembly apparatus for most LHCs in the thylakoid membranes. Interestingly, photosystem I (PSI) proteins were also detected in the Alb3.1 preparations, suggesting that the integration of LHCIs to a PSI core complex to form a PSI-LHCI subcomplex occurs before assembled LHCIs dissociate from the Alb3.1-cpSRP complex.

Fitness Effects of Mutations: An Assessment of PROVEAN Predictions Using Mutation Accumulation Data

Predicting fitness in natural populations is a major challenge in biology. It may be possible to leverage fast-accumulating genomic data sets to infer the fitness effects of mutant alleles, allowing evolutionary questions to be addressed in any organism. In this paper, we investigate the utility of one such tool, called PROVEAN. This program compares a query sequence with existing data to provide an alignment-based score for any protein variant, with scores categorized as neutral or deleterious based on a pre-set threshold. PROVEAN has been used widely in evolutionary studies, for example, to estimate mutation load in natural populations, but has not been formally tested as a predictor of aggregate mutational effects on fitness. Using three large published data sets on the genome sequences of laboratory mutation accumulation lines, we assessed how well PROVEAN predicted the actual fitness patterns observed, relative to other metrics. In most cases, we find that a simple count of the total number of mutant proteins is a better predictor of fitness than the number of proteins with variants scored as deleterious by PROVEAN. We also find that the sum of all mutant protein scores explains variation in fitness better than the number of mutant proteins in one of the data sets. We discuss the implications of these results for studies of populations in the wild.

De Novo Mutation Rate Variation and Its Determinants in Chlamydomonas

De novo mutations are central for evolution, since they provide the raw material for natural selection by regenerating genetic variation. However, studying de novo mutations is challenging and is generally restricted to model species, so we have a limited understanding of the evolution of the mutation rate and spectrum between closely related species. Here, we present a mutation accumulation (MA) experiment to study de novo mutation in the unicellular green alga Chlamydomonas incerta and perform comparative analyses with its closest known relative, Chlamydomonas reinhardtii. Using whole-genome sequencing data, we estimate that the median single nucleotide mutation (SNM) rate in C. incerta is μ = 7.6 × 10-10, and is highly variable between MA lines, ranging from μ = 0.35 × 10-10 to μ = 131.7 × 10-10. The SNM rate is strongly positively correlated with the mutation rate for insertions and deletions between lines (r > 0.97). We infer that the genomic factors associated with variation in the mutation rate are similar to those in C. reinhardtii, allowing for cross-prediction between species. Among these genomic factors, sequence context and complexity are more important than GC content. With the exception of a remarkably high C→T bias, the SNM spectrum differs markedly between the two Chlamydomonas species. Our results suggest that similar genomic and biological characteristics may result in a similar mutation rate in the two species, whereas the SNM spectrum has more freedom to diverge.

Chlamydomonas reinhardtii cellular compartments and their contribution to intracellular calcium signalling

Calcium (Ca2+)-dependent signalling plays a well-characterized role in the response to different environmental stimuli, in both plant and animal cells. In the model organism for green algae, Chlamydomonas reinhardtii, Ca2+ signals were reported to have a crucial role in different physiological processes, such as stress responses, photosynthesis, and flagella functions. Recent reports identified the underlying components of the Ca2+ signalling machinery at the level of specific subcellular compartments and reported in vivo imaging of cytosolic Ca2+ concentration in response to environmental stimuli. The characterization of these Ca2+-related mechanisms and proteins in C. reinhardtii is providing knowledge on how microalgae can perceive and respond to environmental stimuli, but also on how this Ca2+ signalling machinery has evolved. Here, we review current knowledge on the cellular mechanisms underlying the generation, shaping, and decoding of Ca2+ signals in C. reinhardtii, providing an overview of the known and possible molecular players involved in the Ca2+ signalling of its different subcellular compartments. The advanced toolkits recently developed to measure time-resolved Ca2+ signalling in living C. reinhardtii cells are also discussed, suggesting how they can improve the study of the role of Ca2+ signals in the cellular response of microalgae to environmental stimuli.

Evolutionary and biochemical characterization of a Chromochloris zofingiensis MBOAT with wax synthase and diacylglycerol acyltransferase activity

Wax synthase (WS) catalyzes the last step in wax ester biosynthesis in green plants. Two unrelated sub-families of WS, including the bifunctional acyltransferase and plant-like WS have been reported, but the latter is largely uncharacterized in microalgae. Here, we functionally characterized a putative plant-like WS (CzWS1) from the emerging model green microalga Chromochloris zofingiensis. Our results showed that plant-like WS evolved under different selection constraints in plants and microalgae, with positive selection likely contributing to functional divergence. Unlike jojoba with high amounts of wax ester in seeds and a highly active WS enzyme, C. zofingiensis has no detectable wax ester but a high abundance of WS transcripts. Co-expression analysis showed that C. zofingiensis WS has different expression correlation with lipid biosynthetic genes from jojoba, and may have a divergent function. In vitro characterization indicated that CzWS1 had diacylglycerol acyltransferase activity along with WS activity, and overexpression of CzWS1 in yeast and Chlamydomonas reinhardtii affected triacylglycerol accumulation. Moreover, biochemical and bioinformatic analyses revealed the relevance of the C-terminal region of CzWS1 in enzyme function. Taken together, our results indicated a functional divergence of plant-like WS in plants and microalgae, and the importance of its C-terminal region in specialization of enzyme function.

Chlamydomonas POLQ is necessary for CRISPR/Cas9-mediated gene targeting

The use of CRISPR/Cas endonucleases has revolutionized gene editing techniques for research on Chlamydomonas reinhardtii. To better utilize the CRISPR/Cas system, it is essential to develop a more comprehensive understanding of the DNA repair pathways involved in genome editing. In this study, we have analyzed contributions from canonical KU80/KU70-dependent nonhomologous end-joining (cNHEJ) and DNA polymerase theta (POLQ)-mediated end joining on SpCas9-mediated untemplated mutagenesis and homology-directed repair (HDR)/gene inactivation in Chlamydomonas. Using CRISPR/SpCas9 technology, we generated DNA repair-defective mutants ku80, ku70, polQ for gene targeting experiments. Our results show that untemplated repair of SpCas9-induced double strand breaks results in mutation spectra consistent with an involvement of both KU80/KU70 and POLQ. In addition, the inactivation of POLQ was found to negatively affect HDR of the inactivated paromomycin-resistant mut-aphVIII gene when donor single-stranded oligos were used. Nevertheless, mut-aphVIII was still repaired by homologous recombination in these mutants. POLQ inactivation suppressed random integration of transgenes co-transformed with the donor ssDNA. KU80 deficiency did not affect these events but instead was surprisingly found to stimulate HDR/gene inactivation. Our data suggest that in Chlamydomonas, POLQ is the main contributor to CRISPR/Cas-induced HDR and random integration of transgenes, whereas KU80/KU70 potentially plays a secondary role. We expect our results will lead to improvement of genome editing in C. reinhardtii and can be used for future development of algal biotechnology.

The Molecular Function of Plant mTERFs as Key Regulators of Organellar Gene Expression

The protein family of mTERFs (mitochondrial transcription termination factors) was initially studied in mammalian and insect mitochondria before the first Arabidopsis mTERF mutant was characterized. More than 10 years of research on the function of plant mTERFs in the flowering plants Arabidopsis thaliana, Zea mays and the green microalga Chlamydomonas reinhardtii has since highlighted that mTERFs are key regulators of organellar gene expression (OGE) in mitochondria and in chloroplasts. Additional functions to be fulfilled by plant mTERFs (e.g. splicing) and the fact that the expression of two organellar genomes had to be facilitated have led to a massive expansion of the plant mTERF portfolio compared to that found in mammals. Plant mTERFs are implicated in all steps of OGE ranging from the modulation of transcription to the maturation of tRNAs and hence translation. Furthermore, being regulators of OGE, mTERFs are required for a successful long-term acclimation to abiotic stress, retrograde signaling and interorganellar communication. Here, I review the recent progress in the elucidation of molecular mTERF functions.

Characterization of the transient fluorescence wave phenomenon that occurs during H2 production in Chlamydomonas reinhardtii

The redox state of the plastoquinone (PQ) pool in sulfur-deprived, H2-producing Chlamydomonas reinhardtii cells was studied using single flash-induced variable fluorescence decay kinetics. During H2 production, the fluorescence decay kinetics exhibited an unusual post-illumination rise of variable fluorescence, giving a wave-like appearance. The wave showed the transient fluorescence minimum at ~60 ms after the flash, followed by a rise, reaching the transient fluorescence maximum at ~1 s after the flash, before decaying back to the initial fluorescence level. Similar wave-like fluorescence decay kinetics have been reported previously in anaerobically incubated cyanobacteria but not in green algae. From several different electron and proton transfer inhibitors used, polymyxin B, an inhibitor of type II NAD(P)H dehydrogenase (NDA2), had the effect of eliminating the fluorescence wave feature, indicating involvement of NDA2 in this phenomenon. This was further confirmed by the absence of the fluorescence wave in the Δnda2 mutant lacking NDA2. Additionally, Δnda2 mutants have also shown delayed and diminished H2 production (only 23% if compared with the wild type). Our results show that the fluorescence wave phenomenon in C. reinhardtii is observed under highly reducing conditions and is induced by the NDA2-mediated electron flow from the reduced stromal components to the PQ pool. Therefore, the fluorescence wave phenomenon is a sensitive probe for the complex network of redox reactions at the PQ pool level in the thylakoid membrane. It could be used in further characterization and improvement of the electron transfer pathways leading to H2 production in C. reinhardtii.

The pyrenoidal linker protein EPYC1 phase separates with hybrid Arabidopsis-Chlamydomonas Rubisco through interactions with the algal Rubisco small subunit

Photosynthetic efficiencies in plants are restricted by the CO2-fixing enzyme Rubisco but could be enhanced by introducing a CO2-concentrating mechanism (CCM) from green algae, such as Chlamydomonas reinhardtii (hereafter Chlamydomonas). A key feature of the algal CCM is aggregation of Rubisco in the pyrenoid, a liquid-like organelle in the chloroplast. Here we have used a yeast two-hybrid system and higher plants to investigate the protein-protein interaction between Rubisco and essential pyrenoid component 1 (EPYC1), a linker protein required for Rubisco aggregation. We showed that EPYC1 interacts with the small subunit of Rubisco (SSU) from Chlamydomonas and that EPYC1 has at least five SSU interaction sites. Interaction is crucially dependent on the two surface-exposed α-helices of the Chlamydomonas SSU. EPYC1 could be localized to the chloroplast in higher plants and was not detrimental to growth when expressed stably in Arabidopsis with or without a Chlamydomonas SSU. Although EPYC1 interacted with Rubisco in planta, EPYC1 was a target for proteolytic degradation. Plants expressing EPYC1 did not show obvious evidence of Rubisco aggregation. Nevertheless, hybrid Arabidopsis Rubisco containing the Chlamydomonas SSU could phase separate into liquid droplets with purified EPYC1 in vitro, providing the first evidence of pyrenoid-like aggregation for Rubisco derived from a higher plant.

LIP4 Is Involved in Triacylglycerol Degradation in Chlamydomonas reinhardtii

Degradation of the storage compound triacylglycerol (TAG) is a crucial process in response to environmental stimuli. In microalgae, this process is important for re-growth when conditions become favorable after cells have experienced stresses. Mobilization of TAG is initiated by actions of lipases causing the release of glycerol and free fatty acids, which can be further broken down for energy production or recycled to synthesize membrane lipids. Although key enzymes in the process, TAG lipases remain to be characterized in the model green alga Chlamydomonas reinhardtii. Here, we describe the functional analysis of a putative TAG lipase, i.e. LIP4, which shares 44% amino acid identity with the major TAG lipase in Arabidopsis (SUGAR DEPENDENT1-SDP1). The LIP4 transcript level was downregulated during nitrogen deprivation when TAG accumulates, but was upregulated during nitrogen resupply (NR) when TAG was degraded. Both artificial microRNA and insertional mutants showed a delay in TAG mobilization during NR. The difference in TAG degradation was more pronounced when the cultures were incubated without acetate in the dark. Furthermore, the lip4 insertional mutants over-accumulated TAG during optimal growth conditions. Taken together, the results suggest to us that LIP4 likely acts as a TAG lipase and plays a role in TAG homeostasis in Chlamydomonas.

Molecular Genetic Tools and Emerging Synthetic Biology Strategies to Increase Cellular Oil Content in Chlamydomonas reinhardtii

Microalgae constitute a highly diverse group of eukaryotic and photosynthetic microorganisms that have developed extremely efficient systems for harvesting and transforming solar energy into energy-rich molecules such as lipids. Although microalgae are considered to be one of the most promising platforms for the sustainable production of liquid oil, the oil content of these organisms is naturally low, and algal oil production is currently not economically viable. Chlamydomonas reinhardtii (Chlamydomonas) is an established algal model due to its fast growth, high transformation efficiency, and well-understood physiology and to the availability of detailed genome information and versatile molecular tools for this organism. In this review, we summarize recent advances in the development of genetic manipulation tools for Chlamydomonas, from gene delivery methods to state-of-the-art genome-editing technologies and fluorescent dye-based high-throughput mutant screening approaches. Furthermore, we discuss practical strategies and toolkits that enhance transgene expression, such as choice of expression vector and background strain. We then provide examples of how advanced genetic tools have been used to increase oil content in Chlamydomonas. Collectively, the current literature indicates that microalgal oil content can be increased by overexpressing key enzymes that catalyze lipid biosynthesis, blocking lipid degradation, silencing metabolic pathways that compete with lipid biosynthesis and modulating redox state. The tools and knowledge generated through metabolic engineering studies should pave the way for developing a synthetic biological approach to enhance lipid productivity in microalgae.