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

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

JGI Plant Gene Atlas: an updateable transcriptome resource to improve functional gene descriptions across the plant kingdom.

Gene functional descriptions offer a crucial line of evidence for candidate genes underlying trait variation. Conversely, plant responses to environmental cues represent important resources to decipher gene function and subsequently provide molecular targets for plant improvement through gene editing. However, biological roles of large proportions of genes across the plant phylogeny are poorly annotated. Here we describe the Joint Genome Institute (JGI) Plant Gene Atlas, an updateable data resource consisting of transcript abundance assays spanning 18 diverse species. To integrate across these diverse genotypes, we analyzed expression profiles, built gene clusters that exhibited tissue/condition specific expression, and tested for transcriptional response to environmental queues. We discovered extensive phylogenetically constrained and condition-specific expression profiles for genes without any previously documented functional annotation. Such conserved expression patterns and tightly co-expressed gene clusters let us assign expression derived additional biological information to 64 495 genes with otherwise unknown functions. The ever-expanding Gene Atlas resource is available at JGI Plant Gene Atlas (https://plantgeneatlas.jgi.doe.gov) and Phytozome (https://phytozome.jgi.doe.gov/), providing bulk access to data and user-specified queries of gene sets. Combined, these web interfaces let users access differentially expressed genes, track orthologs across the Gene Atlas plants, graphically represent co-expressed genes, and visualize gene ontology and pathway enrichments.

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

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

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.

Widespread polycistronic gene expression in green algae.

Polycistronic gene expression, common in prokaryotes, was thought to be extremely rare in eukaryotes. The development of long-read sequencing of full-length transcript isomers (Iso-Seq) has facilitated a reexamination of that dogma. Using Iso-Seq, we discovered hundreds of examples of polycistronic expression of nuclear genes in two divergent species of green algae: Chlamydomonas reinhardtii and Chromochloris zofingiensis Here, we employ a range of independent approaches to validate that multiple proteins are translated from a common transcript for hundreds of loci. A chromatin immunoprecipitation analysis using trimethylation of lysine 4 on histone H3 marks confirmed that transcription begins exclusively at the upstream gene. Quantification of polyadenylated [poly(A)] tails and poly(A) signal sequences confirmed that transcription ends exclusively after the downstream gene. Coexpression analysis found nearly perfect correlation for open reading frames (ORFs) within polycistronic loci, consistent with expression in a shared transcript. For many polycistronic loci, terminal peptides from both ORFs were identified from proteomics datasets, consistent with independent translation. Synthetic polycistronic gene pairs were transcribed and translated in vitro to recapitulate the production of two distinct proteins from a common transcript. The relative abundance of these two proteins can be modified by altering the Kozak-like sequence of the upstream gene. Replacement of the ORFs with selectable markers or reporters allows production of such heterologous proteins, speaking to utility in synthetic biology approaches. Conservation of a significant number of polycistronic gene pairs between C. reinhardtii, C. zofingiensis, and five other species suggests that this mechanism may be evolutionarily ancient and biologically important in the green algal lineage.

An epigenetic gene silencing pathway selectively acting on transgenic DNA in the green alga Chlamydomonas.

Silencing of exogenous DNA can make transgene expression very inefficient. Genetic screens in the model alga Chlamydomonas have demonstrated that transgene silencing can be overcome by mutations in unknown gene(s), thus producing algal strains that stably express foreign genes to high levels. Here, we show that the silencing mechanism specifically acts on transgenic DNA. Once a permissive chromatin structure has assembled, transgene expression can persist even in the absence of mutations disrupting the silencing pathway. We have identified the gene conferring the silencing and show it to encode a sirtuin-type histone deacetylase. Loss of gene function does not appreciably affect endogenous gene expression. Our data suggest that transgenic DNA is recognized and then quickly inactivated by the assembly of a repressive chromatin structure composed of deacetylated histones. We propose that this mechanism may have evolved to provide protection from potentially harmful types of environmental DNA.

Regulation of Oxygenic Photosynthesis during Trophic Transitions in the Green Alga Chromochloris zofingiensis.

Light and nutrients are critical regulators of photosynthesis and metabolism in plants and algae. Many algae have the metabolic flexibility to grow photoautotrophically, heterotrophically, or mixotrophically. Here, we describe reversible Glc-dependent repression/activation of oxygenic photosynthesis in the unicellular green alga Chromochloris zofingiensis. We observed rapid and reversible changes in photosynthesis, in the photosynthetic apparatus, in thylakoid ultrastructure, and in energy stores including lipids and starch. Following Glc addition in the light, C. zofingiensis shuts off photosynthesis within days and accumulates large amounts of commercially relevant bioproducts, including triacylglycerols and the high-value nutraceutical ketocarotenoid astaxanthin, while increasing culture biomass. RNA sequencing reveals reversible changes in the transcriptome that form the basis of this metabolic regulation. Functional enrichment analyses show that Glc represses photosynthetic pathways while ketocarotenoid biosynthesis and heterotrophic carbon metabolism are upregulated. Because sugars play fundamental regulatory roles in gene expression, physiology, metabolism, and growth in both plants and animals, we have developed a simple algal model system to investigate conserved eukaryotic sugar responses as well as mechanisms of thylakoid breakdown and biogenesis in chloroplasts. Understanding regulation of photosynthesis and metabolism in algae could enable bioengineering to reroute metabolism toward beneficial bioproducts for energy, food, pharmaceuticals, and human health.

Multiomics resolution of molecular events during a day in the life of Chlamydomonas.

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.

A Gelatin Microdroplet Platform for High-Throughput Sorting of Hyperproducing Single-Cell-Derived Microalgal Clones.

Microalgae are an attractive feedstock organism for sustainable production of biofuels, chemicals, and biomaterials, but the ability to rationally engineer microalgae to enhance production has been limited. To enable the evolution-based selection of new hyperproducing variants of microalgae, a method is developed that combines phase-transitioning monodisperse gelatin hydrogel droplets with commercial flow cytometric instruments for high-throughput screening and selection of clonal populations of cells with desirable properties, such as high lipid productivity per time traced over multiple cell cycles. It is found that gelatin microgels enable i) the growth and metabolite (e.g., chlorophyll and lipids) production of single microalgal cells within the compartments, ii) infusion of fluorescent reporter molecules into the hydrogel matrices following a sol-gel transition, iii) selection of high-producing clonal populations of cells using flow cytometry, and iv) cell recovery under mild conditions, enabling regrowth after sorting. This user-friendly method is easily integratable into directed cellular evolution pipelines for strain improvement and can be adopted for other applications that require high-throughput processing, e.g., cellular secretion phenotypes and intercellular interactions.

