Whole-Genome Resequencing Reveals Extensive Natural Variation in the Model Green Alga Chlamydomonas reinhardtii

From nets to networks: tools for deciphering phytoplankton metabolic interactions within communities and their global significance

Our oceans are populated with a wide diversity of planktonic organisms that form complex dynamic communities at the base of marine trophic networks. Within such communities are phytoplankton, unicellular photosynthetic taxa that provide an estimated half of global primary production and support biogeochemical cycles, along with other essential ecosystem services. One of the major challenges for microbial ecologists has been to try to make sense of this complexity. While phytoplankton distributions can be well explained by abiotic factors such as temperature and nutrient availability, there is increasing evidence that their ecological roles are tightly linked to their metabolic interactions with other plankton members through complex mechanisms (e.g. competition and symbiosis). Therefore, unravelling phytoplankton metabolic interactions is the key for inferring their dependency on, or antagonism with, other taxa and better integrating them into the context of carbon and nutrient fluxes in marine trophic networks. In this review, we attempt to summarize the current knowledge brought by ecophysiology, organismal imaging, in silico predictions and co-occurrence networks using 'omics data, highlighting successful combinations of approaches that may be helpful for future investigations of phytoplankton metabolic interactions within their complex communities.This article is part of the theme issue 'Connected interactions: enriching food web research by spatial and social interactions'.

Calredoxin regulates the chloroplast NADPH-dependent thioredoxin reductase in Chlamydomonas reinhardtii

Calredoxin (CRX) is a calcium (Ca2+)-dependent thioredoxin (TRX) in the chloroplast of Chlamydomonas (Chlamydomonas reinhardtii) with a largely unclear physiological role. We elucidated the CRX functionality by performing in-depth quantitative proteomics of wild-type cells compared with a crx insertional mutant (IMcrx), two CRISPR/Cas9 KO mutants, and CRX rescues. These analyses revealed that the chloroplast NADPH-dependent TRX reductase (NTRC) is co-regulated with CRX. Electron transfer measurements revealed that CRX inhibits NADPH-dependent reduction of oxidized chloroplast 2-Cys peroxiredoxin (PRX1) via NTRC and that the function of the NADPH-NTRC complex is under strict control of CRX. Via non-reducing SDS-PAGE assays and mass spectrometry, our data also demonstrated that PRX1 is more oxidized under high light (HL) conditions in the absence of CRX. The redox tuning of PRX1 and control of the NADPH-NTRC complex via CRX interconnect redox control with active photosynthetic electron transport and metabolism, as well as Ca2+ signaling. In this way, an economic use of NADPH for PRX1 reduction is ensured. The finding that the absence of CRX under HL conditions severely inhibited light-driven CO2 fixation underpins the importance of CRX for redox tuning, as well as for efficient photosynthesis.

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

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

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.

Effects of antimony exposure on DNA damage and genome-wide variation in zebrafish (Danio rerio) liver

Antimony (Sb) is a potentially toxic and carcinogenic cumulative contaminant that poses a serious threat to aquatic ecosystems. To better clarify the genotoxicity of Sb and its mechanism of action. In this study, we investigated DNA damage and genome-wide variation in the liver of a model organism, zebrafish (Danio rerio), under subacute Sb exposure and explored its potential toxicological mechanisms. The results showed that medium and high concentrations of Sb significantly reduced the total antioxidant capacity and increased the content of reactive oxygen species in zebrafish liver, and further studies revealed that it increased oxidative DNA damage and DNA-DNA cross-link (DDC), but had little effect on DNA-protein cross-link (DPC). The result of resequencing showed that the mutation sites of the genes with high concentrations of Sb were higher than those with medium concentrations, and the mutation was mainly a single nucleotide. The pathways significantly enriched for nonsynonymous single nucleotide polymorphisms (SNPs) and insertion/deletion mutations (InDels) variant genes in the coding regions of both the medium and high Sb-treated groups were ECM-receptor interactions, and the high Sb-treated group also included lysine degradation, hematopoietic cell lineage, and cytokine-cytokine receptor interactions. This suggests that ECM-receptor interactions play an important role in the mechanism of antimony toxicity to the liver of zebrafish.

Microalgae share key features with human erythrocytes and can safely circulate through the vascular system in mice

As animal cells cannot produce oxygen, erythrocytes are responsible for gas interchange, being able to capture and deliver oxygen upon tissue request. Interestingly, several other cells in nature produce oxygen by photosynthesis, raising the question of whether they could circulate within the vascular networks, acting as an alternative source for oxygen delivery. To address this long-term goal, here some physical and mechanical features of the photosynthetic microalga Chlamydomona reinhardtii were studied and compared with erythrocytes, revealing that both exhibit similar size and rheological properties. Moreover, key biocompatibility aspects of the microalgae were evaluated in vitro and in vivo, showing that C. reinhardtii can be co-cultured with endothelial cells, without affecting each other's morphology and viability. Moreover, short-term systemic perfusion of the microalgae showed a thoroughly intravascular distribution in mice. Finally, the systemic injection of high numbers of microalgae did not trigger deleterious responses in living mice. Altogether, this work provides key scientific insights to support the notion that photosynthetic oxygenation could be achieved by circulating microalgae, representing another important step towards human photosynthesis. KEY POINTS: • C. reinhardtii and endothelial cells are biocompatible in vitro. • C. reinhardtii distribute throughout the entire vasculature after mice perfusion. • C. reinhardtii do not trigger deleterious responses after injection in mice.

