# Whole genome re-sequencing identifies a quantitative trait locus repressing carbon reserve accumulation during optimal growth in *Chlamydomonas reinhardtii*

Hugh Douglas Goold, Hoa Mai Nguyen, Fantao Kong, Audrey Beyly-Adriano, Bertrand Legeret, Emmanuelle Billon, Stéphan Cuiné, Fred Beisson, G.

Peltier, Yonghua Beisson

<figure><img src="image_1.png" /></figure>

## To cite this version:

Hugh Douglas Goold, Hoa Mai Nguyen, Fantao Kong, Audrey Beyly-Adriano, Bertrand Legeret, et al.. Whole genome re-sequencing identifies a quantitative trait locus repressing carbon reserve accumulation during optimal growth in *Chlamydomonas reinhardtii*. Scientific Reports, 2016, 6, pp.25209. 10.1038/srep25209. hal-02065385

HAL Id: hal-02065385

https://hal.science/hal-02065385

Submitted on 10 Jul 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

<figure><img src="image_3.png" /></figure>

OPEN

# Whole Genome Re-Sequencing Identifies a Quantitative Trait Locus Repressing Carbon Reserve Accumulation during Optimal Growth in *Chlamydomonas reinhardtii*

Hugh Douglas Goold<fnref n="1" />,²,³,⁴,*, Hoa Mai Nguyen<fnref n="1" />,²,³,*, Fantao Kong<fnref n="1" />,²,³, Audrey Beyly-Adriano<fnref n="1" />,²,³, Bertrand Légeret<fnref n="1" />,²,³, Emmanuelle Billon<fnref n="1" />,²,³, Stéphan Cuiné<fnref n="1" />,²,³, Fred Beisson<fnref n="1" />,²,³, Gilles Peltier<fnref n="1" />,²,³ & Yonghua Li-Beisson<fnref n="1" />,²,³

Microalgae have emerged as a promising source for biofuel production. Massive oil and starch accumulation in microalgae is possible, but occurs mostly when biomass growth is impaired. The molecular networks underlying the negative correlation between growth and reserve formation are not known. Thus isolation of strains capable of accumulating carbon reserves during optimal growth would be highly desirable. To this end, we screened an insertional mutant library of *Chlamydomonas reinhardtii* for alterations in oil content. A mutant accumulating five times more oil and twice more starch than wild-type during optimal growth was isolated and named constitutive oil accumulator 1 (*coa1*). Growth in photobioreactors under highly controlled conditions revealed that the increase in oil and starch content in *coa1* was dependent on light intensity. Genetic analysis and DNA hybridization pointed to a single insertional event responsible for the phenotype. Whole genome re-sequencing identified in *coa1* a >200 kb deletion on chromosome 14 containing 41 genes. This study demonstrates that, 1), the generation of algal strains accumulating higher reserve amount without compromising biomass accumulation is feasible; 2), light is an important parameter in phenotypic analysis; and 3), a chromosomal region (Quantitative Trait Locus) acts as suppressor of carbon reserve accumulation during optimal growth.

With increasing world population and diminishing fossil fuels, alternative fuel sources that are renewable and do not compete with food production for fresh water and arable land are highly sought after<fnref n="1" />,². Due to their high growth rate and their intrinsic capacity (at least under certain growth conditions) to store large amounts of carbon rich compounds (i.e. starch and oils), microalgae have been investigated intensively as sources for biofuel production (bioethanol from starch; biodiesel from lipids)<fnref n="1" />⁻³. However, most algal species exhibit a maximal accumulation of starch and oil (i.e. triacylglycerols, TAGs) only when exposed to stress conditions, such as a nitrogen starvation, salinity, or high temperatures⁴⁻⁹.

The transition from optimal growth to carbon reserve accumulation requires changes at several systems levels including transcriptome, proteome and metabolome, highlighting the tight control that a cell exerts on the repartition of cellular carbon and energy usage⁸,¹⁰⁻¹⁶. Various genetic approaches have been applied to isolate factors involved in this transition. These reported regulatory proteins of lipid metabolism included a SQUAMOSA

Received: 11 January 2016

Accepted: 13 April 2016

Published: 04 May 2016

<footnote n="1">¹CEA, BIAM, Lab Bioenerget Biotechnol Bacteries & Microalgues, Saint-Paul-lez-Durance, 13108, France. ²CNRS, UMR 7265 Biol Veget & Microbiol Environ, Saint-Paul-lez-Durance, 13108, France. ³Aix Marseille Université, BVME UMR7265, Marseille, 13284, France. ⁴Faculty of Agriculture and the Environment, University of Sydney, Australia.
*These authors contributed equally to this work. Correspondence and requests for materials should be addressed to Y.L.-B. (email: yonghua.li@cea.fr)</footnote>
promoter-binding protein domain transcription factor (the nitrogen-responsive regulator 1, NRR1)¹³, a CXC domain-containing protein (the compromised hydrolysis of triacylglycerols 7, CHT7)¹⁷, and two members of the dual-specificity tyrosine phosphorylation-regulated kinase family, i.e. the triacylglycerol accumulation regulator1 (TAR1)¹⁸ and the starch degradation 1 (STD1)¹⁹. CHT7 is involved in the regulation of oil remobilization, however all three other proteins have been found to control oil amount under nitrogen starvation conditions. These three mutants (*nrr1*, *tar1*, and *std1*) provided regulatory clues as to how cells control oil accumulation in response to nitrogen starvation. Thus, these studies have mostly focused on the transitional changes from a healthy cell to a stressed cell. To our knowledge, no studies have been dedicated to investigate carbon reserve accumulation during optimal growth thus no regulatory gene or protein involved has been reported.

Under environmental conditions favorable for growth, genes of the lipid and starch biosynthetic pathways are expressed at low levels¹⁰, suggesting the occurrence of genetic program that suppresses carbon reserve formation. In this study, we applied a forward genetic screen in the green alga *Chlamydomonas reinhardtii* aiming at identifying new regulatory genes from the screening of mutants with altered cellular oil content during optimal growth. We used *C. reinhardtii* as a model because it has been proven as a suitable alga to study both starch and lipid metabolism²⁰⁻²⁵. We report here the detailed characterization of the mutant *coa1* standing for constitutive oil accumulator 1. This mutant accumulated increased oil and starch relative to the wild-type strain in a light-dependent manner. The isolation of this mutant should provide a molecular basis for our understanding of the control of TAG and starch accumulation during optimal growth.

## Results and Discussion

**Isolation of a Mutant with Higher Oil and Starch Content during Optimal Growth.** As part of our effort to dissect lipid metabolic pathways and identify novel factors involved in carbon reserve formation, an insertional mutant library has been generated as described previously²⁶. The parental line used here is the wild-type CC124. This mutant library was then screened for alterations in oil content during optimal growth conditions, according to the screening method described previously²⁷. Overall, after screening 7,000 independent transformants by Nile red coupled to Flow cytometry, a mutant (initially named 6D8) showing higher level of Nile red fluorescence than WT under normal growth condition was isolated. Higher Nile red fluorescence is generally used as a probe for increased cellular oil content, and indeed, TAG quantification confirmed that the mutant 6D8 made 5 times more oil than WT under nitrogen-replete conditions (Fig. 1A). In addition to an increase in oil content, an increased number of lipid droplet (LDs) was observed in the mutant cells, as revealed by staining with two independent lipophilic dyes i.e. Bodipy (Fig. 1C,D) or Nile red (Fig. 1E,F). Furthermore, the mutant 6D8 accumulated 1.5 times more oil than WT after being starved for nitrogen for 3 days (Fig. 1B). Based on these findings, the mutant 6D8 was re-named as constitutive oil accumulator 1 (*coa1*). No significant difference could be detected in the major membrane lipids (Supplemental Figure S1A,B). Additionally, we also observed that during the recovery phase following a period of nitrogen starvation, the mutant was also impaired in oil remobilization (Supplemental Figure S2). This result could imply that the observed net increase in oil content in the mutant during optimal growth could be, at least partly, due to a defect in oil breakdown. But this and other possibilities should be tested in the future with labelling studies to measure the rate of oil turnover during optimal growth.

