Comparative genomics of Chlamydomonas

Characterization of activating cis-regulatory elements from the histone genes of Chlamydomonas reinhardtii

Regulation of gene expression in eukaryotes is controlled by cis-regulatory modules (CRMs). A major class of CRMs are enhancers which are composed of activating cis-regulatory elements (CREs) responsible for upregulating transcription. To date, most enhancers and activating CREs have been studied in angiosperms; in contrast, our knowledge about these key regulators of gene expression in green algae is limited. In this study, we aimed at characterizing putative activating CREs/CRMs from the histone genes of the unicellular model alga Chlamydomonas reinhardtii. To test the activity of four candidates, reporter constructs consisting of a tetramerized CRE, an established promoter, and a gene for the mCerulean3 fluorescent protein were incorporated into the nuclear genome of C. reinhardtii, and their activity was quantified by flow cytometry. Two tested candidates, Eupstr and Ehist cons, significantly upregulated gene expression and were characterized in detail. Eupstr, which originates from highly expressed genes of C. reinhardtii, is an orientation-independent CRE capable of activating both the RBCS2 and β2-tubulin promoters. Ehist cons, which is a CRM from histone genes of angiosperms, upregulates the β2-tubulin promoter in C. reinhardtii over a distance of at least 1.5 kb. The octamer motif present in Ehist cons was identified in C. reinhardtii and the related green algae Chlamydomonas incerta, Chlamydomonas schloesseri, and Edaphochlamys debaryana, demonstrating its high evolutionary conservation. The results of this investigation expand our knowledge about the regulation of gene expression in green algae. Furthermore, the characterized activating CREs/CRMs can be applied as valuable genetic tools.

Heterochromatin repeat organization at an individual level: Rex1BD and the 14-3-3 protein coordinate to shape the epigenetic landscape within heterochromatin repeats

In eukaryotic cells, heterochromatin is typically composed of tandem DNA repeats and plays crucial roles in gene expression and genome stability. It has been reported that silencing at individual units within tandem heterochromatin repeats exhibits a position-dependent variation. However, how the heterochromatin is organized at an individual repeat level remains poorly understood. Using a novel genetic approach, our recent study identified a conserved protein Rex1BD required for position-dependent silencing within heterochromatin repeats. We further revealed that Rex1BD interacts with the 14-3-3 protein to regulate heterochromatin silencing by linking RNAi and HDAC pathways. In this review, we discuss how Rex1BD and the 14-3-3 protein coordinate to modulate heterochromatin organization at the individual repeat level, and comment on the biological significance of the position-dependent effect in heterochromatin repeats. We also identify the knowledge gaps that still need to be unveiled in the field.

Near telomere-to-telomere genome assemblies of two Chlorella species unveil the composition and evolution of centromeres in green algae

BACKGROUND

Centromeres play a crucial and conserved role in cell division, although their composition and evolutionary history in green algae, the evolutionary ancestors of land plants, remains largely unknown.

RESULTS

We constructed near telomere-to-telomere (T2T) assemblies for two Trebouxiophyceae species, Chlorella sorokiniana NS4-2 and Chlorella pyrenoidosa DBH, with chromosome numbers of 12 and 13, and genome sizes of 58.11 Mb and 53.41 Mb, respectively. We identified and validated their centromere sequences using CENH3 ChIP-seq and found that, similar to humans and higher plants, the centromeric CENH3 signals of green algae display a pattern of hypomethylation. Interestingly, the centromeres of both species largely comprised transposable elements, although they differed significantly in their composition. Species within the Chlorella genus display a more diverse centromere composition, with major constituents including members of the LTR/Copia, LINE/L1, and LINE/RTEX families. This is in contrast to green algae including Chlamydomonas reinhardtii, Coccomyxa subellipsoidea, and Chromochloris zofingiensis, in which centromere composition instead has a pronounced single-element composition. Moreover, we observed significant differences in the composition and structure of centromeres among chromosomes with strong collinearity within the Chlorella genus, suggesting that centromeric sequence evolves more rapidly than sequence in non-centromeric regions.

CONCLUSIONS

This study not only provides high-quality genome data for comparative genomics of green algae but gives insight into the composition and evolutionary history of centromeres in early plants, laying an important foundation for further research on their evolution.

