Dynamic Changes in the Transcriptome and Methylome of Chlamydomonas reinhardtii throughout Its Life Cycle

Protocol to isolate nuclei from Chlamydomonas reinhardtii for ATAC sequencing

The isolation of sufficient amounts of intact nuclei is essential to obtain high-resolution maps of chromatin accessibility via assay for transposase-accessible chromatin using sequencing (ATAC-seq). Here, we present a protocol for tag-free isolation of nuclei from both cell walled and cell wall-deficient strains of the green model alga Chlamydomonas reinhardtii at a suitable quality for ATAC-seq. We describe steps for nuclei isolation, quantification, and downstream ATAC-seq. This protocol is optimized to shorten the time of isolation and quantification of nuclei.

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

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

A gap-free genome assembly of Chlamydomonas reinhardtii and detection of translocations induced by CRISPR-mediated mutagenesis

Genomic assemblies of the unicellular green alga Chlamydomonas reinhardtii have provided important resources for researchers. However, assembly errors, large gaps, and unplaced scaffolds as well as strain-specific variants currently impede many types of analysis. By combining PacBio HiFi and Oxford Nanopore long-read technologies, we generated a de novo genome assembly for strain CC-5816, derived from crosses of strains CC-125 and CC-124. Multiple methods of evaluating genome completeness and base-pair error rate suggest that the final telomere-to-telomere assembly is highly accurate. The CC-5816 assembly enabled previously difficult analyses that include characterization of the 17 centromeres, rDNA arrays on three chromosomes, and 56 insertions of organellar DNA into the nuclear genome. Using Nanopore sequencing, we identified sites of cytosine (CpG) methylation, which are enriched at centromeres. We analyzed CRISPR-Cas9 insertional mutants in the PF23 gene. Two of the three alleles produced progeny that displayed patterns of meiotic inviability that suggested the presence of a chromosomal aberration. Mapping Nanopore reads from pf23-2 and pf23-3 onto the CC-5816 genome showed that these two strains each carry a translocation that was initiated at the PF23 gene locus on chromosome 11 and joined with chromosomes 5 or 3, respectively. The translocations were verified by demonstrating linkage between loci on the two translocated chromosomes in meiotic progeny. The three pf23 alleles display the expected short-cilia phenotype, and immunoblotting showed that pf23-2 lacks the PF23 protein. Our CC-5816 genome assembly will undoubtedly provide an important tool for the Chlamydomonas research community.

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.

DNA Methylation in Algae and Its Impact on Abiotic Stress Responses

Epigenetics, referring to heritable gene regulatory information that is independent of changes in DNA sequences, is an important mechanism involved both in organism development and in the response to environmental events. About the epigenetic marks, DNA methylation is one of the most conserved mechanisms, playing a pivotal role in organism response to several biotic and abiotic stressors. Indeed, stress can induce changes in gene expression through hypo- or hyper-methylation of DNA at specific loci and/or in DNA methylation at the genome-wide level, which has an adaptive significance and can direct genome evolution. Exploring DNA methylation in responses to abiotic stress could have important implications for improving stress tolerance in algae. This article summarises the DNA methylation pattern in algae and its impact on abiotic stress, such as heavy metals, nutrients and temperature. Our discussion provides information for further research in algae for a better comprehension of the epigenetic response under abiotic stress, which could favour important implications to sustain algae growth under abiotic stress conditions, often related to high biosynthesis of interesting metabolites.

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.

The epigenetic origin of life history transitions in plants and algae

Plants and algae have a complex life history that transitions between distinct life forms called the sporophyte and the gametophyte. This phenomenon-called the alternation of generations-has fascinated botanists and phycologists for over 170 years. Despite the mesmerizing array of life histories described in plants and algae, we are only now beginning to learn about the molecular mechanisms controlling them and how they evolved. Epigenetic silencing plays an essential role in regulating gene expression during multicellular development in eukaryotes, raising questions about its impact on the life history strategy of plants and algae. Here, we trace the origin and function of epigenetic mechanisms across the plant kingdom, from unicellular green algae through to angiosperms, and attempt to reconstruct the evolutionary steps that influenced life history transitions during plant evolution. Central to this evolutionary scenario is the adaption of epigenetic silencing from a mechanism of genome defense to the repression and control of alternating generations. We extend our discussion beyond the green lineage and highlight the peculiar case of the brown algae. Unlike their unicellular diatom relatives, brown algae lack epigenetic silencing pathways common to animals and plants yet display complex life histories, hinting at the emergence of novel life history controls during stramenopile evolution.

Deep evolutionary origin of gamete-directed zygote activation by KNOX/BELL transcription factors in green plants

