Sequencing and comparative analysis of three Chlorella genomes provide insights into strain-specific adaptation to wastewater.

Microalgal Chlorella has been demonstrated to process wastewater efficiently from piggery industry, yet optimization through genetic engineering of such a bio-treatment is currently challenging, largely due to the limited data and knowledge in genomics. In this study, we first investigated the differential growth rates among three wastewater-processing Chlorella strains: Chlorella sorokiniana BD09, Chlorella sorokiniana BD08 and Chlorella sp. Dachan, and the previously published Chlorella sorokiniana UTEX 1602, showing us that BD09 maintains the best tolerance in synthetic wastewater. We then performed genome sequencing and analysis, resulting in a high-quality assembly for each genome with scaffold N50 > 2 Mb and genomic completeness ≥91%, as well as genome annotation with 9,668, 10,240, 9,821 high-confidence gene models predicted for BD09, BD08, and Dachan, respectively. Comparative genomics study unravels that metabolic pathways, which are involved in nitrogen and phosphorus assimilation, were enriched in the faster-growing strains. We found that gene structural variation and genomic rearrangement might contribute to differential capabilities in wastewater tolerance among the strains, as indicated by gene copy number variation, domain reshuffling of orthologs involved, as well as a ~1 Mb-length chromosomal inversion we observed in BD08 and Dachan. In addition, we speculated that an associated bacterium, Microbacterium chocolatum, which was identified within Dachan, play a possible role in synergizing nutrient removal. Our three newly sequenced Chlorella genomes provide a fundamental foundation to understand the molecular basis of abiotic stress tolerance in wastewater treatment, which is essential for future genetic engineering and strain improvement.

Efficient coralline algal psbA mini barcoding and High Resolution Melt (HRM) analysis using a simple custom DNA preparation.

Coralline algae form extensive maerl and rhodolith habitats that support a rich biodiversity. Calcium carbonate harvesting as well as trawling activities threatens this ecosystem. Eleven species were recorded so far as maerl-forming in NE Atlantic, but identification based on morphological characters is unreliable. As for most red algae, we now use molecular characters to resolve identification of these taxa. However, obtaining DNA sequences requires time and resource demanding methods. The purpose of our study was to improve methods for achieving simple DNA extraction, amplification, sequencing and sequence analysis to allow robust identification of maerl species and other coralline algae. Our novel and easy DNA preparation method for coralline algae was of sufficient quality for qPCR amplification and sequencing of all 47 tested samples. The new psbA qPCR assay successfully amplified a 350 bp fragment identifying six species and uncovering two new Operational Taxonomic Units. Molecular results were corroborated with anatomical examination using i.e. scanning electron microscopy. Finally, the qPCR assay was coupled with High Resolution Melt analysis that successfully differentiated the closely related species Lithothamnion erinaceum and L. cf. glaciale. This DNA preparation and qPCR technique should vitalize coralline research by reducing time and cost associated with molecular systematics.

Rapid and high efficiency transformation of Chlamydomonas reinhardtii by square-wave electroporation.

Chlamydomonas reinhardtii, the unicellular green algae, is the model organism for studies in various physiological processes and for bioindustrial applications. To explore the molecular mechanisms underlying physiological processes or to establish engineered cell lines, the exogenous DNA needs to be integrated into the genome for the insertional mutagenesis or transgene expression. However, the amount of selected marker DNA is not seriously considered in the existing electroporation methods for mutants library construction. Here, we reported a rapid-and-high-efficiency transformation technique for cell-walled strains using square-wave electroporation system. The final yield with this electroporation method was 2-6 × 103 transformants per µg exogenous DNA for cell-walled strains in a strain-dependent manner. In general, this electroporation technique was the easy and applicable way to build a mutant library for screening phenotypes of interest.

Seasonal succession of ciliate Mesodinium spp. with red, green, or mixed plastids and their association with cryptophyte prey.