The structure, function, and evolution of plant centromeres

Centromeres are essential regions of eukaryotic chromosomes responsible for the formation of kinetochore complexes, which connect to spindle microtubules during cell division. Notably, although centromeres maintain a conserved function in chromosome segregation, the underlying DNA sequences are diverse both within and between species and are predominantly repetitive in nature. The repeat content of centromeres includes high-copy tandem repeats (satellites), and/or specific families of transposons. The functional region of the centromere is defined by loading of a specific histone 3 variant (CENH3), which nucleates the kinetochore and shows dynamic regulation. In many plants, the centromeres are composed of satellite repeat arrays that are densely DNA methylated and invaded by centrophilic retrotransposons. In some cases, the retrotransposons become the sites of CENH3 loading. We review the structure of plant centromeres, including monocentric, holocentric, and metapolycentric architectures, which vary in the number and distribution of kinetochore attachment sites along chromosomes. We discuss how variation in CENH3 loading can drive genome elimination during early cell divisions of plant embryogenesis. We review how epigenetic state may influence centromere identity and discuss evolutionary models that seek to explain the paradoxically rapid change of centromere sequences observed across species, including the potential roles of recombination. We outline putative modes of selection that could act within the centromeres, as well as the role of repeats in driving cycles of centromere evolution. Although our primary focus is on plant genomes, we draw comparisons with animal and fungal centromeres to derive a eukaryote-wide perspective of centromere structure and function.

Optimized transgene expression in the red alga Porphyridium purpureum and efficient recombinant protein secretion into the culture medium

Microalgae represent a promising but yet underexplored production platform for biotechnology. The vast majority of studies on recombinant protein expression in algae have been conducted in a single species, the green alga Chlamydomonas reinhardtii. However, due to epigenetic silencing, transgene expression in Chlamydomonas is often inefficient. Here we have investigated parameters that govern efficient transgene expression in the red microalga Porphyridium purpureum. Porphyridium is unique in that the introduced transformation vectors are episomally maintained as autonomously replicating plasmids in the nucleus. We show that full codon optimization to the preferred codon usage in the Porphyridium genome confers superior transgene expression, not only at the level of protein accumulation, but also at the level of mRNA accumulation, indicating that high translation rates increase mRNA stability. Our optimized expression constructs resulted in YFP accumulation to unprecedented levels of up to 5% of the total soluble protein. We also designed expression cassettes that target foreign proteins to the secretory pathway and lead to efficient protein secretion into the culture medium, thus simplifying recombinant protein harvest and purification. Our study paves the way to the exploration of red microalgae as expression hosts in molecular farming for recombinant proteins and metabolites.

Conserved cobalamin acquisition protein 1 is essential for vitamin B12 uptake in both Chlamydomonas and Phaeodactylum

Microalgae play an essential role in global net primary productivity and global biogeochemical cycling. Despite their phototrophic lifestyle, over half of algal species depend for growth on acquiring an external supply of the corrinoid vitamin B12 (cobalamin), a micronutrient produced only by a subset of prokaryotic organisms. Previous studies have identified protein components involved in vitamin B12 uptake in bacterial species and humans. However, little is known about its uptake in algae. Here, we demonstrate the essential role of a protein, cobalamin acquisition protein 1 (CBA1), in B12 uptake in Phaeodactylum tricornutum using CRISPR-Cas9 to generate targeted knockouts and in Chlamydomonas reinhardtii by insertional mutagenesis. In both cases, CBA1 knockout lines could not take up exogenous vitamin B12. Complementation of the C. reinhardtii mutants with the wild-type CBA1 gene restored B12 uptake, and regulation of CBA1 expression via a riboswitch element enabled control of the phenotype. When visualized by confocal microscopy, a YFP-fusion with C. reinhardtii CBA1 showed association with membranes. Bioinformatics analysis found that CBA1-like sequences are present in all major eukaryotic phyla. In algal taxa, the majority that encoded CBA1 also had genes for B12-dependent enzymes, suggesting CBA1 plays a conserved role. Our results thus provide insight into the molecular basis of algal B12 acquisition, a process that likely underpins many interactions in aquatic microbial communities.

Aberrant light sensing and motility in the green alga Chlamydomonas priscuii from the ice-covered Antarctic Lake Bonney