Furthermore, the *coa1* mutant made almost twice more starch than its parental strain CC124 during normal growth (Fig. 2A), an increase from ~25µg mm⁻³ to over ~50µg mm⁻³. It is well known that the wild-type strains of *C. reinhardtii* store starch as a major form of polysaccharide reserve²⁰, for reasons of simplicity, we have referred to polysaccharide reserves as starch throughout this paper. However it is worth noting here that the additional accumulation of polysaccharides observed in the mutant *coa1* could be starch or any other type of polysaccharide made of α-glucans, or a mixture of both. Under the same cultivation condition, *coa1* mutant contained less chlorophyll than WT (Fig. 2B). The parallel increase in intracellular oil and starch observed in the *coa1* mutant suggests that a common mechanism underlies regulation of these two storage compounds. This regulatory mechanism would be affected in the mutant, resulting in an increased carbon flux towards reserve formation.

**The *coa1* Mutant Forms Palmelloid Structures but Shows Similar Growth Rate as WT.** During cell counting, the *coa1* mutant was found to form aggregates, showing much bigger particles in diameter than WT (Fig. 3A). The aggregated clusters were also observed under a microscope (Fig. 3B). Such a phenotype has been previously described in *Chlamydomonas* mutant strains, and these are referred to as palmelloid colonies²⁸. Formation of palmelloid colonies seem to be common²⁹, and has been suggested as part of microalgal self-defense against stressors³⁰. In an attempt to avoid cell aggregation, strategies were employed to assess whether this phenotype was due to growth conditions, or whether it was intrinsic to the mutant. Colonies were picked and re-plated successively in order to recover a single cellular organization, and also floating cells were recovered from the culture medium by allowing cells to settle at room temperature without shaking for an hour. These isolated cells were transferred to new cultures, thus selecting only those individual cells that are lighter and motile with active flagella, however the growth phenotype persisted. Due to the difficulty in accurately counting clustered cell aggregates, cellular volume has been used for cellular quantification rather than cell numbers. On a total cellular volume basis, similar growth rates were observed between the *coa1* mutant and WT (Fig. 3C).