Mesodinium spp. are commonly found in marine and brackish waters, and several species are known to contain red, green, or both plastids that originate from cryptophyte prey. We observed the seasonal succession of Mesodinium spp. in a Japanese brackish lake, and we analysed the origin and diversity of the various coloured plastids within the cells of Mesodinium spp. using a newly developed primer set that specifically targets the cryptophyte nuclear 18S rRNA gene. Mesodinium rubrum isolated from the lake contained only red plastids originating from cryptophyte Teleaulax amphioxeia. We identified novel Mesodinium sp. that contained only green plastids or both red and green plastids originating from cryptophytes Hemiselmis sp. and Teleaulax acuta. Although the morphology of the newly identified Mesodinium sp. was indistinguishable from that of M. rubrum under normal light microscopy, phylogenetic analysis placed this species between the M. rubrum/major species complex and a well-supported lineage of M. chamaeleon and M. coatsi. Close associations were observed in cryptophyte species composition within cells of Mesodinium spp. and in ambient water samples. The appearance of suitable cryptophyte prey is probably a trigger for succession of Mesodinium spp., and the subsequent abundance of Mesodinium spp. appears to be influenced by water temperature and dissolved inorganic nutrients.

Phylogenetic Analysis Reveals that the 'Radial Centric' Diatom Orthoseira Thwaites (Orthoseiraceae, Bacillariophyta) is a Member of a 'Multipolar' Diatom Lineage.

The diatom genus Orthoseira Thwaites (Bacillariophyta) is a ubiquitous taxon in aerial diatom assemblages, with species found globally. Cylindrical cell shape and radial symmetry of this genus has led to its historical placement in the Coscinodiscophyceae ('radial centric' diatoms), but its systematic relationships have remained uncertain. We present a five-gene phylogeny, based on nuclear (nSSU rDNA) and chloroplast (rbcL, psbC, psbA, and psaB) genes to determine the phylogenetic placement of Orthoseira among the diatoms. The concatenated multi-gene phylogenies and nSSU-only gene tree demonstrate that Orthoseira is deeply embedded within a clade of the Mediophyceae ('multipolar centric' diatoms). Throughout all phylogenetic analyses, Orthoseira was shown to be sister to the genera Terpsinoë and Hydrosera. Through comparison of topologies reflecting competing hypotheses about the placement of Orthoseira, it was determined that the hypothesis that Orthoseira, represented here by O. dendroteres and O. roeseana, is a member of the Melosirales should be rejected. Therefore, lack of morphological similarity between Hydrosera, Orthoseira, and Terpsinoë is hypothesized to be the result of changes in habitat preferences that lead to an ancient divergence event between the Orthoseirales and the Hydrosera, Terpsinoë clade.

The internal transcribed spacer 2 of Jenufa (Chlorophyta, Chlorophyceae) is extraordinarily long: A hypothesis.

Aided by a host of bioinformatics tools, primary and secondary structural analyses of the internal transcribed spacer 2 (ITS2) from the eukaryotic ribosomal RNA repeat have a long and enviable record of service to diversity studies of fungi, plants and protists. Automation of annotation, secondary structure estimation and sequence alignment have become routine for the vast majority of ITS2 sequences. Challenges to the bioinformatics pipeline for ITS2 analysis generally arise in cases where the sequence length lies well outside the norm. These sequences generally defy protocols for annotation and secondary structure prediction. The long ITS2 sequences (ca. 600 nucleotides) from the green alga, Jenufa, offered an opportunity to explore this problem. Custom BLAST parameters revealed the presence of 4-helix structures (200-250 nucleotides) embedded in the 5' portion of several long ITS2 sequences of Jenufa. Of special note is the ITS2 sequence of J. lobulosa where a 4-helix structure was obtained for both the embedded ITS2 and for the complete sequence. Phylogenetic analysis of these typically-sized sequences resolved Golenkinia longispicula as the sister to Jenufa. Our observations indicate that other long ITS2 sequences should be examined for evidence of expansion or duplication. In addition, if the embedded ITS2 sequences are functional, then ribogenesis is almost certainly more diverse than is already apparent from studies of humans and yeast.

A Case of Prototheca zopfii Genotype 1 Infection in a Dog (Canis lupus familiaris).