Palmelloid Formation and Cell Aggregation Are Essential Mechanisms for High Light Tolerance in a Natural Strain of Chlamydomonas reinhardtii

Photosynthetic organisms, such as higher plants and algae, require light to survive. However, an excessive amount of light can be harmful due to the production of reactive oxygen species (ROS), which cause cell damage and, if it is not effectively regulated, cell death. The study of plants' responses to light can aid in the development of methods to improve plants' growth and productivity. Due to the multicellular nature of plants, there may be variations in the results based on plant age and tissue type. Chlamydomonas reinhardtii, a unicellular green alga, has also been used as a model organism to study photosynthesis and photoprotection. Nonetheless, the majority of the research has been conducted with strains that have been consistently utilized in laboratories and originated from the same source. Despite the availability of many field isolates of this species, very few studies have compared the light responses of field isolates. This study examined the responses of two field isolates of Chlamydomonas to high light stress. The light-tolerant strain, CC-4414, managed reactive oxygen species (ROS) slightly better than the sensitive strain, CC-2344, did. The proteomic data of cells subjected to high light revealed cellular modifications of the light-tolerant strain toward membrane proteins. The morphology of cells under light stress revealed that this strain utilized the formation of palmelloid structures and cell aggregation to shield cells from excessive light. As indicated by proteome data, morphological modifications occur simultaneously with the increase in protein degradation and autophagy. By protecting cells from stress, cells are able to continue to upregulate ROS management mechanisms and prevent cell death. This is the first report of palmelloid formation in Chlamydomonas under high light stress.

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.

Pigment Production of Chlamydomonas Strains in Response to Norflurazon and ZnO Nanoparticles

Numerous species of microalgae have been utilized for pigment production. More and more species are gaining popularity due to their ability to accumulate pigments with varying chemical compositions and the fact that some have distinctive byproducts that can be co-produced. Despite the fact that many of the species have unique by-products and traits, they are not being used economically due to high production costs. Utilizing agricultural and industrial wastewater for algae cultivation is one way to lower manufacturing costs. Herbicide-contaminated wastewater can result from agricultural contamination. Norflurazon is a popular pesticide frequently used for weed control. The presence of norflurazon in water renders that water unusable and requires proper treatment. Nanoparticles of ZnO (ZnO NPs), on the other hand, are utilized in a variety of industrial productions of numerous household goods. Water contaminated with ZnO NPs can present potential risks to human health and the environment. In this study, two field isolates of the green microalga Chlamydomonas reinhardtii, a widely used model organism, were examined for their reaction to these two compounds in order to assess the responses of different natural strains to environmental stresses. Norflurazon at 10 µM had a higher inhibitory effect on growth and pigment production than ZnO NPs at 200 mg L−1. Although both norflurazon and ZnO NPs inhibit cell growth and pigmentation, they do so through distinct processes. Norflurazon induces oxidative stress in cells, resulting in photosystem damage. ZnO nanoparticles, on the other hand, did not cause photosystem damage but rather mechanical cell damage and disintegration. In addition, the physiological responses of the two Chlamydomonas strains were distinct, supporting the utilization of natural algal strains for specific types of environmental pollutants.

Studying Whole-Genome Duplication Using Experimental Evolution of Chlamydomonas

In this chapter, we present the use of Chlamydomonas reinhardtii in experiments designed to study the evolutionary impacts of whole genome duplication. We shortly introduce the algal species and depict why it is an excellent model for experimental evolution. Subsequently, we discuss the most relevant steps and methods in the design of a ploidy-related Chlamydomonas experiment. These steps include strain selection, ploidy determination, different methods of making diplo- and polyploid Chlamydomonas cells, replication, culturing conditions, preservation, and the ways to quantify phenotypic and genotypic change.

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.

Morphological, physiological, and transcriptional responses of the freshwater diatom Fragilaria crotonensis to elevated pH conditions

