Whole-Genome Sequencing to Identify Mutants and Polymorphisms in Chlamydomonas reinhardtii

Loss of an extensive ciliary connectome induces proteostasis and cell fate switching in a severe motile ciliopathy

Motile cilia have essential cellular functions in development, reproduction, and homeostasis. Genetic causes for motile ciliopathies have been identified, but the consequences on cellular functions beyond impaired motility remain unknown. Variants in CCDC39 and CCDC40 cause severe disease not explained by loss of motility. Using human cells with pathological variants in these genes, Chlamydomonas genetics, cryo-electron microscopy, single cell RNA transcriptomics, and proteomics, we identified perturbations in multiple cilia-independent pathways. Absence of the axonemal CCDC39/CCDC40 heterodimer results in loss of a connectome of over 90 proteins. The undocked connectome activates cell quality control pathways, switches multiciliated cell fate, impairs microtubule architecture, and creates a defective periciliary barrier. Both cilia-dependent and independent defects are likely responsible for the disease severity. Our findings provide a foundation for reconsidering the broad cellular impact of pathologic variants in ciliopathies and suggest new directions for therapies.

The IFT81-IFT74 complex acts as an unconventional RabL2 GTPase-activating protein during intraflagellar transport

Cilia are important cellular organelles for signaling and motility and are constructed via intraflagellar transport (IFT). RabL2 is a small GTPase that localizes to the basal body of cilia via an interaction with the centriolar protein CEP19 before downstream association with the IFT machinery, which is followed by initiation of IFT. We reconstituted and purified RabL2 with CEP19 or IFT proteins to show that a reconstituted pentameric IFT complex containing IFT81/74 enhances the GTP hydrolysis rate of RabL2. The binding site on IFT81/74 that promotes GTP hydrolysis in RabL2 was mapped to a 70-amino-acid-long coiled-coil region of IFT81/74. We present structural models for RabL2-containing IFT complexes that we validate in vitro and in cellulo and demonstrate that Chlamydomonas IFT81/74 enhances GTP hydrolysis of human RabL2, suggesting an ancient evolutionarily conserved activity. Our results provide an architectural understanding of how RabL2 is incorporated into the IFT complex and a molecular rationale for why RabL2 dissociates from anterograde IFT trains soon after departure from the ciliary base.

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.

The short flagella 1 (SHF1) gene in Chlamydomonas encodes a Crescerin TOG-domain protein required for late stages of flagellar growth

Length control of flagella represents a simple and tractable system to investigate the dynamics of organelle size. Models for flagellar length control in the model organism, Chlamydomonas reinhardtii have focused on the length-dependence of the intraflagellar transport (IFT) system which manages the delivery and removal of axonemal subunits at the tip of the flagella. One of these cargoes, tubulin, is the major axonemal subunit, and its frequency of arrival at the tip plays a central role in size control models. However, the mechanisms determining tubulin dynamics at the tip are still poorly understood. We discovered a loss-of-function mutation that leads to shortened flagella, and found that this was an allele of a previously described gene, SHF1, whose molecular identity had not previously been determined.  We found that SHF1 encodes a Chlamydomonas ortholog of Crescerin, previously identified as a cilia-specific TOG-domain array protein that can bind tubulin via its TOG domains and increase tubulin polymerization rates. In this mutant, flagellar regeneration occurs with the same initial kinetics as wild-type cells, but plateaus at a shorter length. Using a computational model in which the flagellar microtubules are represented by a differential equation for flagellar length combined with a stochastic model for cytoplasmic microtubule dynamics, we found that our experimental results are best described by a model in which Crescerin/SHF1 binds tubulin dimers in the cytoplasm and transports them into the flagellum. We suggest that this TOG-domain protein is necessary to efficiently and preemptively increase intra-flagella tubulin levels to offset decreasing IFT cargo at the tip as flagellar assembly progresses.

Discovery of photosynthesis genes through whole-genome sequencing of acetate-requiring mutants of Chlamydomonas reinhardtii

Large-scale mutant libraries have been indispensable for genetic studies, and the development of next-generation genome sequencing technologies has greatly advanced efforts to analyze mutants. In this work, we sequenced the genomes of 660 Chlamydomonas reinhardtii acetate-requiring mutants, part of a larger photosynthesis mutant collection previously generated by insertional mutagenesis with a linearized plasmid. We identified 554 insertion events from 509 mutants by mapping the plasmid insertion sites through paired-end sequences, in which one end aligned to the plasmid and the other to a chromosomal location. Nearly all (96%) of the events were associated with deletions, duplications, or more complex rearrangements of genomic DNA at the sites of plasmid insertion, and together with deletions that were unassociated with a plasmid insertion, 1470 genes were identified to be affected. Functional annotations of these genes were enriched in those related to photosynthesis, signaling, and tetrapyrrole synthesis as would be expected from a library enriched for photosynthesis mutants. Systematic manual analysis of the disrupted genes for each mutant generated a list of 253 higher-confidence candidate photosynthesis genes, and we experimentally validated two genes that are essential for photoautotrophic growth, CrLPA3 and CrPSBP4. The inventory of candidate genes includes 53 genes from a phylogenomically defined set of conserved genes in green algae and plants. Altogether, 70 candidate genes encode proteins with previously characterized functions in photosynthesis in Chlamydomonas, land plants, and/or cyanobacteria; 14 genes encode proteins previously shown to have functions unrelated to photosynthesis. Among the remaining 169 uncharacterized genes, 38 genes encode proteins without any functional annotation, signifying that our results connect a function related to photosynthesis to these previously unknown proteins. This mutant library, with genome sequences that reveal the molecular extent of the chromosomal lesions and resulting higher-confidence candidate genes, will aid in advancing gene discovery and protein functional analysis in photosynthesis.