The Antarctic green alga Chlamydomonas priscuii is an obligate psychrophile and an emerging model for photosynthetic adaptation to extreme conditions. Endemic to the ice-covered Lake Bonney, this alga thrives at highly unusual light conditions characterized by very low light irradiance (<15 μmol m-2 s-1), a narrow wavelength spectrum enriched in blue light, and an extreme photoperiod. Genome sequencing of C. priscuii exposed an unusually large genome, with hundreds of highly similar gene duplicates and expanded gene families, some of which could be aiding its survival in extreme conditions. In contrast to the described expansion in the genetic repertoire in C. priscuii, here we suggest that the gene family encoding for photoreceptors is reduced when compared to related green algae. This alga also possesses a very small eyespot and exhibits an aberrant phototactic response, compared to the model Chlamydomonas reinhardtii. We also investigated the genome and behavior of the closely related psychrophilic alga Chlamydomonas sp. ICE-MDV, that is found throughout the photic zone of Lake Bonney and is naturally exposed to higher light levels. Our analyses revealed a photoreceptor gene family and a robust phototactic response similar to those in the model Chlamydomonas reinhardtii. These results suggest that the aberrant phototactic response in C. priscuii is a result of life under extreme shading rather than a common feature of all psychrophilic algae. We discuss the implications of these results on the evolution and survival of shade adapted polar algae.

The extracellular matrix of green algae

Green algae display a wide range of extracellular matrix (ECM) components that include various types of cell walls (CW), scales, crystalline glycoprotein coverings, hydrophobic compounds, and complex gels or mucilage. Recently, new information derived from genomic/transcriptomic screening, advanced biochemical analyses, immunocytochemical studies, and ecophysiology has significantly enhanced and refined our understanding of the green algal ECM. In the later diverging charophyte group of green algae, the CW and other ECM components provide insight into the evolution of plants and the ways the ECM modulates during environmental stress. Chlorophytes produce diverse ECM components, many of which have been exploited for various uses in medicine, food, and biofuel production. This review highlights major advances in ECM studies of green algae.

MarFERReT, an open-source, version-controlled reference library of marine microbial eukaryote functional genes

Metatranscriptomics generates large volumes of sequence data about transcribed genes in natural environments. Taxonomic annotation of these datasets depends on availability of curated reference sequences. For marine microbial eukaryotes, current reference libraries are limited by gaps in sequenced organism diversity and barriers to updating libraries with new sequence data, resulting in taxonomic annotation of about half of eukaryotic environmental transcripts. Here, we introduce Marine Functional EukaRyotic Reference Taxa (MarFERReT), a marine microbial eukaryotic sequence library designed for use with taxonomic annotation of eukaryotic metatranscriptomes. We gathered 902 publicly accessible marine eukaryote genomes and transcriptomes and assessed their sequence quality and cross-contamination issues, selecting 800 validated entries for inclusion in MarFERReT. Version 1.1 of MarFERReT contains reference sequences from 800 marine eukaryotic genomes and transcriptomes, covering 453 species- and strain-level taxa, totaling nearly 28 million protein sequences with associated NCBI and PR2 Taxonomy identifiers and Pfam functional annotations. The MarFERReT project repository hosts containerized build scripts, documentation on installation and use case examples, and information on new versions of MarFERReT.

Chloroplast ATP synthase biogenesis requires peripheral stalk subunits AtpF and ATPG and stabilization of atpE mRNA by OPR protein MDE1

Chloroplast ATP synthase contains subunits of plastid and nuclear genetic origin. To investigate the coordinated biogenesis of this complex, we isolated novel ATP synthase mutants in the green alga Chlamydomonas reinhardtii by screening for high light sensitivity. We report here the characterization of mutants affecting the two peripheral stalk subunits b and b', encoded respectively by the atpF and ATPG genes, and of three independent mutants which identify the nuclear factor MDE1, required to stabilize the chloroplast-encoded atpE mRNA. Whole-genome sequencing revealed a transposon insertion in the 3'UTR of ATPG while mass spectrometry shows a small accumulation of functional ATP synthase in this knock-down ATPG mutant. In contrast, knock-out ATPG mutants, obtained by CRISPR-Cas9 gene editing, fully prevent ATP synthase function and accumulation, as also observed in an atpF frame-shift mutant. Crossing ATP synthase mutants with the ftsh1-1 mutant of the major thylakoid protease identifies AtpH as an FTSH substrate, and shows that FTSH significantly contributes to the concerted accumulation of ATP synthase subunits. In mde1 mutants, the absence of atpE transcript fully prevents ATP synthase biogenesis and photosynthesis. Using chimeric atpE genes to rescue atpE transcript accumulation, we demonstrate that MDE1, a novel octotricopeptide repeat (OPR) protein, genetically targets the atpE 5'UTR. In the perspective of the primary endosymbiosis (~1.5 Gy), the recruitment of MDE1 to its atpE target exemplifies a nucleus/chloroplast interplay that evolved rather recently, in the ancestor of the CS clade of Chlorophyceae, ~300 My ago.

