Plant Physiology Preview. Published on October 8, 2015, as DOI:10.1104/pp.15.00861 1 Running head: Chlamydomonas life cycle transcriptome and methylome 2 3 Corresponding author: 4 Matteo Pellegrini 5 Department of Molecular, Cell, and Developmental Biology 6 University of California, Los Angeles 7 Los Angeles, California, 90095, USA 8 Tel: (310) 825-0012 9 [email protected] 10 11 Research area: Genes, Development, and Evolution 12 Secondary research area: Cell Biology 13 1 Downloaded from on June 17, 2017 - Published by www.plantphysiol.org Copyright © 2015 American Society of Plant Biologists. All rights reserved. Copyright 2015 by the American Society of Plant Biologists 14 Dynamic changes in the transcriptome and methylome of Chlamydomonas 15 reinhardtii throughout its life cycle 16 17 David Lopez1,2, Takashi Hamaji3, Janette Kropat4, Peter De Hoff5,6, Marco Morselli2, 18 Liudmilla Rubbi2, Sorel Fitz-Gibbon2,4, Sean D. Gallaher4, Sabeeha S. Merchant4,7, James 19 Umen3, Matteo Pellegrini2,7 20 21 1 22 23 Molecular Biology Institute, University of California, Los Angeles, Los Angeles, California 2 24 Department of Molecular, Cell, and Developmental Biology, University of California, Los Angeles, Los Angeles, California 25 3 Donald Danforth Plant Science Center, St. Louis, Missouri 26 4 Department of Chemistry and Biochemistry, University of California, Los Angeles, Los 27 Angeles, California 28 5 Salk Institute for Biological Studies, La Jolla, California 29 6 Current address: Synthetic Genomics, La Jolla, California 30 7 Institute for Genomics and Proteomics, University of California, Los Angeles, Los 31 Angeles, California 32 33 One-sentence summary: 34 Secretory and algal-specific genes are induced during gamete and zygote development in 35 C. reinhardii, concurrent with a dramatic increase in chloroplast cytosine methylation. 36 37 2 Downloaded from on June 17, 2017 - Published by www.plantphysiol.org Copyright © 2015 American Society of Plant Biologists. All rights reserved. 38 Funding was provided by the National Institutes of Health R24 GM092473, R01 39 GM078376, and by the Office of Science (Biological and Environmental Research), U.S. 40 Department of Energy (grant no. DE–FC02–02ER63421). David Lopez is supported by 41 an NIH T32 Training Fellowship in Genome Analysis (5T32HG002536-13), the Eugene 42 V. Cota-Robles Fellowship, and the Fred Eiserling and Judith Lengyel Doctoral 43 Fellowship. Takashi Hamaji is a JSPS Postdoctral Fellow for Research Abroad (26-495) 44 of the Japan Society for the Promotion of Science. 45 46 47 Author Contributions 48 MP, JU, and SM designed the experiment. JK, MM, LR, and PDH collected the samples 49 and created the sequence libraries in preparation for sequencing. DL and TH aligned the 50 sequence data and performed all bioinformatic analyses, with the exception of the 51 analysis of the genomic data and identification of structural variants, which was carried 52 out by SFG and SG. The manuscript was written by DL, JU, MP, and SM. 53 54 55 3 Downloaded from on June 17, 2017 - Published by www.plantphysiol.org Copyright © 2015 American Society of Plant Biologists. All rights reserved. 56 Abstract 57 The green alga Chlamydomonas reinhardtii undergoes gametogenesis and mating upon 58 nitrogen starvation. While the steps involved in its sexual reproductive cycle have been 59 extensively characterized, the genome-wide transcriptional and epigenetic changes 60 underlying different life cycle stages have yet to be fully described. Here, we performed 61 transcriptome and methylome sequencing to quantify expression and DNA methylation 62 from vegetative and gametic cells of each mating-type and from zygotes. We identified 63 361 gametic genes with mating-type-specific expression patterns, and 627 genes that are 64 specifically induced in zygotes; furthermore, these sex-related gene sets were enriched 65 for secretory pathway and algal-specific genes. We also examined the C. reinhardtii 66 nuclear methylation map with base-level resolution at different life cycle stages. Despite 67 having low global levels of nuclear methylation, we detected 23 hypermethylated loci in 68 gene-poor, repeat-rich regions. We observed mating-type-specific differences in 69 chloroplast DNA methylation levels in plus versus minus mating type gametes followed 70 by chloroplast DNA hypermethylation in zygotes. Lastly, we examined expression of 71 candidate DNA methyltransferases and found three – DMT1a, DMT1b and DMT4 – that 72 are differentially expressed during the life cycle and are candidate DNA methylases. The 73 expression and methylation data we present provide insight into cell-type-specific 74 transcriptional and epigenetic programs during key stages of the C. reinhardtii life cycle. 75 76 4 Downloaded from on June 17, 2017 - Published by www.plantphysiol.org Copyright © 2015 American Society of Plant Biologists. All rights reserved. 77 Introduction 78 Chlamydomonas reinhardtii is a unicellular, biflagellate species of green alga found 79 primarily in fresh water and soil (Harris et al., 2009). C. reinhardtii is an important 80 reference organism for diverse eukaryotic cellular and metabolic processes including 81 photosynthetic biology (Rochaix, 2001), flagellar function and biogenesis (Silflow and 82 Lefebvre, 2001), nutrient homeostasis (Grossman, 2000; Merchant et al., 2006; Glaesener 83 et al., 2013) and sexual cycles (Goodenough et al., 2007). The nuclear and chloroplast 84 genomes of C. reinhardtii have been fully sequenced, enabling genomic and epigenomic 85 analyses (Maul et al., 2002; Merchant et al., 2007). The ~112 Mbp haploid C. reinhardtii 86 nuclear genome comprises 17 chromosomes. The circular chloroplast genome (cpDNA) 87 is 203 kbp and present in 80-100 copies per cell that are organized into 8-10 88 nucleoprotein complexes called nucleoids, which are distributed through the stroma. 89 90 Like many unicellular eukaryotes, C. reinhardtii has a biphasic life cycle where haploid 91 cells can reproduce vegetatively by mitotic division or, alternatively, undergo a sexual 92 cycle. Vegetative cells can propagate indefinitely when provided with nutrients and light. 93 Upon nitrogen starvation, however, cells stop dividing and differentiate into gametes 94 whose mating type (plus or minus) is determined genetically by a ~300 kb mating-type 95 locus (MT) on chromosome 6, with two haplotypes, MT+ and MT− (Umen, 2011; De 96 Hoff et al., 2013). Gametes express a set of mating-related proteins that are different 97 between minus and plus cells and which allow cells of opposite mating-type to recognize 98 each other and fuse to form a quadriflagellate zygote. Upon fertilization, the 99 heterodimeric KNOX/BELL-type homeodomain proteins GSM1 and GSP1 initiate a 100 zygote-specific developmental program which includes flagellar resorption, fusion of 101 organelles including nuclei and chloroplasts, destruction of MT− cpDNA, and secretion 102 of a thick, environmentally-resistant cell wall that protects the zygospore from cold, 103 dessication and other environmental stresses (Cavalier-Smith, 1976; Catt, 1979; Grief et 104 al., 1987; Brawley and Johnson, 1992; Goodenough et al., 2007; Lee et al., 2008). Upon 105 return to favorable conditions of light and nutrients, zygospores undergo meiosis to 106 produce four haploid progeny (2 MT+, 2 MT−) that can reenter the vegetative life cycle. 107 While nuclear loci segregate in a Mendelian pattern of 2:2, both chloroplast and 5 Downloaded from on June 17, 2017 - Published by www.plantphysiol.org Copyright © 2015 American Society of Plant Biologists. All rights reserved. 108 mitochondrial genomes are inherited uniparentally, with cpDNA inherited from the MT+ 109 parent and mitochondrial DNA (mtDNA) from the MT− parent (Nakamura, 2010; 110 Nishimura, 2010). 111 112 While previous high-throughput expression studies have focused on the transcriptional 113 programs underlying processes such as nutrient deprivation (Nguyen et al., 2008; 114 Gonzalez-Ballester et al., 2010; Toepel et al., 2011; Toepel et al., 2013; Schmollinger et 115 al., 2014), environmental responses (Simon et al., 2008; Matsuo et al., 2011; Fang et al., 116 2012; Simon et al., 2013), flagellar biogenesis (Albee et al., 2013), lipid accumulation 117 (Miller et al., 2010; Boyle et al., 2012; Lv et al., 2013), and diurnal rhythms (Idoine et al., 118 2014; Panchy et al., 2014), only a few studies have explored the genome-wide 119 transcriptional and epigenetic changes associated with the sexual cycle (Kubo et al., 120 2008; Ning et al., 2013; Aoyama et al., 2014). Several genes expressed in the early 121 zygote, termed EZY genes, have predicted functions related to cell wall production, 122 vesicular transport, and secretion (Ferris and Goodenough, 1987; Ferris et al., 2002; 123 Kubo et al., 2008). A separate analysis of zygospore transcripts following light-induced 124 germination revealed the upregulation of photosynthetic and methionine synthesis 125 pathways (Aoyama et al., 2014). 126 127 DNA methylation studies have also been conducted on both the nuclear and chloroplast 128 genomes (Hattman et al., 1978; Dyer, 1982). cpDNA methylation has been studied more 129 extensively and shows dramatic changes in 5meC content at different stages of the C. 130 reinhardtii life cycle. Vegetative cells have low levels of 5meC in cpDNA, while 131 gametes show a substantial increase with cpDNA (12% 5meC in MT+ gamete cells and 132 4% in MT–) (Royer and Sager, 1979; Feng and Chiang, 1984). In zygotes, MT− cpDNA 133 is eliminated while MT+ cpDNA becomes hypermethylated. While differential cpDNA 134 methylation was once thought to be part of a restriction-methylation system regulating 135 uniparental inheritance (Burton et al., 1979), this model is unlikely since loss of cpDNA 136 methylation in MT+ cells does not result in its destruction in zygotes (Umen and 137 Goodenough, 2001), and ectopic methylation of MT− cpDNA does not spare it from 138 destruction (Bolen et al., 1982). However, previous studies are consistent with a role for 6 Downloaded from on June 17, 2017 - Published by www.plantphysiol.org Copyright © 2015 American Society of Plant Biologists. All rights reserved. 139 5meC in promoting cpDNA replication upon zygote germination, which can influence the 140 amount of residual MT− cpDNA that is inherited by exceptional progeny (Umen and 141 Goodenough, 2001; Nishiyama et al., 2004). An alternative proposed mechanism 142 involves digestion of MT− cpDNA by differentially-localized or -activated nucleases that 143 are methylation-insensitive early in zygote development before chloroplast fusion 144 (Nishimura et al., 2002). 145 146 Several methyltransferase enzymes that modify chloroplast DNA have been investigated 147 biochemically (Sano et al., 1981) and one candidate chloroplast methyltransfersase gene 148 has been cloned (Nishiyama et al., 2002; Nishiyama et al., 2004). Since this time, the 149 genome sequence of Chlamydomonas reinhardtii has become available (Merchant et al., 150 2007) and extensively annotated (Blaby et al., 2014) so that a comprehensive 151 identification of genes encoding DNA methyltransferases can be undertaken. 152 153 Few studies have focused on the role of nuclear cytosine methylation in Chlamydomonas, 154 but previous work has shown that induced silencing of nuclear transgenes does not 155 correlate with transgene cytosine methylation levels, leaving open the question of what 156 role cytosine methylation plays in chromatin structure and gene expression in C. 157 reinhardtii (Cerutti et al., 1997). To date, the absence of detailed methylation maps has 158 precluded a clear view of methylation patterns in the nuclear and chloroplast genomes of 159 Chlamydomonas. 160 161 Here we have performed RNA-seq based transcriptome analysis and bisulfite-DNA 162 sequencing of C. reinhardtii at different life cycle stages. We identify sex- and mating- 163 related changes in gene expression including genes that are preferentially expressed in 164 gametes of each mating type or in zygotes. We generated a high-resolution map of 165 nuclear and chloroplast cytosine methylation during the life cycle, and identified 166 candidate DNA methyltransferases whose expression profiles correlate with dynamic 167 changes in cpDNA methylation patterns. 168 169 7 Downloaded from on June 17, 2017 - Published by www.plantphysiol.org Copyright © 2015 American Society of Plant Biologists. All rights reserved. 170 171 Results 172 To quantify nuclear gene expression and 5-cytosine DNA methylation throughout the C. 173 reinhardtii life cycle, we collected RNA and DNA samples at various stages for RNA 174 and whole genome bisulfite (WGBS) sequencing (Figure 1). Samples were collected 175 from vegetative and gametic cells of both mating types and zygotic stages to enable 176 multiple comparisons. The design of the experiment was for matched DNA and RNA 177 samples, but the RNA protocol was modified to avoid degradation. We note that for 178 transcriptome studies, our vegetative samples were prepared in a different manner than 179 those used for typical N-starvation studies. We grew all cultures to saturation on solid 180 agar media and then resuspended cells at high density (~2 x 107 mL-1) under illumination 181 in liquid medium (HSM) with (+N) or without (−N) nitrogen for ~3-5 hours. Under these 182 conditions neither culture grew measurably, but the –N cultures expressed the gametic 183 program and efficiently mated at >90% efficiency; while the +N cultures could not mate 184 and did not express gametic marker genes (see below). The advantage of this procedure 185 is that it minimizes the differences between +N and –N cultures that would normally be 186 attributable to growth rate differences and thereby allows more reliable identification of 187 mating-related gene expression. For DNA methylation studies, the vegetative and 188 gametic samples were obtained from growing N-replete and N-starved samples 189 respectively as described in Materials and Methods. Actively growing cultures were used 190 for studies of methylation in vegetative cells since agar plate grown cells are already 191 partially gametic and would require additional rounds of division in order to remove 192 preexisting methylation. 193 194 Sequence polymorphisms between R3 (MT+) and CJU10 (MT–) parental strains 195 As a prelude to transcriptome and methylome analyses, we cataloged the genetic 196 differences (single nucleotide variants, insertions, and deletions) between our two 197 parental strains using genome resequencing (Figure 2A) (Gallaher et al., 2015). In total, 198 the two strains differ by 0.16% in their nuclear genomes, and most of these differences 199 are single nucleotide variants (SNVs). Of all the variants, 98.2% are localized to two 200 regions, one of length 2.2 Mb on chromosome 17, and one of 2 Mb encompassing the 8 Downloaded from on June 17, 2017 - Published by www.plantphysiol.org Copyright © 2015 American Society of Plant Biologists. All rights reserved. 201 mating locus on chromosome 6 (Figure 2B) where mating-type haplotype differences 202 have previously been observed (De Hoff et al., 2013). Additionally, the chloroplast 203 genome of the CJU10 strain contains an insertion adjacent to the atpB gene of a 204 spectinomycin resistance marker (aadA), flanked by the rbcL 5’ promoter and psbC 3’ 205 UTR (Figure 2C) (Goldschmidt-Clermont, 1991; Umen and Goodenough, 2001). 206 207 Gamete and zygote-specific genes 208 Following the quantitation of RNA-seq data for all of our samples, we identified genes 209 with mating type and zygote-specific expression patterns using a series of filters to screen 210 for genes matching expression patterns of known gametic and zygotic genes. We 211 required that mating-type-specific genes be expressed in gametes at least four-fold higher 212 than in vegetative cells of the same mating type, and at least ten fold higher than in 213 gametes or vegetative cells of the opposite mating type. For zygote-specific genes, we 214 required that expression be at least four-fold higher in early zygotes than in any other 215 sample. Using these criteria, we identified 293 and 68 genes whose expression is specific 216 to plus and to minus gametes respectively, and 627 genes whose expression is specific to 217 zygotes (Figure 3A, Supplementary Table 1). Genes whose expression is known to be 218 gamete- or zygote-specific, including GSM1 (minus gametes), SAG1 (plus gametes) and 219 EZY1 (early zygotes) were retained by our filters (Armbrust et al., 1993; Kurvari et al., 220 1998; Ferris et al., 2002; Ferris et al., 2005; Lee et al., 2008). Of the 32 previously 221 described zygotic genes (Matters and Goodenough, 1992; Armbrust et al., 1993; Uchida 222 et al., 1993; Kuriyama et al., 1999; Uchida et al., 1999; Suzuki et al., 2000; Ferris et al., 223 2002; Kubo et al., 2008), 25 were in our zygote data set, while the remaining 7 (EZY6, 224 EZY15, EZY16, EZY17, EZY21, EZY23, and LYZ1A) were found to have significant 225 expression in other samples and were therefore excluded from the zygote-specific set 226 (Supplementary Tables 1 and 2). 227 228 We functionally classified the predicted proteins encoded by genes whose expression was 229 plus-gamete-specific, minus-gamete-specific, or early-zygote-specific. Protein 230 localization prediction revealed the enrichment of putative secretory proteins in plus- 231 gamete (MT+) and zygote-specific sets (Figure 3B). We also assessed conservation and 9 Downloaded from on June 17, 2017 - Published by www.plantphysiol.org Copyright © 2015 American Society of Plant Biologists. All rights reserved. 232 phylogenetic distribution for homologs of each protein in a set of nested phylogenetic 233 domains that encompass different taxonomic levels from cellular organisms (prokarya, 234 archaea and eukarya) to the single species level of C. reinhardtii (Figure 3C). Zygote 235 and gamete expression groups were enriched for Chlamydomonas-specific genes and/or 236 algal specific genes. In addition, gametic genes were under-represented for more widely 237 conserved genes (i.e. those with homologs outside of chlorophyte algae). Consistent with 238 these findings, the gametic and zygotic gene lists were also depleted to various extents 239 for GreenCut2 and CiliaCut genes whose members are associated with conserved 240 photosynthetic and flagella/basal body functions (Merchant et al., 2007; Karpowicz et al., 241 2011; Heinnickel and Grossman, 2013) (Figures 4A and 4B); however, the overall low 242 numbers of genes in these categories precluded obtaining a significant statistical result in 243 all cases but one. Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes 244 (KEGG), and MapMan (Kanehisa and Goto, 2000; Harris et al., 2004; Thimm et al., 245 2004) classification of predicted proteins encoded by gametic and zygotic genes was 246 performed using the Algal Functional Annotation Tool (Lopez et al., 2011), but had 247 limited utility because of the large number of non-conserved proteins in these groups and 248 incomplete annotations. Nonetheless, we found significant enrichment in two MapMan 249 categories for zygotic genes, “Cell Wall” and “Transport”, both of which may relate to 250 the requirement for new cell wall biosynthesis in zygotes (Figures 4C and 4D). Indeed, 251 examination of manually curated early zygotic gene annotations revealed numerous cell- 252 wall related protein coding genes as described below. 253 254 Volvocine cell walls are composed primarily of glycosylated hydroxyproline-rich 255 glycoproteins (HRGPs) that enter the secretory pathway and are exported to the 256 extracellular space where they co-assemble (Woessner et al., 1994). The thick and 257 environmentally-resistant cell walls of zygospores are formed by a specialized set of 258 HRGPs that are synthesized shortly after fertilization (Minami and Goodenough, 1978; 259 Catt, 1979; Grief et al., 1987). Manual annotation and inspection of zygotic up-regulated 260 genes verified the MapMan ontology assignments of cell wall and transport categories as 261 described in Supplementary Table 1. At least 57 zygotic genes are predicted HRGPs or 262 have putative cell wall biogenesis-related functions that include secretion, glycosylation 10 Downloaded from on June 17, 2017 - Published by www.plantphysiol.org Copyright © 2015 American Society of Plant Biologists. All rights reserved. 263 and metabolism of nucleotide sugars (e.g. UDP-glucose 4-epimerase, pyrophosphorylase, 264 dehydrogenase; dTDP-6-deozy-L-lyxo-4-hexulose reductase related, exotosin-like 265 glycosyltransferase), or sugar metabolite transport (e.g. ABC transporter, triose 266 phosphate transporter, UDP-N-acetylglucosamine transporter). Among these were some 267 previously identified early zygotic genes as noted in Supplementary Table 1 (EZY4, 268 EZY11/UPT1/UPTG1, EZY12/UGD1, EZY14/TPT14, EZY22) (Kubo et al., 2008). 269 270 Besides cell wall and secretory pathway genes, we also noted in the zygote gene set 271 predicted functions that may be related to other zygotic processes including elimination 272 of MT– chloroplast DNA, zygotic cpDNA methylation and packaging for long-term 273 dormancy, nuclear fusion (karyogamy), chloroplast fusion, and flagellar resorption 274 (Goodenough et al., 2007). These annotations include predicted chloroplast-targeted 275 DNA binding proteins such as a recA homolog and predicted nuclease 276 (EZY19/Cre07.g314650), chloroplast targeted DNA methyltransferases (DMT1A, 277 DMT1B, DMT4; discussed below), and a chloroplast-targeted dynamin 278 (EZY8/Cre06.g25065) that may be involved in chloroplast fusion (Kubo et al., 2008). 279 Predicted nuclear-targeted zygotic proteins include several DNA binding transcription 280 factors (EZY18/Cre02.g091550, RLS7/Cre14.g617200, ZYS1A/Cre17.g719200) such as 281 RLS7 that contains a SAND-domain (Duncan et al., 2006; Duncan et al., 2007) related to 282 RegA, a repressor of germ cell fate in Volvox carteri (Kirk et al., 1999), and several types 283 of chromatin-related proteins (Cre03.g184900, Cre08.g367000, Cre08.g400200, 284 Cre09.g401812, HON1/Cre13.g567450) that may be involved in nuclear DNA packaging 285 in preparation for zygospore dormancy. Minutes after fertilization and prior to flagellar 286 resorption the four basal bodies and flagella of newly formed zygotes (two from each 287 parent) move to a single apical location via an unknown mechanism. One possible 288 participant in this process could be the striated fiber protein SF-assemblin (SFA1/ 289 Cre07.g332950), a zygote up-regulated gene whose protein product in vegetative cells 290 associates with rootlet microtubules that are proximal to basal bodies and are thought to 291 play a role in organization of the rootlet structure (Lechtreck et al., 2002). Two other 292 cytoskeletal proteins whose genes are up-regulated in zygotes are a flagella-associated 293 protein of unknown function, FAP79 (Cre04.g217908) and the flagella length regulatory 11 Downloaded from on June 17, 2017 - Published by www.plantphysiol.org Copyright © 2015 American Society of Plant Biologists. All rights reserved. 294 protein LF5/FAP279 (Tam et al., 2013) (Cre12.g538300) which may be related to 295 flagellar resorption that begins two or three hours after fertilization. Lastly, we identified 296 a gene for a predicted secreted trypsin-related protease (Cre06.g287750), which could 297 contribute to rapid post-fertilization degradation of gametic plasma membrane surface 298 proteins such as FUS1 and GCS1/HAP2, a process that is thought to restrict polygamy 299 (fusion between >2 gametes) (Liu et al., 2010). 300 301 DNA methylation of the nuclear genome 302 We conducted bisulfite sequencing of vegetative, gametic and zygotic samples to 303 generate DNA methylation profiles. The nuclear genome had an average per-site CG 304 methylation of < 0.75% in all samples, and this level of methylation did not differ 305 significantly between plus and minus strains or at different life cycle stages (Figure 5A). 306 However, CG methylation densities greater than 80% were identified for 23 loci that 307 ranged in size between 10 kb and 22 kb (Figure 5B, Supplementary Table 3). The highly- 308 methylated regions are enriched for repeats, and their overall protein coding gene 309 densities are significantly lower than average, though there are still genes in these regions 310 (Figure 5C). In addition, one example where methylation was strain-specific is shown in 311 Figure 5D, although most hypermethylated sites did not have strain specific methylation 312 patterns. The expression levels of genes overlapping hypermethylated loci are not 313 strongly correlated with degree of methylation (Supplementary Table 4). 314 315 Chloroplast methylation changes during the life cycle 316 In contrast to the relatively stable pattern of cytosine methylation in nuclear DNA, the C. 317 reinhardtii chloroplast genome underwent dynamic changes in cytosine methylation 318 throughout the life cycle (Figure 6A). In the vegetative stage, global per-cytosine 319 methylation was < 2% for both mating types for all cytosine contexts (CG, CHG, CHH). 320 After gametogenesis, cpDNA methylation increased in a mating-type dependent manner. 321 MT+ gametes had an average of ~10% per-site CG methylation while MT− gametes had 322 an average of ~3%. A large increase in 5meC was observed for all sequence contexts 323 during zygote development with 54% (CG) methylation reached by 24 hours. 324 Methylation levels began to drop during germination to an average per-site CG 12 Downloaded from on June 17, 2017 - Published by www.plantphysiol.org Copyright © 2015 American Society of Plant Biologists. All rights reserved. 325 methylation level of 45%. CHG and CHH methylation levels were lower than CG 326 methylation levels in all cases, and this difference was particularly notable in germinated 327 zygotes, where CHH and CHG methylation dropped from their peak zygotic levels much 328 faster than CG methylation did. 329 330 For all life cycle stages and mating types, methylation of cpDNA is uniformly distributed 331 in the chloroplast genome without bias towards genes or other sequence features (Figure 332 6B). With the exception of two large inverted repeats in the chloroplast sequence, where 333 methylation cannot be measured accurately with short reads, average methylation levels 334 tabulated in 1kb intervals rarely deviated more than 5% from the global average, 335 suggesting that no specific regions of the chloroplast genome are targeted for 336 methylation. 337 338 DNA methyltransferases 339 In order to further explore the mechanisms responsible for the dynamic patterns of DNA 340 methylation, we identified candidate cytosine methyltransferases from the predicted C. 341 reinhardtii proteome based on the presence of predicted DNA cytosine methylase (DCM) 342 domains (Materials and Methods). We found a total of six candidate methyltransferases 343 with different domain architectures (Supplementary Figure 1) including a homolog of 344 Volvox carteri MET1, a putative nuclear methyltransferase (Babinger et al., 2007). As 345 described below, three predicted methyltransferases – DMT1a, DMT1b, and DMT4 – had 346 expression patterns and/or predicted localization sequences suggesting a role in 347 chloroplast DNA methylation. Each of these predicted proteins has DNA methylase as 348 well as BAH domains. DMT1a, DMT1b, and DMT4 were all expressed in gametes and 349 showed highest levels of expression in zygotes, coinciding with elevated levels of 350 methylation in the zygotic samples (Figure 7). DMT1a and DMT1b are near the mating 351 locus in its telomere-proximal domain and encode highly similar paralogs. Both genes 352 had higher expression in MT+ than in MT− gametes, a pattern which matches the 353 methylation bias seen in gametic cpDNA from the two mating types (Figure 6A). 354 DMT1a and DMT1b sequences have previously been described as a single gene 355 (Nishiyama et al., 2002), but the published cDNA sequence is a hybrid with 5’ sequences 13 Downloaded from on June 17, 2017 - Published by www.plantphysiol.org Copyright © 2015 American Society of Plant Biologists. All rights reserved. 356 derived from DMT1a and 3’ sequences from DMT1b (NCBI Gene IDs #5722229 and 357 #5722231). Consistent with sub-cellular targeting predictions (Figure 7) the DMT1a pre- 358 sequence directs chloroplast localization (Nishiyama et al., 2002). While the MT+ strain 359 copy of DMT1a appears to be intact, a survey of structural variants derived from genome 360 resequencing data led to the identification of a point insertion in exon 10 of the MT− 361 strain copy of DMT1a leading to a frame shift and premature termination before the 362 methyltransferase domains (Figure 8B). This insertion was observed in MT– 363 transcriptome data. Furthermore, this point insertion is found in the genome of 12 out of 364 13 MT– strains resequenced by Gallaher et al. No variants predicted to be deleterious 365 were found within the DMT1b gene. Targeting predictions of DMT1b suggest it is 366 mitochondria-localized. Although the mitochondrial genome is largely devoid of 367 cytosine methylation (1-2% global per-cytosine methylation), the zygote sample at 24 368 hours shows evidence of methylation (~13% global per-cytosine methylation) 369 (Supplementary Figure 2). DMT4 is also predicted to encode a chloroplast-targeted 370 cytosine methyltransferase and is one of the 361 genes identified with a strong zygotic 371 expression pattern (Figure 3A), consistent with a possible participation of DMT4 in 372 zygotic cpDNA hypermethylation. 373 374 Discussion 375 376 The dynamic changes in gamete- and zygote-specific mRNA abundance and DNA 377 methylation presented in this work provide a framework for understanding cell 378 differentiation during the Chlamydomonas sexual life cycle. A previous study of plus 379 and minus gamete specific genes focused on cell type genes whose expression was up- 380 regulated during the process of mating (Ning et al., 2013). However, direct comparison 381 with the previous data is confounded by differences in annotation between genome 382 versions used to define transcripts as well as in the different clustering criteria used in the 383 two studies. Here we made use of culture conditions that were designed to suppress the 384 differential transcript abundance signal resulting from growing versus non-growing N- 385 starved cultures to help identify gamete specific transcripts. Our method identified 386 known gametic genes whose expression is mating-type limited (expressed preferentially 14 Downloaded from on June 17, 2017 - Published by www.plantphysiol.org Copyright © 2015 American Society of Plant Biologists. All rights reserved. 387 in plus versus minus gametes or vice versa). However, because we used nitrogen 388 resupply of stationary phase cultures to create non-growing ‘vegetative’ samples we may 389 have missed highly stable gametic transcripts that did not turn over after nitrogen 390 addition. Nonetheless, nitrogen resupply for several hours was completely effective at 391 suppressing mating, so the transcriptome differences we were able to identify in our 392 gametic versus vegetative samples are likely tied to gametogenesis-related functions and 393 perhaps less relevant for other N-starvation responses such as neutral lipid accumulation. 394 395 Sex-related genes evolve rapidly and are therefore expected to appear younger and have a 396 more restricted phylogenetic distribution than other genes (Swanson and Vacquier, 2002). 397 Indeed, using phylogenomic profiling we found an enrichment of Chlamydomonas- 398 specific genes belonging to the plus or minus gamete up-regulated categories and zygote 399 up-regulated categories (Figure 3C). This finding underscores the importance of species- 400 specific and clade-specific genes as potential drivers of cell-type specialization related to 401 sex and speciation. Although the functions of these genes are difficult to predict since 402 they have no homologs outside of Chlamydomonas or Volvocine algae, we did find 403 enrichment for secretory pathway targeting signals within the plus gamete and zygote 404 predicted proteins (Figure 3B) which could indicate a role for sex-related proteins in the 405 plasma membrane, cell wall or extracellular space. Manual annotation of zygotic genes 406 showed that their predicted functions match processes such as glycosylation and transport 407 that are associated with cell wall formation. In addition we identified or confirmed 408 expression of zygotic genes that may be associated with other poorly-understood 409 differentiation processes, including chloroplast fusion, MT– cpDNA elimination, DNA 410 methylation, and cytoskeletal remodeling (Supplementary Table 1). Deeper investigation 411 of these genes may yield insights into the cell biology of zygote differentiation that may 412 have parallels in other zygospore-forming algae (Brawley and Johnson, 1992) and even 413 in plants where pollen cells must undergo a similar process of dormancy and reactivation 414 when exposed to appropriate conditions (Brown and Lemmon, 2011). Early zygotes in 415 Chlamydomonas also undergo dramatic changes in their cytoskeleton. Unlike the case in 416 animals where paternal but not maternal centrioles are contributed to zygotes (Avidor- 417 Reiss et al., 2015), in Chlamydomonas both parental basal bodies (structurally similar to 15 Downloaded from on June 17, 2017 - Published by www.plantphysiol.org Copyright © 2015 American Society of Plant Biologists. All rights reserved. 418 centrioles) and flagella are retained initially in zygotes, but eventually are disassembled 419 and rebuilt during germination (Cavalier-Smith, 1976). Nonetheless, investigation of 420 conserved zygotically up-regulated cytoskeletal proteins such as LF5 and SF-assemblin 421 may shed light on more general mechanisms involved in controlling the dynamic 422 behavior of flagella/cilia and basal bodies during cellular remodeling and life cycle 423 transitions. 424 425 DNA methylation changes during the life cycle 426 In contrast to many organisms where a large fraction of the nuclear genome is methylated 427 (Smith and Meissner, 2013; Bestor et al., 2014), we found relatively low levels of 428 cytosine methylation in the nuclear genome of Chlamydomonas, consistent with previous 429 surveys that were not as high resolution as reported here (Hattman et al., 1978; Feng et 430 al., 2010). Interestingly, we did identify 23 hypermethylated loci larger than 10kb 431 (Supplementary Table 3) that tended to occur within gene-poor, repeat-rich regions of the 432 genome. However, the mechanism that leads to the methylation of these loci remains 433 unknown. Previous studies on silencing of transgenes in Chlamydomonas found that 434 inserted transgene repeats were frequently methylated, but reported little correlation 435 between methylation and gene expression (Cerutti et al., 1997). On the other hand, in the 436 related colonial alga Volvox carteri where methylcytosine frequency is slightly higher 437 than in C. reinhardtii (1.1% vs. ~0.75%), nuclear methylation did appear to be associated 438 with gene silencing (Babinger et al., 2001; Babinger et al., 2007). Interestingly, V. 439 carteri also has a more repeat-rich genome than C. reinhardtii with many active 440 transposons (Miller et al., 1993; Ueki and Nishii, 2008; Prochnik et al., 2010), and may 441 have retained or evolved higher levels of nuclear methylation activity than C. reinhardtii 442 to suppress transposon activity. 443 444 The low frequency of cytosine methylation in the nuclear genome contrasts with the 445 dynamic and abundant cytosine methylation in the chloroplast genome. The overall 446 patterns we observed are in agreement with previous findings where vegetative cpDNA 447 from both mating types had low levels of cytosine methylation, and that gametes showed 448 elevated levels with MT+ cpDNA having a two or three fold higher frequency of 16 Downloaded from on June 17, 2017 - Published by www.plantphysiol.org Copyright © 2015 American Society of Plant Biologists. All rights reserved. 449 methylcytosine compared to MT– cpDNA (Royer and Sager, 1979; Dyer, 1982). We 450 observed hypermethylation in zygotes, which has also been seen previously. Our study 451 extended these earlier results by examining genome-wide methylation in cpDNA at 452 single-base resolution. Unlike methylation of nuclear DNA, which was mostly restricted 453 to a few loci, cpDNA cytosine methylation was uniformly distributed across all regions 454 and sequence feature types (genic, intergenic, repeats, exons, introns, etc.). This finding 455 has implications for the function of cpDNA methylation and the enzymes that are 456 responsible for methylating cpDNA (elaborated in next section). Previous models of 457 cpDNA inheritance invoked a methylation-restriction system similar to that in 458 prokaryotes where sequence-specific methyltransferases protect MT+ cpDNA from site- 459 specific restriction endonucleases. However, the non-sequence-specific distribution of 460 methylcytosines we observed and the modest differences between levels of MT+ and 461 MT– cpDNA methylation in gametes are not consistent with a methylation-restriction 462 mechanism. Moreover, if cpDNA methylation were related to methylation and 463 restriction, it would be most effective if it were established at the gametic stage of the life 464 cycle. Instead, the majority of cpDNA methylation occurs in zygotes and is pervasive 465 genome-wide. The purpose of this massive methylation is unknown, but is likely tied to 466 the packaging and protection of cpDNA in zygospores where it may be dormant for many 467 years before germination (Brawley and Johnson, 1992). 468 469 DNA methyltransferases 470 Our survey of candidate methyltransferases revealed two candidates with expression 471 patterns and predicted chloroplast targeting sequences that make them likely responsible 472 for cpDNA methylation—DMT1a and DMT4. Both genes have detectable expression in 473 gametes and zygotes (Figure 7). DMT1a and DMT1b are physically and genetically 474 linked to the mating locus and furthermore, the MT– linked copy of DMT1a has a point 475 mutation that is predicted to generate a truncated, non-functional protein. Although they 476 are linked to the mating locus, the DMT1a/b locus could still recombine with MT or be 477 subject to gene conversion, possibly leading to the creation of strains where levels of 478 gametic cpDNA in MT+ and MT– parents is equivalent. Such strains, if they could be 479 isolated, would be useful for testing the significance of differential gametic cpDNA 17 Downloaded from on June 17, 2017 - Published by www.plantphysiol.org Copyright © 2015 American Society of Plant Biologists. All rights reserved. 480 methylation and the contribution of the DMT1a gene to gametic cpDNA methylation 481 levels. 482 483 Although the predicted DMT1 and DMT4 encoded methyltransferases have accessory 484 BAH domains that are implicated in protein-protein interactions and targeting to specific 485 loci (Callebaut et al., 1999; Yang and Xu, 2013), the non-specific pattern of cpDNA 486 methylation we observed and attribute to the predicted DMT1a and DMT4 proteins 487 suggest that their targeting is non-sequence-specific. The BAH domains in these proteins 488 may serve other functions such as general targeting of the methyltransferase enzymes to 489 chloroplast nucleoids. The lack of sequence specificity predicted for DMT1a and DMT4 490 may also prove useful for biotechnology applications that require non-specific DNA 491 cytosine methylation activity. 492 493 494 18 Downloaded from on June 17, 2017 - Published by www.plantphysiol.org Copyright © 2015 American Society of Plant Biologists. All rights reserved. 495 Materials and Methods 496 Sample generation for bisulfite-treated DNA 497 Chlamydomonas reinhardtii strains R3 (CC-620 R3 NM MT+) and CJU10 (Umen and 498 Goodenough, 2001) were grown in high salt medium (HSM) to concentrations of 3.4 x 499 106 and 3.5 x 106 cells/ml, respectively. 100 ml of each strain culture was used for DNA 500 and RNA isolation as vegetative samples. Cells were collected by centrifugation (3500 501 rpm, 5 min.), resuspended in HSM−N, and after 15 hours, 40 ml of each strain culture 502 was collected as gametic samples. Gametes were mixed and checked for mating 503 efficiency (85% efficient), and after 1 hour, a 40 ml sample was collected, corresponding 504 to the “1 hour zygote” sample. Mixed gametes were split into two flasks: (1) cells in the 505 first were collected after 24 hours, corresponding to the “24 hour zygote sample”; (2) 506 cells in the second were resuspended in water and plated with HSM media, incubated in 507 the light for 24 hours, incubated for 5 days in the dark, and resuspended in Tris-EDTA- 508 NaCl (TEN) buffer with 0.2% NP−40, at which time half of the culture was collected as 6 509 day zygote samples. The remaining half was transferred to TAP medium, incubated in 510 the light for 24 hours, and collected as the germinated zygote sample. 511 512 DNA isolation and purification for bisulfite treatment 513 Cells were resuspended in 4 ml TEN buffer, with the exception of 24 hour and 6 day 514 zygotes. 24 hour and 6 day zygotes were resuspended in 10 ml TEN, mixed for one 515 minute, and repelleted (repeated three times) followed by resuspension in 4 ml TEN. 400 516 μl 20% SDS and 400 μl 20% Sarkosyl were added to each sample. For 24 hour and 6 day 517 zytores, 4 ml zirconium beads was also added followed by vortexing for 5 minutes. 200 518 μl pronase E solution (10 mg/ml) was then added and samples were incubated at 37° C 519 for 30 minutes. 2.5 ml phenol (10 mM Tris-Cl pH 8 saturated) was added, followed by 520 2.5 ml chloroform:isoamylalcohol (24:1). The phases were separated by centrifugation 521 (5000 rpm, 10 min) at 10° C and the upper phase (4.5 ml) was transferred into 9 ml 100% 522 ethanol and incubated overnight at −20° C. Nucleic acids were collected by 523 centrifugation, resuspended in 1 ml 10mM Tris-Cl pH 8 and 100 μl RNase (5 mg/ml) was 524 added. After RNase treatment, samples were extracted again with 1 ml 19 Downloaded from on June 17, 2017 - Published by www.plantphysiol.org Copyright © 2015 American Society of Plant Biologists. All rights reserved. 525 phenol/chloropform/isoamylalcohol (25:24:1) and DNA was precipitated with 0.3M 526 NaAcetate and 70% isopropanol for 30 minutes at room temperature. 527 528 DNA library preparation 529 500 ng of purified genomic DNA was sheared by sonication with Covaris S2 to generate 530 DNA fragments spanning from 100 to 400 bp size range. Library preparation was carried 531 out using NEBNext DNA Library Prep Master Mix (Set for Illumina, NEB, cat #E6040) 532 according to manufacturer instructions with minor modifications. The ligation was 533 performed using Illumina TruSeq Adapters (Illumina, cat # 15025064), DNA size- 534 selection (200-400 bp range) was carried out with AMPure XP beads (Beckman Coulter, 535 Inc.), prior to bisulfite conversion using EZ DNA Methylation-Lightning Kit (Zymo, cat 536 # D5030). The bisulfite treated DNA was amplified using Illumina Primer Cocktail Mix 537 (Illumina, cat #15027084) and MyTaq Mix (Bioline, cat # BIO-25045) according to the 538 following program: 98°C – 2 min; 12 cycles of: 98°C – 15 sec, 60°C – 30 sec, 72°C – 30 539 sec; 72°C – 5 min. 540 541 RNA library preparation 542 Total RNA isolation for RNA-seq analysis was performed as previously described (De 543 Hoff et al., 2013) with additional DNase treatment (2U Roche RNase-free DNase per 110 544 μg total RNA, 37 °C for 20 min.) before final Qiagen RNeasy column purification. 545 546 Genomic variation 547 Genomic reads were aligned using Burrows-Wheeler Aligner v.0.6.2 v.0.6.2 (Li et al., 548 2008), with default parameters, to the version 5.0 assembly of the C. reinhardtii CC-503 549 genome (Merchant et al., 2007). After removing duplicates with Picard MarkDuplicates 550 (http://picard.sourceforge.net), we applied Genome Analyzer Tool Kit (McKenna et al., 551 2010) base quality score recalibration, InDel realignment, and small variant discovery 552 (DePristo et al., 2011). This was followed by hard filtering of variants with extensive 553 manual calibration guided by inspections in IGV (Thorvaldsdottir et al., 2013). 554 555 RNA-sequencing analysis 20 Downloaded from on June 17, 2017 - Published by www.plantphysiol.org Copyright © 2015 American Society of Plant Biologists. All rights reserved. 556 RNA sequencing data was aligned to the C. reinhardtii genome (assembly version 5.5; 557 Phytozome v10.3 gene annotation) with TopHat v2.0.10 using the annotation to guide 558 spliced alignment. Default parameters were kept, with the exception of constraining 559 intron lengths to less than 5kb. Expression levels were quantified using Cufflinks v2.2.1 560 to compute FPKMs. Sequence data has been deposited in NCBI’s Short Read Archive 561 (SRA) under the accessions SRR2048158, SRR2048159, SRR2048160, SRR2048161, 562 SRR2048162, SRR2048163, SRR2048164, and SRR2048165. 563 564 Identification of gamete- and zygote-specific genes 565 Gamete- and zygote-specific genes were identified by applying a series of filters to the 566 FPKM data generated as described above. We defined gamete-specific genes for each 567 mating type as those that have expression in all samples (excluding the zygote) that is 568 less than 10% of the expression level in the gamete of that mating type. Zygote-specific 569 genes were defined as those whose expression in all other samples was less than 25% of 570 the expression level in the zygote. Expression values for samples with replicates were 571 averaged. 572 573 Bisulfite-treated DNA sequencing analysis 574 Raw sequence data was demultiplexed using standard Illumina barcode indices and was 575 checked for quality using FastQC (v0.10.1). Bisulfite-converted sequences were aligned 576 to the Chlamydomonas reinhardtii nuclear genome (Nov. 2011 assembly) and chloroplast 577 genome using the BS-Seeker2 alignment pipeline v2.0.5 (Guo et al., 2013). Default 578 whole-genome bisulfite alignment parameters were chosen with the following 579 exceptions: Bowtie2 was used as the aligner and local alignments were enabled. 580 Methylation levels were called for cytosines covered by at least 4 reads. Sequence data 581 has been deposited in NCBI’s SRA under the accessions SRR2051057, SRR2051058, 582 SRR2051059, SRR2051060, SRR2051061, SRR2051062, SRR2051063, and 583 SRR2051065. 584 585 Identification of candidate methyltransferases 21 Downloaded from on June 17, 2017 - Published by www.plantphysiol.org Copyright © 2015 American Society of Plant Biologists. All rights reserved. 586 To identify candidate DNA methyltransferases, protein sequences of C. reinhardtii 587 (Phytozome v10 gene annotation) were scanned against the Pfam-A database (Release 588 27) using the standalone PfamScan scripts provided by the Wellcome Trust Sanger 589 Institute. The presence of the ‘C-5 cytosine-specific DNA methylase’ domain (PF00145) 590 was used as a criterion for assignment as a cytosine methyltransferase. 591 592 22 Downloaded from on June 17, 2017 - Published by www.plantphysiol.org Copyright © 2015 American Society of Plant Biologists. All rights reserved. 593 Figure Captions 594 595 Figure 1: Chlamydomonas reinhardtii sexual life cycle and sequenced samples 596 597 Vegetative C. reinhardtii cells of each mating type (MT+ and MT−) can be induced to 598 undergo gametogenesis by nitrogen starvation (−N). Gametes of opposite mating type 599 recognize each other through flagellar adhesion and fuse to form a diploid zygote. 600 During zygote maturation MT− chloroplast DNA (cpDNA) is eliminated, flagella are 601 resorbed, and a thick, zygote cell wall forms. Upon return to nitrogen and light 602 zygospores undergo meiosis to form four haploid progeny (two of each mating type, all 603 containing uniparentally inherited MT+ parental cpDNA) that reenter the vegetative 604 cycle. Colored boxes designate samples and material sequenced. Blue and red unfilled 605 circles represent MT+ and MT– chloroplast genomes. Filled pink and light blue circles 606 represent nuclear DNA from MT+ and MT– strains respectively. 607 608 Figure 2: Genetic differences between R3 (MT+) and CJU10 (MT−) parental strains 609 610 (A) Genetic differences, categorized by variant-type (single nucleotide variants [SNVs], 611 insertions, and deletions) are shown with total number of variant loci and total bases. (B) 612 Large-scale view of chromosomes 6 and 17, where 98.5% of the identified variants are 613 located. Locations of sequence differences are shown in red. The region of chromosome 614 6 containing the mating locus is denoted by a blue bar. (C) Schematic representation of 615 the antibiotic resistance transgene insertion in CJU10 chloroplast DNA. The aadA gene, 616 along with the rbcL 5’ promoter and psbC 3’ UTR, is inserted between inverted repeat 617 region “B” (IRB) and the endogenous atpB gene. 618 619 Figure 3: Gamete- and zygote-specific genes 620 621 (A) Box and whisker plots of gene expression profiles for gamete- and zygote-specific 622 genes are plotted. Values are plotted relative to the maximum expression value of each 623 of the genes. 293 and 68 genes were expressed specifically in MT+ and MT– gametes 23 Downloaded from on June 17, 2017 - Published by www.plantphysiol.org Copyright © 2015 American Society of Plant Biologists. All rights reserved. 624 respectively. 627 genes were expressed specifically in zygotes. (B) Localization 625 predictions for proteins encoded by gamete (gam. MT+, gam. MT−) and zygote-specific 626 genes, compared to all predicted proteins in the Chlamydomonas proteome (all). Data are 627 plotted as fraction predicted for each of four compartments with total numbers indicated 628 within each graph portion: secretory in blue; chloroplast (chloro.) in green; mitochondria 629 (mito.) in purple; other in gray. Starred values are significantly different from the total 630 proteome distribution (*: p < 0.05; **: p < 0.01). (C) Predicted protein groups as 631 described in 3B are plotted according to taxonomic distribution of encoded proteins from 632 each category. From bottom (Chlamydomonas-specific proteins, dark green) to top (all 633 cellular organisms, gray) are increasingly broad phylogenetic distributions. Starred 634 samples indicate significant enrichment or depletion in a taxonomic category relative to 635 the distribution of all proteins. 636 637 Figure 4: Functional annotation of proteins encoded by gamete- and zygote-specific 638 genes 639 640 Functional annotations for gamete- and zygote-specific gene lists are shown as bar plots 641 showing the fraction of total proteins in each group with the specified annotation. Total 642 number of genes in each annotation category versus number of genes with the described 643 annotation are shown. Significant differences between all genes and gamete or zygote 644 specific genes were calculated with the hypergeometric test and significant p-values are 645 shown. (A) GreenCut2 genes (Karpowicz et al., 2011; Heinnickel and Grossman, 2013) 646 (B) CiliaCut genes (Merchant et al., 2007). (C) MapMan cell wall related genes (Thimm 647 et al., 2004). (D) MapMan transport related genes. 648 649 Figure 5: Nuclear methylation at different C. reinhardtii life cycle stages 650 651 (A) Bulk averages of 5-methyl cytosine in the nuclear genome are shown for each sample 652 and categorized by the fraction of 5meC in each of three sequence contexts (CG, CHG 653 and CHH). (B) Scaled representations of the 17 C. reinhardtii chromosomes with 23 654 hypermethylated regions larger than 10kb shaded in blue. (C) Plots comparing gene and 24 Downloaded from on June 17, 2017 - Published by www.plantphysiol.org Copyright © 2015 American Society of Plant Biologists. All rights reserved. 655 repeat densities of hypermethylated regions relative to the entire genome. (D) Genome 656 browser display of nuclear cytosine methylation frequency at two representative 657 hypermethylated loci from chromosome 6. Two representative examples of strain- 658 specific hypermethylated loci are shown side-by-side, with a gene track below for 659 genomic context. In the coverage plot, the blue represents unmethylated cytosines 660 (converted by bisulfite treatment), the red represents methylated cytosines (unconverted), 661 and the gray represents the total coverage (including reads from the opposite strand 662 which do not provide methylation information). 663 664 Figure 6: Chloroplast genome methylation at different life cycle stages 665 666 (A) Bulk cytosine methylation frequencies for chloroplast DNA at different life cycle 667 stages are plotted for each methylation context. Mating type of each sample (+, - or 668 diploid +/-) and time of zygote development and germination “(g)” are shown below the 669 graph. (B) Plot of average CG methylation frequency for each sample in 1kb bins across 670 the chloroplast genome. The two inverted repeat regions are shaded in gray. Plots are 671 color coded according to the legend on the right. 672 673 Figure 7: Candidate chloroplast methyltransferases in C. reinhardtii 674 675 The protein domain structure of each candidate methyltransferase is shown schematically 676 (green = DNA methylase domain; red = bromo adjacent homology (BAH) domain). 677 Predicted localizations from PredAlgo and WolfPSORT (Cp=chloroplast, 678 Mt=mitochondria) are shown alongside the log-transformed normalized expression level 679 of each candidate from different RNA-seq samples. 680 681 Figure 8: Single base insertion variant in DMT1a leads to premature stop codon 682 683 (A) The DMT1 locus on chromosome 6 containing DMT1a and DMT1b is shown below 684 its context in the entire chromosome. (B) Upper diagram shows enlarged chromosomal 685 region containing DMT1 locus with genomic coordinates and adjacent genes. The 25 Downloaded from on June 17, 2017 - Published by www.plantphysiol.org Copyright © 2015 American Society of Plant Biologists. All rights reserved. 686 position of the insertion variant in DMT1a is indicated by a yellow triangle. The inset 687 below is an expanded view of exon 10 from DMT1a showing the single-base insertion in 688 the MT– strain (CJU10) and the altered protein coding sequence and premature stop 689 codon caused by the cytosine insertion shown in bold and marked with a yellow triangle. 690 691 Supplementary Figure 1: Candidate methyltransferases in C. reinhardtii 692 693 Candidate methyltransferases were identified by the presence of a DNA 694 methyltransferase domain. Predicted localizations from PredAlgo and WolfPSORT are 695 shown alongside the expression patterns of each candidate methyltransferase, scaled to 696 the maximum and minimum expression values per gene. 697 698 Supplementary Figure 2: Methylation patterns for Chlamydomonas mtDNA 699 700 Cytosine methylation is absent (1-2% per-cytosine methylation) in all samples, with the 701 exception of 24h zygotes (13% per-cytosine methylation). In the coverage plots, blue 702 represents unmethylated cytosines (converted by bisulfite treatment), red represents 703 methylated cytosines (unconverted), and gray represents the total coverage (including 704 reads from the opposite strand which do not provide methylation information). 705 706 26 Downloaded from on June 17, 2017 - Published by www.plantphysiol.org Copyright © 2015 American Society of Plant Biologists. All rights reserved. 707 708 709 710 711 712 713 714 715 716 717 718 719 720 721 722 723 724 725 726 727 728 729 730 731 732 733 734 735 736 737 738 739 740 741 742 743 744 745 746 747 748 749 750 751 Literature Cited Albee, A.J., Kwan, A.L., Lin, H., Granas, D., Stormo, G.D., and Dutcher, S.K. (2013). Identification of cilia genes that affect cell-cycle progression using wholegenome transcriptome analysis in Chlamydomonas reinhardtti. G3 3, 979991. Aoyama, H., Saitoh, S., Kuroiwa, T., and Nakamura, S. (2014). Comparative analysis of zygospore transcripts during early germination in Chlamydomonas reinhardtii. Journal of plant physiology 171, 1685-1692. Armbrust, E.V., Ferris, P.J., and Goodenough, U.W. (1993). A mating type-linked gene cluster expressed in Chlamydomonas zygotes participates in the uniparental inheritance of the chloroplast genome. Cell 74, 801-811. Avidor-Reiss, T., Khire, A., Fishman, E.L., and Jo, K.H. (2015). Atypical centrioles during sexual reproduction. Frontiers in cell and developmental biology 3, 21. Babinger, P., Kobl, I., Mages, W., and Schmitt, R. (2001). A link between DNA methylation and epigenetic silencing in transgenic Volvox carteri. Nucleic acids research 29, 1261-1271. Babinger, P., Volkl, R., Cakstina, I., Maftei, A., and Schmitt, R. (2007). Maintenance DNA methyltransferase (Met1) and silencing of CpG-methylated foreign DNA in Volvox carteri. Plant molecular biology 63, 325-336. Bestor, T.H., Edwards, J.R., and Boulard, M. (2014). Notes on the role of dynamic DNA methylation in mammalian development. Proceedings of the National Academy of Sciences of the United States of America. Blaby, I.K., Blaby-Haas, C.E., Tourasse, N., Hom, E.F., Lopez, D., Aksoy, M., Grossman, A., Umen, J., Dutcher, S., Porter, M., King, S., Witman, G.B., Stanke, M., Harris, E.H., Goodstein, D., Grimwood, J., Schmutz, J., Vallon, O., Merchant, S.S., and Prochnik, S. (2014). The Chlamydomonas genome project: a decade on. Trends in plant science 19, 672-680. Bolen, P.L., Grant, D.M., Swinton, D., Boynton, J.E., and Gillham, N.W. (1982). Extensive methylation of chloroplast DNA by a nuclear gene mutation does not affect chloroplast gene transmission in chlamydomonas. Cell 28, 335343. Boyle, N.R., Page, M.D., Liu, B., Blaby, I.K., Casero, D., Kropat, J., Cokus, S.J., Hong-Hermesdorf, A., Shaw, J., Karpowicz, S.J., Gallaher, S.D., Johnson, S., Benning, C., Pellegrini, M., Grossman, A., and Merchant, S.S. (2012). Three acyltransferases and nitrogen-responsive regulator are implicated in nitrogen starvation-induced triacylglycerol accumulation in Chlamydomonas. The Journal of biological chemistry 287, 15811-15825. Brawley, S.H., and Johnson, L.E. (1992). Gametogenesis, Gametes and Zygotes - an Ecological Perspective on Sexual Reproduction in the Algae. Brit Phycol J 27, 233-252. Brown, R.C., and Lemmon, B.E. (2011). Spores before sporophytes: hypothesizing the origin of sporogenesis at the algal-plant transition. The New phytologist 190, 875-881. 27 Downloaded from on June 17, 2017 - Published by www.plantphysiol.org Copyright © 2015 American Society of Plant Biologists. All rights reserved. 752 753 754 755 756 757 758 759 760 761 762 763 764 765 766 767 768 769 770 771 772 773 774 775 776 777 778 779 780 781 782 783 784 785 786 787 788 789 790 791 792 793 794 795 796 Burton, W.G., Grabowy, C.T., and Sager, R. (1979). Role of methylation in the modification and restriction of chloroplast DNA in Chlamydomonas. Proceedings of the National Academy of Sciences of the United States of America 76, 1390-1394. Callebaut, I., Courvalin, J.C., and Mornon, J.P. (1999). The BAH (bromo-adjacent homology) domain: a link between DNA methylation, replication and transcriptional regulation. FEBS letters 446, 189-193. Catt, J.W. (1979). Isolation and Chemical Composition of the Zygospore Cell-Wall of Chlamydomonas-Reinhardii. Plant Sci Lett 15, 69-74. Cavalier-Smith, T. (1976). Electron microscopy of zygospore formation in Chlamydomonas reinhardii. Protoplasma 87, 297-315. Cerutti, H., Johnson, A.M., Gillham, N.W., and Boynton, J.E. (1997). Epigenetic silencing of a foreign gene in nuclear transformants of Chlamydomonas. The Plant cell 9, 925-945. De Hoff, P.L., Ferris, P., Olson, B.J., Miyagi, A., Geng, S., and Umen, J.G. (2013). Species and population level molecular profiling reveals cryptic recombination and emergent asymmetry in the dimorphic mating locus of C. reinhardtii. PLoS genetics 9, e1003724. DePristo, M.A., Banks, E., Poplin, R., Garimella, K.V., Maguire, J.R., Hartl, C., Philippakis, A.A., del Angel, G., Rivas, M.A., Hanna, M., McKenna, A., Fennell, T.J., Kernytsky, A.M., Sivachenko, A.Y., Cibulskis, K., Gabriel, S.B., Altshuler, D., and Daly, M.J. (2011). A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nature genetics 43, 491-498. Duncan, L., Nishii, I., Howard, A., Kirk, D., and Miller, S.M. (2006). Orthologs and paralogs of regA, a master cell-type regulatory gene in Volvox carteri. Current genetics 50, 61-72. Duncan, L., Nishii, I., Harryman, A., Buckley, S., Howard, A., Friedman, N.R., and Miller, S.M. (2007). The VARL gene family and the evolutionary origins of the master cell-type regulatory gene, regA, in Volvox carteri. Journal of molecular evolution 65, 1-11. Dyer, T.A. (1982). Methylation of chloroplast DNA in Chlamydomonas. Nature 298, 422-423. Fang, W., Si, Y., Douglass, S., Casero, D., Merchant, S.S., Pellegrini, M., Ladunga, I., Liu, P., and Spalding, M.H. (2012). Transcriptome-wide changes in Chlamydomonas reinhardtii gene expression regulated by carbon dioxide and the CO2-concentrating mechanism regulator CIA5/CCM1. The Plant cell 24, 1876-1893. Feng, S., Cokus, S.J., Zhang, X., Chen, P.Y., Bostick, M., Goll, M.G., Hetzel, J., Jain, J., Strauss, S.H., Halpern, M.E., Ukomadu, C., Sadler, K.C., Pradhan, S., Pellegrini, M., and Jacobsen, S.E. (2010). Conservation and divergence of methylation patterning in plants and animals. Proceedings of the National Academy of Sciences of the United States of America 107, 8689-8694. Feng, T.Y., and Chiang, K.S. (1984). The persistence of maternal inheritance in Chlamydomonas despite hypomethylation of chloroplast DNA induced by 28 Downloaded from on June 17, 2017 - Published by www.plantphysiol.org Copyright © 2015 American Society of Plant Biologists. All rights reserved. 797 798 799 800 801 802 803 804 805 806 807 808 809 810 811 812 813 814 815 816 817 818 819 820 821 822 823 824 825 826 827 828 829 830 831 832 833 834 835 836 837 838 839 840 841 842 inhibitors. Proceedings of the National Academy of Sciences of the United States of America 81, 3438-3442. Ferris, P.J., and Goodenough, U.W. (1987). Transcription of novel genes, including a gene linked to the mating-type locus, induced by Chlamydomonas fertilization. Molecular and cellular biology 7, 2360-2366. Ferris, P.J., Armbrust, E.V., and Goodenough, U.W. (2002). Genetic structure of the mating-type locus of Chlamydomonas reinhardtii. Genetics 160, 181-200. Ferris, P.J., Waffenschmidt, S., Umen, J.G., Lin, H., Lee, J.H., Ishida, K., Kubo, T., Lau, J., and Goodenough, U.W. (2005). Plus and minus sexual agglutinins from Chlamydomonas reinhardtii. The Plant cell 17, 597-615. Gallaher, S.D., Fitz-Gibbon, S.T., Glaesener, A.G., Pellegrini, M., and Merchant, S.S. (2015). Chlamydomonas Genome Resource for Laboratory Strains Reveals a Mosaic of Sequence Variation, Identifies True Strain Histories, and Enables Strain-Specific Studies. The Plant cell. Glaesener, A.G., Merchant, S.S., and Blaby-Haas, C.E. (2013). Iron economy in Chlamydomonas reinhardtii. Frontiers in plant science 4, 337. Goldschmidt-Clermont, M. (1991). Transgenic expression of aminoglycoside adenine transferase in the chloroplast: a selectable marker of site-directed transformation of chlamydomonas. Nucleic acids research 19, 4083-4089. Gonzalez-Ballester, D., Casero, D., Cokus, S., Pellegrini, M., Merchant, S.S., and Grossman, A.R. (2010). RNA-seq analysis of sulfur-deprived Chlamydomonas cells reveals aspects of acclimation critical for cell survival. The Plant cell 22, 2058-2084. Goodenough, U., Lin, H., and Lee, J.H. (2007). Sex determination in Chlamydomonas. Seminars in cell & developmental biology 18, 350-361. Grief, C., O'Neill, M.A., and Shaw, P.J. (1987). The zygote cell wall of Chlamydomonas reinhardii: a structural, chemical and immunological approach. Planta 170, 433-445. Grossman, A. (2000). Acclimation of Chlamydomonas reinhardtii to its nutrient environment. Protist 151, 201-224. Guo, W., Fiziev, P., Yan, W., Cokus, S., Sun, X., Zhang, M.Q., Chen, P.Y., and Pellegrini, M. (2013). BS-Seeker2: a versatile aligning pipeline for bisulfite sequencing data. BMC genomics 14, 774. Harris, E.H., Stern, D.B., Witman, G., Harris, E.H., and Harris, E.H. (2009). The Chlamydomonas sourcebook. (Amsterdam ; Boston: Elsevier/Academic Press). Harris, M.A., Clark, J., Ireland, A., Lomax, J., Ashburner, M., Foulger, R., Eilbeck, K., Lewis, S., Marshall, B., Mungall, C., Richter, J., Rubin, G.M., Blake, J.A., Bult, C., Dolan, M., Drabkin, H., Eppig, J.T., Hill, D.P., Ni, L., Ringwald, M., Balakrishnan, R., Cherry, J.M., Christie, K.R., Costanzo, M.C., Dwight, S.S., Engel, S., Fisk, D.G., Hirschman, J.E., Hong, E.L., Nash, R.S., Sethuraman, A., Theesfeld, C.L., Botstein, D., Dolinski, K., Feierbach, B., Berardini, T., Mundodi, S., Rhee, S.Y., Apweiler, R., Barrell, D., Camon, E., Dimmer, E., Lee, V., Chisholm, R., Gaudet, P., Kibbe, W., Kishore, R., Schwarz, E.M., Sternberg, P., Gwinn, M., Hannick, L., Wortman, J., Berriman, M., Wood, V., de la Cruz, N., Tonellato, P., Jaiswal, P., Seigfried, T., White, R., and 29 Downloaded from on June 17, 2017 - Published by www.plantphysiol.org Copyright © 2015 American Society of Plant Biologists. All rights reserved. 843 844 845 846 847 848 849 850 851 852 853 854 855 856 857 858 859 860 861 862 863 864 865 866 867 868 869 870 871 872 873 874 875 876 877 878 879 880 881 882 883 884 885 886 887 Gene Ontology, C. (2004). The Gene Ontology (GO) database and informatics resource. Nucleic acids research 32, D258-261. Hattman, S., Kenny, C., Berger, L., and Pratt, K. (1978). Comparative study of DNA methylation in three unicellular eucaryotes. Journal of bacteriology 135, 1156-1157. Heinnickel, M.L., and Grossman, A.R. (2013). The GreenCut: re-evaluation of physiological role of previously studied proteins and potential novel protein functions. Photosynthesis research 116, 427-436. Idoine, A.D., Boulouis, A., Rupprecht, J., and Bock, R. (2014). The diurnal logic of the expression of the chloroplast genome in Chlamydomonas reinhardtii. PloS one 9, e108760. Kanehisa, M., and Goto, S. (2000). KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic acids research 28, 27-30. Karpowicz, S.J., Prochnik, S.E., Grossman, A.R., and Merchant, S.S. (2011). The GreenCut2 resource, a phylogenomically derived inventory of proteins specific to the plant lineage. The Journal of biological chemistry 286, 2142721439. Kirk, M.M., Stark, K., Miller, S.M., Muller, W., Taillon, B.E., Gruber, H., Schmitt, R., and Kirk, D.L. (1999). regA, a Volvox gene that plays a central role in germ-soma differentiation, encodes a novel regulatory protein. Development 126, 639-647. Kubo, T., Abe, J., Oyamada, T., Ohnishi, M., Fukuzawa, H., Matsuda, Y., and Saito, T. (2008). Characterization of novel genes induced by sexual adhesion and gamete fusion and of their transcriptional regulation in Chlamydomonas reinhardtii. Plant & cell physiology 49, 981-993. Kuriyama, H., Takano, H., Suzuki, L., Uchida, H., Kawano, S., Kuroiwa, H., and Kuroiwa, T. (1999). Characterization of Chlamydomonas reinhardtii zygotespecific cDNAs that encode novel proteins containing ankyrin repeats and WW domains. Plant physiology 119, 873-884. Kurvari, V., Grishin, N.V., and Snell, W.J. (1998). A gamete-specific, sex-limited homeodomain protein in Chlamydomonas. The Journal of cell biology 143, 1971-1980. Lechtreck, K.F., Rostmann, J., and Grunow, A. (2002). Analysis of Chlamydomonas SF-assemblin by GFP tagging and expression of antisense constructs. Journal of cell science 115, 1511-1522. Lee, J.H., Lin, H., Joo, S., and Goodenough, U. (2008). Early sexual origins of homeoprotein heterodimerization and evolution of the plant KNOX/BELL family. Cell 133, 829-840. Liu, Y., Misamore, M.J., and Snell, W.J. (2010). Membrane fusion triggers rapid degradation of two gamete-specific, fusion-essential proteins in a membrane block to polygamy in Chlamydomonas. Development 137, 1473-1481. Lopez, D., Casero, D., Cokus, S.J., Merchant, S.S., and Pellegrini, M. (2011). Algal Functional Annotation Tool: a web-based analysis suite to functionally interpret large gene lists using integrated annotation and expression data. BMC bioinformatics 12, 282. 30 Downloaded from on June 17, 2017 - Published by www.plantphysiol.org Copyright © 2015 American Society of Plant Biologists. All rights reserved. 888 889 890 891 892 893 894 895 896 897 898 899 900 901 902 903 904 905 906 907 908 909 910 911 912 913 914 915 916 917 918 919 920 921 922 923 924 925 926 927 928 929 930 931 932 Lv, H., Qu, G., Qi, X., Lu, L., Tian, C., and Ma, Y. (2013). Transcriptome analysis of Chlamydomonas reinhardtii during the process of lipid accumulation. Genomics 101, 229-237. Matsuo, M., Hachisu, R., Tabata, S., Fukuzawa, H., and Obokata, J. (2011). Transcriptome analysis of respiration-responsive genes in Chlamydomonas reinhardtii: mitochondrial retrograde signaling coordinates the genes for cell proliferation with energy-producing metabolism. Plant & cell physiology 52, 333-343. Matters, G.L., and Goodenough, U.W. (1992). A gene/pseudogene tandem duplication encodes a cysteine-rich protein expressed during zygote development in Chlamydomonas reinhardtii. Molecular & general genetics : MGG 232, 81-88. Maul, J.E., Lilly, J.W., Cui, L., dePamphilis, C.W., Miller, W., Harris, E.H., and Stern, D.B. (2002). The Chlamydomonas reinhardtii plastid chromosome: islands of genes in a sea of repeats. The Plant cell 14, 2659-2679. McKenna, A., Hanna, M., Banks, E., Sivachenko, A., Cibulskis, K., Kernytsky, A., Garimella, K., Altshuler, D., Gabriel, S., Daly, M., and DePristo, M.A. (2010). The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome research 20, 1297-1303. Merchant, S.S., Allen, M.D., Kropat, J., Moseley, J.L., Long, J.C., Tottey, S., and Terauchi, A.M. (2006). Between a rock and a hard place: trace element nutrition in Chlamydomonas. Biochimica et biophysica acta 1763, 578-594. Merchant, S.S., Prochnik, S.E., Vallon, O., Harris, E.H., Karpowicz, S.J., Witman, G.B., Terry, A., Salamov, A., Fritz-Laylin, L.K., Marechal-Drouard, L., Marshall, W.F., Qu, L.H., Nelson, D.R., Sanderfoot, A.A., Spalding, M.H., Kapitonov, V.V., Ren, Q., Ferris, P., Lindquist, E., Shapiro, H., Lucas, S.M., Grimwood, J., Schmutz, J., Cardol, P., Cerutti, H., Chanfreau, G., Chen, C.L., Cognat, V., Croft, M.T., Dent, R., Dutcher, S., Fernandez, E., Fukuzawa, H., Gonzalez-Ballester, D., Gonzalez-Halphen, D., Hallmann, A., Hanikenne, M., Hippler, M., Inwood, W., Jabbari, K., Kalanon, M., Kuras, R., Lefebvre, P.A., Lemaire, S.D., Lobanov, A.V., Lohr, M., Manuell, A., Meier, I., Mets, L., Mittag, M., Mittelmeier, T., Moroney, J.V., Moseley, J., Napoli, C., Nedelcu, A.M., Niyogi, K., Novoselov, S.V., Paulsen, I.T., Pazour, G., Purton, S., Ral, J.P., Riano-Pachon, D.M., Riekhof, W., Rymarquis, L., Schroda, M., Stern, D., Umen, J., Willows, R., Wilson, N., Zimmer, S.L., Allmer, J., Balk, J., Bisova, K., Chen, C.J., Elias, M., Gendler, K., Hauser, C., Lamb, M.R., Ledford, H., Long, J.C., Minagawa, J., Page, M.D., Pan, J., Pootakham, W., Roje, S., Rose, A., Stahlberg, E., Terauchi, A.M., Yang, P., Ball, S., Bowler, C., Dieckmann, C.L., Gladyshev, V.N., Green, P., Jorgensen, R., Mayfield, S., Mueller-Roeber, B., Rajamani, S., Sayre, R.T., Brokstein, P., Dubchak, I., Goodstein, D., Hornick, L., Huang, Y.W., Jhaveri, J., Luo, Y., Martinez, D., Ngau, W.C., Otillar, B., Poliakov, A., Porter, A., Szajkowski, L., Werner, G., Zhou, K., Grigoriev, I.V., Rokhsar, D.S., and Grossman, A.R. (2007). The Chlamydomonas genome reveals the evolution of key animal and plant functions. Science 318, 245-250. 31 Downloaded from on June 17, 2017 - Published by www.plantphysiol.org Copyright © 2015 American Society of Plant Biologists. All rights reserved. 933 934 935 936 937 938 939 940 941 942 943 944 945 946 947 948 949 950 951 952 953 954 955 956 957 958 959 960 961 962 963 964 965 966 967 968 969 970 971 972 973 974 975 976 977 978 Miller, R., Wu, G., Deshpande, R.R., Vieler, A., Gartner, K., Li, X., Moellering, E.R., Zauner, S., Cornish, A.J., Liu, B., Bullard, B., Sears, B.B., Kuo, M.H., Hegg, E.L., Shachar-Hill, Y., Shiu, S.H., and Benning, C. (2010). Changes in transcript abundance in Chlamydomonas reinhardtii following nitrogen deprivation predict diversion of metabolism. Plant physiology 154, 17371752. Miller, S.M., Schmitt, R., and Kirk, D.L. (1993). Jordan, an active Volvox transposable element similar to higher plant transposons. The Plant cell 5, 1125-1138. Minami, S.A., and Goodenough, U.W. (1978). Novel glycopolypeptide synthesis induced by gametic cell fusion in Chlamydomonas reinhardtii. The Journal of cell biology 77, 165-181. Nakamura, S. (2010). Paternal inheritance of mitochondria in Chlamydomonas. Journal of plant research 123, 163-170. Nguyen, A.V., Thomas-Hall, S.R., Malnoe, A., Timmins, M., Mussgnug, J.H., Rupprecht, J., Kruse, O., Hankamer, B., and Schenk, P.M. (2008). Transcriptome for photobiological hydrogen production induced by sulfur deprivation in the green alga Chlamydomonas reinhardtii. Eukaryotic cell 7, 1965-1979. Ning, J., Otto, T.D., Pfander, C., Schwach, F., Brochet, M., Bushell, E., Goulding, D., Sanders, M., Lefebvre, P.A., Pei, J., Grishin, N.V., Vanderlaan, G., Billker, O., and Snell, W.J. (2013). Comparative genomics in Chlamydomonas and Plasmodium identifies an ancient nuclear envelope protein family essential for sexual reproduction in protists, fungi, plants, and vertebrates. Genes & development 27, 1198-1215. Nishimura, Y. (2010). Uniparental inheritance of cpDNA and the genetic control of sexual differentiation in Chlamydomonas reinhardtii. Journal of plant research 123, 149-162. Nishimura, Y., Misumi, O., Kato, K., Inada, N., Higashiyama, T., Momoyama, Y., and Kuroiwa, T. (2002). An mt(+) gamete-specific nuclease that targets mt() chloroplasts during sexual reproduction in C. reinhardtii. Genes & development 16, 1116-1128. Nishiyama, R., Ito, M., Yamaguchi, Y., Koizumi, N., and Sano, H. (2002). A chloroplast-resident DNA methyltransferase is responsible for hypermethylation of chloroplast genes in Chlamydomonas maternal gametes. Proceedings of the National Academy of Sciences of the United States of America 99, 5925-5930. Nishiyama, R., Wada, Y., Mibu, M., Yamaguchi, Y., Shimogawara, K., and Sano, H. (2004). Role of a nonselective de novo DNA methyltransferase in maternal inheritance of chloroplast genes in the green alga, Chlamydomonas reinhardtii. Genetics 168, 809-816. Panchy, N., Wu, G., Newton, L., Tsai, C.H., Chen, J., Benning, C., Farre, E.M., and Shiu, S.H. (2014). Prevalence, Evolution, and cis-Regulation of Diel Transcription in Chlamydomonas reinhardtii. G3 4, 2461-2471. Prochnik, S.E., Umen, J., Nedelcu, A.M., Hallmann, A., Miller, S.M., Nishii, I., Ferris, P., Kuo, A., Mitros, T., Fritz-Laylin, L.K., Hellsten, U., Chapman, J., 32 Downloaded from on June 17, 2017 - Published by www.plantphysiol.org Copyright © 2015 American Society of Plant Biologists. All rights reserved. 979 980 981 982 983 984 985 986 987 988 989 990 991 992 993 994 995 996 997 998 999 1000 1001 1002 1003 1004 1005 1006 1007 1008 1009 1010 1011 1012 1013 1014 1015 1016 1017 1018 1019 1020 1021 1022 1023 1024 Simakov, O., Rensing, S.A., Terry, A., Pangilinan, J., Kapitonov, V., Jurka, J., Salamov, A., Shapiro, H., Schmutz, J., Grimwood, J., Lindquist, E., Lucas, S., Grigoriev, I.V., Schmitt, R., Kirk, D., and Rokhsar, D.S. (2010). Genomic analysis of organismal complexity in the multicellular green alga Volvox carteri. Science 329, 223-226. Rochaix, J.D. (2001). Assembly,function, and dynamics of the photosynthetic machinery in Chlamydomonas reinhardtii. Plant physiology 127, 1394-1398. Royer, H.D., and Sager, R. (1979). Methylation of chloroplast DNAs in the life cycle of Chlamydomonas. Proceedings of the National Academy of Sciences of the United States of America 76, 5794-5798. Sano, H., Grabowy, C., and Sager, R. (1981). Differential activity of DNA methyltransferase in the life cycle of Chlamydomonas reinhardi. Proceedings of the National Academy of Sciences of the United States of America 78, 3118-3122. Schmollinger, S., Muhlhaus, T., Boyle, N.R., Blaby, I.K., Casero, D., Mettler, T., Moseley, J.L., Kropat, J., Sommer, F., Strenkert, D., Hemme, D., Pellegrini, M., Grossman, A.R., Stitt, M., Schroda, M., and Merchant, S.S. (2014). Nitrogen-Sparing Mechanisms in Chlamydomonas Affect the Transcriptome, the Proteome, and Photosynthetic Metabolism. The Plant cell 26, 1410-1435. Silflow, C.D., and Lefebvre, P.A. (2001). Assembly and motility of eukaryotic cilia and flagella. Lessons from Chlamydomonas reinhardtii. Plant physiology 127, 1500-1507. Simon, D.F., Descombes, P., Zerges, W., and Wilkinson, K.J. (2008). Global expression profiling of Chlamydomonas reinhardtii exposed to trace levels of free cadmium. Environmental toxicology and chemistry / SETAC 27, 16681675. Simon, D.F., Domingos, R.F., Hauser, C., Hutchins, C.M., Zerges, W., and Wilkinson, K.J. (2013). Transcriptome sequencing (RNA-seq) analysis of the effects of metal nanoparticle exposure on the transcriptome of Chlamydomonas reinhardtii. Applied and environmental microbiology 79, 4774-4785. Smith, Z.D., and Meissner, A. (2013). DNA methylation: roles in mammalian development. Nature reviews. Genetics 14, 204-220. Suzuki, L., Woessner, J.P., Uchida, H., Kuroiwa, H., Yuasa, Y., Waffenschmidt, S., Goodenough, U.W., and Kuroiwa, T. (2000). Zygote-specific protein with hydroxyproline-rich glycoprotein domains and lectin-like domains involved in the assembly of the cell wall of Chlamydomonas reinhardtii (Chlorophyta). J Phycol 36, 571-583. Swanson, W.J., and Vacquier, V.D. (2002). The rapid evolution of reproductive proteins. Nature reviews. Genetics 3, 137-144. Tam, L.W., Ranum, P.T., and Lefebvre, P.A. (2013). CDKL5 regulates flagellar length and localizes to the base of the flagella in Chlamydomonas. Molecular biology of the cell 24, 588-600. Thimm, O., Blasing, O., Gibon, Y., Nagel, A., Meyer, S., Kruger, P., Selbig, J., Muller, L.A., Rhee, S.Y., and Stitt, M. (2004). MAPMAN: a user-driven tool to display genomics data sets onto diagrams of metabolic pathways and other 33 Downloaded from on June 17, 2017 - Published by www.plantphysiol.org Copyright © 2015 American Society of Plant Biologists. All rights reserved. 1025 1026 1027 1028 1029 1030 1031 1032 1033 1034 1035 1036 1037 1038 1039 1040 1041 1042 1043 1044 1045 1046 1047 1048 1049 1050 1051 1052 1053 1054 1055 1056 1057 1058 1059 1060 biological processes. The Plant journal : for cell and molecular biology 37, 914-939. Thorvaldsdottir, H., Robinson, J.T., and Mesirov, J.P. (2013). Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration. Briefings in bioinformatics 14, 178-192. Toepel, J., Albaum, S.P., Arvidsson, S., Goesmann, A., la Russa, M., Rogge, K., and Kruse, O. (2011). Construction and evaluation of a whole genome microarray of Chlamydomonas reinhardtii. BMC genomics 12, 579. Toepel, J., Illmer-Kephalides, M., Jaenicke, S., Straube, J., May, P., Goesmann, A., and Kruse, O. (2013). New insights into Chlamydomonas reinhardtii hydrogen production processes by combined microarray/RNA-seq transcriptomics. Plant biotechnology journal 11, 717-733. Uchida, H., Kawano, S., Sato, N., and Kuroiwa, T. (1993). Isolation and characterization of novel genes which are expressed during the very early stage of zygote formation in Chlamydomonas reinhardtii. Current genetics 24, 296-300. Uchida, H., Suzuki, L., Anai, T., Doi, K., Takano, H., Yamashita, H., Oka, T., Kawano, S., Tomizawa, K.I., Kawazu, T., Kuroiwa, H., and Kuroiwa, T. (1999). A pair of invertedly repeated genes in Chlamydomonas reinhardtii encodes a zygote-specific protein whose expression is UV-sensitive. Current genetics 36, 232-240. Ueki, N., and Nishii, I. (2008). Idaten is a new cold-inducible transposon of Volvox carteri that can be used for tagging developmentally important genes. Genetics 180, 1343-1353. Umen, J.G. (2011). Evolution of sex and mating loci: an expanded view from Volvocine algae. Current opinion in microbiology 14, 634-641. Umen, J.G., and Goodenough, U.W. (2001). Chloroplast DNA methylation and inheritance in Chlamydomonas. Genes & development 15, 2585-2597. Woessner, J.P., Molendijk, A.J., van Egmond, P., Klis, F.M., Goodenough, U.W., and Haring, M.A. (1994). Domain conservation in several volvocalean cell wall proteins. Plant molecular biology 26, 947-960. Yang, N., and Xu, R.M. (2013). Structure and function of the BAH domain in chromatin biology. Critical reviews in biochemistry and molecular biology 48, 211-221. 34 Downloaded from on June 17, 2017 - Published by www.plantphysiol.org Copyright © 2015 American Society of Plant Biologists. All rights reserved. CJU10 (mt–) R3 (mt+) vegetative cells Genome (Whole) Transcriptome (RNA-seq) 1 2 1 2 Methylome (BS-Seq) gametogenesis 1 2 Replicates mating gametic cells 1 2 fertilization cell fusion zygote development !! ! ! ! ! ! !! ! !! ! ! new zygote pooled (0.5, 1, 2 hrs) !! !! ! ! ! ! ! ! ! ! !! ! ! ! ! degradation of mt– cpDNA nuclear & organellar fusion 1h 2-4 h !! ! ! ! zygospore maturation 24 h meiosis & cpDNA replication daughter cells 6d Figure 1. Chlamydomonas reinhardtii sexual life cycle and sequenced samples Vegetative C. reinhardtii cells of each mating type (MT+ and MT−) can be induced to undergo gametogenesis by nitrogen starvation (−N). Gametes of opposite mating type recognize each other through flagellar adhesion and fuse to form a diploid zygote. During zygote maturation MT− chloroplast DNA (cpDNA) is eliminated, flagella are resorbed, and a thick, zygote cell wall forms. Upon return to nitrogen and light zygospores undergo meiosis to form four haploid progeny (two of each mating type, all containing uniparentally inherited MT+ parentalDownloaded cpDNA) that reenter the vegetative cycle. Colored boxes designate samples and material sequenced. Blue from on June 17, 2017 - Published by www.plantphysiol.org and red unfilled circlesCopyright represent MT+ and MT– chloroplast genomes. Filled pink and light blue circles represent nuclear DNA from © 2015 American Society of Plant Biologists. All rights reserved. MT+ and MT– strains respectively. A B Genetic Differences R3 (MT+) è CJU (MT–) Category mating locus Loci Total Bases 101,452 101,452 Insertions 7,510 35,999 Deletions 8,113 41,209 SNVs Chr. 6 C R3 (mt+) IRB Chr. 17 atpB Insertion (cpDNA) CJU (mt–) IRB rbcL (promoter) aadA psbC atpB (UTR) Figure 2. Genetic differences between R3 (MT+) and CJU10 (MT–) parental strains (A) Genetic differences, categorized by variant-type (single nucleotide variants [SNVs], insertions, and deletions) are shown with total number of variant loci and total bases. (B) Large-scale view of chromosomes 6 and 17, where 98.5% of the identified variants are located. Locations of sequence differences are shown in red. The region of chromosome 6 containing the mating locus is denoted by a blue bar. (C) Schematic representation of the antibiotic resistance transgene insertion in CJU10 chloroplast DNA. The aadA gene, along with the rbcL 5’ promoter and psbC 3’ UTR, is inserted between inverted repeat region “B” (IRB) and the endogenous atpB gene. Downloaded from on June 17, 2017 - Published by www.plantphysiol.org Copyright © 2015 American Society of Plant Biologists. All rights reserved. gamete-specific (mt+) 293 genes gamete-specific (mt–) 68 genes ValsOfGenesInNewListsInNewData −− JU ChlReiV5 −− Tophat_to_Cufflinks −− gamPlus ValsOfGenesInNewListsInNewData −− JU ChlReiV5 −− Tophat_to_Cufflinks −− gamMinus 1.0 0.8 ● gamPlus gamMinus zygote 100% 2628 14 78 64 * ** ** ** ** 80% 10 3361 72 45 70% 60% 131 102 19 126 102 10 17 2836 71 55 95 73 9 50% 0.0 vegPlus vegMinus gamPlus C 187 129 * 216 275 ** ** 60% ** 40% 50% 50% gamPlus 7 97 93 ** 17 18 17 * 13 15 14 12 6058 6058 50 Secretory Chloroplast ** 2712 20% 20% 20% 50 68 3 2 11* *3 8164 81 ** ** ** 179 179 ** 153** 45 45 46 35 ** 46 37 78 ** ** * 78 * 73 34 ** 2712 143 ** 2712 34 143 191 ** 38 ** ** ** 198 ** 198 131 ** 3605 10% 10% zygote * 1 5 7 5 2 1300 Mitochondria 1300 1300 878 Others 40% 40% 878 878 30% 30% gamMinus **52 ** 13 14 70% 70% 60% 60% vegMinus 21 ** 21 ** 29 80% 80% 6058 10% 10% vegPlus 52 80% 50% Other 20% 43 30 veg. veg. gam. gam. zyg. (mt+) (mt–) (mt+) (mt–) 3191 90% 90% 90%3191 3191 30% 30% zygote 100% 100% 100% 70% 40% 8916 gamMinus veg. veg. gam. gam. (mt+) (mt–) (mt+) (mt–) Secretory vegMinus ● ● ● ● ● ● ● ● ● ● ● ● 0.0 0.0 ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● veg. veg. gam. gam. (mt+) (mt–) (mt+) (mt–) vegPlus ● ● ● ● ● Chloro. ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● 0.2 0.2 ● ● ● ● ● ● ● ● ● ● ● ● ● 0.2 0.4 0.4 0.4 0.6 0.6 0.6 0.4 90% ValsOfGenesInNewListsInNewData −− JU ChlReiV5 −− Tophat_to_Cufflinks −− zygote 1.0 0.8 0.6 Min B zygote-specific 627 genes 0.8 0.8 0.2 Gene Expression 1.0 Max Mito. A ** 3605 3605 0% 0% all all + – gam. mt minus gam. mt plus plus gamete gamete minus gamete gamete zygote zygote zygote 0% all 0% plus allall gam. mt+ plus minus gam. mt– minus zygote zygote zygote Figure 3. Gamete- and zygote-specific genes (A) Box and whisker plots of gene expression profiles for gamete- and zygote-specific genes are plotted. Values are plotted relative to the maximum expression value of each of the genes. 293 and 68 genes were expressed specifically in MT+ and MT– gametes respectively. 627 genes were expressed specifically in zygotes. (B) Localization predictions for proteins encoded by gamete (gam. MT+, gam. MT−) and zygote-specific genes, compared to all predicted proteins in the Chlamydomonas proteome (all). Data are plotted as fraction predicted for each of four compartments with total numbers indicated within each graph portion: secretory in blue; chloroplast (chloro.) in green; mitochondria (mito.) in purple; other in gray. Starred values are significantly different from the total proteome distribution (*: p < 0.05; **: p < 0.01). (C) Predicted protein groups as described in 3B are plotted according to taxonomic distribution of encoded proteins from each category. From bottom (Chlamydomonas-specific proteins, dark green) to top (all cellular organisms, gray) are increasingly broad phylogenetic distributions. Starred samples indicate significant enrichment or depletion in a taxonomic category relative to the distribution of all proteins. A B GreenCut2 Genes 6% 4% 2.0% p = 0.0098 p = 0.0018 1.6% Genes In List Genes In List 5% 555/ 17,741 3% 2% 9/627 1% 1/293 0% all C CiliaCut Genes gam. (mt+) 1.2% 1/68 161/ 17,741 0.8% 3/627 1/293 0.4% 0/68 gam. (mt–) 0.0% zygote MapMan: "Cell Wall" all D gam. (mt+) gam. (mt–) zygote MapMan: "Transport" 2.0% 1.5% 7/559 1.0% 0.5% p = 1.386 x 10-3 6.0% Genes In List Genes In List p = 7.276 x 10-5 24/559 4.0% 368/ 17,302 2.0% 33/ 17,302 0.0% 0.0% all zygote all zygote Figure 4. Functional annotation of proteins encoded by gamete- and zygote-specific genes Functional annotations for gamete- and zygote-specific gene lists are shown as bar plots showing the fraction of total proteins in each group with the specified annotation. Total number of genes in each annotation category versus number of genes with the described annotation are shown. Significant differences between all genes and gamete or zygote specific genes were calculated with the hypergeometric test and significant p-values are shown. (A) GreenCut2 genes (Karpowicz et al., 2011; Heinnickel and Grossman, 2013) (B) CiliaCut genes (Merchant et al., 2007). (C) MapMan cell wall related genes (Thimm et al., 2004). (D) MapMan transport related genes. A B C D veg.(mt+) veg.(mt–) gam. (mt+) gam. (mt–) zyg. (1h) zyg. (24h) zyg. (6d) zyg. (germ.) Genes Figure 5. Nuclear methylation at different C. reinhardtii life cycle stages (A) Bulk averages of 5-methyl cytosine in the nuclear genome are shown for each sample and categorized by the fraction of 5meC in each of three sequence contexts (CG, CHG and CHH). (B) Scaled representations of the 17 C. reinhardtii chromosomes with 23 hypermethylated regions larger than 10kb shaded in blue. (C) Plots comparing gene and repeat densities of hypermethylated regions relative to the entire genome. (D) Genome browser display of nuclear cytosine methylation frequency at two representative hypermethylated loci from chromosome 6. Two representative examples of strain-specific hypermethylated loci are shown side-by-side, with a gene track below for genomic context. In the coverage plots, blue represents unmethylated cytosines (converted by bisulfite treatment), red represents methylated cytosines (unconverted), and gray represents the total coverage (including reads from the opposite strand which do not provide methylation information). 60 CG CHG CHH 0 10 % Methylation 20 30 40 50 A mt S1 + – S3 + S5 S7 – +/– +/– +/– +/– 1h 24h 6d (g) time vegetative gametes zygotes B Figure 6. Chloroplast genome methylation at different life cycle stages (A) Bulk cytosine methylation frequencies for chloroplast DNA at different life cycle stages are plotted for each methylation context. Mating type of each sample (+, - or diploid +/-) and time of zygote development and germination “(g)” are shown below the graph. (B) Plot of average CG methylation frequency for each sample in 1kb bins across the chloroplast genome. The two inverted repeat regions are shaded in gray. Plots are color coded according to the legend on the right. Protein Architecture Predicted Localization Color Key 20 40 60 Life Cycle Expression 80 100 Value Cre06.g249350 (dmt1a) Cp Cre06.g249500 (dmt1b) Mt Cp 80 1 60 2 40 veg. (mt–) gam. (mt+) gam. (mt–) 20 V5 veg. (mt+) V4 PredAlgo WolfPSORT V3 DNA methylase domain Bromo adjacent homology (BAH) domain 3 Cp V2 Cp V1 Cre06.g286650 (dmt4) 0 zyg. Figure 7. Candidate chloroplast methyltransferases in C. reinhardtii The protein domain structure of each candidate methyltransferase is shown schematically (green = DNA methylase domain; red = bromo adjacent homology (BAH) domain). Predicted localizations from PredAlgo and WolfPSORT (Cp=chloroplast, Mt=mitochondria) are shown alongside the normalized expression level of each candidate from different RNA-seq samples. FPKM A D. Cre06.g249350.t1.1 dmt1a Cre06.g249400.t1.1 Cre06.g249500.t1.1 dmt1b Cre06.g249450.t2.1 B ………… Exon 10 R3 (mt+) ATTCCAG GTCCGCCGTTTCCTGCGGCCCGAGGACCTGTGTCCGACCTGGCC P CJU (mt–) P F P A A R G P V S D L A ATTCCAG GTCCGCCGTTTCCTGCGGCCCGAGGACCTGTGTCCCGACCTGGC P P F P A A R G P V S R P G ……… Exon 10 GGCCAGGAGGAGGAGGCTGAGTGGGCG G Q E E E A E W A CGGCCAGGAGGAGGAGGCTGAGTGGGC T R P G G G * Figure 8. Single base insertion variant in DMT1a leads to premature stop codon (A) The DMT1 locus on chromosome 6 containing DMT1a and DMT1b is shown below its context in the entire chromosome. (B) Upper diagram shows enlarged chromosomal region containing DMT1 locus with genomic coordinates and adjacent genes. The position of the insertion variant in DMT1a is indicated by a yellow triangle. The inset below is an expanded view of exon 10 from DMT1a showing the single-base insertion in the MT– strain (CJU10) and the altered protein coding sequence and premature stop codon caused by the cytosine insertion shown in bold and marked with a yellow triangle. Parsed Citations Albee, A.J., Kwan, A.L., Lin, H., Granas, D., Stormo, G.D., and Dutcher, S.K. (2013). Identification of cilia genes that affect cell-cycle progression using whole-genome transcriptome analysis in Chlamydomonas reinhardtti. G3 3, 979-991. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Aoyama, H., Saitoh, S., Kuroiwa, T., and Nakamura, S. (2014). Comparative analysis of zygospore transcripts during early germination in Chlamydomonas reinhardtii. Journal of plant physiology 171, 1685-1692. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Armbrust, E.V., Ferris, P.J., and Goodenough, U.W. (1993). A mating type-linked gene cluster expressed in Chlamydomonas zygotes participates in the uniparental inheritance of the chloroplast genome. Cell 74, 801-811. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Avidor-Reiss, T., Khire, A., Fishman, E.L., and Jo, K.H. (2015). Atypical centrioles during sexual reproduction. Frontiers in cell and developmental biology 3, 21. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Babinger, P., Kobl, I., Mages, W., and Schmitt, R. (2001). A link between DNA methylation and epigenetic silencing in transgenic Volvox carteri. Nucleic acids research 29, 1261-1271. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Babinger, P., Volkl, R., Cakstina, I., Maftei, A., and Schmitt, R. (2007). Maintenance DNA methyltransferase (Met1) and silencing of CpG-methylated foreign DNA in Volvox carteri. Plant molecular biology 63, 325-336. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Bestor, T.H., Edwards, J.R., and Boulard, M. (2014). Notes on the role of dynamic DNA methylation in mammalian development. Proceedings of the National Academy of Sciences of the United States of America. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Blaby, I.K., Blaby-Haas, C.E., Tourasse, N., Hom, E.F., Lopez, D., Aksoy, M., Grossman, A., Umen, J., Dutcher, S., Porter, M., King, S., Witman, G.B., Stanke, M., Harris, E.H., Goodstein, D., Grimwood, J., Schmutz, J., Vallon, O., Merchant, S.S., and Prochnik, S. (2014). The Chlamydomonas genome project: a decade on. Trends in plant science 19, 672-680. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Bolen, P.L., Grant, D.M., Swinton, D., Boynton, J.E., and Gillham, N.W. (1982). Extensive methylation of chloroplast DNA by a nuclear gene mutation does not affect chloroplast gene transmission in chlamydomonas. Cell 28, 335-343. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Boyle, N.R., Page, M.D., Liu, B., Blaby, I.K., Casero, D., Kropat, J., Cokus, S.J., Hong-Hermesdorf, A., Shaw, J., Karpowicz, S.J., Gallaher, S.D., Johnson, S., Benning, C., Pellegrini, M., Grossman, A., and Merchant, S.S. (2012). Three acyltransferases and nitrogen-responsive regulator are implicated in nitrogen starvation-induced triacylglycerol accumulation in Chlamydomonas. The Journal of biological chemistry 287, 15811-15825. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Brawley, S.H., and Johnson, L.E. (1992). Gametogenesis, Gametes and Zygotes - an Ecological Perspective on Sexual Reproduction in the Algae. Brit Phycol J 27, 233-252. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Brown, R.C., and Lemmon, B.E. (2011). Spores before sporophytes: hypothesizing the origin of sporogenesis at the algal-plant transition. The New phytologist 190, 875-881. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Burton, W.G., Grabowy, C.T., and Sager, R. (1979). Role of methylation in the modification and restriction of chloroplast DNA in Chlamydomonas. Proceedings of the National Academy of Sciences of the United States of America 76, 1390-1394. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Callebaut, I., Courvalin, J.C., and Mornon, J.P. (1999). The BAH (bromo-adjacent homology) domain: a link between DNA methylation, replication and transcriptional regulation. FEBS letters 446, 189-193. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Catt, J.W. (1979). Isolation and Chemical Composition of the Zygospore Cell-Wall of Chlamydomonas-Reinhardii. Plant Sci Lett 15, 69-74. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Cavalier-Smith, T. (1976). Electron microscopy of zygospore formation in Chlamydomonas reinhardii. Protoplasma 87, 297-315. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Cerutti, H., Johnson, A.M., Gillham, N.W., and Boynton, J.E. (1997). Epigenetic silencing of a foreign gene in nuclear transformants of Chlamydomonas. The Plant cell 9, 925-945. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title De Hoff, P.L., Ferris, P., Olson, B.J., Miyagi, A., Geng, S., and Umen, J.G. (2013). Species and population level molecular profiling reveals cryptic recombination and emergent asymmetry in the dimorphic mating locus of C. reinhardtii. PLoS genetics 9, e1003724. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title DePristo, M.A., Banks, E., Poplin, R., Garimella, K.V., Maguire, J.R., Hartl, C., Philippakis, A.A., del Angel, G., Rivas, M.A., Hanna, M., McKenna, A., Fennell, T.J., Kernytsky, A.M., Sivachenko, A.Y., Cibulskis, K., Gabriel, S.B., Altshuler, D., and Daly, M.J. (2011). A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nature genetics 43, 491-498. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Duncan, L., Nishii, I., Howard, A., Kirk, D., and Miller, S.M. (2006). Orthologs and paralogs of regA, a master cell-type regulatory gene in Volvox carteri. Current genetics 50, 61-72. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Duncan, L., Nishii, I., Harryman, A., Buckley, S., Howard, A., Friedman, N.R., and Miller, S.M. (2007). The VARL gene family and the evolutionary origins of the master cell-type regulatory gene, regA, in Volvox carteri. Journal of molecular evolution 65, 1-11. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Dyer, T.A. (1982). Methylation of chloroplast DNA in Chlamydomonas. Nature 298, 422-423. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Fang, W., Si, Y., Douglass, S., Casero, D., Merchant, S.S., Pellegrini, M., Ladunga, I., Liu, P., and Spalding, M.H. (2012). Transcriptome-wide changes in Chlamydomonas reinhardtii gene expression regulated by carbon dioxide and the CO2concentrating mechanism regulator CIA5/CCM1. The Plant cell 24, 1876-1893. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Feng, S., Cokus, S.J., Zhang, X., Chen, P.Y., Bostick, M., Goll, M.G., Hetzel, J., Jain, J., Strauss, S.H., Halpern, M.E., Ukomadu, C., Sadler, K.C., Pradhan, S., Pellegrini, M., and Jacobsen, S.E. (2010). Conservation and divergence of methylation patterning in plants and animals. Proceedings of the National Academy of Sciences of the United States of America 107, 8689-8694. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Feng, T.Y., and Chiang, K.S. (1984). The persistence of maternal inheritance in Chlamydomonas despite hypomethylation of chloroplast DNA induced by inhibitors. Proceedings of the National Academy of Sciences of the United States of America 81, 34383442. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Ferris, P.J., and Goodenough, U.W. (1987). Transcription of novel genes, including a gene linked to the mating-type locus, induced by Chlamydomonas fertilization. Molecular and cellular biology 7, 2360-2366. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Ferris, P.J., Armbrust, E.V., and Goodenough, U.W. (2002). Genetic structure of the mating-type locus of Chlamydomonas reinhardtii. Genetics 160, 181-200. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Ferris, P.J., Waffenschmidt, S., Umen, J.G., Lin, H., Lee, J.H., Ishida, K., Kubo, T., Lau, J., and Goodenough, U.W. (2005). Plus and minus sexual agglutinins from Chlamydomonas reinhardtii. The Plant cell 17, 597-615. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Gallaher, S.D., Fitz-Gibbon, S.T., Glaesener, A.G., Pellegrini, M., and Merchant, S.S. (2015). Chlamydomonas Genome Resource for Laboratory Strains Reveals a Mosaic of Sequence Variation, Identifies True Strain Histories, and Enables Strain-Specific Studies. The Plant cell. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Glaesener, A.G., Merchant, S.S., and Blaby-Haas, C.E. (2013). Iron economy in Chlamydomonas reinhardtii. Frontiers in plant science 4, 337. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Goldschmidt-Clermont, M. (1991). Transgenic expression of aminoglycoside adenine transferase in the chloroplast: a selectable marker of site-directed transformation of chlamydomonas. Nucleic acids research 19, 4083-4089. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Gonzalez-Ballester, D., Casero, D., Cokus, S., Pellegrini, M., Merchant, S.S., and Grossman, A.R. (2010). RNA-seq analysis of sulfurdeprived Chlamydomonas cells reveals aspects of acclimation critical for cell survival. The Plant cell 22, 2058-2084. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Goodenough, U., Lin, H., and Lee, J.H. (2007). Sex determination in Chlamydomonas. Seminars in cell & developmental biology 18, 350-361. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Grief, C., O'Neill, M.A., and Shaw, P.J. (1987). The zygote cell wall of Chlamydomonas reinhardii: a structural, chemical and immunological approach. Planta 170, 433-445. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Grossman, A. (2000). Acclimation of Chlamydomonas reinhardtii to its nutrient environment. Protist 151, 201-224. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Guo, W., Fiziev, P., Yan, W., Cokus, S., Sun, X., Zhang, M.Q., Chen, P.Y., and Pellegrini, M. (2013). BS-Seeker2: a versatile aligning pipeline for bisulfite sequencing data. BMC genomics 14, 774. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Harris, E.H., Stern, D.B., Witman, G., Harris, E.H., and Harris, E.H. (2009). The Chlamydomonas sourcebook. (Amsterdam ; Boston: Elsevier/Academic Press). Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Harris, M.A., Clark, J., Ireland, A., Lomax, J., Ashburner, M., Foulger, R., Eilbeck, K., Lewis, S., Marshall, B., Mungall, C., Richter, J., Rubin, G.M., Blake, J.A., Bult, C., Dolan, M., Drabkin, H., Eppig, J.T., Hill, D.P., Ni, L., Ringwald, M., Balakrishnan, R., Cherry, J.M., Christie, K.R., Costanzo, M.C., Dwight, S.S., Engel, S., Fisk, D.G., Hirschman, J.E., Hong, E.L., Nash, R.S., Sethuraman, A., Theesfeld, C.L., Botstein, D., Dolinski, K., Feierbach, B., Berardini, T., Mundodi, S., Rhee, S.Y., Apweiler, R., Barrell, D., Camon, E., Dimmer, E., Lee, V., Chisholm, R., Gaudet, P., Kibbe, W., Kishore, R., Schwarz, E.M., Sternberg, P., Gwinn, M., Hannick, L., Wortman, J., Berriman, M., Wood, V., de la Cruz, N., Tonellato, P., Jaiswal, P., Seigfried, T., White, R., and Gene Ontology, C. (2004). The Gene Ontology (GO) database and informatics resource. Nucleic acids research 32, D258-261. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Hattman, S., Kenny, C., Berger, L., and Pratt, K. (1978). Comparative study of DNA methylation in three unicellular eucaryotes. Journal of bacteriology 135, 1156-1157. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Heinnickel, M.L., and Grossman, A.R. (2013). The GreenCut: re-evaluation of physiological role of previously studied proteins and potential novel protein functions. Photosynthesis research 116, 427-436. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Idoine, A.D., Boulouis, A., Rupprecht, J., and Bock, R. (2014). The diurnal logic of the expression of the chloroplast genome in Chlamydomonas reinhardtii. PloS one 9, e108760. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Kanehisa, M., and Goto, S. (2000). KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic acids research 28, 27-30. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Karpowicz, S.J., Prochnik, S.E., Grossman, A.R., and Merchant, S.S. (2011). The GreenCut2 resource, a phylogenomically derived inventory of proteins specific to the plant lineage. The Journal of biological chemistry 286, 21427-21439. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Kirk, M.M., Stark, K., Miller, S.M., Muller, W., Taillon, B.E., Gruber, H., Schmitt, R., and Kirk, D.L. (1999). regA, a Volvox gene that plays a central role in germ-soma differentiation, encodes a novel regulatory protein. Development 126, 639-647. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Kubo, T., Abe, J., Oyamada, T., Ohnishi, M., Fukuzawa, H., Matsuda, Y., and Saito, T. (2008). Characterization of novel genes induced by sexual adhesion and gamete fusion and of their transcriptional regulation in Chlamydomonas reinhardtii. Plant & cell physiology 49, 981-993. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Kuriyama, H., Takano, H., Suzuki, L., Uchida, H., Kawano, S., Kuroiwa, H., and Kuroiwa, T. (1999). Characterization of Chlamydomonas reinhardtii zygote-specific cDNAs that encode novel proteins containing ankyrin repeats and WW domains. Plant physiology 119, 873-884. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Kurvari, V., Grishin, N.V., and Snell, W.J. (1998). A gamete-specific, sex-limited homeodomain protein in Chlamydomonas. The Journal of cell biology 143, 1971-1980. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Lechtreck, K.F., Rostmann, J., and Grunow, A. (2002). Analysis of Chlamydomonas SF-assemblin by GFP tagging and expression of antisense constructs. Journal of cell science 115, 1511-1522. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Lee, J.H., Lin, H., Joo, S., and Goodenough, U. (2008). Early sexual origins of homeoprotein heterodimerization and evolution of the plant KNOX/BELL family. Cell 133, 829-840. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Liu, Y., Misamore, M.J., and Snell, W.J. (2010). Membrane fusion triggers rapid degradation of two gamete-specific, fusionessential proteins in a membrane block to polygamy in Chlamydomonas. Development 137, 1473-1481. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Lopez, D., Casero, D., Cokus, S.J., Merchant, S.S., and Pellegrini, M. (2011). Algal Functional Annotation Tool: a web-based analysis suite to functionally interpret large gene lists using integrated annotation and expression data. BMC bioinformatics 12, 282. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Lv, H., Qu, G., Qi, X., Lu, L., Tian, C., and Ma, Y. (2013). Transcriptome analysis of Chlamydomonas reinhardtii during the process of lipid accumulation. Genomics 101, 229-237. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Matsuo, M., Hachisu, R., Tabata, S., Fukuzawa, H., and Obokata, J. (2011). Transcriptome analysis of respiration-responsive genes in Chlamydomonas reinhardtii: mitochondrial retrograde signaling coordinates the genes for cell proliferation with energyproducing metabolism. Plant & cell physiology 52, 333-343. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Matters, G.L., and Goodenough, U.W. (1992). A gene/pseudogene tandem duplication encodes a cysteine-rich protein expressed during zygote development in Chlamydomonas reinhardtii. Molecular & general genetics : MGG 232, 81-88. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Maul, J.E., Lilly, J.W., Cui, L., dePamphilis, C.W., Miller, W., Harris, E.H., and Stern, D.B. (2002). The Chlamydomonas reinhardtii plastid chromosome: islands of genes in a sea of repeats. The Plant cell 14, 2659-2679. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title McKenna, A., Hanna, M., Banks, E., Sivachenko, A., Cibulskis, K., Kernytsky, A., Garimella, K., Altshuler, D., Gabriel, S., Daly, M., and DePristo, M.A. (2010). The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome research 20, 1297-1303. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Merchant, S.S., Allen, M.D., Kropat, J., Moseley, J.L., Long, J.C., Tottey, S., and Terauchi, A.M. (2006). Between a rock and a hard place: trace element nutrition in Chlamydomonas. Biochimica et biophysica acta 1763, 578-594. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Merchant, S.S., Prochnik, S.E., Vallon, O., Harris, E.H., Karpowicz, S.J., Witman, G.B., Terry, A., Salamov, A., Fritz-Laylin, L.K., Marechal-Drouard, L., Marshall, W.F., Qu, L.H., Nelson, D.R., Sanderfoot, A.A., Spalding, M.H., Kapitonov, V.V., Ren, Q., Ferris, P., Lindquist, E., Shapiro, H., Lucas, S.M., Grimwood, J., Schmutz, J., Cardol, P., Cerutti, H., Chanfreau, G., Chen, C.L., Cognat, V., Croft, M.T., Dent, R., Dutcher, S., Fernandez, E., Fukuzawa, H., Gonzalez-Ballester, D., Gonzalez-Halphen, D., Hallmann, A., Hanikenne, M., Hippler, M., Inwood, W., Jabbari, K., Kalanon, M., Kuras, R., Lefebvre, P.A., Lemaire, S.D., Lobanov, A.V., Lohr, M., Manuell, A., Meier, I., Mets, L., Mittag, M., Mittelmeier, T., Moroney, J.V., Moseley, J., Napoli, C., Nedelcu, A.M., Niyogi, K., Novoselov, S.V., Paulsen, I.T., Pazour, G., Purton, S., Ral, J.P., Riano-Pachon, D.M., Riekhof, W., Rymarquis, L., Schroda, M., Stern, D., Umen, J., Willows, R., Wilson, N., Zimmer, S.L., Allmer, J., Balk, J., Bisova, K., Chen, C.J., Elias, M., Gendler, K., Hauser, C., Lamb, M.R., Ledford, H., Long, J.C., Minagawa, J., Page, M.D., Pan, J., Pootakham, W., Roje, S., Rose, A., Stahlberg, E., Terauchi, A.M., Yang, P., Ball, S., Bowler, C., Dieckmann, C.L., Gladyshev, V.N., Green, P., Jorgensen, R., Mayfield, S., Mueller-Roeber, B., Rajamani, S., Sayre, R.T., Brokstein, P., Dubchak, I., Goodstein, D., Hornick, L., Huang, Y.W., Jhaveri, J., Luo, Y., Martinez, D., Ngau, W.C., Otillar, B., Poliakov, A., Porter, A., Szajkowski, L., Werner, G., Zhou, K., Grigoriev, I.V., Rokhsar, D.S., and Grossman, A.R. (2007). The Chlamydomonas genome reveals the evolution of key animal and plant functions. Science 318, 245-250. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Miller, R., Wu, G., Deshpande, R.R., Vieler, A., Gartner, K., Li, X., Moellering, E.R., Zauner, S., Cornish, A.J., Liu, B., Bullard, B., Sears, B.B., Kuo, M.H., Hegg, E.L., Shachar-Hill, Y., Shiu, S.H., and Benning, C. (2010). Changes in transcript abundance in Chlamydomonas reinhardtii following nitrogen deprivation predict diversion of metabolism. Plant physiology 154, 1737-1752. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Miller, S.M., Schmitt, R., and Kirk, D.L. (1993). Jordan, an active Volvox transposable element similar to higher plant transposons. The Plant cell 5, 1125-1138. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Minami, S.A., and Goodenough, U.W. (1978). Novel glycopolypeptide synthesis induced by gametic cell fusion in Chlamydomonas reinhardtii. The Journal of cell biology 77, 165-181. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Nakamura, S. (2010). Paternal inheritance of mitochondria in Chlamydomonas. Journal of plant research 123, 163-170. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Nguyen, A.V., Thomas-Hall, S.R., Malnoe, A., Timmins, M., Mussgnug, J.H., Rupprecht, J., Kruse, O., Hankamer, B., and Schenk, P.M. (2008). Transcriptome for photobiological hydrogen production induced by sulfur deprivation in the green alga Chlamydomonas reinhardtii. Eukaryotic cell 7, 1965-1979. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Ning, J., Otto, T.D., Pfander, C., Schwach, F., Brochet, M., Bushell, E., Goulding, D., Sanders, M., Lefebvre, P.A., Pei, J., Grishin, N.V., Vanderlaan, G., Billker, O., and Snell, W.J. (2013). Comparative genomics in Chlamydomonas and Plasmodium identifies an ancient nuclear envelope protein family essential for sexual reproduction in protists, fungi, plants, and vertebrates. Genes & development 27, 1198-1215. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Nishimura, Y. (2010). Uniparental inheritance of cpDNA and the genetic control of sexual differentiation in Chlamydomonas reinhardtii. Journal of plant research 123, 149-162. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Nishimura, Y., Misumi, O., Kato, K., Inada, N., Higashiyama, T., Momoyama, Y., and Kuroiwa, T. (2002). An mt(+) gamete-specific nuclease that targets mt(-) chloroplasts during sexual reproduction in C. reinhardtii. Genes & development 16, 1116-1128. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Nishiyama, R., Ito, M., Yamaguchi, Y., Koizumi, N., and Sano, H. (2002). A chloroplast-resident DNA methyltransferase is responsible for hypermethylation of chloroplast genes in Chlamydomonas maternal gametes. Proceedings of the National Academy of Sciences of the United States of America 99, 5925-5930. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Nishiyama, R., Wada, Y., Mibu, M., Yamaguchi, Y., Shimogawara, K., and Sano, H. (2004). Role of a nonselective de novo DNA methyltransferase in maternal inheritance of chloroplast genes in the green alga, Chlamydomonas reinhardtii. Genetics 168, 809816. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Panchy, N., Wu, G., Newton, L., Tsai, C.H., Chen, J., Benning, C., Farre, E.M., and Shiu, S.H. (2014). Prevalence, Evolution, and cisRegulation of Diel Transcription in Chlamydomonas reinhardtii. G3 4, 2461-2471. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Prochnik, S.E., Umen, J., Nedelcu, A.M., Hallmann, A., Miller, S.M., Nishii, I., Ferris, P., Kuo, A., Mitros, T., Fritz-Laylin, L.K., Hellsten, U., Chapman, J., Simakov, O., Rensing, S.A., Terry, A., Pangilinan, J., Kapitonov, V., Jurka, J., Salamov, A., Shapiro, H., Schmutz, J., Grimwood, J., Lindquist, E., Lucas, S., Grigoriev, I.V., Schmitt, R., Kirk, D., and Rokhsar, D.S. (2010). Genomic analysis of organismal complexity in the multicellular green alga Volvox carteri. Science 329, 223-226. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Rochaix, J.D. (2001). Assembly,function, and dynamics of the photosynthetic machinery in Chlamydomonas reinhardtii. Plant physiology 127, 1394-1398. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Royer, H.D., and Sager, R. (1979). Methylation of chloroplast DNAs in the life cycle of Chlamydomonas. Proceedings of the National Academy of Sciences of the United States of America 76, 5794-5798. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Sano, H., Grabowy, C., and Sager, R. (1981). Differential activity of DNA methyltransferase in the life cycle of Chlamydomonas reinhardi. Proceedings of the National Academy of Sciences of the United States of America 78, 3118-3122. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Schmollinger, S., Muhlhaus, T., Boyle, N.R., Blaby, I.K., Casero, D., Mettler, T., Moseley, J.L., Kropat, J., Sommer, F., Strenkert, D., Hemme, D., Pellegrini, M., Grossman, A.R., Stitt, M., Schroda, M., and Merchant, S.S. (2014). Nitrogen-Sparing Mechanisms in Chlamydomonas Affect the Transcriptome, the Proteome, and Photosynthetic Metabolism. The Plant cell 26, 1410-1435. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Silflow, C.D., and Lefebvre, P.A. (2001). Assembly and motility of eukaryotic cilia and flagella. Lessons from Chlamydomonas reinhardtii. Plant physiology 127, 1500-1507. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Simon, D.F., Descombes, P., Zerges, W., and Wilkinson, K.J. (2008). Global expression profiling of Chlamydomonas reinhardtii exposed to trace levels of free cadmium. Environmental toxicology and chemistry / SETAC 27, 1668-1675. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Simon, D.F., Domingos, R.F., Hauser, C., Hutchins, C.M., Zerges, W., and Wilkinson, K.J. (2013). Transcriptome sequencing (RNAseq) analysis of the effects of metal nanoparticle exposure on the transcriptome of Chlamydomonas reinhardtii. Applied and environmental microbiology 79, 4774-4785. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Smith, Z.D., and Meissner, A. (2013). DNA methylation: roles in mammalian development. Nature reviews. Genetics 14, 204-220. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Suzuki, L., Woessner, J.P., Uchida, H., Kuroiwa, H., Yuasa, Y., Waffenschmidt, S., Goodenough, U.W., and Kuroiwa, T. (2000). Zygote-specific protein with hydroxyproline-rich glycoprotein domains and lectin-like domains involved in the assembly of the cell wall of Chlamydomonas reinhardtii (Chlorophyta). J Phycol 36, 571-583. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Swanson, W.J., and Vacquier, V.D. (2002). The rapid evolution of reproductive proteins. Nature reviews. Genetics 3, 137-144. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Tam, L.W., Ranum, P.T., and Lefebvre, P.A. (2013). CDKL5 regulates flagellar length and localizes to the base of the flagella in Chlamydomonas. Molecular biology of the cell 24, 588-600. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Thimm, O., Blasing, O., Gibon, Y., Nagel, A., Meyer, S., Kruger, P., Selbig, J., Muller, L.A., Rhee, S.Y., and Stitt, M. (2004). MAPMAN: a user-driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes. The Plant journal : for cell and molecular biology 37, 914-939. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Thorvaldsdottir, H., Robinson, J.T., and Mesirov, J.P. (2013). Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration. Briefings in bioinformatics 14, 178-192. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Toepel, J., Albaum, S.P., Arvidsson, S., Goesmann, A., la Russa, M., Rogge, K., and Kruse, O. (2011). Construction and evaluation of a whole genome microarray of Chlamydomonas reinhardtii. BMC genomics 12, 579. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Toepel, J., Illmer-Kephalides, M., Jaenicke, S., Straube, J., May, P., Goesmann, A., and Kruse, O. (2013). New insights into Chlamydomonas reinhardtii hydrogen production processes by combined microarray/RNA-seq transcriptomics. Plant biotechnology journal 11, 717-733. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Uchida, H., Kawano, S., Sato, N., and Kuroiwa, T. (1993). Isolation and characterization of novel genes which are expressed during the very early stage of zygote formation in Chlamydomonas reinhardtii. Current genetics 24, 296-300. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Uchida, H., Suzuki, L., Anai, T., Doi, K., Takano, H., Yamashita, H., Oka, T., Kawano, S., Tomizawa, K.I., Kawazu, T., Kuroiwa, H., and Kuroiwa, T. (1999). A pair of invertedly repeated genes in Chlamydomonas reinhardtii encodes a zygote-specific protein whose expression is UV-sensitive. Current genetics 36, 232-240. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Ueki, N., and Nishii, I. (2008). Idaten is a new cold-inducible transposon of Volvox carteri that can be used for tagging developmentally important genes. Genetics 180, 1343-1353. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Umen, J.G. (2011). Evolution of sex and mating loci: an expanded view from Volvocine algae. Current opinion in microbiology 14, 634-641. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Umen, J.G., and Goodenough, U.W. (2001). Chloroplast DNA methylation and inheritance in Chlamydomonas. Genes & development 15, 2585-2597. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Woessner, J.P., Molendijk, A.J., van Egmond, P., Klis, F.M., Goodenough, U.W., and Haring, M.A. (1994). Domain conservation in several volvocalean cell wall proteins. Plant molecular biology 26, 947-960. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Yang, N., and Xu, R.M. (2013). Structure and function of the BAH domain in chromatin biology. Critical reviews in biochemistry and molecular biology 48, 211-221. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title
© Copyright 2026 Paperzz