Harmful algal blooms (HABs) caused by the toxin-producing cyanobacteria Microcystis spp., can increase water column pH. While the effect(s) of these basified conditions on the bloom formers are a high research priority, how these pH shifts affect other biota remains understudied. Recently, it was shown these high pH levels decrease growth and Si deposition rates in the freshwater diatom Fragilaria crotonensis and natural Lake Erie (Canada-US) diatom populations. However, the physiological mechanisms and transcriptional responses of diatoms associated with these observations remain to be documented. Here, we examined F. crotonensis with a set of morphological, physiological, and transcriptomic tools to identify cellular responses to high pH. We suggest 2 potential mechanisms that may contribute to morphological and physiological pH effects observed in F. crotonensis. Moreover, we identified a significant upregulation of mobile genetic elements in the F. crotonensis genome which appear to be an extreme transcriptional response to this abiotic stress to enhance cellular evolution rates-a process we have termed "genomic roulette." We discuss the ecological and biogeochemical effects high pH conditions impose on fresh waters and suggest a means by which freshwater diatoms such as F. crotonensis may evade high pH stress to survive in a "basified" future.

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.

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.

Light chain 2 is a Tctex-type related axonemal dynein light chain that regulates directional ciliary motility in Trypanosoma brucei

Flagellar motility is essential for the cell morphology, viability, and virulence of pathogenic kinetoplastids. Trypanosoma brucei flagella beat with a bending wave that propagates from the flagellum's tip to its base, rather than base-to-tip as in other eukaryotes. Thousands of dynein motor proteins coordinate their activity to drive ciliary bending wave propagation. Dynein-associated light and intermediate chains regulate the biophysical mechanisms of axonemal dynein. Tctex-type outer arm dynein light chain 2 (LC2) regulates flagellar bending wave propagation direction, amplitude, and frequency in Chlamydomonas reinhardtii. However, the role of Tctex-type light chains in regulating T. brucei motility is unknown. Here, we used a combination of bioinformatics, in-situ molecular tagging, and immunofluorescence microscopy to identify a Tctex-type light chain in the procyclic form of T. brucei (TbLC2). We knocked down TbLC2 expression using RNAi in both wild-type and FLAM3, a flagellar attachment zone protein, knockdown cells and quantified TbLC2's effects on trypanosome cell biology and biophysics. We found that TbLC2 knockdown reduced the directional persistence of trypanosome cell swimming, induced an asymmetric ciliary bending waveform, modulated the bias between the base-to-tip and tip-to-base beating modes, and increased the beating frequency. Together, our findings are consistent with a model of TbLC2 as a down-regulator of axonemal dynein activity that stabilizes the forward tip-to-base beating ciliary waveform characteristic of trypanosome cells. Our work sheds light on axonemal dynein regulation mechanisms that contribute to pathogenic kinetoplastids' unique tip-to-base ciliary beating nature and how those mechanisms underlie dynein-driven ciliary motility more generally.

A YAK1-type protein kinase, triacylglycerol accumulation regulator 1, in the green alga Chlamydomonas reinhardtii is a potential regulator of cell division and differentiation into gametes during photoautotrophic nitrogen deficiency

Yet another kinase (YAK) 1 is a conserved eukaryotic protein kinase coordinating growth and development. We previously isolated a mutant of Chlamydomonas reinhardtii defective in the YAK1 ortholog triacylglycerol (TAG) accumulation regulator 1 (TAR1). The mutant tar1-1 displayed higher levels of chlorophyll, starch, TAG, and biomass than the parental strain C9 (renamed as C9-3) in photoautotrophic nitrogen (N)-deficient conditions. However, we found that the parental C9-3 showed faster chlorosis upon N-deficiency than the original C9 (C9-1) freshly recovered from cryopreservation, suggesting that C9-3 had acquired particular characteristics during long-term subculturing. To exclude phenotypes dependent on a particular parental strain, we newly created tar1 mutants from two wild-types, C9-1 and CC 125. Like tar1-1, the new tar1 mutants showed higher levels of chlorophyll and TAG/starch than the parental strain. Upon removal of N, Chlamydomonas cells divide once before ceasing further division. Previously, the single division after N-removal was arrested in tar1-1 in photomixotrophic conditions, but this phenotype was not observed in photoautotrophic conditions because of the particular characteristics of the parental C9-3. However, using C9- 1 and CC-125 as parental strains, we showed that cell division after N-removal was impaired in new tar1 mutants in photoautotrophic conditions. Consistent with the view that the division under N-deficiency is necessary for gametic differentiation, new tar1 mutants showed lower mating efficiency than the parental strains. Taken together, TAR1 was suggested to promote differentiation into gametes through the regulation of cell division in response to N-deficiency.

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.

Nonsense Mutations in Eukaryotes

Nonsense mutations are a type of mutations which results in a premature termination codon occurrence. In general, these mutations have been considered to be among the most harmful ones which lead to premature protein translation termination and result in shortened nonfunctional polypeptide. However, there is evidence that not all nonsense mutations are harmful as well as some molecular mechanisms exist which allow to avoid pathogenic effects of these mutations. This review addresses relevant information on nonsense mutations in eukaryotic genomes, characteristics of these mutations, and different molecular mechanisms preventing or mitigating harmful effects thereof.

Comparative Phenotyping of Two Commonly Used Chlamydomonas reinhardtii Background Strains: CC-1690 (21gr) and CC-5325 (The CLiP Mutant Library Background)