Predictions of DNA mechanical properties at a genomic scale reveal potentially new functional roles of DNA flexibility

Mechanical properties of DNA have been implied to influence many of its biological functions. Recently, a new high-throughput method, called loop-seq, which allows measuring the intrinsic bendability of DNA fragments, has been developed. Using loop-seq data, we created a deep learning model to explore the biological significance of local DNA flexibility in a range of different species from different kingdoms. Consistently, we observed a characteristic and largely dinucleotide-composition-driven change of local flexibility near transcription start sites. In the presence of a TATA-box, a pronounced peak of high flexibility can be observed. Furthermore, depending on the transcription factor investigated, flanking-sequence-dependent DNA flexibility was identified as a potential factor influencing DNA binding. Compared to randomized genomic sequences, depending on species and taxa, actual genomic sequences were observed both with increased and lowered flexibility. Furthermore, in Arabidopsis thaliana, mutation rates, both de novo and fixed, were found to be associated with relatively rigid sequence regions. Our study presents a range of significant correlations between characteristic DNA mechanical properties and genomic features, the significance of which with regard to detailed molecular relevance awaits further theoretical and experimental exploration.

Telomerase-independent survival leads to a mosaic of complex subtelomere rearrangements in Chlamydomonas reinhardtii

Telomeres and subtelomeres, the genomic regions located at chromosome extremities, are essential for genome stability in eukaryotes. In the absence of the canonical maintenance mechanism provided by telomerase, telomere shortening induces genome instability. The landscape of the ensuing genome rearrangements is not accessible by short-read sequencing. Here, we leverage Oxford Nanopore Technologies long-read sequencing to survey the extensive repertoire of genome rearrangements in telomerase mutants of the model green microalga Chlamydomonas reinhardtii In telomerase-mutant strains grown for hundreds of generations, most chromosome extremities were capped by short telomere sequences that were either recruited de novo from other loci or maintained in a telomerase-independent manner. Other extremities did not end with telomeres but only with repeated subtelomeric sequences. The subtelomeric elements, including rDNA, were massively rearranged and involved in breakage-fusion-bridge cycles, translocations, recombinations, and chromosome circularization. These events were established progressively over time and displayed heterogeneity at the subpopulation level. New telomere-capped extremities composed of sequences originating from more internal genomic regions were associated with high DNA methylation, suggesting that de novo heterochromatin formation contributes to the restoration of chromosome end stability in C. reinhardtii The diversity of alternative strategies present in the same organism to maintain chromosome integrity and the variety of rearrangements found in telomerase mutants are remarkable, and illustrate genome plasticity at short timescales.

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.

Small RNAs >26 nt in length associate with AGO1 and are upregulated by nutrient deprivation in the alga Chlamydomonas

Small RNAs (sRNAs) associate with ARGONAUTE (AGO) proteins forming effector complexes with key roles in gene regulation and defense responses against molecular parasites. In multicellular eukaryotes, extensive duplication and diversification of RNA interference (RNAi) components have resulted in intricate pathways for epigenetic control of gene expression. The unicellular alga Chlamydomonas reinhardtii also has a complex RNAi machinery, including 3 AGOs and 3 DICER-like proteins. However, little is known about the biogenesis and function of most endogenous sRNAs. We demonstrate here that Chlamydomonas contains uncommonly long (>26 nt) sRNAs that associate preferentially with AGO1. Somewhat reminiscent of animal PIWI-interacting RNAs, these >26 nt sRNAs are derived from moderately repetitive genomic clusters and their biogenesis is DICER-independent. Interestingly, the sequences generating these >26-nt sRNAs have been conserved and amplified in several Chlamydomonas species. Moreover, expression of these longer sRNAs increases substantially under nitrogen or sulfur deprivation, concurrently with the downregulation of predicted target transcripts. We hypothesize that the transposon-like sequences from which >26-nt sRNAs are produced might have been ancestrally targeted for silencing by the RNAi machinery but, during evolution, certain sRNAs might have fortuitously acquired endogenous target genes and become integrated into gene regulatory networks.

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 origin and early evolution of plants