Protothecosis is a rare disease caused by environmental algae of the genus Prototheca. These are saprophytic, non-photosynthetic, aerobic, colorless algae that belong to the Chlorellaceae family. Seven different species have been described. Prototheca zopfii genotype 2 and P. wickerhamii are most commonly involved in pathogenic infections in humans and animals. The objective of this work is to describe, for the first time, a case of protothecosis caused by P. zopfii genotype 1 in a dog. The dog, a 4-year-old mix bred male, was presented to a veterinary clinic in Montevideo, Uruguay, with multiple skin nodules, one of which was excised by surgical biopsy. The sample was examined histologically and processed by PCR, DNA sequencing, and restriction fragments length polymorphisms for the detection and genotyping of P. zopfii. In addition, transmission electron microscopy and scanning electron microscopy were performed. Histology showed severe ulcerative granulomatous dermatitis and panniculitis with myriads of pleomorphic algae. Algal cells were 4-17 µm in size, with an amphophilic, 2-4-µm-thick wall frequently surrounded by a clear halo, contained flocculant material and a deeply basophilic nucleus, and internal septae with daughter cells (endospores) consistent with endosporulation. Ultrastructurally, algal cells/endospores at different stages of development were found within parasitophorous vacuoles in macrophages. Prototheca zopfii genotype 1 was identified by molecular testing, confirming the etiologic diagnosis of protothecosis.

Chloropicophyceae, a new class of picophytoplanktonic prasinophytes.

Prasinophytes are a paraphyletic group of nine lineages of green microalgae that are currently classified either at the class or order level or as clades without formal taxonomic description. Prasinophyte clade VII comprises picoplanktonic algae that are important components of marine phytoplankton communities, particularly in moderately oligotrophic waters. Despite first being cultured in the 1960s, this clade has yet to be formally described. Previous phylogenetic analyses using the 18S rRNA gene divided prasinophyte clade VII into three lineages, termed A, B and C, the latter formed by a single species, Picocystis salinarum, that to date has only been found in saline lakes. Strains from lineages A and B cannot be distinguished by light microscopy and have very similar photosynthetic pigment profiles corresponding to the prasino-2A pigment group. We obtained phenotypic and genetic data on a large set of prasinophyte clade VII culture strains that allowed us to clarify the taxonomy of this important marine group. We describe two novel classes, the Picocystophyceae and the Chloropicophyceae, the latter containing two novel genera, Chloropicon and Chloroparvula, and eight new species of marine picoplanktonic green algae.

Ecology and pathology of Phytophthora ITS clade 3 species in forests in western Oregon, USA.

Phytophthora species are widespread and diverse in forest ecosystems, but little is known about their ecology. We explore ecological attributes of the closely related clade 3 species that occur sympatrically in western North American forests. We address the population structure, pathology, and epidemiology of P. ilicis, P. nemorosa, P. pluvialis, P. pseudosyringae, and P. psychrophila. Phytophthora species were isolated from plant tissues, rainwater falling through the forest canopy, streams, and soils in forests in western Oregon. Species identifications were based on morphology in culture with molecular confirmation using COX spacer and internal transcribed spacer (ITS) sequences. All five clade 3 Phytophthora species are present in western Oregon forests, although P. ilicis (only 1 forest isolate) and P. psychrophila (only 12 isolates) are apparently rare. P. ilicis is known only from holly in horticultural situations and once from a naturalized seedling in an urban forest. The known distribution of P. nemorosa in forest settings coincides with the ranges of its principle hosts, tanoak and myrtlewood, in Oregon and California. Although it is regularly identified from streams within the tanoak range, it has not been recovered from streams beyond that range. P. pluvialis is primarily associated with Douglas-fir canopies. It was identified from scattered locations throughout western Oregon in rain traps beneath Douglas-fir plantations and from diseased needles. P. pseudosyringae is also isolated from tanoak and myrtlewood in southwest Oregon and California, but its distribution, in streams at least, extends throughout much of western Oregon. P. psychrophila in Oregon is known only from rain traps beneath tanoak trees. Little intraspecific variation was detected in the nuclear rDNA ITS of clade 3 species. Variation in the mitochondrial COX spacer region was more frequent, with 2 to 10 haplotypes identified in the clade 3 species, for which we had multiple isolates.

Use of Highly Specific Molecular Markers Reveals Positive Correlation between Abundances of Mesodinium cf. major and Its Preferred Prey, Teleaulax amphioxeia, During Red Water Blooms in the Columbia River Estuary.