The unicellular green alga Chlamydomonas reinhardtii is an excellent model organism to investigate many essential cellular processes in photosynthetic eukaryotes. Two commonly used background strains of Chlamydomonas are CC-1690 and CC-5325. CC-1690, also called 21gr, has been used for the Chlamydomonas genome project and several transcriptome analyses. CC-5325 is the background strain for the Chlamydomonas Library Project (CLiP). Photosynthetic performance in CC-5325 has not been evaluated in comparison with CC-1690. Additionally, CC-5325 is often considered to be cell-wall deficient, although detailed analysis is missing. The circadian rhythms in CC-5325 are also unclear. To fill these knowledge gaps and facilitate the use of the CLiP mutant library for various screens, we performed phenotypic comparisons between CC-1690 and CC-5325. Our results showed that CC-5325 grew faster heterotrophically in dark and equally well in mixotrophic liquid medium as compared to CC-1690. CC-5325 had lower photosynthetic efficiency and was more heat-sensitive than CC-1690. Furthermore, CC-5325 had an intact cell wall which had comparable integrity to that in CC-1690 but appeared to have reduced thickness. Additionally, CC-5325 could perform phototaxis, but could not maintain a sustained circadian rhythm of phototaxis as CC1690 did. Finally, in comparison to CC-1690, CC-5325 had longer cilia in the medium with acetate but slower swimming speed in the medium without nitrogen and acetate. Our results will be useful for researchers in the Chlamydomonas community to choose suitable background strains for mutant analysis and employ the CLiP mutant library for genome-wide mutant screens under appropriate conditions, especially in the areas of photosynthesis, thermotolerance, cell wall, and circadian rhythms.

Fourth generation biofuel from genetically modified algal biomass: Challenges and future directions

Genetic engineering applications in the field of biofuel are rapidly expanding due to their potential to boost biomass productivity while lowering its cost and enhancing its quality. Recently, fourth-generation biofuel (FGB), which is biofuel obtained from genetically modified (GM) algae biomass, has gained considerable attention from academic and industrial communities. However, replacing fossil resources with FGB is still beset with many challenges. Most notably, technical aspects of genetic modification operations need to be more fully articulated and elaborated. However, relatively little attention has been paid to GM algal biomass. There is a limited number of reviews on the progress and challenges faced in the algal genetics of FGB. Therefore, the present review aims to fill this gap in the literature by recapitulating the findings of recent studies and achievements on safe and efficient genetic manipulation in the production of FGB. Then, the essential issues and parameters related to genome editing in algal strains are highlighted. Finally, the main challenges to FGB pertaining to the diffusion risk and regulatory frameworks are addressed. This review concluded that the technical and biosafety aspects of FGB, as well as the complexity and diversity of the related regulations, legitimacy concerns, and health and environmental risks, are among the most important challenges that require a strong commitment at the national/international levels to reach a global consensus.

Large-scale genome sequencing reveals the driving forces of viruses in microalgal evolution

Being integral primary producers in diverse ecosystems, microalgal genomes could be mined for ecological insights, but representative genome sequences are lacking for many phyla. We cultured and sequenced 107 microalgae species from 11 different phyla indigenous to varied geographies and climates. This collection was used to resolve genomic differences between saltwater and freshwater microalgae. Freshwater species showed domain-centric ontology enrichment for nuclear and nuclear membrane functions, while saltwater species were enriched in organellar and cellular membrane functions. Further, marine species contained significantly more viral families in their genomes (p = 8e-4). Sequences from Chlorovirus, Coccolithovirus, Pandoravirus, Marseillevirus, Tupanvirus, and other viruses were found integrated into the genomes of algal from marine environments. These viral-origin sequences were found to be expressed and code for a wide variety of functions. Together, this study comprehensively defines the expanse of protein-coding and viral elements in microalgal genomes and posits a unified adaptive strategy for algal halotolerance.

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

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

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

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

Recombination Rate Variation and Infrequent Sex Influence Genetic Diversity in Chlamydomonas reinhardtii

Recombination confers a major evolutionary advantage by breaking up linkage disequilibrium between harmful and beneficial mutations, thereby facilitating selection. However, in species that are only periodically sexual, such as many microbial eukaryotes, the realized rate of recombination is also affected by the frequency of sex, meaning that infrequent sex can increase the effects of selection at linked sites despite high recombination rates. Despite this, the rate of sex of most facultatively sexual species is unknown. Here, we use genomewide patterns of linkage disequilibrium to infer fine-scale recombination rate variation in the genome of the facultatively sexual green alga Chlamydomonas reinhardtii. We observe recombination rate variation of up to two orders of magnitude and find evidence of recombination hotspots across the genome. Recombination rate is highest flanking genes, consistent with trends observed in other nonmammalian organisms, though intergenic recombination rates vary by intergenic tract length. We also find a positive relationship between nucleotide diversity and physical recombination rate, suggesting a widespread influence of selection at linked sites in the genome. Finally, we use estimates of the effective rate of recombination to calculate the rate of sex that occurs in natural populations, estimating a sexual cycle roughly every 840 generations. We argue that the relatively infrequent rate of sex and large effective population size creates a population genetic environment that increases the influence of selection on linked sites across the genome.