Plant (archaeplastid) evolution has transformed the biosphere, but we are only now beginning to learn how this took place through comparative genomics, phylogenetics, and the fossil record. This has illuminated the phylogeny of Archaeplastida, Viridiplantae, and Streptophyta, and has resolved the evolution of key characters, genes, and genomes - revealing that many key innovations evolved long before the clades with which they have been casually associated. Molecular clock analyses estimate that Streptophyta and Viridiplantae emerged in the late Mesoproterozoic to late Neoproterozoic, whereas Archaeplastida emerged in the late-mid Palaeoproterozoic. Together, these insights inform on the coevolution of plants and the Earth system that transformed ecology and global biogeochemical cycles, increased weathering, and precipitated snowball Earth events, during which they would have been key to oxygen production and net primary productivity (NPP).

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.

A novel method for identifying Chlamydomonas reinhardtii (Chlorophyta) and closely related species from nature

Here, we introduce a new method for efficiently sampling Chlamydomonas reinhardtii and closely related species using a colony PCR-based screen with novel primer sets designed to specifically detect these important model microalgae. To demonstrate the utility of our new method, we collected 130 soil samples from a wide range of habitats in Ontario, Canada and identified 33 candidate algae, which were barcoded by sequencing a region of the rbcL plastid gene. For select isolates, 18S rRNA gene and YPT4 nuclear markers were also sequenced. Based on phylogenetic and haplotype network analyses of these three loci, seven novel isolates were identified as C. reinhardtii, and one additional isolate appeared to be more closely related to C. reinhardtii than any other known species. All seven new C. reinhardtii strains were interfertile with previously collected C. reinhardtii field isolates, validating the effectiveness of our molecular screen.

Toxicology and digestibility of Chlamydomonas debaryana green algal biomass

There is an economic interest, both for food security and for the non-meat-eating population, in the development of novel, sustainable sources of high-quality protein. The green algae Chlamydomonas reinhardtii has already been developed for this purpose, and the closely related species, Chlamydomonas debaryana, is a complementary source that also presents some additional advantages, such as reduced production cost. To determine whether C. debaryana may have a similar safety profile to that of C. reinhardtii, a wild type strain was obtained, designated TS04 after confirmation of its identity, and subjected to a battery of preclinical studies. Genetic toxicity was evaluated using a bacterial reverse mutation test, an in vitro mammalian chromosomal aberration test, and an in vivo mammalian micronucleus test in a mouse model. No genotoxic potential (e.g., mutagenicity and clastogenicity) was observed in these tests under the employed conditions up to maximum recommended concentrations or doses. To assess general toxicity, a 90-day repeated-dose oral toxicity study was conducted in rats. No mortality or adverse effects were observed, and no target organs were identified up to the maximum feasible dose, due to solubility, of 4,000 mg/kg bw/day. The no-observed-adverse-effect level was determined as the highest dose tested. A digestibility study in simulated gastric fluid was conducted and determined that TS04 has low allergenic potential, exhibiting rapid digestion of proteins. Due to the negative results of our evaluation, it is reasonable to proceed with further development and additional investigations to contribute towards a safety assessment of the proposed use in food for human consumption.

Rates and spectra of de novo structural mutations in Chlamydomonas reinhardtii