In a previous study, Teleaulax amphioxeia-the preferred prey of Mesodinium in the Columbia River estuary-were undetectable within intense annual blooms, suggesting blooms are prey-limited or prey are acquired outside of bloom patches. We used a novel molecular approach specifically targeting the prey (i.e., Unique Sequence Element [USE] within the ribosomal RNA 28S D2 regions of T. amphioxeia nucleus and nucleomorph) in estuarine water samples acquired autonomously with an Environmental Sample Processor integrated within a monitoring network (ESP-SATURN). This new approach allowed for both more specific detection of the prey and better constraint of sample variability. A positive correlation was observed between abundances of M. cf. major and T. amphioxeia during bloom periods. The correlation was stronger at depth (> 8.2 m) and weak or nonexistent in the surface, suggesting that predator-prey dynamics become uncoupled when stratification is strong. We confirmed exclusive selectivity for T. amphioxeia by M. cf. major and observed the incorporation of the prey nucleus into a 4-nuclei complex, where it remained functionally active. The specific biomarker for T. amphioxeia was also recovered in M. cf. major samples from a Namibian coastal bloom, suggesting that a specific predator-prey relationship might be widespread between M. cf. major and T. amphioxeia.

Microbialite response to an anthropogenic salinity gradient in Great Salt Lake, Utah.

A railroad causeway across Great Salt Lake, Utah (GSL), has restricted water flow since its construction in 1959, resulting in a more saline North Arm (NA; 24%-31% salinity) and a less saline South Arm (SA; 11%-14% salinity). Here, we characterized microbial carbonates collected from the SA and the NA to evaluate the effect of increased salinity on community composition and abundance and to determine whether the communities present in the NA are still actively precipitating carbonate or if they are remnant features from prior to causeway construction. SSU rRNA gene abundances associated with the NA microbialite were three orders of magnitude lower than those associated with the SA microbialite, indicating that the latter community is more productive. SSU rRNA gene sequencing and functional gene microarray analyses indicated that SA and NA microbialite communities are distinct. In particular, abundant sequences affiliated with photoautotrophic taxa including cyanobacteria and diatoms that may drive carbonate precipitation and thus still actively form microbialites were identified in the SA microbialite; sequences affiliated with photoautotrophic taxa were in low abundance in the NA microbialite. SA and NA microbialites comprise smooth prismatic aragonite crystals. However, the SA microbialite also contained micritic aragonite, which can be formed as a result of biological activity. Collectively, these observations suggest that NA microbialites are likely to be remnant features from prior to causeway construction and indicate a strong decrease in the ability of NA microbialite communities to actively precipitate carbonate minerals. Moreover, the results suggest a role for cyanobacteria and diatoms in carbonate precipitation and microbialite formation in the SA of GSL.

Primary productivity of snow algae communities on stratovolcanoes of the Pacific Northwest.

The majority of geomicrobiological research conducted on glacial systems to date has focused on glaciers that override primarily carbonate or granitic bedrock types, with little known of the processes that support microbial life in glacial systems overriding volcanic terrains (e.g., basalt or andesite). To better constrain the role of the supraglacial ecosystems in the carbon and nitrogen cycles, to gain insight into microbiome composition and function in alpine glacial systems overriding volcanic terrains, and to constrain potential elemental sequestration or release through weathering processes associated with snow algae communities, we examined the microbial community structure and primary productivity of snow algae communities on stratovolcanoes in the Cascade Range of the Pacific Northwest. Here, we present the first published values for carbon fixation rates of snow algae communities on glaciers in the Pacific Northwest. We observed varying levels of light-dependent carbon fixation on supraglacial and periglacial snowfields at Mt. Hood, Mt. Adams, and North Sister. Recovery of abundant 18S rRNA transcripts affiliated with photoautotrophs and 16S rRNA transcripts affiliated with heterotrophic bacteria is consistent with previous studies indicating the majority of primary productivity on snow and ice can be attributed to photoautotrophs. In contrast to previous observations of glacial ecosystems, our geochemical, isotopic, and microcosm data suggest these assemblages are not limited by phosphorus or fixed nitrogen availability. Furthermore, our data indicate these snow algae communities actively sequester Fe, Mn, and P leached from minerals sourced from the local rocks. Our observations of light-dependent primary productivity on snow are consistent with similar studies in polar ecosystems; however, our data may suggest that DIC may be a limiting nutrient in contrast to phosphorus or fixed nitrogen as has been observed in other glacial ecosystems. Our data underscore the need for similar studies on glacier surfaces and seasonal snowfields to better constrain the role of local bedrock and nutrient delivery on carbon fixation and biogeochemical cycling in these ecosystems.

Species concept and nomenclatural changes within the genera Elliptochloris and Pseudochlorella (Trebouxiophyceae) based on an integrative approach.