KNOX and BELL transcription factors regulate distinct steps of diploid development in plants. In the green alga Chlamydomonas reinhardtii, KNOX and BELL proteins are inherited by gametes of the opposite mating types and heterodimerize in zygotes to activate diploid development. By contrast, in land plants such as Physcomitrium patens and Arabidopsis thaliana, KNOX and BELL proteins function in sporophyte and spore formation, meristem maintenance and organogenesis during the later stages of diploid development. However, whether the contrasting functions of KNOX and BELL were acquired independently in algae and land plants is currently unknown. Here, we show that in the basal land plant species Marchantia polymorpha, gamete-expressed KNOX and BELL are required to initiate zygotic development by promoting nuclear fusion in a manner strikingly similar to that in C. reinhardtii. Our results indicate that zygote activation is the ancestral role of KNOX/BELL transcription factors, which shifted toward meristem maintenance as land plants evolved.

Genome-Wide Mapping of Cytosine Methylation Revealed Dynamic DNA Methylation Patterns Associated with Sporophyte Development of Saccharina japonica

Cytosine methylation plays vital roles in regulating gene expression and plant development. However, the function of DNA methylation in the development of macroalgae remains unclear. Through the genome-wide bisulfite sequencing of cytosine methylation in holdfast, stipe and blade, we obtained the complete 5-mC methylation landscape of Saccharina japonica sporophyte. Our results revealed that the total DNA methylation level of sporophyte was less than 0.9%, and the content of CHH contexts was dominant. Moreover, the distribution of CHH methylation within the genes exhibited exon-enriched characteristics. Profiling of DNA methylation in three parts revealed the diverse methylation pattern of sporophyte development. These pivotal DMRs were involved in cell motility, cell cycle and cell wall/membrane biogenesis. In comparison with stipe and blade, hypermethylation of mannuronate C5-epimerase in holdfast decreased the transcript abundance, which affected the synthesis of alginate, the key component of cell walls. Additionally, 5-mC modification participated in the regulation of blade and holdfast development by the glutamate content respectively via glutamine synthetase and amidophosphoribosyl transferase, which may act as the epigenetic regulation signal. Overall, our study revealed the global methylation characteristics of the well-defined holdfast, stipe and blade, and provided evidence for epigenetic regulation of sporophyte development in brown macroalgae.

Architecture and evolution of subtelomeres in the unicellular green alga Chlamydomonas reinhardtii

In most eukaryotes, subtelomeres are dynamic genomic regions populated by multi-copy sequences of different origins, which can promote segmental duplications and chromosomal rearrangements. However, their repetitive nature has complicated the efforts to sequence them, analyse their structure and infer how they evolved. Here, we use recent genome assemblies of Chlamydomonas reinhardtii based on long-read sequencing to comprehensively describe the subtelomere architecture of the 17 chromosomes of this model unicellular green alga. We identify three main repeated elements present at subtelomeres, which we call Sultan, Subtile and Suber, alongside three chromosome extremities with ribosomal DNA as the only identified component of their subtelomeres. The most common architecture, present in 27 out of 34 subtelomeres, is a heterochromatic array of Sultan elements adjacent to the telomere, followed by a transcribed Spacer sequence, a G-rich microsatellite and transposable elements. Sequence similarity analyses suggest that Sultan elements underwent segmental duplications within each subtelomere and rearranged between subtelomeres at a much lower frequency. Analysis of other green algae reveals species-specific repeated elements that are shared across subtelomeres, with an overall organization similar to C. reinhardtii. This work uncovers the complexity and evolution of subtelomere architecture in green algae.

Architecture and evolution of subtelomeres in the unicellular green alga Chlamydomonas reinhardtii

In most eukaryotes, subtelomeres are dynamic genomic regions populated by multi-copy sequences of different origins, which can promote segmental duplications and chromosomal rearrangements. However, their repetitive nature has complicated the efforts to sequence them, analyse their structure and infer how they evolved. Here, we use recent genome assemblies of Chlamydomonas reinhardtii based on long-read sequencing to comprehensively describe the subtelomere architecture of the 17 chromosomes of this model unicellular green alga. We identify three main repeated elements present at subtelomeres, which we call Sultan, Subtile and Suber, alongside three chromosome extremities with ribosomal DNA as the only identified component of their subtelomeres. The most common architecture, present in 27 out of 34 subtelomeres, is a heterochromatic array of Sultan elements adjacent to the telomere, followed by a transcribed Spacer sequence, a G-rich microsatellite and transposable elements. Sequence similarity analyses suggest that Sultan elements underwent segmental duplications within each subtelomere and rearranged between subtelomeres at a much lower frequency. Analysis of other green algae reveals species-specific repeated elements that are shared across subtelomeres, with an overall organization similar to C. reinhardtii. This work uncovers the complexity and evolution of subtelomere architecture in green algae.