Genetic variation originates from several types of spontaneous mutation, including single-nucleotide substitutions, short insertions and deletions (indels), and larger structural changes. Structural mutations (SMs) drive genome evolution and are thought to play major roles in evolutionary adaptation, speciation, and genetic disease, including cancers. Sequencing of mutation accumulation (MA) lines has provided estimates of rates and spectra of single-nucleotide and indel mutations in many species, yet the rate of new SMs is largely unknown. Here, we use long-read sequencing to determine the full mutation spectrum in MA lines derived from two strains (CC-1952 and CC-2931) of the green alga Chlamydomonas reinhardtii The SM rate is highly variable between strains and between MA lines, and SMs represent a substantial proportion of all mutations in both strains (CC-1952 6%; CC-2931 12%). The SM spectra differ considerably between the two strains, with almost all inversions and translocations occurring in CC-2931 MA lines. This variation is associated with heterogeneity in the number and type of active transposable elements (TEs), which comprise major proportions of SMs in both strains (CC-1952 22%; CC-2931 38%). In CC-2931, a Crypton and a previously undescribed type of DNA element have caused 71% of chromosomal rearrangements, whereas in CC-1952, a Dualen LINE is associated with 87% of duplications. Other SMs, notably large duplications in CC-2931, are likely products of various double-strand break repair pathways. Our results show that diverse types of SMs occur at substantial rates, and support prominent roles for SMs and TEs in evolution.

Trait drift in microalgae and applications for strain improvement

Microalgae are increasingly used to generate a wide range of commercial products, and there is growing evidence that microalgae-based products can be produced sustainably. However, industrial production of microalgal biomass is not as developed as other biomanufacturing platform technologies. In addition, results of bench-scale research often fail to translate to large-scale or mass production systems. This disconnect may result from trait drift and evolution occurring, through time, in response to unique drivers in each environment, such as cultivation regimes, weather, and pests. Moreover, outdoor and indoor cultivation of microalgae has the potential to impose negative selection pressures, which makes the maintenance of desired traits a challenge. In this context, this review sheds the light on our current understanding of trait drift and evolution in microalgae. We delineate the basics of phenotype plasticity and evolution, with a focus on how microalgae respond under various conditions. In addition, we review techniques that exploit phenotypic plasticity and evolution for strain improvement in view of industrial commercial applications, highlighting associated advantages and shortcomings. Finally, we suggest future research directions and recommendations to overcome unwanted trait drift and evolution in microalgae cultivation.

Proteomics unveil a central role for peroxisomes in butyrate assimilation of the heterotrophic Chlorophyte alga Polytomella sp

Volatile fatty acids found in effluents of the dark fermentation of biowastes can be used for mixotrophic growth of microalgae, improving productivity and reducing the cost of the feedstock. Microalgae can use the acetate in the effluents very well, but butyrate is poorly assimilated and can inhibit growth above 1 gC.L-1. The non-photosynthetic chlorophyte alga Polytomella sp. SAG 198.80 was found to be able to assimilate butyrate fast. To decipher the metabolic pathways implicated in butyrate assimilation, quantitative proteomics study was developed comparing Polytomella sp. cells grown on acetate and butyrate at 1 gC.L-1. After statistical analysis, a total of 1772 proteins were retained, of which 119 proteins were found to be overaccumulated on butyrate vs. only 46 on acetate, indicating that butyrate assimilation necessitates additional metabolic steps. The data show that butyrate assimilation occurs in the peroxisome via the β-oxidation pathway to produce acetyl-CoA and further tri/dicarboxylic acids in the glyoxylate cycle. Concomitantly, reactive oxygen species defense enzymes as well as the branched amino acid degradation pathway were strongly induced. Although no clear dedicated butyrate transport mechanism could be inferred, several membrane transporters induced on butyrate are identified as potential condidates. Metabolic responses correspond globally to the increased needs for central cofactors NAD, ATP and CoA, especially in the peroxisome and the cytosol.

Endogenous giant viruses contribute to intraspecies genomic variability in the model green alga Chlamydomonas reinhardtii

Chlamydomonas reinhardtii is a unicellular eukaryotic alga that has been studied as a model organism for decades. Despite an extensive history as a model system, phylogenetic and genetic characteristics of viruses infecting this alga have remained elusive. We analyzed high-throughput genome sequence data of C. reinhardtii field isolates, and in six we discovered sequences belonging to endogenous giant viruses that reach up to several 100 kb in length. In addition, we have also discovered the entire genome of a closely related giant virus that is endogenized within the genome of Chlamydomonas incerta, the closest sequenced relative of C. reinhardtii. Endogenous giant viruses add hundreds of new gene families to the host strains, highlighting their contribution to the pangenome dynamics and interstrain genomic variability of C. reinhardtii. Our findings suggest that the endogenization of giant viruses may have important implications for structuring the population dynamics and ecology of protists in the environment.