Co-expression networks in Chlamydomonas reveal significant rhythmicity in batch cultures and empower gene function discovery

The unicellular green alga Chlamydomonas reinhardtii is a choice reference system for the study of photosynthesis and chloroplast metabolism, cilium assembly and function, lipid and starch metabolism, and metal homeostasis. Despite decades of research, the functions of thousands of genes remain largely unknown, and new approaches are needed to categorically assign genes to cellular pathways. Growing collections of transcriptome and proteome data now allow a systematic approach based on integrative co-expression analysis. We used a dataset comprising 518 deep transcriptome samples derived from 58 independent experiments to identify potential co-expression relationships between genes. We visualized co-expression potential with the R package corrplot, to easily assess co-expression and anti-correlation between genes. We extracted several hundred high-confidence genes at the intersection of multiple curated lists involved in cilia, cell division, and photosynthesis, illustrating the power of our method. Surprisingly, Chlamydomonas experiments retained a significant rhythmic component across the transcriptome, suggesting an underappreciated variable during sample collection, even in samples collected in constant light. Our results therefore document substantial residual synchronization in batch cultures, contrary to assumptions of asynchrony. We provide step-by-step protocols for the analysis of co-expression across transcriptome data sets from Chlamydomonas and other species to help foster gene function discovery.

Identifying RNA splicing factors using IFT genes in Chlamydomonas reinhardtii

Intraflagellar transport moves proteins in and out of flagella/cilia and it is essential for the assembly of these organelles. Using whole-genome sequencing, we identified splice site mutations in two IFT genes, IFT81 (fla9) and IFT121 (ift121-2), which lead to flagellar assembly defects in the unicellular green alga Chlamydomonas reinhardtii. The splicing defects in these ift mutants are partially corrected by mutations in two conserved spliceosome proteins, DGR14 and FRA10. We identified a dgr14 deletion mutant, which suppresses the 3′ splice site mutation in IFT81, and a frameshift mutant of FRA10, which suppresses the 5′ splice site mutation in IFT121. Surprisingly, we found dgr14-1 and fra10 mutations suppress both splice site mutations. We suggest these two proteins are involved in facilitating splice site recognition/interaction; in their absence some splice site mutations are tolerated. Nonsense mutations in SMG1, which is involved in nonsense-mediated decay, lead to accumulation of aberrant transcripts and partial restoration of flagellar assembly in the ift mutants. The high density of introns and the conservation of noncore splicing factors, together with the ease of scoring the ift mutant phenotype, make Chlamydomonas an attractive organism to identify new proteins involved in splicing through suppressor screening.

MAPINS, a Highly Efficient Detection Method That Identifies Insertional Mutations and Complex DNA Rearrangements

Insertional mutagenesis, in which a piece of exogenous DNA is integrated randomly into the genomic DNA of the recipient cell, is a useful method to generate new mutants with phenotypes of interest. The unicellular green alga Chlamydomonas reinhardtii is an outstanding model for studying many biological processes. We developed a new computational algorithm, MAPINS (mapping insertions), to efficiently identify insertion sites created by the integration of an APHVIII (aminoglycoside 3'-phosphotransferase VIII) cassette that confers paromomycin resistance. Using whole-genome sequencing data, this method eliminates the need for genomic DNA manipulation and retains all the sequencing information provided by paired-end sequencing. We experimentally verified 38 insertion sites out of 41 sites (93%) identified by MAPINS from 20 paromomycin-resistant strains. Using meiotic analysis of 18 of these strains, we identified insertion sites that completely cosegregate with paromomycin resistance. In six of the seven strains with a mutant phenotype, we demonstrated complete cosegregation of the mutant phenotype and the verified insertion site. In addition, we provide direct evidence of complex rearrangements of genomic DNA in five strains, three of which involve the APHVIII insertion site. We suggest that strains obtained by insertional mutagenesis are more complicated than expected from previous analyses in Chlamydomonas To map the locations of some complex insertions, we designed 49 molecular markers based on differences identified via whole-genome sequencing between wild-type strains CC-124 and CC-125. Overall, MAPINS provides a low-cost, efficient method to characterize insertional mutants in Chlamydomonas.

Building a multipurpose insertional mutant library for forward and reverse genetics in Chlamydomonas

BACKGROUND

The unicellular green alga, Chlamydomonas reinhardtii, is a classic model for studying flagella and biofuel. However, precise gene editing, such as Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) and CRISPR-associated protein (Cas9) system, is not widely used in this organism. Screening of random insertional mutant libraries by polymerase chain reaction provides an alternate strategy to obtain null mutants of individual gene. But building, screening, and maintaining such a library was time-consuming and expensive.

RESULTS

By selecting a suitable parental strain, keeping individual mutants using the agar plate, and designing an insertion cassette-specific primer for library screening, we successfully generated and maintained ~150,000 insertional mutants of Chlamydomonas, which was used for both reverse and forward genetics analysis. We obtained 26 individual mutants corresponding to 20 genes and identified 967 motility-defect mutants including 10 mutants with defective accumulation of intraflagellar transport complex at the basal body. We also obtained 929 mutants defective in oil droplet assembly after nitrogen deprivation. Furthermore, a new insertion cassette with splicing donor sequences at both ends was also constructed, which increased the efficiency of gene interruption.

CONCLUSION

In summary, this library provides a multifunctional platform both for obtaining mutants of interested genes and for screening of mutants with specific phenotype.

Isolation, phenotypic characterization and genome wide analysis of a Chlamydomonas reinhardtii strain naturally modified under laboratory conditions: towards enhanced microalgal biomass and lipid production for biofuels

BACKGROUND

Microalgal strain development through genetic engineering has received much attention as a way to improve the traits of microalgae suitable for biofuel production. However, there are still some limitations in application of genetically modified organisms. In this regard, there has been recent interest in the isolation and characterization of superior strains naturally modified and/or adapted under a certain condition and on the interpretation of phenotypic changes through the whole genome sequencing.