Highly Contiguous Nanopore Genome Assembly of Chlamydomonas reinhardtii CC-1690

The current Chlamydomonas reinhardtii reference genome remains fragmented due to gaps stemming from large repetitive regions. To overcome the vast majority of these gaps, publicly available Oxford Nanopore Technology data were used to create a new reference-quality de novo genome assembly containing only 21 contigs, 30/34 telomeric ends, and a genome size of 111 Mb.

The current Chlamydomonas reinhardtii reference genome remains fragmented due to gaps stemming from large repetitive regions. To overcome the vast majority of these gaps, publicly available Oxford Nanopore Technology data were used to create a new reference-quality de novo genome assembly containing only 21 contigs, 30/34 telomeric ends, and a genome size of 111 Mb.

Single Nucleotide Polymorphism Charting of P. patens Reveals Accumulation of Somatic Mutations During in vitro Culture on the Scale of Natural Variation by Selfing

Introduction: Physcomitrium patens (Hedw.) Mitten (previously known as Physcomitrella patens) was collected by H.L.K. Whitehouse in Gransden Wood (Huntingdonshire, United Kingdom) in 1962 and distributed across the globe starting in 1974. Hence, the Gransden accession has been cultured in vitro in laboratories for half a century. Today, there are more than 13 different pedigrees derived from the original accession. Additionally, accessions from other sites worldwide were collected during the last decades. Methods and Results: In this study, 250 high throughput RNA sequencing (RNA-seq) samples and 25 gDNA samples were used to detect single nucleotide polymorphisms (SNPs). Analyses were performed using five different P. patens accessions and 13 different Gransden pedigrees. SNPs were overlaid with metadata and known phenotypic variations. Unique SNPs defining Gransden pedigrees and accessions were identified and experimentally confirmed. They can be successfully employed for PCR-based identification. Conclusion: We show independent mutations in different Gransden laboratory pedigrees, demonstrating that somatic mutations occur and accumulate during in vitro culture. The frequency of such mutations is similar to those observed in naturally occurring populations. We present evidence that vegetative propagation leads to accumulation of deleterious mutations, and that sexual reproduction purges those. Unique SNP sets for five different P. patens accessions were isolated and can be used to determine individual accessions as well as Gransden pedigrees. Based on that, laboratory methods to easily determine P. patens accessions and Gransden pedigrees are presented.

The Seminavis robusta genome provides insights into the evolutionary adaptations of benthic diatoms

Benthic diatoms are the main primary producers in shallow freshwater and coastal environments, fulfilling important ecological functions such as nutrient cycling and sediment stabilization. However, little is known about their evolutionary adaptations to these highly structured but heterogeneous environments. Here, we report a reference genome for the marine biofilm-forming diatom Seminavis robusta, showing that gene family expansions are responsible for a quarter of all 36,254 protein-coding genes. Tandem duplications play a key role in extending the repertoire of specific gene functions, including light and oxygen sensing, which are probably central for its adaptation to benthic habitats. Genes differentially expressed during interactions with bacteria are strongly conserved in other benthic diatoms while many species-specific genes are strongly upregulated during sexual reproduction. Combined with re-sequencing data from 48 strains, our results offer insights into the genetic diversity and gene functions in benthic diatoms.

Convergent Loss of an EDS1/PAD4 Signaling Pathway in Several Plant Lineages Reveals Coevolved Components of Plant Immunity and Drought Response

Plant innate immunity relies on nucleotide binding leucine-rich repeat receptors (NLRs) that recognize pathogen-derived molecules and activate downstream signaling pathways. We analyzed the variation in NLR gene copy number and identified plants with a low number of NLR genes relative to sister species. We specifically focused on four plants from two distinct lineages, one monocot lineage (Alismatales) and one eudicot lineage (Lentibulariaceae). In these lineages, the loss of NLR genes coincides with loss of the well-known downstream immune signaling complex ENHANCED DISEASE SUSCEPTIBILITY 1 (EDS1)/PHYTOALEXIN DEFICIENT 4 (PAD4). We expanded our analysis across whole proteomes and found that other characterized immune genes were absent only in Lentibulariaceae and Alismatales. Additionally, we identified genes of unknown function that were convergently lost together with EDS1/PAD4 in five plant species. Gene expression analyses in Arabidopsis (Arabidopsis thaliana) and Oryza sativa revealed that several homologs of the candidates are differentially expressed during pathogen infection, drought, and abscisic acid treatment. Our analysis provides evolutionary evidence for the rewiring of plant immunity in some plant lineages, as well as the coevolution of the EDS1/PAD4 pathway and drought responses.