Single-cell RNA sequencing of batch Chlamydomonas cultures reveals heterogeneity in their diurnal cycle phase

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.

Integrative analysis of chloroplast DNA methylation in a marine alga-Saccharina japonica

KEY MESSAGE

We applied an integrative approach using multiple methods to verify cytosine methylation in the chloroplast DNA of the multicellular brown alga Saccharina japonica. Cytosine DNA methylation is a heritable process which plays important roles in regulating development throughout the life cycle of an organism. Although methylation of nuclear DNA has been studied extensively, little is known about the state and role of DNA methylation in chloroplast genomes, especially in marine algae. Here, we have applied an integrated approach encompassing whole-genome bisulfite sequencing, methylated DNA immunoprecipitation, gene co-expression networks and photophysiological analyses to provide evidence for the role of chloroplast DNA methylation in a marine alga, the multicellular brown alga Saccharina japonica. Although the overall methylation level was relatively low in the chloroplast genome of S. japonica, gametophytes exhibited higher methylation levels than sporophytes. Gene-specific bisulfite-cloning sequencing provided additional evidence for the methylation of key photosynthetic genes. Many of them were highly expressed in sporophytes whereas genes involved in transcription, translation and biosynthesis were strongly expressed in gametophytes. Nucleus-encoded photosynthesis genes were co-expressed with their chloroplast-encoded counterparts potentially contributing to the higher photosynthetic performance in sporophytes compared to gametophytes where these co-expression networks were less pronounced. A nucleus-encoded DNA methyltransferase of the DNMT2 family is assumed to be responsible for the methylation of the chloroplast genome because it is predicted to possess a plastid transit peptide.

The small RNA locus map for Chlamydomonas reinhardtii

Small (s)RNAs play crucial roles in the regulation of gene expression and genome stability across eukaryotes where they direct epigenetic modifications, post-transcriptional gene silencing, and defense against both endogenous and exogenous viruses. It is known that Chlamydomonas reinhardtii, a well-studied unicellular green algae species, possesses sRNA-based mechanisms that are distinct from those of land plants. However, definition of sRNA loci and further systematic classification is not yet available for this or any other algae. Here, using data-driven machine learning approaches including Multiple Correspondence Analysis (MCA) and clustering, we have generated a comprehensively annotated and classified sRNA locus map for C. reinhardtii. This map shows some common characteristics with higher plants and animals, but it also reveals distinct features. These results are consistent with the idea that there was diversification in sRNA mechanisms after the evolutionary divergence of algae from higher plant lineages.

Inhibition of DNA Methylation in Picochlorum soloecismus Alters Algae Productivity

Eukaryotic organisms regulate the organization, structure, and accessibility of their genomes through chromatin remodeling that can be inherited as epigenetic modifications. These DNA and histone protein modifications are ultimately responsible for an organism's molecular adaptation to the environment, resulting in distinctive phenotypes. Epigenetic manipulation of algae holds yet untapped potential for the optimization of biofuel production and bioproduct formation; however, epigenetic machinery and modes-of-action have not been well characterized in algae. We sought to determine the extent to which the biofuel platform species Picochlorum soloecismus utilizes DNA methylation to regulate its genome. We found candidate genes with domains for DNA methylation in the P. soloecismus genome. Whole-genome bisulfite sequencing revealed DNA methylation in all three cytosine contexts (CpG, CHH, and CHG). While global DNA methylation is low overall (∼1.15%), it occurs in appreciable quantities (12.1%) in CpG dinucleotides in a bimodal distribution in all genomic contexts, though terminators contain the greatest number of CpG sites per kilobase. The P. soloecismus genome becomes hypomethylated during the growth cycle in response to nitrogen starvation. Algae cultures were treated daily across the growth cycle with 20 μM 5-aza-2'-deoxycytidine (5AZA) to inhibit propagation of DNA methylation in daughter cells. 5AZA treatment significantly increased optical density and forward and side scatter of cells across the growth cycle (16 days). This increase in cell size and complexity correlated with a significant increase (∼66%) in lipid accumulation. Site specific CpG DNA methylation was significantly altered with 5AZA treatment over the time course, though nitrogen starvation itself induced significant hypomethylation in CpG contexts. Genes involved in several biological processes, including fatty acid synthesis, had altered methylation ratios in response to 5AZA; we hypothesize that these changes are potentially responsible for the phenotype of early induction of carbon storage as lipids. This is the first report to utilize epigenetic manipulation strategies to alter algal physiology and phenotype. Collectively, these data suggest these strategies can be utilized to fine-tune metabolic responses, alter growth, and enhance environmental adaption of microalgae for desired outcomes.