RESULTS

In this study, we isolated and characterized a novel derivative of C. reinhardtii, whose phenotypic traits diverged significantly from its ancestral strain, C. reinhardtii CC-124. This strain, designated as CC-124H, displayed cell population containing increased numbers of larger cells, which resulted in an increased biomass productivity compared to its ancestor CC-124. CC-124H was further compared with the CC-124 wild-type strain which underwent long-term storage under low light condition, designated as CC-124L. In an effort to evaluate the potential of CC-124H for biofuel production, we also found that CC-124H accumulated 116 and 66% greater lipids than that of the CC-124L, after 4 days under nitrogen and sulfur depleted conditions, respectively. Taken together, our results revealed that CC-124H had significantly increased fatty acid methyl ester (FAME) yields that were 2.66 and 1.98 times higher than that of the CC-124L at 4 days after the onset of cultivation under N and S depleted conditions, respectively, and these higher FAME yields were still maintained by day 8. We next analyzed single nucleotide polymorphisms (SNPs) and insertion/deletions (indels) based on the whole genome sequencing. The result revealed that of the 44 CDS region alterations, 34 resulted in non-synonymous substitutions within 33 genes which may mostly be involved in cell cycle, division or proliferation.

CONCLUSION

Our phenotypic analysis, which emphasized lipid productivity, clearly revealed that CC-124H had a dramatically enhanced biomass and lipid content compared to the CC-124L. Moreover, SNPs and indels analysis enabled us to identify 34 of non-synonymous substitutions which may result in phenotypic changes of CC-124H. All of these results suggest that the concept of adaptive evolution combined with genome wide analysis can be applied to microalgal strain development for biofuel production.

Volvox: A simple algal model for embryogenesis, morphogenesis and cellular differentiation

Patterning of a multicellular body plan involves a coordinated set of developmental processes that includes cell division, morphogenesis, and cellular differentiation. These processes have been most intensively studied in animals and land plants; however, deep insight can also be gained by studying development in simpler multicellular organisms. The multicellular green alga Volvox carteri (Volvox) is an excellent model for the investigation of developmental mechanisms and their evolutionary origins. Volvox has a streamlined body plan that contains only a few thousand cells and two distinct cell types: reproductive germ cells and terminally differentiated somatic cells. Patterning of the Volvox body plan is achieved through a stereotyped developmental program that includes embryonic cleavage with asymmetric cell division, morphogenesis, and cell-type differentiation. In this review we provide an overview of how these three developmental processes give rise to the adult form in Volvox and how developmental mutants have provided insights into the mechanisms behind these events. We highlight the accessibility and tractability of Volvox and its relatives that provide a unique opportunity for studying development.

Mapping Challenging Mutations by Whole-Genome Sequencing

Whole-genome sequencing provides a rapid and powerful method for identifying mutations on a global scale, and has spurred a renewed enthusiasm for classical genetic screens in model organisms. The most commonly characterized category of mutation consists of monogenic, recessive traits, due to their genetic tractability. Therefore, most of the mapping methods for mutation identification by whole-genome sequencing are directed toward alleles that fulfill those criteria (i.e., single-gene, homozygous variants). However, such approaches are not entirely suitable for the characterization of a variety of more challenging mutations, such as dominant and semidominant alleles or multigenic traits. Therefore, we have developed strategies for the identification of those classes of mutations, using polymorphism mapping in Caenorhabditis elegans as our model for validation. We also report an alternative approach for mutation identification from traditional recombinant crosses, and a solution to the technical challenge of sequencing sterile or terminally arrested strains where population size is limiting. The methods described herein extend the applicability of whole-genome sequencing to a broader spectrum of mutations, including classes that are difficult to map by traditional means.

Evidence That an Unconventional Actin Can Provide Essential F-Actin Function and That a Surveillance System Monitors F-Actin Integrity in Chlamydomonas

Actin is one of the most conserved eukaryotic proteins. It is thought to have multiple essential cellular roles and to function primarily or exclusively as filaments ("F-actin"). Chlamydomonas has been an enigma, because a null mutation (ida5-1) in its single gene for conventional actin does not affect growth. A highly divergent actin gene, NAP1, is upregulated in ida5-1 cells, but it has been unclear whether NAP1 can form filaments or provide actin function. Here, we used the actin-depolymerizing drug latrunculin B (LatB), the F-actin-specific probe Lifeact-Venus, and genetic and molecular methods to resolve these issues. LatB-treated wild-type cells continue to proliferate; they initially lose Lifeact-stained structures but recover them concomitant with upregulation of NAP1. Thirty-nine LatB-sensitive mutants fell into four genes (NAP1 and LAT1-LAT3) in which we identified the causative mutations using a novel combinatorial pool-sequencing strategy. LAT1-LAT3 are required for NAP1 upregulation upon LatB treatment, and ectopic expression of NAP1 largely rescues the LatB sensitivity of the lat1-lat3 mutants, suggesting that the LAT gene products comprise a regulatory hierarchy with NAP1 expression as the major functional output. Selection of LatB-resistant revertants of a nap1 mutant yielded dominant IDA5 mutations that presumably render F-IDA5 resistant to LatB, and nap1 and lat mutations are synthetically lethal with ida5-1 in the absence of LatB. We conclude that both IDA5 and the divergent NAP1 can form filaments and redundantly provide essential F-actin functions and that a novel surveillance system, probably responding to a loss of F-actin, triggers NAP1 expression and perhaps other compensatory responses.

A NIMA-Related Kinase Suppresses the Flagellar Instability Associated with the Loss of Multiple Axonemal Structures