The genera Elliptochloris and Pseudochlorella were erected for Chlorella-like green algae producing two types of autospores and cell packages, respectively. Both genera are widely distributed in different soil habitats, either as free living or as photobionts of lichens. The species of these genera are often difficult to identify because of the high phenotypic plasticity and occasional lack of characteristic features. The taxonomic and nomenclatural status of these species, therefore, remains unclear. In this study, 34 strains were investigated using an integrative approach. Phylogenetic analyses demonstrated that the isolates belong to two independent lineages of the Trebouxiophyceae (Elliptochloris and Prasiola clades) and confirmed that the genera are not closely related. The comparison of morphology, molecular phylogeny, and analyses of secondary structures of SSU and ITS rDNA sequences revealed that all of the strains belong to three genera: Elliptochloris, Pseudochlorella, and Edaphochlorella. As a consequence of the taxonomic revisions, we propose two new combinations (Elliptochloris antarctica and Pseudochlorella signiensis) and validate Elliptochloris reniformis, which is invalidly described according to the International Code for Nomenclature (ICN), by designating a holotype. To reflect the high phenotypic plasticity of P. signiensis, two new varieties were described: P. signiensis var. magna and P. signiensis var. communis. Chlorella mirabilis was not closely related to any of these genera and was, therefore, transferred to the new genus Edaphochlorella. All of the taxonomic changes were highly supported by all phylogenetic analyses and were confirmed by the ITS-2 Barcodes using the ITS-2/CBC approach.

Subtype Distribution of Blastocystis in Communities along the Chao Phraya River, Thailand.

Blastocystis is a common zoonotic enteric protozoan that has been classified into 17 distinct subtypes (STs). A cross-sectional study was conducted to determine the prevalence and subtype distributions of Blastocystis in villagers living along the Chao Phraya River, Ayutthaya Province, Thailand, and to assess the risk of zoonotic infection. In total, 220 stool samples were collected, and DNA was extracted. PCR and sequencing were performed with primers targeting the small-subunit ribosomal RNA (SSU rRNA) genes. Blastocystis was present in 5.9% (13/220) of samples, and ST3 (5.0%; 11/220) was the predominant subtype, followed by ST2 (0.45%; 1/220) and ST6 (0.45%; 1/220). Phylogenetic trees were constructed with the maximum-likelihood method based on the Hasegawa-Kishino-Yano + G + I model, neighbor-joining, and maximum parsimony methods. The percentage of bootstrapped trees in which the associated taxa clustered together was relatively high. All the sequences of the Blastocystis-positive samples (KU051524-KU051536) were closely related to those from animals (pig, cattle, and chicken), indicating a zoonotic risk. Therefore, the villagers require proper health education, especially regarding the prevention of parasitic infection, to improve their personal hygiene and community health. Further studies are required to investigate the Blastocystis STs in the animals living in these villages.

Phylogenetic diversity and biogeography of the Mamiellophyceae lineage of eukaryotic phytoplankton across the oceans.

High-throughput diversity amplicon sequencing of marine microbial samples has revealed that members of the Mamiellophyceae lineage are successful phytoplankton in many oceanic habitats. Indeed, these eukaryotic green algae can dominate the picoplanktonic biomass, however, given the broad expanses of the oceans, their geographical distributions and the phylogenetic diversity of some groups remain poorly characterized. As these algae play a foundational role in marine food webs, it is crucial to assess their global distribution in order to better predict potential changes in abundance and community structure. To this end, we analyzed the V9-18S small subunit rDNA sequences deposited from the Tara Oceans expedition to evaluate the diversity and biogeography of these phytoplankton. Our results show that the phylogenetic composition of Mamiellophyceae communities is in part determined by geographical provenance, and do not appear to be influenced - in the samples recovered - by water depth, at least at the resolution possible with the V9-18S. Phylogenetic classification of Mamiellophyceae sequences revealed that the Dolichomastigales order encompasses more sequence diversity than other orders in this lineage. These results indicate that a large fraction of the Mamiellophyceae diversity has been hitherto overlooked, likely because of a combination of size fraction, sequencing and geographical limitations.

Complete genome sequence of mitochondrial DNA (mtDNA) of Chlorella sorokiniana.

The complete sequence of mitochondrial genome of the Chlorella sorokiniana strain (SAG 111-8 k) is presented in this work. Within the Chlorella genus, it represents the second species with a complete sequenced and annotated mitochondrial genome (GenBank accession no. KM241869). The genome consists of circular chromosomes of 52,528 bp and encodes a total of 31 protein coding genes, 3 rRNAs and 26 tRNAs. The overall AT contents of the C. sorokiniana mtDNA is 70.89%, while the coding sequence is of 97.4%.