Responses of a Newly Evolved Auxotroph of Chlamydomonas to B12 Deprivation

The corrinoid B₁₂ is synthesized only by prokaryotes yet is widely required by eukaryotes as an enzyme cofactor. Microalgae have evolved B₁₂ dependence on multiple occasions, and we previously demonstrated that experimental evolution of the non-B₁₂-requiring alga Chlamydomonas reinhardtii in media supplemented with B₁₂ generated a B₁₂-dependent mutant (hereafter metE7). This clone provides a unique opportunity to study the physiology of a nascent B₁₂ auxotroph. Our analyses demonstrate that B₁₂ deprivation of metE7 disrupts C1 metabolism, causes an accumulation of starch and triacylglycerides, and leads to a decrease in photosynthetic pigments, proteins, and free amino acids. B₁₂ deprivation also caused a substantial increase in reactive oxygen species, which preceded rapid cell death. Survival could be improved without compromising growth by simultaneously depriving the cells of nitrogen, suggesting a type of cross protection. Significantly, we found further improvements in survival under B₁₂ limitation and an increase in B₁₂ use efficiency after metE7 underwent a further period of experimental evolution, this time in coculture with a B₁₂-producing bacterium. Therefore, although an early B₁₂-dependent alga would likely be poorly adapted to coping with B₁₂ deprivation, association with B₁₂-producers can ensure long-term survival whilst also providing a suitable environment for evolving mechanisms to tolerate B₁₂ limitation better.

Diatom Molecular Research Comes of Age: Model Species for Studying Phytoplankton Biology and Diversity

Diatoms are the world's most diverse group of algae, comprising at least 100,000 species. Contributing ∼20% of annual global carbon fixation, they underpin major aquatic food webs and drive global biogeochemical cycles. Over the past two decades, Thalassiosira pseudonana and Phaeodactylum tricornutum have become the most important model systems for diatom molecular research, ranging from cell biology to ecophysiology, due to their rapid growth rates, small genomes, and the cumulative wealth of associated genetic resources. To explore the evolutionary divergence of diatoms, additional model species are emerging, such as Fragilariopsis cylindrus and Pseudo-nitzschia multistriata Here, we describe how functional genomics and reverse genetics have contributed to our understanding of this important class of microalgae in the context of evolution, cell biology, and metabolic adaptations. Our review will also highlight promising areas of investigation into the diversity of these photosynthetic organisms, including the discovery of new molecular pathways governing the life of secondary plastid-bearing organisms in aquatic environments.

Single-base methylome profiling of the giant kelp Saccharina japonica reveals significant differences in DNA methylation to microalgae and plants

Brown algae have convergently evolved plant-like body plans and reproductive cycles, which in plants are controlled by differential DNA methylation. This contribution provides the first single-base methylome profiles of haploid gametophytes and diploid sporophytes of a multicellular alga. Although only c. 1.4% of cytosines in Saccharina japonica were methylated mainly at CHH sites and characterized by 5-methylcytosine (5mC), there were significant differences between life-cycle stages. DNA methyltransferase 2 (DNMT2), known to efficiently catalyze tRNA methylation, is assumed to methylate the genome of S. japonica in the structural context of tRNAs as the genome does not encode any other DNA methyltransferases. Circular and long noncoding RNA genes were the most strongly methylated regulatory elements in S. japonica. Differential expression of genes was negatively correlated with DNA methylation with the highest methylation levels measured in both haploid gametophytes. Hypomethylated and highly expressed genes in diploid sporophytes included genes involved in morphogenesis and halogen metabolism. The data herein provide evidence that cytosine methylation, although occurring at a low level, is significantly contributing to the formation of different life-cycle stages, tissue differentiation and metabolism in brown algae.

When Unity Is Strength: The Strategies Used by Chlamydomonas to Survive Environmental Stresses

The unicellular green alga Chlamydomonas reinhardtii is a valuable model system to study a wide spectrum of scientific fields, including responses to environmental conditions. Most studies are performed under optimal growth conditions or under mild stress. However, when environmental conditions become harsher, the behavior of this unicellular alga is less well known. In this review we will show that despite being a unicellular organism, Chlamydomonas can survive very severe environmental conditions. To do so, and depending on the intensity of the stress, the strategies used by Chlamydomonas can range from acclimation to the formation of multicellular structures, or involve programmed cell death.

A Series of Fortunate Events: Introducing Chlamydomonas as a Reference Organism