CCDC39 and CCDC40 were first identified as causative mutations in primary ciliary dyskinesia patients; cilia from patients show disorganized microtubules, and they are missing both N-DRC and inner dynein arms proteins. In Chlamydomonas, we used immunoblots and microtubule sliding assays to show that mutants in CCDC40 (PF7) and CCDC39 (PF8) fail to assemble N-DRC, several inner dynein arms, tektin, and CCDC39. Enrichment screens for suppression of pf7; pf8 cells led to the isolation of five independent extragenic suppressors defined by four different mutations in a NIMA-related kinase, CNK11. These alleles partially rescue the flagellar length defect, but not the motility defect. The suppressor does not restore the missing N-DRC and inner dynein arm proteins. In addition, the cnk11 mutations partially suppress the short flagella phenotype of N-DRC and axonemal dynein mutants, but do not suppress the motility defects. The tpg1 mutation in TTLL9, a tubulin polyglutamylase, partially suppresses the length phenotype in the same axonemal dynein mutants. In contrast to cnk11, tpg1 does not suppress the short flagella phenotype of pf7. The polyglutamylated tubulin in the proximal region that remains in the tpg1 mutant is reduced further in the pf7; tpg1 double mutant by immunofluorescence. CCDC40, which is needed for docking multiple other axonemal complexes, is needed for tubulin polyglutamylation in the proximal end of the flagella. The CCDC39 and CCDC40 proteins are likely to be involved in recruiting another tubulin glutamylase(s) to the flagella. Another difference between cnk11-1 and tpg1 mutants is that cnk11-1 cells show a faster turnover rate of tubulin at the flagellar tip than in wild-type flagella and tpg1 flagella show a slower rate. The double mutant shows a turnover rate similar to tpg1, which suggests the faster turnover rate in cnk11-1 flagella requires polyglutamylation. Thus, we hypothesize that many short flagella mutants in Chlamydomonas have increased instability of axonemal microtubules. Both CNK11 and tubulin polyglutamylation play roles in regulating the stability of axonemal microtubules.

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

We performed whole-genome resequencing of 12 field isolates and eight commonly studied laboratory strains of the model organism Chlamydomonas reinhardtii to characterize genomic diversity and provide a resource for studies of natural variation. Our data support previous observations that Chlamydomonas is among the most diverse eukaryotic species. Nucleotide diversity is ∼3% and is geographically structured in North America with some evidence of admixture among sampling locales. Examination of predicted loss-of-function mutations in field isolates indicates conservation of genes associated with core cellular functions, while genes in large gene families and poorly characterized genes show a greater incidence of major effect mutations. De novo assembly of unmapped reads recovered genes in the field isolates that are absent from the CC-503 assembly. The laboratory reference strains show a genomic pattern of polymorphism consistent with their origin as the recombinant progeny of a diploid zygospore. Large duplications or amplifications are a prominent feature of laboratory strains and appear to have originated under laboratory culture. Extensive natural variation offers a new source of genetic diversity for studies of Chlamydomonas, including naturally occurring alleles that may prove useful in studies of gene function and the dissection of quantitative genetic traits.

The Octatricopeptide Repeat Protein Raa8 Is Required for Chloroplast trans Splicing

The mRNA maturation of the tripartite chloroplast psaA gene from the green alga Chlamydomonas reinhardtii depends on various nucleus-encoded factors that participate in trans splicing of two group II introns. Recently, a multiprotein complex was identified that is involved in processing the psaA precursor mRNA. Using coupled tandem affinity purification (TAP) and mass spectrometry analyses with the trans-splicing factor Raa4 as a bait protein, we recently identified a multisubunit ribonucleoprotein (RNP) complex comprising the previously characterized trans-splicing factors Raa1, Raa3, Raa4, and Rat2 plus novel components. Raa1 and Rat2 share a structural motif, an octatricopeptide repeat (OPR), that presumably functions as an RNA interaction module. Two of the novel RNP complex components also exhibit a predicted OPR motif and were therefore considered potential trans-splicing factors. In this study, we selected bacterial artificial chromosome (BAC) clones encoding these OPR proteins and conducted functional complementation assays using previously generated trans-splicing mutants. Our assay revealed that the trans-splicing defect of mutant F19 was restored by a new factor we named RAA8; molecular characterization of complemented strains verified that Raa8 participates in splicing of the first psaA group II intron. Three of six OPR motifs are located in the C-terminal end of Raa8, which was shown to be essential for restoring psaA mRNA trans splicing. Our results support the important role played by OPR proteins in chloroplast RNA metabolism and also demonstrate that combining TAP and mass spectrometry with functional complementation studies represents a vigorous tool for identifying trans-splicing factors.

Genetic tools and techniques for Chlamydomonas reinhardtii

The development of tools has always been a major driving force for the advancement of science. Optical microscopes were the first instruments that allowed discovery and descriptive studies of the subcellular features of microorganisms. Although optical and electron microscopes remained at the forefront of microbiological research tools since their inventions, the advent of molecular genetics brought about questions which had to be addressed with new "genetic tools". The unicellular green microalgal genus Chlamydomonas, especially the most prominent species C. reinhardtii, has become a frequently used model organism for many diverse fields of research and molecular genetic analyses of C. reinhardtii, as well as the available genetic tools and techniques, have become increasingly sophisticated throughout the last decades. The aim of this review is to provide an overview of the molecular key features of C. reinhardtii and summarize the progress related to the development of tools and techniques for genetic engineering of this organism, from pioneering DNA transformation experiments to state-of-the-art techniques for targeted nuclear genome editing and high-throughput screening approaches.

Molecular techniques to interrogate and edit the Chlamydomonas nuclear genome