**The High-Oil and High-Starch Phenotype of the *coa1* Mutant is Dependent on Light.** When comparing oil contents of wild-type and mutant cells grown in flasks, we sometimes observed big variations in measurements despite the fact that cells were apparently cultivated under similar conditions. We then hypothesized that such variations may be related to a light effect, as the light absorbed by algal cells is greatly reduced during the time course of a batch culture performed in flasks, due to a shading effect as cell density gets higher³¹. In order to test this hypothesis, we cultivated the mutant and WT in well-controlled conditions of illumination, by using photobioreactors (PBRs) operated as turbidostats under photoautotrophic conditions. In this experimental
<figure><img src="image_4.png" /><figcaption>**Figure 1.** The constitutive *accumulator 1* (*coa1*) mutant of *Chlamydomonas reinhardtii* over-accumulates oil during optimal growth in shake-flask cultures. (A) Triacylglycerol content as quantified by thin layer chromatography (TLC). Data are means of five biological replicates with 95% confidence intervals shown. Stars denote significant increase (Student's t-test; *p* < 0.05). (B) Oil content in WT and *coa1* (6D8) mutant strains after nitrogen depletion for 3 days. Data are means of three biological replicates and two technical replicates, with 95% confidence intervals shown. Stars denote significant increase (Student's t-test; *p* < 0.05). (C,D) Images of Bodipy-stained cells of WT and the mutant *coa1*. Scale bars = 50 µm. (E,F) Magnification of lipid droplets stained by Nile red. Scale bars = 5 µm. Cells were cultivated in shake flasks in standard TAP medium at a light density of 100 µmol photons m⁻² s⁻¹. Because the *coa1* mutant forms cell clusters, biochemical quantification is made based on total cellular volume instead of cell numbers. It is worth noting that for non-aggregated cells of *Chlamydomonas*, 1 mm⁻³ cellular volumes are equivalent to ~5 million cells.</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>**Figure 2.** The *coa1* mutant accumulates twice more starch than WT during optimal growth in shake-flask cultures. (A) Starch content (B) Chlorophyll content. Data are means of six replicates with 95% confidence intervals shown. Stars denote significant increase (Student's t-test; *p* < 0.05). Cells were cultivated in shake flasks in standard TAP medium at a light density of 100 µmol photons m⁻² s⁻¹.</figcaption></figure>

set-up, cell density was assessed by measuring the turbidity (OD₈八十nm) and was maintained at a constant value (OD₈八十nm = 0.4 corresponding to 2 million cells mL⁻¹) by step-wise injection of fresh culture medium. Using this system, photoautotrophic cultures can be maintained at physiological steady states under different light regimes. It should be mentioned here that because PBRs are radially illuminated and cell density maintained at a relatively low level, the light required in this set-up to observe physiological effects of the same magnitude as in flasks are much lower (a factor of 4 is observed: i.e. 40 µmol photons m⁻² s⁻¹ in the PBR set-up being approximately equivalent to ~160 µmol photons m⁻² s⁻¹ in shake flask culture). Increasing light intensity from 40 to 120 µmol photons m⁻² s⁻¹ had very little effect on intracellular oil levels in the WT strain, but led to an increase in oil content from 0.5 to 2.5 µg mm⁻³ in mutant cells (Fig. 4A).

Intracellular starch increased in both WT and the *coa1* mutant with increasing light intensity, but the increase was more pronounced in the mutant (Fig. 4B). Such an effect of light on starch accumulation has been previously reported for different species, including duckweed³², *Scenedesmus* sp.³³ and *Chlorella* sp.³¹, and also in *Chlamydomonas reinhardtii*³⁴, however a physiological explanation has not been put forward. When cells were grown under low light, the chlorophyll content was about 50% lower in the *coa1* mutant than in the WT (Fig. 4C). However, the efficiency of the photosynthetic machinery was only partially affected in *coa1*, as shown by photosystem II yields measured under various light intensities by means of chlorophyll fluorescence (Supplemental Figures S3). As generally observed during high light acclimation³⁵, the chlorophyll content of both strains decreased (Fig. 4C), reaching a similar value under high light. It is well known that when algae or plants are grown under high light (HL) intensity, photosystems develop smaller antennae sizes and contain lower amounts of Chlorophyll *b* (Chl *b*) thus resulting in an increase in the chlorophyll *a/b* (Chl *a/b*) ratio³⁶. An increase in the Chl *a/b* was observed in both strains under HL, the Chl *a/b* ratio being however higher in the *coa1* mutant than in the WT (Fig. 5A). It is important to note that when cells were grown in the dark both strains displayed similar Chl *a/b* ratios (Supplemental Figure S4) showing that this effect is related to light acclimation. Altogether, these data suggest that the *coa1* mutant is affected in its ability to regulate the antenna size and the chlorophyll content in response to HL conditions. This interpretation is supported by the immunoblot analysis which shows lower amounts of the light harvesting complex (LHCII) antenna protein LHCB5 in the *coa1* in comparison to the WT, with LHCB5 decreasing in response to HL acclimation in both strains (Supplemental Figure S5).

The increased starch and the decreased chlorophyll content observed in the mutant together suggest that the *coa1* mutant is somehow more sensitive to light as we have observed repeatedly, i.e. behaving under lower light similarly to WT under high light. Moreover, light is found to trigger higher oil accumulation in the mutant under increasing light intensity, whereas oil content increased only slightly in the WT under the light range tested. The *coa1* mutant therefore provides a means to probe the link between light intensity and reserve formation, a phenomenon which has sometimes been observed in the literature³⁷ but has rarely been explored at a mechanistic level.

**The Phenotype of the *coa1* Mutant Results from a Single Insertional Event.** Insertional mutagenesis can lead to insertion of several copies of DNA cassettes randomly in a genome³⁸,³⁹. As a first step to determine whether the *coa1* phenotype is genetically linked to the insertion of the paromomycin cassette, we performed a DNA/DNA hybridization analysis (Southern blot). After digesting the *coa1* genomic DNA with several restriction enzymes independently, a single hybridizing band was observed in all cases using a probe specific to the DNA cassette harboring paromomycin gene *AphVIII*. This suggested that there is only one insertion of *AphVIII* cassette in the *coa1* genome (Fig. 6A).
A.

<figure><img src="image_6.png" /></figure>

B.

<figure><img src="image_7.png" /></figure>

C.

<figure><img src="image_8.png" /></figure>

**Figure 3.** Growth characteristics of the *coa1* mutant and WT. (A) Population of WT and *coa1* mutant cells was monitored using a cell counter. The mutant cells appear much bigger (~9–10µm) than WT cells (~5–6µm) under the same cultivation conditions. (B) Clusters of the *coa1* mutant cells were observed under light microscope. Bars = 10µm. (C) Cell growth for WT and *coa1* mutant cells were monitored based either on cell number (solid lines) or on cellular volume (dashed lines). Cells were grown in shake flasks in standard TAP medium at a light density of 100µmol photons m⁻² s⁻¹. Error bars represent standard deviation based on three biological replicates. Abbreviations: TAP, Tris-Acetate Phosphate.

As a second step, we performed genetic backcrosses between *coa1* (*mt*⁻) and the wild-type strain CC125 (*mt*⁺) which is the *mt*⁻ strain mostly closely related to the *coa1*'s progenitor (CC124). From the analysis of 3 complete tetrads (Fig. 6B) and 9 incomplete tetrads (Supplemental Figure S6A), we found that all paromomycin-resistant strains accumulated higher oil amount than paromomycin-sensitive strains and showed persistent palmelloid colony formation. Chlorophyll content in the three complete tetrads also showed the same segregation as cellular oil content, i.e. those progenies sensitive to paromomycin also contained less chlorophyll, which was already obvious from visual inspection of the color of the culture (Supplemental Figure S6B). No strain issued from this cross harbored the constitutive TAG accumulating phenotype in the absence of paromomycin resistance, thus indicating a genetic link between the oil phenotype and the single insertion of the antibiotic cassette.
<figure><img src="image_9.png" /><figcaption>**Figure 4.** The *coa1* phenotype is positively linked to higher light intensity in photobioreactors operated as turbidostat. (A) TAG content in cells cultivated under increasing light irradiance. (B) Starch content in cells cultivated under increasing light irradiance. (C) Chlorophyll content in cells cultivated under increasing light irradiance. It is important to note that 40 µmol photons m⁻² s⁻¹ in the PBR set-up being approximately equivalent to 160 µmol photons m⁻² s⁻¹ in the flask batch set-up, roughly a 4-fold change. Data are means of three replicates, and error bars denote 95% confidence intervals. Cells were cultivated in PBRs under strict photoautotrophic conditions maintaining a constant OD₈₈₀nm = 0.4 (eq. = 2 million cells mL⁻¹ for WT).</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>**Figure 5.** Alterations in the ratio of chlorophyll a versus chlorophyll b (Chl a/b). (A) Difference in Chl a/b ratio and its evolution in response to increasing light intensity. (B) The chlorophyll cycle and the requirement for a chlorophyll b reductase. Data are means of three biological replicates with two technical replicates each. Error bars represent standard deviations. Chl: chlorophyll.</figcaption></figure>

**Whole Genome Re-Sequencing Identifies a Large Deletion in Chromosome 14.** Although Southern blot analysis revealed a single insertional event occurred in the *coa1* mutant (Fig. 6A), two distant sites of insertion were identified on chromosome 14 by sequencing flanking regions of the insertion (one at 69,200 bp and the other at 315,218 bp) using Genome Walker technique. In order to understand how a single insertion event could lead to the identification of two distant sites, we re-sequenced the whole genome of *Chlamydomonas*. The genome sequencing generated 34257132 clean reads, with 90 bp per read and has a calculated coverage of 23.7x the whole *Chlamydomonas* genome (Table 1).
<figure><img src="image_11.png" /><figcaption>**Figure 6.** Genetic analyses confirm a single insertion in the *coa1* genome. (A) DNA blotting of digested genomic DNA of WT and *coa1* mutant with three different enzymes and then hybridized with the *AphVIII* probe. (B) Oil content for the progenies issued from tetrad analyses of *coa1* and CC124. (C) Whole genome re-sequencing identifies >200 kb deletion in chromosome 14 between coordinates 69,266 and 317,163. T: tetrad.</figcaption></figure>

<table><caption>**Table 1.** Summary of whole genome re-sequencing results. Q20 is considered a high quality data commonly used by sequencing facilities. It indicates a 99% certainty that the base has been called correctly. Sequencing was performed by the Beijing Genome Institute (BGI, China).</caption><thead><tr><td>Insert size</td><td>Read length</td><td>Clean reads</td><td>Clean bases</td><td>Q20(%)</td><td>GC(%)</td></tr></thead><tbody><tr><td>500 bp</td><td>90 bp</td><td>34257132</td><td>3083141880</td><td>92.04</td><td>60.36</td></tr></tbody></table>

Genomic reads were paired using FASTQ joiner⁴⁰ and mapped to the *C. reinhardtii* genome⁴¹ (version 4.3) using the BAM tools package available on the Galaxy project website (https://main.g2.bx.psu.edu/). Global and local coverage was determined by analysis of BAM output by GATK (Genome Analysis Tool Kit). The genomic region between 69,200 bp and 315,218 bp compared to the reference *Chlamydomonas* genome version 4.3 revealed low read coverage in the region of the insertion/deletion (mean coverage depth of 1.6, compared to a mean coverage depth of 25.4 for the entire chromosome 14). Peaks of high read coverage in the region between 69,200 bp and 317,163 bp on the reference genome were identified as highly repetitive motifs found elsewhere in the *Chlamydomonas* genome. Furthermore, reads mapped to the 3' and 5' end of the gene cassette, and reads
paired to those which mapped to the cassette revealed that the genomic loci adjacent to the gene cassette had sequence identity to 69,200 bp and 315,218 bp on chromosome 14.

Taken together, the whole genome re-sequencing analysis and flanking genomic DNA sequencing by Genome Walker revealed that a single complete antisense *AphVIII* cassette was inserted in a 3' –>5' orientation between the 69,200 bp and 317,163 bp on chromosome 14, resulting in a substantial genomic deletion of >200 kb in chromosome 14 between the coordinates 69,266 and 317,163 (Fig. 6C). Chromosomal deletions caused by insertional mutagenesis in *C. reinhardtii* have previously been reported¹⁷,³⁵. Within the >200 kb region, 41 putative protein-coding genes have been found (Table 2). The phenotype associated with *coa1* could thus be attributed to either the result of one of the genes removed in the deleted locus of the genome, or a combinatorial effect of multiple genes. This and other possibilities are discussed in the next section.

**Suppressors of Carbon Reserve Accumulation Within Chromosome_14:69,266.317,163.** The link between the lesion on chromosome 14 and the observed increased level of TAG and starch under high light conditions suggest that carbon reserve formation is repressed in optimal growth in wild-type strains. In the absence of this chromosomic region (the *coa1* mutant), the cellular homeostasis of carbon reserves is disturbed, thus creating strains that over-accumulate oil and starch simultaneous to growth. Furthermore, it is found that one or more of the genes encoded on the same region may be key in the orchestration of the light response in *C. reinhardtii*. The reason *Chlamydomonas* exhibits an increased proportion of TAG and starch during high light growth remains poorly defined, and may remain difficult to analyze without the access to genome editing in *Chlamydomonas*, however the generation of an indexed genome-wide mutant library for *C. reinhardtii* (http://jonikaslab.dpb.carnegiescience.edu/chlamy-mutant-library )⁴² may aid in more rapidly elucidating which of these genes encode key functions in regulation of TAG homeostasis during high light.

While there is little information known about the proteins in the locus of interest (Table 2), some proteins do deserve particular note due to their potential involvement in the observed phenotype. For example, the presence of a lipase-like protein, a putative lipid transfer protein or a Coenzyme A binding protein in this region could perhaps contribute to the observed oil content phenotype in *coa1* mutant. The presence of a Myb domain-containing putative DNA binding protein is another example. Homologs of these proteins have been identified as potential key players in the carbon concentration mechanism (CCM) as transcription factors involved in the response to the physiological pressure of low carbon growth⁴³. A gene encoding a chlorophyllide *b* reductase (Cre14.g608800) is also located in this chromosomic region. Chlorophyllide *b* reductase catalyzes the first step in the degradation of chlorophyll *b* in the chlorophyll cycle⁴⁴. Its deletion in the *coa1* mutant may have perturbed the chlorophyll cycle thus explaining the altered chlorophyll content and chlorophyll *a/b* ratio observed in the mutant (Fig. 5B). Increase in chlorophyll *a/b* ratio is an integral feature of cells' acclimation to high light conditions⁴⁵. Although molecular mechanisms remain to be elucidated, the mutant *coa1* displays several phenotypic features (such as decreased chlorophyll content, increased chlorophyll *a/b* ratio and increased intracellular starch), indicating a higher sensitivity to light that may contribute to an increased overall flux in cellular carbon and an increase in oil content.

**Outlook.** One of the major obstacles to economically feasible production of algal fuel is the requirement of nitrogen starvation or other stress conditions for higher reserve accumulation². Here we isolated a mutant of *Chlamydomonas*, *coa1*, which accumulates five times more oil and twice more starch than the WT during optimal growth. To our knowledge, this is the first report of such a mutant accumulating higher oil amount in parallel to optimal growth for the green alga *C. reinhardtii*. Through the study of the mutant *coa1*, we also pointed out a missing link underlying the regulation of reserve accumulation under high light conditions. Molecular characterization of the *coa1* mutant identified a QTL responsible for the observed phenotype. Detailed examination of the genes, via cross-references to published transcriptomic datasets (Table 2)³⁴, encoded in the missing region should help understanding the molecular mechanisms that led to increased amount of carbon reserves in the mutant, and should provide molecular tools for uncoupling lipid accumulation from impairment in cell growth. This report in conjunction with the availability of genome-wide indexed mutant libraries (http://jonikaslab.dpb.carnegiescience.edu/chlamy-mutant-library )⁴², and systems analysis of cells' response to various stresses¹⁰,¹³,³⁴, should make it faster to assign the *coa1* phenotype to particular gene(s) located within this region on chromosome 14. Thus the isolation of the *coa1* mutant demonstrates that mutants over-accumulating carbon reserves under non-stress conditions can be isolated and further provides a list of candidate genes involved in suppressing carbon reserve accumulation during optimal growth.

**Methods**

**Strains, Growth Conditions and Generation of Insertional Mutants.** *Chlamydomonas reinhardtii* wild-type strain CC124 (*mt*⁻ *nit1 nit2*), obtained from the Chlamydomonas resource center (http://chlamycol-lection.org), was used to generate the insertional mutant library in our laboratory as described previously in Nguyen *et al*.²⁶. Briefly, the antibiotic resistance gene *AphVIII* was inserted into the genome of *C. reinhardtii* via electroporation⁴⁶. Paromomycin-resistant clones were screened for alterations in oil content by Nile red staining coupled to Flow Cytometry as described in detail in Cagnon *et al*.²⁷. Unless otherwise stated, all strains were cultivated at 25°C in Tris-acetate-phosphate media (TAP) in shake flasks kept in incubators (Multitron, Infors HT, Switzerland) under continuous illumination (120 µmol photons m⁻² s⁻¹) with shaking (120 rpm)²⁸. For nitrogen (N) starvation experiments, log-phase grown cells were centrifuged for 3 min at 1000 g, cell pellets were washed twice in TAP-N media, then resuspended into TAP-N for 3 days.

**Cell Counting and Microscopy.** Cell numbers, sizes and total cellular volumes were quantified by the use of a Beckman Coulter Multisizer 4 (Multisizer™ Coulter Counter, Beckman Coulter, USA).
<table><caption>**Table 2.** List of the putative genes encoded in the chromosome region deleted in *coal1* and their relative expression under high light or nitrogen starvation. The listed putative genes were encoded between the coordinates 69,266 and 317,163 in chromosome 14. Values are log2-transformed average transcript reads of biological replicates for the light experiment (LL: 41 µmol photons m⁻² s⁻¹ versus HL: 145 µmol photons m⁻² s⁻¹)³⁴, and are normalized mRNA abundance in RPKM (stands for: reads per kilobase of transcript per million reads mapped) for nitrogen (N) starvation experiment¹³. Abbreviations: N, nitrogen; h, hour; HL, high light; ‘-’ no annotation or no known domain identified; N.D. refers to those genes whose expression were not detected.</caption><thead><tr><th rowspan="2">Gene ID<br>(Phytozome)</th><th rowspan="2">Annotated name or putative functions</th><th colspan="2">LL6h</th><th colspan="2">+N</th><th colspan="2">-N8h</th></tr><tr><th colspan="2">Based on that of<sup>34</sup></th><th colspan="2">Based on that of<sup>13</sup></th><th colspan="2"></th></tr></thead><tbody><tr><td>Cre14.g608200</td><td>Actin binding protein</td><td>11.09</td><td>10.12</td><td>1.04</td><td>1.64</td><td></td><td></td></tr><tr><td>Cre14.g608250</td><td>Type II secretion system protein</td><td>4.86</td><td>4.19</td><td>0.54</td><td>0.31</td><td></td><td></td></tr><tr><td>Cre14.g608275</td><td>-</td><td></td><td></td><td>1.92</td><td>2.89</td><td></td><td></td></tr><tr><td>Cre14.g608300</td><td>Proline rich extensin family</td><td></td><td></td><td>2.33</td><td>2.61</td><td></td><td></td></tr><tr><td>Cre14.g608350</td><td>-</td><td>4.14</td><td>4.93</td><td>0.64</td><td>N.D.</td><td></td><td></td></tr><tr><td>Cre14.g608400</td><td>-</td><td>6.94</td><td>6.71</td><td>0.60</td><td>0.67</td><td></td><td></td></tr><tr><td>Cre14.g608452</td><td>-</td><td></td><td></td><td>7.73</td><td>4.63</td><td></td><td></td></tr><tr><td>Cre14.g608500</td><td>Carbonic anhydrase</td><td></td><td></td><td>1.17</td><td>0.84</td><td></td><td></td></tr><tr><td>Cre14.g608550</td><td>-</td><td></td><td></td><td>1.37</td><td>2.81</td><td></td><td></td></tr><tr><td>Cre14.g608600</td><td>-</td><td></td><td></td><td>1.85</td><td>3.28</td><td></td><td></td></tr><tr><td>Cre14.g608652</td><td>-</td><td></td><td></td><td>N.D.</td><td></td><td></td><td></td></tr><tr><td>Cre14.g608700</td><td>2-amino adipate aminotransferase</td><td></td><td></td><td>2.54</td><td>1.60</td><td></td><td></td></tr><tr><td>Cre14g608750</td><td>Gametogenesis family protein</td><td>3.08</td><td>3.46</td><td>1.76</td><td>1.71</td><td></td><td></td></tr><tr><td>Cre14.g608788</td><td>Arabinogalactan protein</td><td></td><td></td><td>1.73</td><td>1.70</td><td></td><td></td></tr><tr><td>Cre14.g608800</td><td>Chlorophyll(ide) b reductase</td><td></td><td></td><td>10.26</td><td>4.36</td><td></td><td></td></tr><tr><td>Cre14.g608850</td><td>COP9 signalosome</td><td></td><td></td><td>6.50</td><td>3.20</td><td></td><td></td></tr><tr><td>Cre14.g608900</td><td>Lipid transfer protein</td><td></td><td></td><td>0.38</td><td>1.33</td><td></td><td></td></tr><tr><td>Cre14.g608950</td><td>-</td><td>1.76</td><td>1.57</td><td>5.97</td><td>4.76</td><td></td><td></td></tr><tr><td>Cre14.g608970</td><td>Ribosomal protein L34</td><td></td><td></td><td>10.67</td><td>5.09</td><td></td><td></td></tr><tr><td>Cre14.g608990</td><td>-</td><td></td><td></td><td>4.51</td><td>5.55</td><td></td><td></td></tr><tr><td>Cre14.g609010</td><td>LRR (leucine repeat) family protein</td><td></td><td></td><td>1.45</td><td>0.33</td><td></td><td></td></tr><tr><td>Cre14.g609030</td><td>COP9 signalosome</td><td></td><td></td><td>10.48</td><td>9.58</td><td></td><td></td></tr><tr><td>Cre14.g609050</td><td>-</td><td></td><td></td><td>12.82</td><td>3.50</td><td></td><td></td></tr><tr><td>Cre14.g609100</td><td>LRR (leucine repeat) family</td><td>6.46</td><td>6.95</td><td>3.68</td><td>2.49</td><td></td><td></td></tr><tr><td>Cre14.g609150</td><td>-</td><td>9.18</td><td>9.38</td><td>1.76</td><td>2.83</td><td></td><td></td></tr><tr><td>Cre14.g609202</td><td>Matrix metalloproteinase</td><td></td><td></td><td>1.18</td><td>2.35</td><td></td><td></td></tr><tr><td>Cre14.g609250</td><td>Lipase/acetylhydralase</td><td></td><td></td><td>0.48</td><td>0.85</td><td></td><td></td></tr><tr><td>Cre14.g609300</td><td>N-acetylglucosamine-6-phosphate deacetylase</td><td></td><td></td><td>4.71</td><td>0.65</td><td></td><td></td></tr><tr><td>Cre14.g609350</td><td>Myb-like DNA binding domain</td><td></td><td></td><td>6.09</td><td>4.08</td><td></td><td></td></tr><tr><td>Cre14.g609400</td><td>Peptide methionine sulfoxide reducatse</td><td>12.01</td><td>11.78</td><td>12.27</td><td>4.79</td><td></td><td></td></tr><tr><td>Cre14.g609450</td><td>Predicted PWR protein</td><td></td><td></td><td>12.66</td><td>4.11</td><td></td><td></td></tr><tr><td>Cre14.g609500</td><td>Protein of unknown function</td><td></td><td></td><td>3.11</td><td>1.25</td><td></td><td></td></tr><tr><td>Cre14.g609551</td><td>-</td><td></td><td></td><td>2.11</td><td>1.55</td><td></td><td></td></tr><tr><td>Cre14.g609600</td><td>-</td><td>5.88</td><td>5.69</td><td>6.93</td><td>5.64</td><td></td><td></td></tr><tr><td>Cre14.g609650</td><td>CoA binding protein</td><td>10.27</td><td>10.71</td><td>6.06</td><td>9.92</td><td></td><td></td></tr><tr><td>Cre14.g609700</td><td>DOMON-domain containing protein</td><td></td><td></td><td>6.04</td><td>3.95</td><td></td><td></td></tr><tr><td>Cre14.g609750</td><td>Cyclin-dependent kinase inhibitor</td><td>1.69</td><td>1.51</td><td>1.56</td><td>0.39</td><td></td><td></td></tr><tr><td>Cre14.g609800</td><td>-</td><td>1.42</td><td>1.31</td><td></td><td></td><td></td><td></td></tr><tr><td>Cre14.g609850</td><td>-</td><td>9.08</td><td>9.36</td><td>4.34</td><td>2.33</td><td></td><td></td></tr><tr><td>Cre14.g609900</td><td>Predicted membrane protein</td><td></td><td></td><td>1.16</td><td>3.39</td><td></td><td></td></tr><tr><td>Cre14.g609926</td><td>-</td><td></td><td></td><td>19.25</td><td>21.37</td><td></td><td></td></tr></tbody></table>

To observe lipid droplets, cells were first stained with a solution of Nile red (at a final concentration of 0.1 µg mL⁻¹ in methanol; Sigma, Saint Louis, USA)⁴⁷ or Bodipy (Bodipy 515, Invitrogen) at a final concentration of 1.5 µg mL⁻¹ from a stock solution of 1 mg mL⁻¹ in DMSO⁴⁸. After incubation at room temperature for 20 min, cells were visualized under a Leica DMRXA epifluorescent microscope (Leica Microsystems, Germany). Flow cytometry measurement of Nile red fluorescence has previously been described in detail²⁷.

**Cultivation in Photobioreactors (PBRs).** PBR cultures were cultivated in minimal medium (MM, i.e. without buffer neither carbon source)²⁸ in automated ‘Biostat A Plus’ photobioreactor systems (Sartorius Stedim
Biotech) equipped with a biomass probe (Excell probe, Exner; measuring OD<sub>880 nm</sub> with a 2-cm light path) and operated as turbidostats. Cells were maintained at a constant density (OD<sub>880nm</sub> = 0.4) throughout experimentation by automated injection of fresh media (Stepdos FEM03TT18RC; KNF) and ejection of excess biomass through a waste tube. Culture pH was maintained constant by addition of 0.2 N HCl or 0.2 N KOH using the Biostat module. A gas mixture containing air plus 1.8% CO<sub>2</sub> was generated using two mass flow meters (EL flow; Bronkhorst) and bubbled into the PBRs at a flow rate of 0.5 L min<sup>-1</sup>. Eight fluorescent bulbs (Osram Dulux L 18 W) placed around the PBR were used to supply light in a range of light intensities up to 400 µmol photons m<sup>-2</sup> s<sup>-1</sup>.

**Lipid Extractions and Lipid Quantifications by Thin Layer Chromatography (TLC).** A modified Bligh and Dyer method⁴⁹ was used during this work to isolate total cellular lipids as described previously in Siaut *et al.*⁴. Briefly, samples were vortexed in quenching solution (1 mM EDTA, 0.15 M acetic acid) in a glass tube, then to which 3 mL of methanol:chloroform (2:1 v/v) was added. Samples were then vortexed for 10 min, 1 mL of chloroform and 0.8 mL of KCl (0.8% v/v) was added, and samples were vortexed briefly. Phases were separated by centrifugation at 600 g for 2 min at 4 °C in a bench top centrifuge. The lower phase was then isolated and 1 mL of hexane was added to re-extract the lipids, and samples were vortexed again for 10 min and phases were separated. The upper phase containing lipids was then transferred to the chloroform extract using a glass Pasteur pipette, and the rest was discarded. The combined chloroform and hexane extracts were then dried under a flow of nitrogen gas and re-suspended in an appropriate amount of solvent for subsequent analyses. Cellular content of lipid classes were quantified with Thin Layer Chromatography (TLC) using an automated "High-Performance Thin-Layer Chromatography" platform (CAMAG). These procedures were carried out as described previously⁴.

Starch and Chlorophyll Quantification.

**Starch and Chlorophyll Quantification.** A given number of cells (~2 million cells) were harvested by centrifugation (13000 g for 10 min) and cell pellets were resuspended in 1 mL methanol. The mixture was kept at -20 °C for at least 1 h before the organic phase was taken for chlorophyll quantification using spectrophotometer⁵⁰. The pellet was air dried under a fume hood, and resuspended in MilliQ water, and autoclaved for 20 min at 121 °C. Amyloglucosidase (0.2 U, Roche) was added and samples were incubated for 1 h at 60 °C. Samples were briefly centrifuged, and measured for glucose concentration using a YSI2700 select sugar analyzer (YSI Life Sciences, Yellow Springs, USA)⁴.

**Mating and Tetrad Dissection.** The *coal* mutant was back-crossed with the wild-type CC125 (*mt*<sup>+</sup>). Genetic crosses were performed following the protocol described previously⁵¹. Briefly, *Chlamydomonas* strains of *mt*<sup>-</sup> and *mt*<sup>+</sup> were cultivated in TAP medium until mid-log phase. Cells were gently centrifuged at 360 g for 3 min, washed twice with TAP-N medium and resuspended in the same medium. Cultures were kept in TAP-N for 24 h at 25 °C under shaking. Cells with opposite mating types were mixed together in a glass tube (1 mL: 1 mL) and then the tubes were kept in an incline position in a laminar hood for 3 hr. The mixture (1 mL) was then spread into center of TAP plates containing 3% agar (w/v), the plates were dried and kept under low light overnight. Then the plates were wrapped in aluminum foil, kept in the dark for 5-10 days to induce zygote germination upon induction by light. All the vegetative cells were scraped from the plate with a dull scalpel. To further eliminate the vegetative cells, the plates were exposed to chloroform vapor for 20-45 sec. The plates were then placed under low light for 24 h. The tetrads were dissected using a micromanipulator needle connected to a light microscopy under a laminar hood.

**DNA Extractions.** Total DNA was extracted and purified from exponentially grown *Chlamydomonas* cells using a Qiagen Plant DNeasy Max kit (with QIAshredder). Briefly, a 200-mL culture of exponentially grown cells in TAP media was centrifuged at 600 g for 5 min at 4 °C. Cells were snap frozen in liquid nitrogen, pellets were briefly pulverised with a sterile mortar and pestle. The rest of the procedure was following that described by the manufacture. Extracted DNA was quantified using a spectrophotometer, and DNA purity was calculated relative to the ratio of the absorbance of A260 nm/A280 nm.

**Southern Blot.** A fragment of the *AphVIII* gene was cloned by PCR (with primer AphVIII-F: (5'-CGAAGCATGGACGATGCGTT-3' and AphVIII-R: 5'-CGAGACTGCGATCGAACGGACA-3') and labelled with Digoxigenin-dUTP (DIG, Roche) to be used as a hybridization probe. Before hybridization, the probe was denatured by heating at 100 °C for 10 min, and then kept on ice for 5 min. Genomic DNA (~8 µg) was digested overnight by *Not*I, *Stu*I, or *Apa*I (NEB Biolabs). Digested DNA was loaded onto a 0.8% agarose gel, and each lane representing DNA cut by one specific enzyme and the gel was run slowly at 50 V to better separate fragmented DNA. Migration, transfer, and hybridization were described in detail in⁵². Detection was revealed with addition of a chemiluminescent substrate CSPD® (disodium 3-(4-methoxyspiro-3,2'-(5'-chloro)tricyclo(3.3.1.13,7)decan)-4-yl)phenyl phosphate ) (Roche), images were taken by G:BOX (Syngene). The number of hybridized bands in each lane indicates the potential number of *AphVIII* insertions in the genome.

**Genome Walker, Genome Re-Sequencing and Data Analyses.** To identify the position of cassette insertion in the genome, genome walking PCR was carried out using a modified protocol of the commercial Genome Walker™ kit (Clontech). In brief, DNA (1 µg) was digested using restriction enzymes *Fsp*I, *Nae*I, *Pml*I, *Pvu*II, or *Sma*I (New-England Labs). Digested DNA was then ligated to adapter sequences (5'-GTAATACGACTCACTATAGGGCACGCGTGGACGGCCCGGGCTGGT-3'). Nested PCRs were carried out using two sets of oligonucleotides, i.e. GWaph8fw1 (5'-CTGGTGCTGCGCGAGCTGGCCACGAGGAG-3') or GWaph8rev1 (5'-CCAGCGCGAGATCGGAG TGCCGGTCCG-3') with AP1 (5'-GTAATACGACTCACTATAGGGC-3'), and GWaph8fw2 (5'-TGGTTCGGGCCGGAGTGTTCCGCGGCGTT-3') or GWaph8rev2 (5'-CGAGACTGCGATC GAACGGACACCGC-3') with AP2 (5'-ACTATAGGGCACGCGTGGT-3'), to the cassette, and the adapter
sequence in a final volume of 50 µL. PCR product was then cloned and sequenced. The obtained sequence was then put through a BLAST against the genome of *C. reinhardtii* to Phytozome⁴¹ to identify the flanking region.

Highly pure DNA (~15 µg per sample at a concentration of 30 ng µL⁻¹; and two biological replicates were prepared; OD₂₆₀/₂₈₀ = 1.66) was sent to Beijing Genome Institute (BGI, China) for whole genome re-sequencing. Two Genomic read files in FASTQ (fq) format obtained from the Beijing Genome Institute (BGI, China) were uploaded to Galaxy (https://main.g2.bx.psu.edu/)⁵³ and renamed "data 1" and "data 2". Data1.fq and data2.fq were then processed using "FASTQ Groomer", found in the NGS Toolbox Beta software package. Data were then filtered using fastQ quality filter. Processed fastq files were then joined using FASTQ joiner, also found in the NGS Toolbox Beta suite of software⁴⁰. The *Chlamydomonas* genome v5 was uploaded as a template genome, the genome was downloaded as *C. reinhardtii*.110311.tar.gz from http://www.chlamy.org. Joined reads were then mapped using Bowtie with a gap length of 500 bp, and visualized with Trackster⁵³,⁵⁴.

## References

1. Radakovits, R., Jinkerson, R. E., Darzins, A. & Posewitz, M. C. Genetic engineering of algae for enhanced biofuel production. *Eukaryot. Cell* **9**, 486–501, doi: 10.1128/ec.00364-09 (2010).

2. Wijffels, R. H. & Barbosa, M. J. An outlook on microalgal biofuels. *Science* **329**, 796–799, doi: 10.1126/science.1189003 (2010).

3. Work, V. H., D'Adamo, S., Radakovits, R., Jinkerson, R. E. & Posewitz, M. C. Improving photosynthesis and metabolic networks for the competitive production of phototroph-derived biofuels. *Current Opinion in Biotechnology* **23**, 290–297 (2012).

4. Siaut, M. et al. Oil accumulation in the model green alga *Chlamydomonas reinhardtii*: characterization, variability between common laboratory strains and relationship with starch reserves. *BMC Biotechnology* **11**, 7 (2011).

5. Wang, Z. T., Ullrich, N., Joo, S., Waffenschmidt, S. & Goodenough, U. Algal lipid bodies: stress induction, purification, and biochemical characterization in wild-type and starchless *Chlamydomonas reinhardtii*. *Eukaryot Cell* **8**, 1856–1868 (2009).

6. Moellering, E. R. & Benning, C. RNA interference silencing of a major lipid droplet protein affects lipid droplet size in *Chlamydomonas reinhardtii*. *Eukaryot. Cell* **9**, 97–106, doi: 10.1128/ec.00203-09 (2010).

7. Légeret, B. et al. Lipidomic and transcriptomic analyses of *Chlamydomonas reinhardtii* under heat stress unveil a direct route for the conversion of membrane lipids into storage lipids. *Plant, Cell & Environment* n/a-n/a, doi: 10.1111/pce.12656 (2016).

8. Hemme, D. et al. Systems-wide analysis of acclimation responses to long-term heat stress and recovery in the photosynthetic model organism *Chlamydomonas reinhardtii*. *The Plant Cell Online* **26**, 4270–4297, doi: 10.1105/tpc.114.130997 (2014).

9. Takeshita, T. et al. Starch and lipid accumulation in eight strains of six *Chlorella* species under comparatively high light intensity and aeration culture conditions. *Bioresour. Technol.* **158**, 127–134, doi: 10.1016/j.biortech.2014.01.135 (2014).

10. Miller, R. et al. Changes in transcript abundance in *Chlamydomonas reinhardtii* following nitrogen deprivation predict diversion of metabolism. *Plant Physiology* **154**, 1737–1752, doi: 10.1104/pp.110.165159 (2010).

11. Schmollinger, S. et al. Nitrogen-sparing mechanisms in *Chlamydomonas* affect the transcriptome, the proteome, and photosynthetic metabolism. *Plant Cell* **26**, 1410–1435, doi: 10.1105/tpc.113.122523 (2014).

12. Goodenough, U. et al. The path to triacylglyceride obesity in the *sta6* strain of *Chlamydomonas reinhardtii*. *Eukaryot. Cell* **13**, 591–613, doi: 10.1128/ec.00013-14 (2014).

13. Boyle, N. R. et al. Three acyltransferases and nitrogen-responsive regulator are implicated in nitrogen starvation-induced triacylglycerol accumulation in *Chlamydomonas*. *Journal of Biological Chemistry* **287**, 15811–15825, doi: 10.1074/jbc.M111.334052 (2012).

14. Jinkerson, R. E., Radakovits, R. & Posewitz, M. C. Genomic insights from the oleaginous model alga *Nannochloropsis gaditana*. *Bioengineered* **4**, 0–1 (2013).

15. Park, J.-J. et al. The response of *Chlamydomonas reinhardtii* to nitrogen deprivation: a systems biology analysis. *The Plant Journal* **81**, 611–624, doi: 10.1111/tpj.12747 (2015).

16. Nguyen, H. M. et al. Proteomic profiling of oil bodies isolated from the unicellular green microalga *Chlamydomonas reinhardtii*: with focus on proteins involved in lipid metabolism. *Proteomics* **11**, 4266–4273, doi: 10.1002/pmic.201100114 (2011).

17. Tsai, C.-H. et al. The protein Compromised Hydrolysis of Triacylglycerols 7 (CHT7) acts as a repressor of cellular quiescence in *Chlamydomonas*. *Proceedings of the National Academy of Sciences* **111**, 15833–15838, doi: 10.1073/pnas.1414567111 (2014).

18. Kajikawa, M. et al. Algal dual-specificity tyrosine phosphorylation-regulated kinase, triacylglycerol accumulation regulator1, regulates accumulation of triacylglycerol in nitrogen or sulfur deficiency. *Plant Physiology* **168**, 752–764, doi: 10.1104/pp.15.00319 (2015).

19. Schulz-Raffelt, M. et al. Hyper-accumulation of starch and oil in a *Chlamydomonas* mutant affected in a plant-specific DYRK kinase. *Biotechnology for Biofuels* **9**, 55, doi: 10.1186/s13068-016-0469-2 (2016).

20. Tunçay, H. et al. A forward genetic approach in *Chlamydomonas reinhardtii* as a strategy for exploring starch catabolism. *PLoS ONE* **8**, e74763, doi: 10.1371/journal.pone.0074763 (2013).

21. Ball, S. G. The intricate pathway of starch biosynthesis and degradation in the monocellular alga *Chlamydomonas reinhardtii*. *Australian Journal of Chemistry* **55**, 49–59, doi: 10.1071/ch02052 (2002).

22. Hicks, G. et al. When simpler is better. Unicellular green algae for discovering new genes and functions in carbohydrate metabolism. *Plant Physiol* **127**, 1334–1338 (2001).

23. Li-Beisson, Y., Beisson, F. & Riekhof, W. Metabolism of acyl-lipids in *Chlamydomonas reinhardtii*. *Plant J.* **82**, 504–522, doi: 10.1111/tpj.12787 (2015).

24. Liu, B. & Benning, C. Lipid metabolism in microalgae distinguishes itself. *Current Opinion in Biotechnology* **24**, 300–309, 10.1016/j.copbio.2012.08.008 (2013).

25. Merchant, S. S., Kropat, J., Liu, B., Shaw, J. & Warakanont, J. TAG, You're it! *Chlamydomonas* as a reference organism for understanding algal triacylglycerol accumulation. *Current Opinion in Biotechnology* **23**, 352–363 (2012).

26. Nguyen, H. M. et al. The green microalga *Chlamydomonas reinhardtii* has a single ω-3 fatty acid desaturase that localizes to the chloroplast and impacts both plastidic and extraplastidic membrane lipids. *Plant Physiology* **163**, 914–928, doi: 10.1104/pp.113.223941 (2013).

27. Cagnon, C. et al. Development of a forward genetic screen to isolate oil mutants in the green microalga *Chlamydomonas reinhardtii*. *Biotechnology for Biofuels* **6**, 178 (2013).

28. Harris, E. H. *The Chlamydomonas sourcebook: introduction to Chlamydomonas and its laboratory use*. 2nd edn, (Elsevier, 2009).

29. Matsuo, T. et al. A systematic forward genetic analysis identified components of the *Chlamydomonas* circadian system. *Genes & Development* **22**, 918–930, doi: 10.1101/gad.1650408 (2008).

30. Lurling, M. & Beekman, W. Palmelloids formation in *Chlamydomonas reinhardtii*: defence against rotifer predators? *Ann. Limnol. - Int. J. Lim.* **42**, 65–72 (2006).

31. Branyikova, I. et al. Microalgae-novel highly efficient starch producers. *Biotechnology and Bioengineering* **108**, 766–776, doi: 10.1002/bit.23016 (2011).

32. Yin, Y. H. et al. The influence of light intensity and photoperiod on duckweed biomass and starch accumulation for bioethanol production. *Bioresour. Technol.* **187**, 84–90, doi: 10.1016/j.biortech.2015.03.097 (2015).
33. Ho, S. H., Chen, C. Y. & Chang, J. S. Effect of light intensity and nitrogen starvation on CO2 fixation and lipid/carbohydrate production of an indigenous microalga *Scenedesmus obliquus* CNW-N. *Bioresour. Technol.* **113**, 244–252, doi: 10.1016/j.biortech.2011.11.133 (2012).

34. Mettler, T. et al. Systems analysis of the response of photosynthesis, metabolism, and growth to an increase in irradiance in the photosynthetic model organism *Chlamydomonas reinhardtii*. *The Plant Cell* **26**, 2310–2350, doi: 10.1105/tpc.114.124537 (2014).

35. Peers, G. et al. An ancient light-harvesting protein is critical for the regulation of algal photosynthesis. *Nature* **462**, 518–521, doi: 10.1038/nature08587 (2009).

36. Melis, A. In *Photostasis and Related Phenomena* (eds Williams, Theodore P. & Thistle, Anne B.) 207–221 (Springer US, 1998).

37. Mussgnug, J. H. et al. Engineering photosynthetic light capture: impacts on improved solar energy to biomass conversion. *Plant Biotechnology Journal* **5**, 802–814, doi: 10.1111/j.1467-7652.2007.00285.x (2007).

38. Dent, R. M., Haglund, C. M., Chin, B. L., Kobayashi, M. C. & Niyogi, K. K. Functional genomics of eukaryotic photosynthesis using insertional mutagenesis of *Chlamydomonas reinhardtii*. *Plant Physiology* **137**, 545–556, doi: 10.1104/pp.104.055244 (2005).

39. González-Ballester, D., de Montaigu, A., Higuera, J. J., Galván, A. & Fernández, E. Functional genomics of the regulation of the nitrate assimilation pathway in *Chlamydomonas*. *Plant Physiology* **137**, 522–533, doi: 10.1104/pp.104.050914 (2005).

40. Blankenberg, D. et al. Manipulation of FASTQ data with Galaxy. *Bioinformatics* **26**, 1783–1785, doi: 10.1093/bioinformatics/btq281 (2010).

41. Merchant, S. S. et al. The Chlamydomonas genome reveals the evolution of key animal and plant functions. *Science* **318**, 245–250, doi: 10.1126/science.1143609 (2007).

42. Li, X. et al. An indexed, mapped mutant library enables reverse genetics studies of biological processes in *Chlamydomonas reinhardtii*. *The Plant Cell* **28**, 367–387, doi: 10.1105/tpc.15.00465 (2016).

43. Moroney, J. V. & Ynalvez, R. A. Proposed carbon dioxide concentrating mechanism in *Chlamydomonas reinhardtii*. *Eukaryot. Cell* **6**, 1251–1259, doi: 10.1128/ec.00064-07 (2007).

44. Sato, R., Ito, H. & Tanaka, A. Chlorophyll *b* degradation by chlorophyll *b* reductase under high-light conditions. *Photosynth Res* **126**, 249–259, doi: 10.1007/s11120-015-0145-6 (2015).

45. Kitajima, K. & Hogan, K. P. Increases of chlorophyll a/b ratios during acclimation of tropical woody seedlings to nitrogen limitation and high light. *Plant, Cell & Environment* **26**, 857–865, doi: 10.1046/j.1365-3040.2003.01017.x (2003).

46. Shimogawara, K., Fujiwara, S., Grossman, A. & Usuda, H. High-efficiency transformation of *Chlamydomonas reinhardtii* by electroporation. *Genetics* **148**, 1821–1828 (1998).

47. Greenspan, P., Mayer, E. & Fowler, S. Nile red - a selective fluorescent stain for intracellular lipid droplets. *J Cell Biol* **100**, 965–973 (1985).

48. Cooper, M. S., Hardin, W. R., Petersen, T. W. & Cattolico, R. A. Visualizing “green oil” in live algal cells. *Journal of Bioscience and Bioengineering* **109**, 198–201, doi: 10.1016/j.jbiosc.2009.08.004 (2010).

49. Bligh, E. & Dyer, W. A rapid method of total lipid extraction and purification. *Can J Biochem Physiol* **37**, 911–917 (1959).

50. Lichtenthaler, H. Chlorophylls and carotenoids - pigments of photosynthetic biomembranes. *Method Enzymol* **148**, 350–382 (1987).

51. Jiang, X. & Stern, D. Mating and tetrad separation of *Chlamydomonas reinhardtii* for genetic analysis. *e1274*, doi: 10.3791/1274 (2009).

52. Green, M. R. & Sambrook, J. *Molecular cloning a laboratory manual*. 4th edition Vol. 1 (eds Green, M. R. & Sambrook, J.) Ch. 2, 133–140 (Cold Spring Harbor Laboratory press, 2012).

53. Giardine, B. et al. Galaxy: A platform for interactive large-scale genome analysis. *Genome Research* **15**, 1451–1455, doi: 10.1101/gr.4086505 (2005).

54. Goecks, J. et al. Web-based visual analysis for high-throughput genomics. *BMC Genomics* **14**, 397 (2013).

## Acknowledgements

Work in the authors' laboratory is supported by the French Agence Nationale pour la Recherche (ANR-Diesalg) and by the project Oseo "Exploitation Industrielle des Micro-Algues". Hugh Douglas Goold acknowledges The University of Sydney for an Australian Postgraduate Award, and the Paul Finlay scholarship. Hoa Mai Nguyen thanks the program "Irtelis, CEA' for a PhD studentship. We thank Dr. Xenie Johnson for helpful discussions and Isabelle Te for technical assistance.

## Author Contributions

H.D.G., H.M.N., F.K., A.B.-A., B.L., E.B. and S.C. performed research. H.D.G., H.M.N., F.B., S.C., G.P. and Y.L.-B. designed research. All authors analyzed data and wrote the paper.

## Additional Information

Supplementary information accompanies this paper at http://www.nature.com/srep

Competing financial interests: The authors declare no competing financial interests.

**How to cite this article:** Goold, H. D. *et al.* Whole Genome Re-Sequencing Identifies a Quantitative Trait Locus Repressing Carbon Reserve Accumulation during Optimal Growth in *Chlamydomonas reinhardtii*. *Sci. Rep.* **6**, 25209; doi: 10.1038/srep25209 (2016).

<figure><img src="image_12.png" /></figure>

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/