Identification and analysis of smORFs in Chlamydomonas reinhardtii

Small open reading frames (smORFs) have been acknowledged as an important partner in organism functions ranging from bacteria to higher eukaryotes. However, lack of investigation of smORFs in green algae, despite their importance in ecology and evolution. We applied bioinformatic analysis, ribosome profiling, and small peptide proteomics to provide a genome-wide and high-confident smORF database in the model green alga Chlamydomonas reinhardtii. The whole genome was screened first to mine potential coding smORFs. Then conservative analysis, ribosome profiling, and proteomics data were processed to identify conserved smORFs and generate translation evidence. The combination of procedures resulted in 2014 smORFs that might exist in the C. reinhardtii genome. The expression of smORFs in Cd treatment suggested that two smORFs might participate in redox reaction, three in inorganic phosphate transport, and one in DNA repair under stress. Our study built a genome-widely database in C. reinhardtii, providing target smORFs for further research.

The distribution of fitness effects of spontaneous mutations in Chlamydomonas reinhardtii inferred using frequency changes under experimental evolution

The distribution of fitness effects (DFE) for new mutations is fundamental for many aspects of population and quantitative genetics. In this study, we have inferred the DFE in the single-celled alga Chlamydomonas reinhardtii by estimating changes in the frequencies of 254 spontaneous mutations under experimental evolution and equating the frequency changes of linked mutations with their selection coefficients. We generated seven populations of recombinant haplotypes by crossing seven independently derived mutation accumulation lines carrying an average of 36 mutations in the haploid state to a mutation-free strain of the same genotype. We then allowed the populations to evolve under natural selection in the laboratory by serial transfer in liquid culture. We observed substantial and repeatable changes in the frequencies of many groups of linked mutations, and, surprisingly, as many mutations were observed to increase as decrease in frequency. Mutation frequencies were highly repeatable among replicates, suggesting that selection was the cause of the observed allele frequency changes. We developed a Bayesian Monte Carlo Markov Chain method to infer the DFE. This computes the likelihood of the observed distribution of changes of frequency, and obtains the posterior distribution of the selective effects of individual mutations, while assuming a two-sided gamma distribution of effects. We infer that the DFE is a highly leptokurtic distribution, and that approximately equal proportions of mutations have positive and negative effects on fitness. This result is consistent with what we have observed in previous work on a different C. reinhardtii strain, and suggests that a high fraction of new spontaneously arisen mutations are advantageous in a simple laboratory environment.

Sex-linked deubiquitinase establishes uniparental transmission of chloroplast DNA

Most sexual organisms inherit organelles from one parent, commonly by excluding organelles from the smaller gametes. However, post-mating elimination of organelles derived from one gamete ensures uniparental inheritance, where the underlying mechanisms to distinguish organelles by their origin remain obscure. Mating in Chlamydomonas reinhardtii combines isomorphic plus and minus gametes, but chloroplast DNA from minus gametes is selectively degraded in zygotes. Here, we identify OTU2p (otubain protein 2), encoded in the plus mating-type locus MT+, as the protector of plus chloroplast. Otu2p is an otubain-like deubiquitinase, which prevents proteasome-mediated degradation of the preprotein translocase of the outer chloroplast membrane (TOC) during gametogenesis. Using OTU2p-knockouts and proteasome inhibitor treatment, we successfully redirect selective DNA degradation in chloroplasts with reduced TOC levels regardless of mating type, demonstrating that plus-specific Otu2p establishes uniparental chloroplast DNA inheritance. Our work documents that a sex-linked organelle quality control mechanism drives the uniparental organelle inheritance without dimorphic gametes.

Most sexual organisms ensure that organelles are inherited from a single parent. Here, the authors describe OTU2p, a Chlamydomonas deubiquitinase that drives uniparental organelle inheritance without gametic dimorphism by preventing proteasome-mediated degradation exclusively in gametes of the plus mating type.

Green Algal Models for Multicellularity