The success of the green alga Chlamydomonas reinhardtii as a model organism is to a large extent due to the wide range of molecular techniques that are available for its characterization. Here, we review some of the techniques currently used to modify and interrogate the C. reinhardtii nuclear genome and explore several technologies under development. Nuclear mutants can be generated with ultraviolet (UV) light and chemical mutagens, or by insertional mutagenesis. Nuclear transformation methods include biolistic delivery, agitation with glass beads, and electroporation. Transforming DNA integrates into the genome at random sites, and multiple strategies exist for mapping insertion sites. A limited number of studies have demonstrated targeted modification of the nuclear genome by approaches such as zinc-finger nucleases and homologous recombination. RNA interference is widely used to knock down expression levels of nuclear genes. A wide assortment of transgenes has been successfully expressed in the Chlamydomonas nuclear genome, including transformation markers, fluorescent proteins, reporter genes, epitope tagged proteins, and even therapeutic proteins. Optimized expression constructs and strains help transgene expression. Emerging technologies such as the CRISPR/Cas9 system, high-throughput mutant identification, and a whole-genome knockout library are being developed for this organism. We discuss how these advances will propel future investigations.

Metabolism of acyl-lipids in Chlamydomonas reinhardtii

Microalgae are emerging platforms for production of a suite of compounds targeting several markets, including food, nutraceuticals, green chemicals, and biofuels. Many of these products, such as biodiesel or polyunsaturated fatty acids (PUFAs), derive from lipid metabolism. A general picture of lipid metabolism in microalgae has been deduced from well characterized pathways of fungi and land plants, but recent advances in molecular and genetic analyses of microalgae have uncovered unique features, pointing out the necessity to study lipid metabolism in microalgae themselves. In the past 10 years, in addition to its traditional role as a model for photosynthetic and flagellar motility processes, Chlamydomonas reinhardtii has emerged as a model organism to study lipid metabolism in green microalgae. Here, after summarizing data on total fatty acid composition, distribution of acyl-lipid classes, and major acyl-lipid molecular species found in C. reinhardtii, we review the current knowledge on the known or putative steps for fatty acid synthesis, glycerolipid desaturation and assembly, membrane lipid turnover, and oil remobilization. A list of characterized or putative enzymes for the major steps of acyl-lipid metabolism in C. reinhardtii is included, and subcellular localizations and phenotypes of associated mutants are discussed. Biogenesis and composition of Chlamydomonas lipid droplets and the potential importance of lipolytic processes in increasing cellular oil content are also highlighted.

MLT1 links cytoskeletal asymmetry to organelle placement in chlamydomonas

Asymmetric placement of the photosensory eyespot organelle in Chlamydomonas is patterned by mother-daughter differences between the two basal bodies, which template the anterior flagella. Each basal body is associated with two bundled microtubule rootlets, one with two microtubules and one with four, forming a cruciate pattern. In wild-type cells, the single eyespot is positioned at the equator in close proximity to the plus end of the daughter rootlet comprising four microtubules, the D4. Here we identify mutations in two linked loci, MLT1 and MLT2, which cause multiple eyespots. Antiserum raised against MLT1 localized the protein along the D4 rootlet microtubules, from the basal bodies to the eyespot. MLT1 associates immediately with the new D4 as it extends during cell division, before microtubule acetylation. MLT1 is a low-complexity protein of over 300,000 Daltons. The expression or stability of MLT1 is dependent on MLT2, predicted to encode a second large, low-complexity protein. MLT1 was not restricted to the D4 rootlet in cells with the vfl2-220 mutation in the gene encoding the basal body-associated protein centrin. The cumulative data highlight the role of mother-daughter basal body differences in establishing asymmetry in associated rootlets, and suggest that eyespot components are directed to the correct location by MLT1 on the D4 microtubules.

Large-scale insertional mutagenesis of Chlamydomonas supports phylogenomic functional prediction of photosynthetic genes and analysis of classical acetate-requiring mutants

Chlamydomonas reinhardtii is a unicellular green alga that is a key model organism in the study of photosynthesis and oxidative stress. Here we describe the large-scale generation of a population of insertional mutants that have been screened for phenotypes related to photosynthesis and the isolation of 459 flanking sequence tags from 439 mutants. Recent phylogenomic analysis has identified a core set of genes, named GreenCut2, that are conserved in green algae and plants. Many of these genes are likely to be central to the process of photosynthesis, and they are over-represented by sixfold among the screened insertional mutants, with insertion events isolated in or adjacent to 68 of 597 GreenCut2 genes. This enrichment thus provides experimental support for functional assignments based on previous bioinformatic analysis. To illustrate one of the uses of the population, a candidate gene approach based on genome position of the flanking sequence of the insertional mutant CAL027_01_20 was used to identify the molecular basis of the classical C. reinhardtii mutation ac17. These mutations were shown to affect the gene PDH2, which encodes a subunit of the plastid pyruvate dehydrogenase complex. The mutants and associated flanking sequence data described here are publicly available to the research community, and they represent one of the largest phenotyped collections of algal insertional mutants to date.

Genetic and genomic approaches to identify genes involved in flagellar assembly in Chlamydomonas reinhardtii

Flagellar assembly requires intraflagellar transport of components from the cell body to the flagellar tip for assembly. The understanding of flagellar assembly has been aided by the ease of biochemistry and the availability of mutants in the unicellular green alga, Chlamydomonas reinhardtii. In this chapter, we discuss means to identify genes involved in these processes using forward and reverse genetics. In particular, the ease and low cost of whole genome sequencing (WGS) will help to make gene identification easier and promote the understanding of this important process.

Fast forward genetics to identify mutations causing a high light tolerant phenotype in Chlamydomonas reinhardtii by whole-genome-sequencing

BACKGROUND

