Spontaneous mutation rate as a source of diversity for improving desirable traits in cultured microalgae

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.

Salt stress alters the spectrum of de novo mutation available to selection during experimental adaptation of Chlamydomonas reinhardtii

The genetic basis of adaptation is driven by both selection and the spectrum of available mutations. Given that the rate of mutation is not uniformly distributed across the genome and varies depending on the environment, understanding the signatures of selection across the genome is aided by first establishing what the expectations of genetic change are from mutation. To determine the interaction between salt stress, selection, and mutation across the genome, we compared mutations observed in a selection experiment for salt tolerance in Chlamydomonas reinhardtii to those observed in mutation accumulation (MA) experiments with and without salt exposure. MA lines evolved under salt stress had a single-nucleotide mutation rate of 1.1 × 10 - 9 $1.1 \times 10^{-9}$ , similar to that of MA lines under standard conditions ( 9.6 × 10 - 10 $9.6 \times 10^{-10}$ ). However, we found that salt stress led to an increased rate of indel mutations, but that many of these mutations were removed under selection. Finally, lines adapted to salt also showed excess clustering of mutations in the genome and the co-expression network, suggesting a role for positive selection in retaining mutations in particular compartments of the genome during the evolution of salt tolerance. Our study shows that characterizing mutation rates and spectra expected under stress helps disentangle the effects of environment and selection during adaptation.

Random Mutagenesis as a Promising Tool for Microalgal Strain Improvement towards Industrial Production

Microalgae have become a promising novel and sustainable feedstock for meeting the rising demand for food and feed. However, microalgae-based products are currently hindered by high production costs. One major reason for this is that commonly cultivated wildtype strains do not possess the robustness and productivity required for successful industrial production. Several strain improvement technologies have been developed towards creating more stress tolerant and productive strains. While classical methods of forward genetics have been extensively used to determine gene function of randomly generated mutants, reverse genetics has been explored to generate specific mutations and target phenotypes. Site-directed mutagenesis can be accomplished by employing different gene editing tools, which enable the generation of tailor-made genotypes. Nevertheless, strategies promoting the selection of randomly generated mutants avoid the introduction of foreign genetic material. In this paper, we review different microalgal strain improvement approaches and their applications, with a primary focus on random mutagenesis. Current challenges hampering strain improvement, selection, and commercialization will be discussed. The combination of these approaches with high-throughput technologies, such as fluorescence-activated cell sorting, as tools to select the most promising mutants, will also be discussed.

Strain Development in Microalgal Biotechnology-Random Mutagenesis Techniques

Microalgal biomass and metabolites can be used as a renewable source of nutrition, pharmaceuticals and energy to maintain or improve the quality of human life. Microalgae’s high volumetric productivity and low impact on the environment make them a promising raw material in terms of both ecology and economics. To optimize biotechnological processes with microalgae, improving the productivity and robustness of the cell factories is a major step towards economically viable bioprocesses. This review provides an overview of random mutagenesis techniques that are applied to microalgal cell factories, with a particular focus on physical and chemical mutagens, mutagenesis conditions and mutant characteristics.

Codon Usage Bias in Phytoplankton

Non-random usage of synonymous codons, known as “codon bias”, has been described in many organisms, from bacteria to Drosophila, but little is known about it in phytoplankton. This phenomenon is thought to be driven by selection for translational efficiency. As the efficacy of selection is proportional to the effective population size, species with large population sizes, such as phytoplankton, are expected to have strong codon bias. To test this, we measured codon bias in 215 strains from Haptophyta, Chlorophyta, Ochrophyta (except diatoms that were studied previously), Dinophyta, Cryptophyta, Ciliophora, unicellular Rhodophyta and Chlorarachniophyta. Codon bias is modest in most groups, despite the astronomically large population sizes of marine phytoplankton. The strength of the codon bias, measured with the effective number of codons, is the strongest in Haptophyta and the weakest in Chlorarachniophyta. The optimal codons are GC-ending in most cases, but several shifts to AT-ending codons were observed (mainly in Ochrophyta and Ciliophora). As it takes a long time to reach a new equilibrium after such shifts, species having AT-ending codons show a lower frequency of optimal codons compared to other species. Genetic diversity, calculated for species with more than three strains sequenced, is modest, indicating that the effective population sizes are many orders of magnitude lower than the astronomically large census population sizes, which helps to explain the modest codon bias in marine phytoplankton. This study represents the first comparative analysis of codon bias across multiple major phytoplankton groups.

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.

Evolution of Mutation Rate in Astronomically Large Phytoplankton Populations