The repeated evolution of multicellularity across the tree of life has profoundly affected the ecology and evolution of nearly all life on Earth. Many of these origins were in different groups of photosynthetic eukaryotes, or algae. Here, we review the evolution and genetics of multicellularity in several groups of green algae, which include the closest relatives of land plants. These include millimeter-scale, motile spheroids of up to 50,000 cells in the volvocine algae; decimeter-scale seaweeds in the genus Ulva (sea lettuce); and very plantlike, meter-scale freshwater algae in the genus Chara (stoneworts). We also describe algae in the genus Caulerpa, which are giant, multinucleate, morphologically complex single cells. In each case, we review the life cycle, phylogeny, and genetics of traits relevant to the evolution of multicellularity, and genetic and genomic resources available for the group in question. Finally, we suggest routes toward developing these groups as model organisms for the evolution of multicellularity. Expected final online publication date for the Annual Review of Genetics, Volume 55 is November 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

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

BACKGROUND

The volvocine algae, which include the single-celled species Chlamydomonas reinhardtii and the colonial species Volvox carteri, serve as a model in which to study the evolution of multicellularity and cellular differentiation. Studies reconstructing the history of this group have by and large relied on datasets of one to a few genes for phylogenetic inference and ancestral character state reconstruction. As a result, volvocine phylogenies lack concordance depending on the number and/or type of genes (i.e., chloroplast vs nuclear) chosen for phylogenetic inference. While multiple studies suggest that multicellularity evolved only once in the volvocine algae, that each of its three colonial families is monophyletic, and that there have been at least three independent origins of cellular differentiation in the group, other studies call into question one or more of these conclusions. An accurate assessment of the evolutionary history of the volvocine algae requires inference of a more robust phylogeny.

RESULTS

We performed RNA sequencing (RNA-seq) on 55 strains representing 47 volvocine algal species and obtained similar data from curated databases on 13 additional strains. We then compiled a dataset consisting of transcripts for 40 single-copy, protein-coding, nuclear genes and subjected the predicted amino acid sequences of these genes to maximum likelihood, Bayesian inference, and coalescent-based analyses. These analyses show that multicellularity independently evolved at least twice in the volvocine algae and that the colonial family Goniaceae is not monophyletic. Our data further indicate that cellular differentiation arose independently at least four, and possibly as many as six times, within the volvocine algae.

CONCLUSIONS

Altogether, our results demonstrate that multicellularity and cellular differentiation are evolutionarily labile in the volvocine algae, affirming the importance of this group as a model system for the study of major transitions in the history of life.

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.

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.

An ancient clade of Penelope-like retroelements with permuted domains is present in the green lineage and protists, and dominates many invertebrate genomes

Penelope-like elements (PLEs) are an enigmatic clade of retrotransposons whose reverse transcriptases (RTs) share a most recent common ancestor with telomerase RTs. The single ORF of canonical endonuclease (EN)+ PLEs encodes RT and a C-terminal GIY–YIG EN that enables intrachromosomal integration, whereas EN− PLEs lack EN and are generally restricted to chromosome termini. EN+ PLEs have only been found in animals, except for one case of horizontal transfer to conifers, whereas EN− PLEs occur in several kingdoms. Here, we report a new, deep-branching PLE clade with a permuted domain order, whereby an N-terminal GIY–YIG EN is linked to a C-terminal RT by a short domain with a characteristic CxC motif. These N-terminal EN+ PLEs share a structural organization, including pseudo-LTRs and complex tandem/inverted insertions, with canonical EN+ PLEs from Penelope/Poseidon, Neptune, and Nematis clades, and show insertion bias for microsatellites, but lack canonical hammerhead ribozyme motifs. However, their phylogenetic distribution is much broader. The Naiads, found in numerous invertebrate phyla, can reach tens of thousands of copies per genome. In spiders and clams, Naiads independently evolved to encode selenoproteins containing multiple selenocysteines. Chlamys, which lack the CCHH motif universal to PLE ENs, occur in green algae, spike mosses (targeting ribosomal DNA), and slime molds. Unlike canonical PLEs, RTs of N-terminal EN+ PLEs contain the insertion-in-fingers domain (IFD), strengthening the link between PLEs and telomerases. Additionally, we describe Hydra, a novel metazoan C-terminal EN+ clade. Overall, we conclude that PLE diversity, taxonomic distribution, and abundance are comparable with non-LTR and LTR-retrotransposons.

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.