Ubiquity or not ubiquity: That is the question

The nature of population structure in eukaryotic microbes has been the subject of intense debate, but until recently the tools to test these hypotheses were either problematic (e.g., allozymes that cannot detect all genetic changes) or beyond financial and technological limits of most laboratories (e.g., high throughput sequencing). In a recent issue of Molecular Ecology, Craig et al. (2019) use a genomic approach to investigate the population structure of a model alga, the chlorophyte Chlamydomonas reinhardtii (Figure 1). Using high throughput sequencing, read mapping, and variant calling, they detected strong signals of differentiation at a continental scale, while local patterns of admixture were complex. Population genomic techniques such as these have not been used extensively in studies of microbial eukaryotes and the fields of conservation genetics and evolution stand to benefit vastly from the adoption of these techniques to studies of diverse protist lineages.

Genomic diversity analysis and identification of novel SSR markers in four tobacco varieties by high-throughput resequencing

Genome resequencing was carried out on two varieties of flue-cured tobacco (LY1306 and Qinyan 96), one variety of sun-cured tobacco (Wanmao 3), and one variety of air-cured Maryland tobacco (Wufeng 1), for a comparative analysis of genomic variation across the four varieties. Single nucleotide polymorphisms (SNPs), insertions and deletions (InDels), structural variations (SVs), and copy-number variations (CNVs) were then identified in each tobacco variety. Furthermore, a functional analysis of mutated genes was carried out. Through in-depth comparative analysis of genomes of different tobacco varieties, we identified genome variations in a number of SNPs, InDels, SVs, and CNVs, respectively. Computational analysis to predict the function of mutated genes containing these differential SNPs, InDels, SVs, and CNVs showed that they were mainly involved in different functions, such as carbohydrate metabolism and secondary metabolites biosynthesis. We mainly focused on genes that were involved in the biosynthesis of secondary metabolites and nicotine metabolism. In addition, we identified five simple sequence repeat (SSR)-based markers and verified them by PCR amplification in 10 tobacco varieties. Taken together, our study increases the understanding of genetic differences between tobacco types or varieties and identifies five SSR markers to classify tobacco varieties or types.

A genomics approach reveals the global genetic polymorphism, structure, and functional diversity of ten accessions of the marine model diatom Phaeodactylum tricornutum

Diatoms emerged in the Mesozoic period and presently constitute one of the main primary producers in the world's ocean and are of a major economic importance. In the current study, using whole genome sequencing of ten accessions of the model diatom Phaeodactylum tricornutum, sampled at broad geospatial and temporal scales, we draw a comprehensive landscape of the genomic diversity within the species. We describe strong genetic subdivisions of the accessions into four genetic clades (A-D) with constituent populations of each clade possessing a conserved genetic and functional makeup, likely a consequence of the limited dispersal of P. tricornutum in the open ocean. We further suggest dominance of asexual reproduction across all the populations, as implied by high linkage disequilibrium. Finally, we show limited yet compelling signatures of genetic and functional convergence inducing changes in the selection pressure on many genes and metabolic pathways. We propose these findings to have significant implications for understanding the genetic structure of diatom populations in nature and provide a framework to assess the genomic underpinnings of their ecological success and impact on aquatic ecosystems where they play a major role. Our work provides valuable resources for functional genomics and for exploiting the biotechnological potential of this model diatom species.

Consequences of recombination for the evolution of the mating type locus in Chlamydomonas reinhardtii

Recombination suppression in sex chromosomes and mating type loci can lead to degeneration as a result of reduced selection efficacy and Muller's ratchet effects. However, genetic exchange in the form of noncrossover gene conversions may still take place within crossover-suppressed regions. Recent work has found evidence that gene conversion may explain the low degrees of allelic differentiation in the dimorphic mating-type locus (MT) of the isogamous alga Chlamydomonas reinhardtii. However, no one has tested whether gene conversion is sufficient to avoid the degeneration of functional sequence within MT. Here, we calculate degree of linkage disequilibrium (LD) across MT as a proxy for recombination rate and investigate its relationship to patterns of population genetic variation and the efficacy of selection in the region. We find that degree of LD predicts selection efficacy across MT, and that purifying selection is stronger in shared genes than in MT-limited genes to the point of being equivalent to that of autosomal genes. We argue that while crossover suppression is needed in the mating-type loci of many isogamous systems, these loci are less likely to experience selection to differentiate further. Thus, recombination can act in these regions and prevent degeneration caused by Hill-Robertson effects.

Genome-wide association mapping of date palm fruit traits