The unicellular alga Chlamydomonas reinhardtii is a classical reference organism for studying photosynthesis, chloroplast biology, cell cycle control, and cilia structure and function. It is also an emerging model for studying sensory cilia, the production of high-value bioproducts, and in situ structural determination. Much of the early appeal of Chlamydomonas was rooted in its promise as a genetic system, but like other classic model organisms, this rise to prominence predated the discovery of the structure of DNA, whole-genome sequences, and molecular techniques for gene manipulation. The haploid genome of C. reinhardtii facilitates genetic analyses and offers many of the advantages of microbial systems applied to a photosynthetic organism. C. reinhardtii has contributed to our understanding of chloroplast-based photosynthesis and cilia biology. Despite pervasive transgene silencing, technological advances have allowed researchers to address outstanding lines of inquiry in algal research. The most thoroughly studied unicellular alga, C. reinhardtii, is the current standard for algal research, and although genome editing is still far from efficient and routine, it nevertheless serves as a template for other algae. We present a historical retrospective of the rise of C. reinhardtii to illuminate its past and present. We also present resources for current and future scientists who may wish to expand their studies to the realm of microalgae.

A vitamin-C-derived DNA modification catalysed by an algal TET homologue

Methylation of cytosine to 5-methylcytosine (5mC) is a prevalent DNA modification found in many organisms. Sequential oxidation of 5mC by ten-eleven translocation (TET) dioxygenases results in a cascade of additional epigenetic marks and promotes demethylation of DNA in mammals1,2. However, the enzymatic activity and function of TET homologues in other eukaryotes remains largely unexplored. Here we show that the green alga Chlamydomonas reinhardtii contains a 5mC-modifying enzyme (CMD1) that is a TET homologue and catalyses the conjugation of a glyceryl moiety to the methyl group of 5mC through a carbon-carbon bond, resulting in two stereoisomeric nucleobase products. The catalytic activity of CMD1 requires Fe(II) and the integrity of its binding motif His-X-Asp, which is conserved in Fe-dependent dioxygenases3. However, unlike previously described TET enzymes, which use 2-oxoglutarate as a co-substrate4, CMD1 uses L-ascorbic acid (vitamin C) as an essential co-substrate. Vitamin C donates the glyceryl moiety to 5mC with concurrent formation of glyoxylic acid and CO2. The vitamin-C-derived DNA modification is present in the genome of wild-type C. reinhardtii but at a substantially lower level in a CMD1 mutant strain. The fitness of CMD1 mutant cells during exposure to high light levels is reduced. LHCSR3, a gene that is critical for the protection of C. reinhardtii from photo-oxidative damage under high light conditions, is hypermethylated and downregulated in CMD1 mutant cells compared to wild-type cells, causing a reduced capacity for photoprotective non-photochemical quenching. Our study thus identifies a eukaryotic DNA base modification that is catalysed by a divergent TET homologue and unexpectedly derived from vitamin C, and describes its role as a potential epigenetic mark that may counteract DNA methylation in the regulation of photosynthesis.

The dynamin-like protein Fzl promotes thylakoid fusion and resistance to light stress in Chlamydomonas reinhardtii

Large GTPases of the Dynamin Related Proteins (DRP) family shape lipid bilayers through membrane fission or fusion processes. Despite the highly organized photosynthetic membranes of thylakoids, a single DRP is known to be targeted inside the chloroplast. Fzl from the land plant Arabidopsis thaliana is inserted in the inner envelope and thylakoid membranes to regulate their morphology. Fzl may promote the fusion of thylakoids but this remains to be proven. Moreover, the physiological requirement for fusing thylakoids is currently unknown. Here, we find that the unicellular microalga Chlamydomonas reinhardtii encodes an Fzl ortholog (CrFzl) that is localized in the chloroplast where it is soluble. To explore its function, the CRISPR/Cas9 technology was employed to generate multiple CrFzl knock out strains. Phenotypic analyzes revealed a specific requirement of CrFzl for survival upon light stress. Consistent with this, strong irradiance lead to increased photoinhibition of photosynthesis in mutant cells. Fluorescence and electron microscopy analysis demonstrated that upon exposure to high light, CrFzl mutants show defects in chloroplast morphology but also large cytosolic vacuoles in close contact with the plastid. We further observe that strong irradiance induces an increased recruitment of the DRP to thylakoid membranes. Most importantly, we show that CrFzl is required for the fusion of thylakoids during mating. Together, our results suggest that thylakoids fusion may be necessary for resistance to light stress.

All eukaryotic cells are composed of compartments with defined functions. Among those, mitochondria generate the main source of energy in human and animal cells. Their capacity to generate and diffuse energy in the cell is regulated by fusion and fragmentation processes.

Together with mitochondria that produce energy from oxygen, plant cells include an additional compartment called the chloroplast that produces energy from light. The machinery that converts light into energy is more precisely located inside the chloroplast within stacks of membranes called the thylakoids. Here, we elucidate the function of CrFzl, a previously uncharacterized protein encoded by the genome of the unicellular alga Chlamydomonas reinhardtii. Algal cells that do not contain CrFzl are impaired in their capacity to grow when they receive too much light and our results indicate that CrFzl promotes the fusion of thylakoids during mating. These results suggest that membrane fusion is an essential tool for energy production in stress conditions by living organisms.