RNA Purification from the Unicellular Green Alga, Chromochloris zofingiensis.

Chromochloris zofingiensis is a unicellular green alga that is an emerging model species for studies in fields such as biofuel production, ketocarotenoid biosynthesis and metabolism. The recent availability of a high-quality genome assembly facilitates systems-level analysis, such as RNA-Seq. However, cells of this alga have a tough cell wall, which is a challenge for RNA purification. This protocol was designed specifically to breach the cell wall and isolate high-quality RNA suitable for RNA-Seq studies.

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.

Bilin-Dependent Photoacclimation in Chlamydomonas reinhardtii.

In land plants, linear tetrapyrrole (bilin)-based phytochrome photosensors optimize photosynthetic light capture by mediating massive reprogramming of gene expression. But, surprisingly, many green algal genomes lack phytochrome genes. Studies of the heme oxygenase mutant (hmox1) of the green alga Chlamydomonas reinhardtii suggest that bilin biosynthesis in plastids is essential for proper regulation of a nuclear gene network implicated in oxygen detoxification during dark-to-light transitions. hmox1 cannot grow photoautotrophically and photoacclimates poorly to increased illumination. We show that these phenotypes are due to reduced accumulation of photosystem I (PSI) reaction centers, the PSI electron acceptors 5'-monohydroxyphylloquinone and phylloquinone, and the loss of PSI and photosystem II antennae complexes during photoacclimation. The hmox1 mutant resembles chlorophyll biosynthesis mutants phenotypically, but can be rescued by exogenous biliverdin IXα, the bilin produced by HMOX1. This rescue is independent of photosynthesis and is strongly dependent on blue light. RNA-seq comparisons of hmox1, genetically complemented hmox1, and chemically rescued hmox1 reveal that tetrapyrrole biosynthesis and known photoreceptor and photosynthesis-related genes are not impacted in the hmox1 mutant at the transcript level. We propose that a bilin-based, blue-light-sensing system within plastids evolved together with a bilin-based retrograde signaling pathway to ensure that a robust photosynthetic apparatus is sustained in light-grown Chlamydomonas.

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.

Mechanisms of Groucho-mediated repression revealed by genome-wide analysis of Groucho binding and activity.