Date palms (Phoenix dactylifera) are an important fruit crop of arid regions of the Middle East and North Africa. Despite its importance, few genomic resources exist for date palms, hampering evolutionary genomic studies of this perennial species. Here we report an improved long-read genome assembly for P. dactylifera that is 772.3 Mb in length, with contig N50 of 897.2 Kb, and use this to perform genome-wide association studies (GWAS) of the sex determining region and 21 fruit traits. We find a fruit color GWAS at the R2R3-MYB transcription factor VIRESCENS gene and identify functional alleles that include a retrotransposon insertion and start codon mutation. We also find a GWAS peak for sugar composition spanning deletion polymorphisms in multiple linked invertase genes. MYB transcription factors and invertase are implicated in fruit color and sugar composition in other crops, demonstrating the importance of parallel evolution in the evolutionary diversification of domesticated species.

Date palm is an important fruit crop in the Middle East and North Africa. Here, the authors report an improved genome assembly of this species and perform GWAS mapping of sex determining region and 21 fruit traits using high density SNP data generated from re-sequencing of the mapping population.

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

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

From molecular manipulation of domesticated Chlamydomonas reinhardtii to survival in nature

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

Environmental fluctuations accelerate molecular evolution of thermal tolerance in a marine diatom

Diatoms contribute roughly 20% of global primary production, but the factors determining their ability to adapt to global warming are unknown. Here we quantify the capacity for adaptation to warming in the marine diatom Thalassiosira pseudonana. We find that evolutionary rescue under severe (32 °C) warming is slow, but adaptation to more realistic scenarios where temperature increases are moderate (26 °C) or fluctuate between benign and severe conditions is rapid and linked to phenotypic changes in metabolic traits and elemental composition. Whole-genome re-sequencing identifies genetic divergence among populations selected in the different warming regimes and between the evolved and ancestral lineages. Consistent with the phenotypic changes, the most rapidly evolving genes are associated with transcriptional regulation, cellular responses to oxidative stress and redox homeostasis. These results demonstrate that the evolution of thermal tolerance in marine diatoms can be rapid, particularly in fluctuating environments, and is underpinned by major genomic and phenotypic change.

Diff-seq: A high throughput sequencing-based mismatch detection assay for DNA variant enrichment and discovery

Much of the within species genetic variation is in the form of single nucleotide polymorphisms (SNPs), typically detected by whole genome sequencing (WGS) or microarray-based technologies. However, WGS produces mostly uninformative reads that perfectly match the reference, while microarrays require genome-specific reagents. We have developed Diff-seq, a sequencing-based mismatch detection assay for SNP discovery without the requirement for specialized nucleic-acid reagents. Diff-seq leverages the Surveyor endonuclease to cleave mismatched DNA molecules that are generated after cross-annealing of a complex pool of DNA fragments. Sequencing libraries enriched for Surveyor-cleaved molecules result in increased coverage at the variant sites. Diff-seq detected all mismatches present in an initial test substrate, with specific enrichment dependent on the identity and context of the variation. Application to viral sequences resulted in increased observation of variant alleles in a biologically relevant context. Diff-Seq has the potential to increase the sensitivity and efficiency of high-throughput sequencing in the detection of variation.

Are Nonsense Alleles of Drosophila melanogaster Genes under Any Selection?

A gene which carries a bona fide loss-of-function mutation effectively becomes a functionless pseudogene, free from selective constraint. However, there is a number of molecular mechanisms that may lead to at least a partial preservation of the function of genes carrying even drastic alleles. We performed a direct measurement of the strength of negative selection acting on nonsense alleles of protein-coding genes in the Zambian population of Drosophila melanogaster. Within those exons that carry nonsense mutations, negative selection, assayed by the ratio of missense over synonymous nucleotide diversity levels, appears to be absent, consistent with total loss of function. In other exons of nonsense alleles, negative selection was deeply relaxed but likely not completely absent, and the per site number of missense alleles declined significantly with the distance from the premature stop codon. This pattern may be due to alternative splicing which preserves function of some isoforms of nonsense alleles of genes.

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.

Fitness change in relation to mutation number in spontaneous mutation accumulation lines of Chlamydomonas reinhardtii

Although all genetic variation ultimately stems from mutations, their properties are difficult to study directly. Here, we used multiple mutation accumulation (MA) lines derived from five genetic backgrounds of the green algae Chlamydomonas reinhardtii that have been previously subjected to whole genome sequencing to investigate the relationship between the number of spontaneous mutations and change in fitness from a nonevolved ancestor. MA lines were on average less fit than their ancestors and we detected a significantly negative correlation between the change in fitness and the total number of accumulated mutations in the genome. Likewise, the number of mutations located within coding regions significantly and negatively impacted MA line fitness. We used the fitness data to parameterize a maximum likelihood model to estimate discrete categories of mutational effects, and found that models containing one to two mutational effect categories (one neutral and one deleterious category) fitted the data best. However, the best‐fitting mutational effects models were highly dependent on the genetic background of the ancestral strain.

Revision of the Genus Micromonas Manton et Parke (Chlorophyta, Mamiellophyceae), of the Type Species M. pusilla (Butcher) Manton & Parke and of the Species M. commoda van Baren, Bachy and Worden and Description of Two New Species Based on the Genetic and Phenotypic Characterization of Cultured Isolates