From molecular manipulation of domesticated Chlamydomonas reinhardtii to survival in nature

In the mid-20th century, the unicellular and genetically tractable green alga Chlamydomonas reinhardtii was first developed as a model organism to elucidate fundamental cellular processes such as photosynthesis, light perception and the structure, function and biogenesis of cilia. Various studies of C. reinhardtii have profoundly advanced plant and cell biology, and have also impacted algal biotechnology and our understanding of human disease. However, the 'real' life of C. reinhardtii in the natural environment has largely been neglected. To extend our understanding of the biology of C. reinhardtii, it will be rewarding to explore its behavior in its natural habitats, learning more about its abundance and life cycle, its genetic and physiological diversity, and its biotic and abiotic interactions.

Chloroplast nucleoids as a transformable network revealed by live imaging with a microfluidic device

Chloroplast DNA is organized into DNA–protein conglomerates called chloroplast nucleoids, which are replicated, transcribed, and inherited. We applied live-imaging technology with a microfluidic device to examine the nature of chloroplast nucleoids in Chlamydomonas reinhardtii. We observed the dynamic and reversible dispersion of globular chloroplast nucleoids into a network structure in dividing chloroplasts. In the monokaryotic chloroplast (moc) mutant, in which chloroplast nucleoids are unequally distributed following chloroplast division due to a defect in MOC1, the early stages of chloroplast nucleoid formation occurred mainly in the proximal area. This suggests the chloroplast nucleoid transformable network consists of a highly compact core with proximal areas associated with cpDNA replication and nucleoid formation.

Yoshitaka Kamimura and colleagues combine live-imaging technology with microfluidics to examine chloroplast DNA organization in nucleoids. They find that these structures form a network structure in dividing chloroplasts, and propose a mechanism for their inheritance in organelle replication.

Surprisal analysis of genome-wide transcript profiling identifies differentially expressed genes and pathways associated with four growth conditions in the microalga Chlamydomonas

The usual cultivation mode of the green microalga Chlamydomonas is liquid medium and light. However, the microalga can also be grown on agar plates and in darkness. Our aim is to analyze and compare gene expression of cells cultivated in these different conditions. For that purpose, RNA-seq data are obtained from Chlamydomonas samples of two different labs grown in four environmental conditions (agar@light, agar@dark, liquid@light, liquid@dark). The RNA seq data are analyzed by surprisal analysis, which allows the simultaneous meta-analysis of all the samples. First we identify a balance state, which defines a state where the expression levels are similar in all the samples irrespectively of their growth conditions, or lab origin. In addition our analysis identifies additional constraints needed to quantify the deviation with respect to the balance state. The first constraint differentiates the agar samples versus the liquid ones; the second constraint the dark samples versus the light ones. The two constraints are almost of equal importance. Pathways involved in stress responses are found in the agar phenotype while the liquid phenotype comprises ATP and NADH production pathways. Remodeling of membrane is suggested in the dark phenotype while photosynthetic pathways characterize the light phenotype. The same trends are also present when performing purely statistical analysis such as K-means clustering and differentially expressed genes.

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

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

Cell-Type Transcriptomes of the Multicellular Green Alga Volvox carteri Yield Insights into the Evolutionary Origins of Germ and Somatic Differentiation Programs

Germ-soma differentiation is a hallmark of complex multicellular organisms, yet its origins are not well understood. Volvox carteri is a simple multicellular green alga that has recently evolved a simple germ-soma dichotomy with only two cell-types: large germ cells called gonidia and small terminally differentiated somatic cells. Here, we provide a comprehensive characterization of the gonidial and somatic transcriptomes of V. carteri to uncover fundamental differences between the molecular and metabolic programming of these cell-types. We found extensive transcriptome differentiation between cell-types, with somatic cells expressing a more specialized program overrepresented in younger, lineage-specific genes, and gonidial cells expressing a more generalist program overrepresented in more ancient genes that shared striking overlap with stem cell-specific genes from animals and land plants. Directed analyses of different pathways revealed a strong dichotomy between cell-types with gonidial cells expressing growth-related genes and somatic cells expressing an altruistic metabolic program geared toward the assembly of flagella, which support organismal motility, and the conversion of storage carbon to sugars, which act as donors for production of extracellular matrix (ECM) glycoproteins whose secretion enables massive organismal expansion. V. carteri orthologs of diurnally controlled genes from C. reinhardtii, a single-celled relative, were analyzed for cell-type distribution and found to be strongly partitioned, with expression of dark-phase genes overrepresented in somatic cells and light-phase genes overrepresented in gonidial cells- a result that is consistent with cell-type programs in V. carteri arising by cooption of temporal regulons in a unicellular ancestor. Together, our findings reveal fundamental molecular, metabolic, and evolutionary mechanisms that underlie the origins of germ-soma differentiation in V. carteri and provide a template for understanding the acquisition of germ-soma differentiation in other multicellular lineages.