High light tolerance of microalgae is a desired phenotype for efficient cultivation in large scale production systems under fluctuating outdoor conditions. Outdoor cultivation requires the use of either wild-type or non-GMO derived mutant strains due to safety concerns. The identification and molecular characterization of such mutants derived from untagged forward genetics approaches was limited previously by the tedious and time-consuming methods involving techniques such as classical meiotic mapping. The combination of mapping with next generation sequencing technologies offers alternative strategies to identify genes involved in high light adaptation in untagged mutants.

RESULTS

We used the model alga Chlamydomonas reinhardtii in a non-GMO mutation strategy without any preceding crossing step or pooled progeny to identify genes involved in the regulatory processes of high light adaptation. To generate high light tolerant mutants, wildtype cells were mutagenized only to a low extent, followed by a stringent selection. We performed whole-genome sequencing of two independent mutants hit1 and hit2 and the parental wildtype. The availability of a reference genome sequence and the removal of shared bakground variants between the wildtype strain and each mutant, enabled us to identify two single nucleotide polymorphisms within the same gene Cre02.g085050, hereafter called LRS1 (putative Light Response Signaling protein 1). These two independent single amino acid exchanges are both located in the putative WD40 propeller domain of the corresponding protein LRS1. Both mutants exhibited an increased rate of non-photochemical-quenching (NPQ) and an improved resistance against chemically induced reactive oxygen species. In silico analyses revealed homology of LRS1 to the photoregulatory protein COP1 in plants.

CONCLUSIONS

In this work we identified the nuclear encoded gene LRS1 as an essential factor for high light adaptation in C. reinhardtii. The causative random mutation within this gene was identified by a rapid and efficient method, avoiding any preceding crossing step, meiotic mapping, or pooled progeny. Our results open up new insights into mechanisms of high light adaptation in microalgae and at the same time provide a simplified strategy for non-GMO forward genetics, a crucial precondition that could result in the identification of key factors for economically relevant biological processes within algae.

The Chlamydomonas mutant pf27 reveals novel features of ciliary radial spoke assembly

To address the mechanisms of ciliary radial spoke assembly, we took advantage of the Chlamydomonas pf27 mutant. The radial spokes that assemble in pf27 are localized to the proximal quarter of the axoneme, but otherwise are fully assembled into 20S radial spoke complexes competent to bind spokeless axonemes in vitro. Thus, pf27 is not defective in radial spoke assembly or docking to the axoneme. Rather, our results suggest that pf27 is defective in the transport of spoke complexes. During ciliary regeneration in pf27, radial spoke assembly occurs asynchronously from other axonemal components. In contrast, during ciliary regeneration in wild-type Chlamydomonas, radial spokes and other axonemal components assemble concurrently as the axoneme grows. Complementation in temporary dikaryons between wild-type and pf27 reveals rescue of radial spoke assembly that begins at the distal tip, allowing further assembly to proceed from tip to base of the axoneme. Notably, rescued assembly of radial spokes occurred independently of the established proximal radial spokes in pf27 axonemes in dikaryons. These results reveal that 20S radial spokes can assemble proximally in the pf27 cilium but as the cilium lengthens, spoke assembly requires transport. We postulate that PF27 encodes an adaptor or modifier protein required for radial spoke–IFT interaction.

A high-throughput fatty acid profiling screen reveals novel variations in fatty acid biosynthesis in Chlamydomonas reinhardtii and related algae

Analysis of fatty acid methyl esters (FAMEs) by gas chromatography (GC) is a common technique for the quantitative and qualitative analysis of acyl lipids. Methods for FAME preparation are typically time-consuming and labor-intensive and require multiple transfers of reagents and products between reaction tubes and autosampler vials. In order to increase throughput and lower the time and materials costs required for FAME preparation prior to GC analysis, we have developed a method in which 10-to-20-mg samples of microbial biomass are transferred to standard GC autosampler vials, transesterified using an emulsion of methanolic trimethylsulfonium hydroxide and hexane, and analyzed directly by GC without further sample handling. This method gives results that are essentially identical to those obtained by the more labor- and material-intensive FAME preparation methods, such as transmethylation with methanolic HCl. We applied this method to the screening of laboratory and environmental isolates of the green alga Chlamydomonas for variations in fatty acid composition. This screening method facilitated two novel discoveries. First, we identified a common laboratory strain of C. reinhardtii, CC-620, completely lacking all ω-3 fatty acids normally found in this organism and showed that this strain contains an inactivating mutation in the CrFAD7 gene, encoding the sole ω-3 desaturase activity in this organism. Second, we showed that some species of Chlamydomonas make Δ6-unsaturated polyunsaturated fatty acids (PUFA) rather than the Δ5 species normally made by the previously characterized laboratory strains of Chlamydomonas, suggesting that there is species-specific variation in the regiospecificity and substrate selectivity of front-end desaturases in this algal genus.

Using next-generation sequencing to isolate mutant genes from forward genetic screens

The long-lasting success of forward genetic screens relies on the simple molecular basis of the characterized phenotypes, which are typically caused by mutations in single genes. Mapping the location of causal mutations using genetic crosses has traditionally been a complex, multistep procedure, but next-generation sequencing now allows the rapid identification of causal mutations at single-nucleotide resolution even in complex genetic backgrounds. Recent advances of this mapping-by-sequencing approach include methods that are independent of reference genome sequences, genetic crosses and any kind of linkage information, which make forward genetics amenable for species that have not been considered for forward genetic screens so far.

Development of a forward genetic screen to isolate oil mutants in the green microalga Chlamydomonas reinhardtii

BACKGROUND

Oils produced by microalgae are precursors to biodiesel. To achieve a profitable production of biodiesel from microalgae, identification of factors governing oil synthesis and turnover is desirable. The green microalga Chlamydomonas reinhardtii is amenable to genetic analyses and has recently emerged as a model to study oil metabolism. However, a detailed method to isolate various types of oil mutants that is adapted to Chlamydomonas has not been reported.

RESULTS