The green picoalgal genus Micromonas is broadly distributed in estuaries, coastal marine habitats and open oceans, from the equator to the poles. Phylogenetic, ecological and genomic analyses of culture strains and natural populations have suggested that this cosmopolitan genus is composed of several cryptic species corresponding to genetic lineages. We performed a detailed analysis of variations in morphology, pigment content, and sequences of the nuclear-encoded small-subunit rRNA gene and the second internal transcribed spacer (ITS2) from strains isolated worldwide. A new morphological feature of the genus, the presence of tip hairs at the extremity of the hair point, was discovered and subtle differences in hair point length were detected between clades. Clear non-homoplasious synapomorphies were identified in the small-subunit rRNA gene and ITS2 spacer sequences of five genetic lineages. These findings lead us to provide emended descriptions of the genus Micromonas, of the type species M. pusilla, and of the recently described species M. commoda, as well as to describe 2 new species, M. bravo and M. polaris. By clarifying the status of the genetic lineages identified within Micromonas, these formal descriptions will facilitate further interpretations of large-scale analyses investigating ecological trends in time and space for this widespread picoplankter.

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

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

kWIP: The k-mer weighted inner product, a de novo estimator of genetic similarity

Modern genomics techniques generate overwhelming quantities of data. Extracting population genetic variation demands computationally efficient methods to determine genetic relatedness between individuals (or "samples") in an unbiased manner, preferably de novo. Rapid estimation of genetic relatedness directly from sequencing data has the potential to overcome reference genome bias, and to verify that individuals belong to the correct genetic lineage before conclusions are drawn using mislabelled, or misidentified samples. We present the k-mer Weighted Inner Product (kWIP), an assembly-, and alignment-free estimator of genetic similarity. kWIP combines a probabilistic data structure with a novel metric, the weighted inner product (WIP), to efficiently calculate pairwise similarity between sequencing runs from their k-mer counts. It produces a distance matrix, which can then be further analysed and visualised. Our method does not require prior knowledge of the underlying genomes and applications include establishing sample identity and detecting mix-up, non-obvious genomic variation, and population structure. We show that kWIP can reconstruct the true relatedness between samples from simulated populations. By re-analysing several published datasets we show that our results are consistent with marker-based analyses. kWIP is written in C++, licensed under the GNU GPL, and is available from https://github.com/kdmurray91/kwip.

Systems level analysis of the Chlamydomonas reinhardtii metabolic network reveals variability in evolutionary co-conservation

Metabolic networks, which are mathematical representations of organismal metabolism, are reconstructed to provide computational platforms to guide metabolic engineering experiments and explore fundamental questions on metabolism. Systems level analyses, such as interrogation of phylogenetic relationships within the network, can provide further guidance on the modification of metabolic circuitries. Chlamydomonas reinhardtii, a biofuel relevant green alga that has retained key genes with plant, animal, and protist affinities, serves as an ideal model organism to investigate the interplay between gene function and phylogenetic affinities at multiple organizational levels. Here, using detailed topological and functional analyses, coupled with transcriptomics studies on a metabolic network that we have reconstructed for C. reinhardtii, we show that network connectivity has a significant concordance with the co-conservation of genes; however, a distinction between topological and functional relationships is observable within the network. Dynamic and static modes of co-conservation were defined and observed in a subset of gene-pairs across the network topologically. In contrast, genes with predicted synthetic interactions, or genes involved in coupled reactions, show significant enrichment for both shorter and longer phylogenetic distances. Based on our results, we propose that the metabolic network of C. reinhardtii is assembled with an architecture to minimize phylogenetic profile distances topologically, while it includes an expansion of such distances for functionally interacting genes. This arrangement may increase the robustness of C. reinhardtii's network in dealing with varied environmental challenges that the species may face. The defined evolutionary constraints within the network, which identify important pairings of genes in metabolism, may offer guidance on synthetic biology approaches to optimize the production of desirable metabolites.

The effect of sex on the repeatability of evolution in different environments

The adaptive function of sex has been extensively studied, while less consideration has been given to the potential downstream consequences of sex on evolution. Here, we investigate one such potential consequence, the effect of sex on the repeatability of evolution. By affecting the repeatability of evolution, sex could have important implications for biodiversity, and for our ability to make predictions about the outcome of environmental change. We allowed asexual and sexual populations of Chlamydomonas reinhardtii to evolve in novel environments and monitored both their change in fitness and variance in fitness after evolution. Sex affected the repeatability of evolution by changing the importance of the effect of selection, chance, and ancestral constraints on the outcome of the evolutionary process. In particular, the effects of sex were highly dependent on the initial genetic composition of the population and on the environment. Given the lack of a consistent effect of sex on repeatability across the environments used here, further studies to dissect in more detail the underlying reasons for these differences as well as studies in additional environments are required if we are to have a general understanding of the effects of sex on the repeatability of evolution.