High-throughput sequencing of the chloroplast and mitochondrion of Chlamydomonas reinhardtii to generate improved de novo assemblies, analyze expression patterns and transcript speciation, and evaluate diversity among laboratory strains and wild isolates

Chlamydomonas reinhardtii is a unicellular chlorophyte alga that is widely studied as a reference organism for understanding photosynthesis, sensory and motile cilia, and for development of an algal-based platform for producing biofuels and bio-products. Its highly repetitive, ~205-kbp circular chloroplast genome and ~15.8-kbp linear mitochondrial genome were sequenced prior to the advent of high-throughput sequencing technologies. Here, high coverage shotgun sequencing was used to assemble both organellar genomes de novo. These new genomes correct dozens of errors in the prior genome sequences and annotations. Genome sequencing coverage indicates that each cell contains on average 83 copies of the chloroplast genome and 130 copies of the mitochondrial genome. Using protocols and analyses optimized for organellar transcripts, RNA-Seq was used to quantify their relative abundances across 12 different growth conditions. Forty-six percent of total cellular mRNA is attributable to high expression from a few dozen chloroplast genes. RNA-Seq data were used to guide gene annotation, to demonstrate polycistronic gene expression, and to quantify splicing of psaA and psbA introns. In contrast to a conclusion from a recent study, we found that chloroplast transcripts are not edited. Unexpectedly, cytosine-rich polynucleotide tails were observed at the 3'-end of all mitochondrial transcripts. A comparative genomics analysis of eight laboratory strains and 11 wild isolates of C. reinhardtii identified 2658 variants in the organellar genomes, which is 1/10th as much genetic diversity as is found in the nucleus.

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

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

Gene Regulatory Networks for the Haploid-to-Diploid Transition of Chlamydomonas reinhardtii

The sexual cycle of the unicellular Chlamydomonas reinhardtii culminates in the formation of diploid zygotes that differentiate into dormant spores that eventually undergo meiosis. Mating between gametes induces rapid cell wall shedding via the enzyme g-lysin; cell fusion is followed by heterodimerization of sex-specific homeobox transcription factors, GSM1 and GSP1, and initiation of zygote-specific gene expression. To investigate the genetic underpinnings of the zygote developmental pathway, we performed comparative transcriptome analysis of both pre- and post-fertilization samples. We identified 253 transcripts specifically enriched in early zygotes, 82% of which were not up-regulated in gsp1 null zygotes. We also found that the GSM1/GSP1 heterodimer negatively regulates the vegetative wall program at the posttranscriptional level, enabling prompt transition from vegetative wall to zygotic wall assembly. Annotation of the g-lysin-induced and early zygote genes reveals distinct vegetative and zygotic wall programs, supported by concerted up-regulation of genes encoding cell wall-modifying enzymes and proteins involved in nucleotide-sugar metabolism. The haploid-to-diploid transition in Chlamydomonas is masterfully controlled by the GSM1/GSP1 heterodimer, translating fertilization and gamete coalescence into a bona fide differentiation program. The fertilization-triggered integration of genes required to make related, but structurally and functionally distinct organelles-the vegetative versus zygote cell wall-presents a likely scenario for the evolution of complex developmental gene regulatory networks.

Chromosome-level genome assembly and transcriptome of the green alga Chromochloris zofingiensis illuminates astaxanthin production

Microalgae have potential to help meet energy and food demands without exacerbating environmental problems. There is interest in the unicellular green alga Chromochloris zofingiensis, because it produces lipids for biofuels and a highly valuable carotenoid nutraceutical, astaxanthin. To advance understanding of its biology and facilitate commercial development, we present a C. zofingiensis chromosome-level nuclear genome, organelle genomes, and transcriptome from diverse growth conditions. The assembly, derived from a combination of short- and long-read sequencing in conjunction with optical mapping, revealed a compact genome of ∼58 Mbp distributed over 19 chromosomes containing 15,274 predicted protein-coding genes. The genome has uniform gene density over chromosomes, low repetitive sequence content (∼6%), and a high fraction of protein-coding sequence (∼39%) with relatively long coding exons and few coding introns. Functional annotation of gene models identified orthologous families for the majority (∼73%) of genes. Synteny analysis uncovered localized but scrambled blocks of genes in putative orthologous relationships with other green algae. Two genes encoding beta-ketolase (BKT), the key enzyme synthesizing astaxanthin, were found in the genome, and both were up-regulated by high light. Isolation and molecular analysis of astaxanthin-deficient mutants showed that BKT1 is required for the production of astaxanthin. Moreover, the transcriptome under high light exposure revealed candidate genes that could be involved in critical yet missing steps of astaxanthin biosynthesis, including ABC transporters, cytochrome P450 enzymes, and an acyltransferase. The high-quality genome and transcriptome provide insight into the green algal lineage and carotenoid production.