Genetic diversity is expected to be proportional to population size, yet, there is a well-known, but unexplained lack of genetic diversity in large populations-the "Lewontin's paradox." Larger populations are expected to evolve lower mutation rates, which may help to explain this paradox. Here, we test this conjecture by measuring the spontaneous mutation rate in a ubiquitous unicellular marine phytoplankton species Emiliania huxleyi (Haptophyta) that has modest genetic diversity despite an astronomically large population size. Genome sequencing of E. huxleyi mutation accumulation lines revealed 455 mutations, with an unusual GC-biased mutation spectrum. This yielded an estimate of the per site mutation rate µ = 5.55×10-10 (CI 95%: 5.05×10-10 - 6.09×10-10), which corresponds to an effective population size Ne ∼ 2.7×106. Such a modest Ne is surprising for a ubiquitous and abundant species that accounts for up to 10% of global primary productivity in the oceans. Our results indicate that even exceptionally large populations do not evolve mutation rates lower than ∼10-10 per nucleotide per cell division. Consequently, the extreme disparity between modest genetic diversity and astronomically large population size in the plankton species cannot be explained by an unusually low mutation rate.

Light spectra as triggers for sorting improved strains of Tisochrysis lutea

It is known that microalgae respond to different light colors, but not at single-cell level. This work aimed to assess if different light colors could be used as triggers to sort over-producing cells. Six light spectra were used: red + green + blue (RGBL), blue (BL), red (RL), green (GL), blue + red (BRL) and blue + green (BGL). Fluorescence-activated cell sorting method was used to analyse single-cell fluorescence and sort cells. BGL and RGBL lead to the highest fucoxanthin production, while RL showed the lowest. Therefore, it was hypothesized that hyper-producing cells can be isolated efficiently under the adverse condition (RL). After exposure to all light colors for 14 days, the top 1% fucoxanthin producing cells were sorted. A sorted strain from RL showed higher (16-19%) growth rate and fucoxanthin productivity. This study showed how light spectra affected single-cell fucoxanthin and lipid contents and productivities. Also, it supplied an approach to sort for high-fucoxanthin or high-lipid cells.

First Estimation of the Spontaneous Mutation Rate in Diatoms

Mutations are the origin of genetic diversity, and the mutation rate is a fundamental parameter to understand all aspects of molecular evolution. The combination of mutation-accumulation experiments and high-throughput sequencing enabled the estimation of mutation rates in most model organisms, but several major eukaryotic lineages remain unexplored. Here, we report the first estimation of the spontaneous mutation rate in a model unicellular eukaryote from the Stramenopile kingdom, the diatom Phaeodactylum tricornutum (strain RCC2967). We sequenced 36 mutation accumulation lines for an average of 181 generations per line and identified 156 de novo mutations. The base substitution mutation rate per site per generation is μbs = 4.77 × 10-10 and the insertion-deletion mutation rate is μid = 1.58 × 10-11. The mutation rate varies as a function of the nucleotide context and is biased toward an excess of mutations from GC to AT, consistent with previous observations in other species. Interestingly, the mutation rates between the genomes of organelles and the nucleus differ, with a significantly higher mutation rate in the mitochondria. This confirms previous claims based on indirect estimations of the mutation rate in mitochondria of photosynthetic eukaryotes that acquired their plastid through a secondary endosymbiosis. This novel estimate enables us to infer the effective population size of P. tricornutum to be Ne∼8.72 × 106.

Evolution of Codon Usage Bias in Diatoms

Codon usage bias (CUB)-preferential use of one of the synonymous codons, has been described in a wide range of organisms from bacteria to mammals, but it has not yet been studied in marine phytoplankton. CUB is thought to be caused by weak selection for translational accuracy and efficiency. Weak selection can overpower genetic drift only in species with large effective population sizes, such as Drosophila that has relatively strong CUB, while organisms with smaller population sizes (e.g., mammals) have weak CUB. Marine plankton species tend to have extremely large populations, suggesting that CUB should be very strong. Here we test this prediction and describe the patterns of codon usage in a wide range of diatom species belonging to 35 genera from 4 classes. We report that most of the diatom species studied have surprisingly modest CUB (mean Effective Number of Codons, ENC = 56), with some exceptions showing stronger codon bias (ENC = 44). Modest codon bias in most studied diatom species may reflect extreme disparity between astronomically large census and modest effective population size (N_e), with fluctuations in population size and linked selection limiting long-term N_e and rendering selection for optimal codons less efficient. For example, genetic diversity (pi ~0.02 at silent sites) in Skeletonema marinoi corresponds to N_e of about 10 million individuals, which is likely many orders of magnitude lower than its census size. Still, N_e ~10⁷ should be large enough to make selection for optimal codons efficient. Thus, we propose that an alternative process-frequent changes of preferred codons, may be a more plausible reason for low CUB despite highly efficient selection for preferred codons in diatom populations. The shifts in the set of optimal codons should result in the changes of the direction of selection for codon usage, so the actual codon usage never catches up with the moving target of the optimal set of codons and the species never develop strong CUB. Indeed, we detected strong shifts in preferential codon usage within some diatom genera, with switches between preferentially GC-rich and AT-rich 3^rd codon positions (GC3). For example, GC3 ranges from 0.6 to 1 in most Chaetoceros species, while for Chaetoceros dichaeta GC3 = 0.1. Both variation in selection intensity and mutation spectrum may drive such shifts in codon usage and limit the observed CUB. Our study represents the first genome-wide analysis of CUB in diatoms and the first such analysis for a major phytoplankton group.