Delineating a New Heterothallic Species of Volvox (Volvocaceae, Chlorophyceae) Using New Strains of "Volvox africanus".

The volvocine algae represent an excellent model lineage in which to study evolution of female and male genders based on comparative analyses of related species. Among these species, Volvox carteri has been extensively studied as a model of an oogamous and complex organism. However, it may have unique derived features that are not present in other species of Volvox. Therefore, information regarding the characteristics of sexual reproduction of other species of Volvox is also important. In 1971, Starr studied four types of sexuality in several global strains identified as Volvox africanus; however, further taxonomic studies of these strains have been lacking, and strains of three of the four sexual types are not available. Here, we studied the morphology, sexual reproduction, and taxonomy of two V. africanus-like species isolated recently from Lake Biwa, Japan. These two species were very similar to two sexual types described by Starr in 1971: one producing dioecious sexual spheroids in heterothallic strains and the other forming both male spheroids and monoecious spheroids in a single strain. The former species produced zygotes with a reticulate cell wall, whereas a smooth zygote wall was observed in the latter species as in V. africanus previously reported from various localities around the world. Our multigene phylogenetic analysis demonstrated that these are sister species to each other. However, the presence of a compensatory base change in the most conserved region of the secondary structure of nuclear ribosomal DNA internal transcribed spacer-2, hybrid inviability demonstrated by intercrossing experiments, and morphological differences in the density of abutment between the gelatinous material of adjacent cells (individual sheaths) in the spheroid supported the recognition of the two species, V. africanus having a smooth zygote wall and V. reticuliferus Nozaki sp. nov. having a reticulate zygote wall.

Distribution patterns of haplotypes for symbionts from Umbilicaria esculenta and U. muehlenbergii reflect the importance of reproductive strategy in shaping population genetic structure.

BACKGROUND: The diversity of lichen fungal components and their photosynthetic partners reflects both ecological and evolutionary factors. In present study, molecular investigations of the internal transcribed spacer of the nuclear ribosomal DNA (ITS nrDNA) region were conducted to analyze the genetic diversity of Umbilicaria esculenta and U. muehlenbergii together with their associated green algae. RESULT: It was here demonstrated that the reproductive strategy is a principal reason for fungal selectivity to algae. U. muehlenbergii, which disperses via sexual spores, exhibits lower selectivity to its photosynthetic partners than U. esculenta, which has a vegetative reproductive strategy. The difference of genotypic diversity (both fungal and algal) between these two Umbilicaria species is low, although their nucleotide diversity can vary greatly. CONCLUSIONS: The present study illustrates that lichen-forming fungi with sexual reproductive strategies are less selective with respect to their photobionts; and reveals that both sexual and vegetative reproduction allow lichens to generate similar amounts of diversity to adapt to the environments. The current study will be helpful for elucidating how lichens with different reproductive strategies adapt to changing environments.

Dynamic Changes in the Transcriptome and Methylome of Chlamydomonas reinhardtii throughout Its Life Cycle.

The green alga Chlamydomonas reinhardtii undergoes gametogenesis and mating upon nitrogen starvation. While the steps involved in its sexual reproductive cycle have been extensively characterized, the genome-wide transcriptional and epigenetic changes underlying different life cycle stages have yet to be fully described. Here, we performed transcriptome and methylome sequencing to quantify expression and DNA methylation from vegetative and gametic cells of each mating type and from zygotes. We identified 361 gametic genes with mating type-specific expression patterns and 627 genes that are specifically induced in zygotes; furthermore, these sex-related gene sets were enriched for secretory pathway and alga-specific genes. We also examined the C. reinhardtii nuclear methylation map with base-level resolution at different life cycle stages. Despite having low global levels of nuclear methylation, we detected 23 hypermethylated loci in gene-poor, repeat-rich regions. We observed mating type-specific differences in chloroplast DNA methylation levels in plus versus minus mating type gametes followed by chloroplast DNA hypermethylation in zygotes. Lastly, we examined the expression of candidate DNA methyltransferases and found three, DMT1a, DMT1b, and DMT4, that are differentially expressed during the life cycle and are candidate DNA methylases. The expression and methylation data we present provide insight into cell type-specific transcriptional and epigenetic programs during key stages of the C. reinhardtii life cycle.