BACKGROUND: The transcriptional corepressor Groucho (Gro) is required for the function of many developmentally regulated DNA binding repressors, thus helping to define the gene expression profile of each cell during development. The ability of Gro to repress transcription at a distance together with its ability to oligomerize and bind to histones has led to the suggestion that Gro may spread along chromatin. However, much is unknown about the mechanism of Gro-mediated repression and about the dynamics of Gro targeting. RESULTS: Our chromatin immunoprecipitation sequencing analysis of temporally staged Drosophila embryos shows that Gro binds in a highly dynamic manner primarily to clusters of discrete (<1 kb) segments. Consistent with the idea that Gro may facilitate communication between silencers and promoters, Gro binding is enriched at both cis-regulatory modules, as well as within the promotors of potential target genes. While this Gro-recruitment is required for repression, our data show that it is not sufficient for repression. Integration of Gro binding data with transcriptomic analysis suggests that, contrary to what has been observed for another Gro family member, Drosophila Gro is probably a dedicated repressor. This analysis also allows us to define a set of high confidence Gro repression targets. Using publically available data regarding the physical and genetic interactions between these targets, we are able to place them in the regulatory network controlling development. Through analysis of chromatin associated pre-mRNA levels at these targets, we find that genes regulated by Gro in the embryo are enriched for characteristics of promoter proximal paused RNA polymerase II. CONCLUSIONS: Our findings are inconsistent with a one-dimensional spreading model for long-range repression and suggest that Gro-mediated repression must be regulated at a post-recruitment step. They also show that Gro is likely a dedicated repressor that sits at a prominent highly interconnected regulatory hub in the developmental network. Furthermore, our findings suggest a role for RNA polymerase II pausing in Gro-mediated repression.

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.

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

The green alga Chlamydomonas reinhardtii undergoes gametogenesis and mating upon nitrogen starvation. While the steps involved in its sexual reproductive cycle have been extensively characterized, the genome-wide transcriptional and epigenetic changes underlying different life cycle stages have yet to be fully described. Here, we performed transcriptome and methylome sequencing to quantify expression and DNA methylation from vegetative and gametic cells of each mating type and from zygotes. We identified 361 gametic genes with mating type-specific expression patterns and 627 genes that are specifically induced in zygotes; furthermore, these sex-related gene sets were enriched for secretory pathway and alga-specific genes. We also examined the C. reinhardtii nuclear methylation map with base-level resolution at different life cycle stages. Despite having low global levels of nuclear methylation, we detected 23 hypermethylated loci in gene-poor, repeat-rich regions. We observed mating type-specific differences in chloroplast DNA methylation levels in plus versus minus mating type gametes followed by chloroplast DNA hypermethylation in zygotes. Lastly, we examined the expression of candidate DNA methyltransferases and found three, DMT1a, DMT1b, and DMT4, that are differentially expressed during the life cycle and are candidate DNA methylases. The expression and methylation data we present provide insight into cell type-specific transcriptional and epigenetic programs during key stages of the C. reinhardtii life cycle.

Chlamydomonas Genome Resource for Laboratory Strains Reveals a Mosaic of Sequence Variation, Identifies True Strain Histories, and Enables Strain-Specific Studies.

Chlamydomonas reinhardtii is a widely used reference organism in studies of photosynthesis, cilia, and biofuels. Most research in this field uses a few dozen standard laboratory strains that are reported to share a common ancestry, but exhibit substantial phenotypic differences. In order to facilitate ongoing Chlamydomonas research and explain the phenotypic variation, we mapped the genetic diversity within these strains using whole-genome resequencing. We identified 524,640 single nucleotide variants and 4812 structural variants among 39 commonly used laboratory strains. Nearly all (98.2%) of the total observed genetic diversity was attributable to the presence of two, previously unrecognized, alternate haplotypes that are distributed in a mosaic pattern among the extant laboratory strains. We propose that these two haplotypes are the remnants of an ancestral cross between two strains with ∼2% relative divergence. These haplotype patterns create a fingerprint for each strain that facilitates the positive identification of that strain and reveals its relatedness to other strains. The presence of these alternate haplotype regions affects phenotype scoring and gene expression measurements. Here, we present a rich set of genetic differences as a community resource to allow researchers to more accurately conduct and interpret their experiments with Chlamydomonas.

Copper economy in Chlamydomonas: prioritized allocation and reallocation of copper to respiration vs. photosynthesis.

Inorganic elements, although required only in trace amounts, permit life and primary productivity because of their functions in catalysis. Every organism has a minimal requirement of each metal based on the intracellular abundance of proteins that use inorganic cofactors, but elemental sparing mechanisms can reduce this quota. A well-studied copper-sparing mechanism that operates in microalgae faced with copper deficiency is the replacement of the abundant copper protein plastocyanin with a heme-containing substitute, cytochrome (Cyt) c6. This switch, which is dependent on a copper-sensing transcription factor, copper response regulator 1 (CRR1), dramatically reduces the copper quota. We show here that in a situation of marginal copper availability, copper is preferentially allocated from plastocyanin, whose function is dispensable, to other more critical copper-dependent enzymes like Cyt oxidase and a ferroxidase. In the absence of an extracellular source, copper allocation to Cyt oxidase includes CRR1-dependent proteolysis of plastocyanin and quantitative recycling of the copper cofactor from plastocyanin to Cyt oxidase. Transcriptome profiling identifies a gene encoding a Zn-metalloprotease, as a candidate effecting copper recycling. One reason for the retention of genes encoding both plastocyanin and Cyt c6 in algal and cyanobacterial genomes might be because plastocyanin provides a competitive advantage in copper-depleted environments as a ready source of copper.