We describe here a forward genetic approach to isolate mutants altered in oil synthesis and turnover from C. reinhardtii. It consists of a three-step screening procedure: a primary screen by flow cytometry of Nile red stained transformants grown in 96-deep-well plates under three sequential conditions (presence of nitrogen, then absence of nitrogen, followed by oil remobilization); a confirmation step using Nile red stained biological triplicates; and a validation step consisting of the quantification by thin layer chromatography of oil content of selected strains. Thirty-one mutants were isolated by screening 1,800 transformants generated by random insertional mutagenesis (1.7%). Five showed increased oil accumulation under the nitrogen-replete condition and 13 had altered oil content under nitrogen-depletion. All mutants were affected in oil remobilization.

CONCLUSION

This study demonstrates that various types of oil mutants can be isolated in Chlamydomonas based on the method set-up here, including mutants accumulating oil under optimal biomass growth. The strategy conceived and the protocol set-up should be applicable to other microalgal species such as Nannochloropsis and Chlorella, thus serving as a useful tool in Chlamydomonas oil research and algal biotechnology.

New mutations in flagellar motors identified by whole genome sequencing in Chlamydomonas

Background

The building of a cilium or flagellum requires molecular motors and associated proteins that allow the relocation of proteins from the cell body to the distal end and the return of proteins to the cell body in a process termed intraflagellar transport (IFT). IFT trains are carried out by kinesin and back to the cell body by dynein.

Methods

We used whole genome sequencing to identify the causative mutations for two temperature-sensitive flagellar assembly mutants in Chlamydomonas and validated the changes using reversion analysis. We examined the effect of these mutations on the localization of IFT81, an IFT complex B protein, the cytoplasmic dynein heavy chain (DHC1b), and the dynein light intermediate chain (D1bLIC).

Results

The strains, fla18 and fla24, have mutations in kinesin-2 and cytoplasmic dynein, respectively. The fla18 mutation alters the same glutamic acid (E₂₄G) mutated in the fla10-14 allele (E₂₄K). The fla18 strain loses flagella at 32?C more rapidly than the E₂₄K allele but less rapidly than the fla10-1 allele. The fla18 mutant loses its flagella by detachment rather than by shortening. The fla24 mutation falls in cytoplasmic dynein and changes a completely conserved amino acid (L₃₂₄₃P) in an alpha helix in the AAA5 domain. The fla24 mutant loses its flagella by shortening within 6 hours at 32?C. DHC1b protein is reduced by 18-fold and D1bLIC is reduced by 16-fold at 21?C compared to wild-type cells. We identified two pseudorevertants (L₃₂₄₃S and L₃₂₄₃R), which remain flagellated at 32?C. Although fla24 cells assemble full-length flagella at 21?C, IFT81 protein localization is dramatically altered. Instead of localizing at the basal body and along the flagella, IFT81 is concentrated at the proximal end of the flagella. The pseudorevertants show wild-type IFT81 localization at 21?C, but proximal end localization of IFT81 at 32?C.

Conclusions

The change in the AAA5 domain of the cytoplasmic dynein in fla24 may block the recycling of IFT trains after retrograde transport. It is clear that different alleles in the flagellar motors reveal different functions and roles. Multiple alleles will be important for understanding structure-function relationships.

Whole genome sequencing identifies a deletion in protein phosphatase 2A that affects its stability and localization in Chlamydomonas reinhardtii

Whole genome sequencing is a powerful tool in the discovery of single nucleotide polymorphisms (SNPs) and small insertions/deletions (indels) among mutant strains, which simplifies forward genetics approaches. However, identification of the causative mutation among a large number of non-causative SNPs in a mutant strain remains a big challenge. In the unicellular biflagellate green alga Chlamydomonas reinhardtii, we generated a SNP/indel library that contains over 2 million polymorphisms from four wild-type strains, one highly polymorphic strain that is frequently used in meiotic mapping, ten mutant strains that have flagellar assembly or motility defects, and one mutant strain, imp3, which has a mating defect. A comparison of polymorphisms in the imp3 strain and the other 15 strains allowed us to identify a deletion of the last three amino acids, Y313F314L315, in a protein phosphatase 2A catalytic subunit (PP2A3) in the imp3 strain. Introduction of a wild-type HA-tagged PP2A3 rescues the mutant phenotype, but mutant HA-PP2A3 at Y313 or L315 fail to rescue. Our immunoprecipitation results indicate that the Y313, L315, or YFLΔ mutations do not affect the binding of PP2A3 to the scaffold subunit, PP2A-2r. In contrast, the Y313, L315, or YFLΔ mutations affect both the stability and the localization of PP2A3. The PP2A3 protein is less abundant in these mutants and fails to accumulate in the basal body area as observed in transformants with either wild-type HA-PP2A3 or a HA-PP2A3 with a V310T change. The accumulation of HA-PP2A3 in the basal body region disappears in mated dikaryons, which suggests that the localization of PP2A3 may be essential to the mating process. Overall, our results demonstrate that the terminal YFL tail of PP2A3 is important in the regulation on Chlamydomonas mating.

IFT80 is essential for chondrocyte differentiation by regulating Hedgehog and Wnt signaling pathways