A Plant Cryptochrome Controls Key Features of the Chlamydomonas Circadian Clock and Its Life Cycle

Cryptochromes are flavin-binding proteins that act as blue light receptors in bacteria, fungi, plants, and insects and are components of the circadian oscillator in mammals. Animal and plant cryptochromes are evolutionarily divergent, although the unicellular alga Chlamydomonas reinhardtii (Chlamydomonas throughout) has both an animal-like cryptochrome and a plant cryptochrome (pCRY; formerly designated CPH1). Here, we show that the pCRY protein accumulates at night as part of a complex. Functional characterization of pCRY was performed based on an insertional mutant that expresses only 11% of the wild-type pCRY level. The pcry mutant is defective for central properties of the circadian clock. In the mutant, the period is lengthened significantly, ultimately resulting in arrhythmicity, while blue light-based phase shifts show large deviations from what is observed in wild-type cells. We also show that pCRY is involved in gametogenesis in Chlamydomonas pCRY is down-regulated in pregametes and gametes, and in the pcry mutant, there is altered transcript accumulation under blue light of the strictly light-dependent, gamete-specific gene GAS28 pCRY acts as a negative regulator for the induction of mating ability in the light and for the loss of mating ability in the dark. Moreover, pCRY is necessary for light-dependent germination, during which the zygote undergoes meiosis that gives rise to four vegetative cells. In sum, our data demonstrate that pCRY is a key blue light receptor in Chlamydomonas that is involved in both circadian timing and life cycle progression.

Identification and Characterization of a cis-Regulatory Element for Zygotic Gene Expression in Chlamydomonas reinhardtii

Upon fertilization Chlamydomonas reinhardtii zygotes undergo a program of differentiation into a diploid zygospore that is accompanied by transcription of hundreds of zygote-specific genes. We identified a distinct sequence motif we term a zygotic response element (ZYRE) that is highly enriched in promoter regions of C reinhardtii early zygotic genes. A luciferase reporter assay was used to show that native ZYRE motifs within the promoter of zygotic gene ZYS3 or intron of zygotic gene DMT4 are necessary for zygotic induction. A synthetic luciferase reporter with a minimal promoter was used to show that ZYRE motifs introduced upstream are sufficient to confer zygotic upregulation, and that ZYRE-controlled zygotic transcription is dependent on the homeodomain transcription factor GSP1. We predict that ZYRE motifs will correspond to binding sites for the homeodomain proteins GSP1-GSM1 that heterodimerize and activate zygotic gene expression in early zygotes.

Genome and methylome of the oleaginous diatom Cyclotella cryptica reveal genetic flexibility toward a high lipid phenotype

BACKGROUND

Improvement in the performance of eukaryotic microalgae for biofuel and bioproduct production is largely dependent on characterization of metabolic mechanisms within the cell. The marine diatom Cyclotella cryptica, which was originally identified in the Aquatic Species Program, is a promising strain of microalgae for large-scale production of biofuel and bioproducts, such as omega-3 fatty acids.

RESULTS

We sequenced the nuclear genome and methylome of this oleaginous diatom to identify the genetic traits that enable substantial accumulation of triacylglycerol. The genome is comprised of highly methylated repetitive sequence, which does not significantly change under silicon starved lipid induction, and data further suggests the primary role of DNA methylation is to suppress DNA transposition. Annotation of pivotal glycolytic, lipid metabolism, and carbohydrate degradation processes reveal an expanded enzyme repertoire in C. cryptica that would allow for an increased metabolic capacity toward triacylglycerol production. Identification of previously unidentified genes, including those involved in carbon transport and chitin metabolism, provide potential targets for genetic manipulation of carbon flux to further increase its lipid phenotype. New genetic tools were developed, bringing this organism on a par with other microalgae in terms of genetic manipulation and characterization approaches.

CONCLUSIONS

Functional annotation and detailed cross-species comparison of key carbon rich processes in C. cryptica highlights the importance of enzymatic subcellular compartmentation for regulation of carbon flux, which is often overlooked in photosynthetic microeukaryotes. The availability of the genome sequence, as well as advanced genetic manipulation tools enable further development of this organism for deployment in large-scale production systems.