Partial mutation of intraflagellar transport 80 (IFT80) in humans causes Jeune asphyxiating thoracic dystrophy (JATD) and short-rib polydactyly (SRP) syndrome type III. These diseases are autosomal recessive chondrodysplasias that share clinical similarities, including shortened long bones and constricted thoracic cage. However, the role and mechanism of IFT80 in the regulation of chondrocyte differentiation and function remain largely unknown. We hypothesize that IFT80 is required for the formation and function of cilia and plays a critical role in chondrogenic differentiation by regulating Hedgehog (Hh) and Wingless (Wnt) signaling pathways. To test this hypothesis, we first analyzed the IFT80 expression pattern and found that IFT80 was predominantly expressed in growth plate chondrocytes and during chondrogenic differentiation. Silencing IFT80 impaired cilia formation and chondrogenic differentiation in mouse bone marrow derived stromal cells (BMSCs), and decreased the expression of chondrocyte marker genes--collagen II and aggrecan. Additionally, silencing IFT80 down-regulated Hh signaling activity whereas up-regulated Wnt signaling activity. The overexpression of Gli2 in IFT80-silenced cells promoted chondrogenesis and recovered the chondrogenic deficiency from IFT80 silencing. Overall, our results demonstrate that IFT80 is essential for chondrocyte differentiation by regulating the Hh and Wnt signaling pathways.

Katanin localization requires triplet microtubules in Chlamydomonas reinhardtii

Centrioles and basal bodies are essential for a variety of cellular processes that include the recruitment of proteins to these structures for both centrosomal and ciliary function. This recruitment is compromised when centriole/basal body assembly is defective. Mutations that cause basal body assembly defects confer supersensitivity to Taxol. These include bld2, bld10, bld12, uni3, vfl1, vfl2, and vfl3. Flagellar motility mutants do not confer sensitivity with the exception of mutations in the p60 (pf19) and p80 (pf15) subunits of the microtubule severing protein katanin. We have identified additional pf15 and bld2 (ε-tubulin) alleles in screens for Taxol sensitivity. Null pf15 and bld2 alleles are viable and are not essential genes in Chlamydomonas. Analysis of double mutant strains with the pf15-3 and bld2-6 null alleles suggests that basal bodies in Chlamydomonas may recruit additional proteins beyond katanin that affect spindle microtubule stability. The bld2-5 allele is a hypomorphic allele and its phenotype is modulated by nutritional cues. Basal bodies in bld2-5 cells are missing proximal ends. The basal body mutants show aberrant localization of an epitope-tagged p80 subunit of katanin. Unlike IFT proteins, katanin p80 does not localize to the transition fibers of the basal bodies based on an analysis of the uni1 mutant as well as the lack of colocalization of katanin p80 with IFT74. We suggest that the triplet microtubules are likely to play a key role in katanin p80 recruitment to the basal body of Chlamydomonas rather than the transition fibers that are needed for IFT localization.

Dissecting the sequential assembly and localization of intraflagellar transport particle complex B in Chlamydomonas

Intraflagellar transport (IFT), the key mechanism for ciliogenesis, involves large protein particles moving bi-directionally along the entire ciliary length. IFT particles contain two large protein complexes, A and B, which are constructed with proteins in a core and several peripheral proteins. Prior studies have shown that in Chlamydomonas reinhardtii, IFT46, IFT52, and IFT88 directly interact with each other and are in a subcomplex of the IFT B core. However, ift46, bld1, and ift88 mutants differ in phenotype as ift46 mutants are able to form short flagella, while the other two lack flagella completely. In this study, we investigated the functional differences of these individual IFT proteins contributing to complex B assembly, stability, and basal body localization. We found that complex B is completely disrupted in bld1 mutant, indicating an essential role of IFT52 for complex B core assembly. Ift46 mutant cells are capable of assembling a relatively intact complex B, but such complex is highly unstable and prone to degradation. In contrast, in ift88 mutant cells the complex B core still assembles and remains stable, but the peripheral proteins no longer attach to the B core. Moreover, in ift88 mutant cells, while complex A and the anterograde IFT motor FLA10 are localized normally to the transition fibers, complex B proteins instead are accumulated at the proximal ends of the basal bodies. In addition, in bld2 mutant, the IFT complex B proteins still localize to the proximal ends of defective centrioles which completely lack transition fibers. Taken together, these results revealed a step-wise assembly process for complex B, and showed that the complex first localizes to the proximal end of the centrioles and then translocates onto the transition fibers via an IFT88-dependent mechanism.

Genome-wide characterization of genetic variation in the unicellular, green alga Chlamydomonas reinhardtii

Chlamydomonas reinhardtii is a model system for studying cilia, photosynthesis, and other core features of eukaryotes, and is also an emerging source of biofuels. Despite its importance to basic and applied biological research, the level and pattern of genetic variation in this haploid green alga has yet to be characterized on a genome-wide scale. To improve understanding of C. reinhardtii's genetic variability, we generated low coverage whole genome resequencing data for nearly all of the available isolates of this species, which were sampled from a number of sites in North America over the past ∼70 years. Based on the analysis of more than 62,000 single nucleotide polymorphisms, we identified two groups of isolates that represent geographical subpopulations of the species. We also found that measurements of genetic diversity were highly variable throughout the genome, in part due to technical factors. We studied the level and pattern of linkage disequilibrium (LD), and observed one chromosome that exhibits elevated LD. Furthermore, we detected widespread evidence of recombination across the genome, which implies that outcrossing occurs in natural populations of this species. In summary, our study provides multiple insights into the sequence diversity of C. reinhardtii that will be useful to future studies of natural genetic variation in this organism.

Genomics of Volvocine Algae

Volvocine algae are a group of chlorophytes that together comprise a unique model for evolutionary and developmental biology. The species Chlamydomonas reinhardtii and Volvox carteri represent extremes in morphological diversity within the Volvocine clade. Chlamydomonas is unicellular and reflects the ancestral state of the group, while Volvox is multicellular and has evolved numerous innovations including germ-soma differentiation, sexual dimorphism, and complex morphogenetic patterning. The Chlamydomonas genome sequence has shed light on several areas of eukaryotic cell biology, metabolism and evolution, while the Volvox genome sequence has enabled a comparison with Chlamydomonas that reveals some of the underlying changes that enabled its transition to multicellularity, but also underscores the subtlety of this transition. Many of the tools and resources are in place to further develop Volvocine algae as a model for evolutionary genomics.