Genome-wide assessment of genetic diversity and transcript variations in 17 accessions of the model diatom Phaeodactylum tricornutum.

Diatoms, a prominent group of phytoplankton, have a significant impact on both the oceanic food chain and carbon sequestration, thereby playing a crucial role in regulating the climate. These highly diverse organisms show a wide geographic distribution across various latitudes. In addition to their ecological significance, diatoms represent a vital source of bioactive compounds that are widely used in biotechnology applications. In the present study, we investigated the genetic and transcriptomic diversity of 17 accessions of the model diatom Phaeodactylum tricornutum including those sampled a century ago as well as more recently collected accessions. The analysis of the data reveals a higher genetic diversity and the emergence of novel clades, indicating an increasing diversity within the P. tricornutum population structure, compared to the previous study and a persistent long-term balancing selection of genes in old and newly sampled accessions. However, the study did not establish a clear link between the year of sampling and genetic diversity, thereby, rejecting the hypothesis of loss of heterozygoty in cultured strains. Transcript analysis identified novel transcript including noncoding RNA and other categories of small RNA such as PiwiRNAs. Additionally, transcripts analysis using differential expression as well as Weighted Gene Correlation Network Analysis has provided evidence that the suppression or downregulation of genes cannot be solely attributed to loss-of-function mutations. This implies that other contributing factors, such as epigenetic modifications, may play a crucial role in regulating gene expression. Our study provides novel genetic resources, which are now accessible through the platform PhaeoEpiview (https://PhaeoEpiView.univ-nantes.fr), that offer both ease of use and advanced tools to further investigate microalgae biology and ecology, consequently enriching our current understanding of these organisms.

The Polycomb repressive complex 2 deposits H3K27me3 and represses transposable elements in a broad range of eukaryotes.

The mobility of transposable elements (TEs) contributes to evolution of genomes. Their uncontrolled activity causes genomic instability; therefore, expression of TEs is silenced by host genomes. TEs are marked with DNA and H3K9 methylation, which are associated with silencing in flowering plants, animals, and fungi. However, in distantly related groups of eukaryotes, TEs are marked by H3K27me3 deposited by the Polycomb repressive complex 2 (PRC2), an epigenetic mark associated with gene silencing in flowering plants and animals. The direct silencing of TEs by PRC2 has so far only been shown in one species of ciliates. To test if PRC2 silences TEs in a broader range of eukaryotes, we generated mutants with reduced PRC2 activity and analyzed the role of PRC2 in extant species along the lineage of Archaeplastida and in the diatom P. tricornutum. In this diatom and the red alga C. merolae, a greater proportion of TEs than genes were repressed by PRC2, whereas a greater proportion of genes than TEs were repressed by PRC2 in bryophytes. In flowering plants, TEs contained potential cis-elements recognized by transcription factors and associated with neighbor genes as transcriptional units repressed by PRC2. Thus, silencing of TEs by PRC2 is observed not only in Archaeplastida but also in diatoms and ciliates, suggesting that PRC2 deposited H3K27me3 to silence TEs in the last common ancestor of eukaryotes. We hypothesize that during the evolution of Archaeplastida, TE fragments marked with H3K27me3 were selected to shape transcriptional regulation, controlling networks of genes regulated by PRC2.

Chromosome-Wide Distribution and Characterization of H3K36me3 and H3K27Ac in the Marine Model Diatom Phaeodactylum tricornutum.

Histone methylation and acetylation play a crucial role in response to developmental cues and environmental changes. Previously, we employed mass spectrometry to identify histone modifications such as H3K27ac and H3K36me3 in the model diatom Phaeodactylum tricornutum, which have been shown to be important for transcriptional activation in animal and plant species. To further investigate their evolutionary implications, we utilized chromatin immunoprecipitation followed by deep sequencing (ChIP-Seq) and explored their genome-wide distribution in P. tricornutum. Our study aimed to determine their role in transcriptional regulation of genes and transposable elements (TEs) and their co-occurrence with other histone marks. Our results revealed that H3K27ac and H3K36me3 were predominantly localized in promoters and genic regions indicating a high conservation pattern with studies of the same marks in plants and animals. Furthermore, we report the diversity of genes encoding H3 lysine 36 (H3K36) trimethylation-specific methyltransferase in microalgae leveraging diverse sequencing resources including the Marine Microbial Eukaryote Transcriptome Sequencing Project database (MMETSP). Our study expands the repertoire of epigenetic marks in a model microalga and provides valuable insights into the evolutionary context of epigenetic-mediated gene regulation. These findings shed light on the intricate interplay between histone modifications and gene expression in microalgae, contributing to our understanding of the broader epigenetic landscape in eukaryotic organisms.

PhaeoEpiView: an epigenome browser of the newly assembled genome of the model diatom Phaeodactylum tricornutum.

Recent advances in DNA sequencing technologies particularly long-read sequencing, greatly improved genomes assembly. However, this has created discrepancies between published annotations and epigenome tracks, which have not been updated to keep pace with the new assemblies. Here, we used the latest improved telomere-to-telomere assembly of the model pennate diatom Phaeodactylum tricornutum to lift over the gene models from Phatr3, a previously annotated reference genome. We used the lifted genes annotation and newly published transposable elements to map the epigenome landscape, namely DNA methylation and post-translational modifications of histones. This provides the community with PhaeoEpiView, a browser that allows the visualization of epigenome data and transcripts on an updated and contiguous reference genome, to better understand the biological significance of the mapped data. We updated previously published histone marks with a more accurate peak calling using mono instead of poly(clonal) antibodies and deeper sequencing. PhaeoEpiView ( https://PhaeoEpiView.univ-nantes.fr ) will be continuously updated with the newly published epigenomic data, making it the largest and richest epigenome browser of any stramenopile. In the upcoming era of molecular environmental studies, where epigenetics plays a significant role, we anticipate that PhaeoEpiView will become a widely used tool.

The model diatom Phaeodactylum tricornutum provides insights into the diversity and function of microeukaryotic DNA methyltransferases.

Cytosine methylation is an important epigenetic mark involved in the transcriptional control of transposable elements in mammals, plants and fungi. The Stramenopiles-Alveolate-Rhizaria (SAR) lineages are a major group of ecologically important marine microeukaryotes, including the phytoplankton groups diatoms and dinoflagellates. However, little is known about their DNA methyltransferase diversity. Here, we performed an in-silico analysis of DNA methyltransferases found in marine microeukaryotes and showed that they encode divergent DNMT3, DNMT4, DNMT5 and DNMT6 enzymes. Furthermore, we found three classes of enzymes within the DNMT5 family. Using a CRISPR/Cas9 strategy we demonstrated that the loss of the DNMT5a gene correlates with a global depletion of DNA methylation and overexpression of young transposable elements in the model diatom Phaeodactylum tricornutum. The study provides a view of the structure and function of a DNMT family in the SAR supergroup using an attractive model species.

Characterization of a Marine Diatom Chitin Synthase Using a Combination of Meta-Omics, Genomics, and Heterologous Expression Approaches.

ß-Chitin has important ecological and physiological roles and potential for widespread applications, but the characterization of chitin-related enzymes from ß-chitin producers was rarely reported. Querying against the Tara Oceans Gene Atlas, 4,939 chitin-related unique sequences from 12 Pfam accessions were found in Bacillariophyta metatranscriptomes. Putative chitin synthase (CHS) sequences are decreasingly present in Crustacea (39%), Stramenopiles (16%) and Insecta (14%) from the Marine Atlas of Tara Oceans Unigenes version 1 Metatranscriptomes (MATOUv1+T) database. A CHS gene from the model diatom Thalassiosira pseudonana (Thaps3_J4413, designated TpCHS1) was identified. Homology analysis of TpCHS1 in Marine Microbial Eukaryote Transcriptome Sequencing Project (MMETSP), PhycoCosm, and the PLAZA diatom omics data set showed that Mediophyceae and Thalassionemales species were potential new ß-chitin producers besides Thalassiosirales. TpCHS1 was overexpressed in Saccharomyces cerevisiae and Phaeodactylum tricornutum. In transgenic P. tricornutum lines, TpCHS1-eGFP localizes to the Golgi apparatus and plasma membrane and predominantly accumulates in the cleavage furrow during cell division. Enhanced TpCHS1 expression could induce abnormal cell morphology and reduce growth rates in P. tricornutum, which might be ascribed to the inhibition of the G2/M phase. S. cerevisiae was proved to be a better system for expressing large amounts of active TpCHS1, which effectively incorporates UDP-N-acetylglucosamine in radiometric in vitro assays. Our study expands the knowledge on chitin synthase taxonomic distribution in marine eukaryotic microbes, and is the first to collectively characterize an active marine diatom CHS which may play an important role during cell division. IMPORTANCE As the most abundant biopolymer in the oceans, the significance of chitin and its biosynthesis is rarely demonstrated in diatoms, which are the main contributors to the primary productivity of the oceans, ascribed to their huge biomass and efficient photosynthesis. We retrieved genes involved in chitin-based metabolism against the Tara Oceans Gene Atlas to expand our knowledge about their diversity and distribution in the marine environment. Potential new producers of chitin were found from the analysis of various algal transcriptome and genome databases. Heterologous expression confirms that Thalassiosira pseudonana contains an active chitin synthase (CHS) which may play an important role in the cell division process of diatoms. This study provides new insight into CHS geographic and taxonomic distribution in marine eukaryotic microbes, as well as into a new CHS functioning in the biosynthesis of ß-chitin in diatoms.

An Efficient Chromatin Immunoprecipitation Protocol for the Analysis of Histone Modification Distributions in the Brown Alga Ectocarpus.

The brown algae are an important but understudied group of multicellular marine organisms. A number of genetic and genomic tools have been developed for the model brown alga Ectocarpus; this includes, most recently, chromatin immunoprecipitation methodology, which allows genome-wide detection and analysis of histone post-translational modifications. Post-translational modifications of histone molecules have been shown to play an important role in gene regulation in organisms from other major eukaryotic lineages, and this methodology will therefore be a very useful tool to investigate genome function in the brown algae. This article provides a detailed, step-by-step description of the Ectocarpus ChIP protocol, which effectively addresses the difficult problem of efficiently extracting chromatin from cells protected by a highly resistant cell wall. The protocol described here will be an essential tool for the future application of chromatin analysis methodologies in brown algal research.

Chromatin landscape associated with sexual differentiation in a UV sex determination system.

In many eukaryotes, such as dioicous mosses and many algae, sex is determined by UV sex chromosomes and is expressed during the haploid phase of the life cycle. In these species, the male and female developmental programs are initiated by the presence of the U- or V-specific regions of the sex chromosomes but, as in XY and ZW systems, sexual differentiation is largely driven by autosomal sex-biased gene expression. The mechanisms underlying the regulation of sex-biased expression of genes during sexual differentiation remain elusive. Here, we investigated the extent and nature of epigenomic changes associated with UV sexual differentiation in the brown alga Ectocarpus, a model UV system. Six histone modifications were quantified in near-isogenic lines, leading to the identification of 16 chromatin signatures across the genome. Chromatin signatures correlated with levels of gene expression and histone PTMs changes in males versus females occurred preferentially at genes involved in sex-specific pathways. Despite the absence of chromosome scale dosage compensation and the fact that UV sex chromosomes recombine across most of their length, the chromatin landscape of these chromosomes was remarkably different to that of autosomes. Hotspots of evolutionary young genes in the pseudoautosomal regions appear to drive the exceptional chromatin features of UV sex chromosomes.

Genome wide natural variation of H3K27me3 selectively marks genes predicted to be important for cell differentiation in Phaeodactylum tricornutum.

In multicellular organisms, Polycomb Repressive Complex2 (PRC2) is known to deposit tri-methylation of lysine 27 of histone H3 (H3K27me3) to establish and maintain gene silencing, critical for developmentally regulated processes. The PRC2 complex is absent in both widely studied model yeasts, which initially suggested that PRC2 arose with the emergence of multicellularity. However, its discovery in several unicellular species including microalgae questions its role in unicellular eukaryotes. Here, we use Phaeodactylum tricornutum enhancer of zeste E(z) knockouts and show that P. tricornutum E(z) is responsible for di- and tri-methylation of lysine 27 of histone H3. H3K27me3 depletion abolishes cell morphology in P. tricornutum providing evidence for its role in cell differentiation. Genome-wide profiling of H3K27me3 in fusiform and triradiate cells further revealed genes that may specify cell identity. These results suggest a role for PRC2 and its associated mark in cell differentiation in unicellular species, and highlight their ancestral function in a broader evolutionary context than currently is appreciated.

Genome-wide analysis of allele-specific expression of genes in the model diatom Phaeodactylum tricornutum.

Recent advances in next generation sequencing technologies have allowed the discovery of widespread autosomal allele-specific expression (aASE) in mammals and plants with potential phenotypic effects. Extensive numbers of genes with allele-specific expression have been described in the diatom Fragilariopsis cylindrus in association with adaptation to external cues, as well as in Fistulifera solaris in the context of natural hybridization. However, the role of aASE and its extent in diatoms remain elusive. In this study, we investigate allele-specific expression in the model diatom Phaeodactylum tricornutum by the re-analysis of previously published whole genome RNA sequencing data and polymorphism calling. We found that 22% of P. tricornutum genes show moderate bias in allelic expression while 1% show nearly complete monoallelic expression. Biallelic expression associates with genes encoding components of protein metabolism while moderately biased genes associate with functions in catabolism and protein transport. We validated candidate genes by pyrosequencing and found that moderate biases in allelic expression were less stable than monoallelically expressed genes that showed consistent bias upon experimental validations at the population level and in subcloning experiments. Our approach provides the basis for the analysis of aASE in P. tricornutum and could be routinely implemented to test for variations in allele expression under different environmental conditions.

Histone modifications during the life cycle of the brown alga Ectocarpus.

BACKGROUND: Brown algae evolved complex multicellularity independently of the animal and land plant lineages and are the third most developmentally complex phylogenetic group on the planet. An understanding of developmental processes in this group is expected to provide important insights into the evolutionary events necessary for the emergence of complex multicellularity. Here, we focus on mechanisms of epigenetic regulation involving post-translational modifications of histone proteins. RESULTS: A total of 47 histone post-translational modifications are identified, including a novel mark H2AZR38me1, but Ectocarpus lacks both H3K27me3 and the major polycomb complexes. ChIP-seq identifies modifications associated with transcription start sites and gene bodies of active genes and with transposons. H3K79me2 exhibits an unusual pattern, often marking large genomic regions spanning several genes. Transcription start sites of closely spaced, divergently transcribed gene pairs share a common nucleosome-depleted region and exhibit shared histone modification peaks. Overall, patterns of histone modifications are stable through the life cycle. Analysis of histone modifications at generation-biased genes identifies a correlation between the presence of specific chromatin marks and the level of gene expression. CONCLUSIONS: The overview of histone post-translational modifications in the brown alga presented here will provide a foundation for future studies aimed at understanding the role of chromatin modifications in the regulation of brown algal genomes.

Paleo-diatom composition from Santa Barbara Basin deep-sea sediments: a comparison of 18S-V9 and diat-rbcL metabarcoding vs shotgun metagenomics.

Sedimentary ancient DNA (sedaDNA) analyses are increasingly used to reconstruct marine ecosystems. The majority of marine sedaDNA studies use a metabarcoding approach (extraction and analysis of specific DNA fragments of a defined length), targeting short taxonomic marker genes. Promising examples are 18S-V9 rRNA (~121-130 base pairs, bp) and diat-rbcL (76 bp), targeting eukaryotes and diatoms, respectively. However, it remains unknown how 18S-V9 and diat-rbcL derived compositional profiles compare to metagenomic shotgun data, the preferred method for ancient DNA analyses as amplification biases are minimised. We extracted DNA from five Santa Barbara Basin sediment samples (up to ~11 000 years old) and applied both a metabarcoding (18S-V9 rRNA, diat-rbcL) and a metagenomic shotgun approach to (i) compare eukaryote, especially diatom, composition, and (ii) assess sequence length and database related biases. Eukaryote composition differed considerably between shotgun and metabarcoding data, which was related to differences in read lengths (~112 and ~161 bp, respectively), and overamplification of short reads in metabarcoding data. Diatom composition was influenced by reference bias that was exacerbated in metabarcoding data and characterised by increased representation of Chaetoceros, Thalassiosira and Pseudo-nitzschia. Our results are relevant to sedaDNA studies aiming to accurately characterise paleo-ecosystems from either metabarcoding or metagenomic data.

PhaeoNet: A Holistic RNAseq-Based Portrait of Transcriptional Coordination in the Model Diatom Phaeodactylum tricornutum.

Transcriptional coordination is a fundamental component of prokaryotic and eukaryotic cell biology, underpinning the cell cycle, physiological transitions, and facilitating holistic responses to environmental stress, but its overall dynamics in eukaryotic algae remain poorly understood. Better understanding of transcriptional partitioning may provide key insights into the primary metabolism pathways of eukaryotic algae, which frequently depend on intricate metabolic associations between the chloroplasts and mitochondria that are not found in plants. Here, we exploit 187 publically available RNAseq datasets generated under varying nitrogen, iron and phosphate growth conditions to understand the co-regulatory principles underpinning transcription in the model diatom Phaeodactylum tricornutum. Using WGCNA (Weighted Gene Correlation Network Analysis), we identify 28 merged modules of co-expressed genes in the P. tricornutum genome, which show high connectivity and correlate well with previous microarray-based surveys of gene co-regulation in this species. We use combined functional, subcellular localization and evolutionary annotations to reveal the fundamental principles underpinning the transcriptional co-regulation of genes implicated in P. tricornutum chloroplast and mitochondrial metabolism, as well as the functions of diverse transcription factors underpinning this co-regulation. The resource is publically available as PhaeoNet, an advanced tool to understand diatom gene co-regulation.

Genome-enabled phylogenetic and functional reconstruction of an araphid pennate diatom Plagiostriata sp. CCMP470, previously assigned as a radial centric diatom, and its bacterial commensal.

Diatoms are an ecologically fundamental and highly diverse group of algae, dominating marine primary production in both open-water and coastal communities. The diatoms include both centric species, which may have radial or polar symmetry, and the pennates, which include raphid and araphid species and arose within the centric lineage. Here, we use combined microscopic and molecular information to reclassify a diatom strain CCMP470, previously annotated as a radial centric species related to Leptocylindrus danicus, as an araphid pennate species in the staurosiroid lineage, within the genus Plagiostriata. CCMP470 shares key ultrastructural features with Plagiostriata taxa, such as the presence of a sternum with parallel striae, and the presence of a highly reduced labiate process on its valve; and this evolutionary position is robustly supported by multigene phylogenetic analysis. We additionally present a draft genome of CCMP470, which is the first genome available for a staurosiroid lineage. 270 Pfams (19%) found in the CCMP470 genome are not known in other diatom genomes, which otherwise does not hold big novelties compared to genomes of non-staurosiroid diatoms. Notably, our DNA library contains the genome of a bacterium within the Rhodobacterales, an alpha-proteobacterial lineage known frequently to associate with algae. We demonstrate the presence of commensal alpha-proteobacterial sequences in other published algal genome and transcriptome datasets, which may indicate widespread and persistent co-occurrence.

Single-base methylome profiling of the giant kelp Saccharina japonica reveals significant differences in DNA methylation to microalgae and plants.

Brown algae have convergently evolved plant-like body plans and reproductive cycles, which in plants are controlled by differential DNA methylation. This contribution provides the first single-base methylome profiles of haploid gametophytes and diploid sporophytes of a multicellular alga. Although only c. 1.4% of cytosines in Saccharina japonica were methylated mainly at CHH sites and characterized by 5-methylcytosine (5mC), there were significant differences between life-cycle stages. DNA methyltransferase 2 (DNMT2), known to efficiently catalyze tRNA methylation, is assumed to methylate the genome of S. japonica in the structural context of tRNAs as the genome does not encode any other DNA methyltransferases. Circular and long noncoding RNA genes were the most strongly methylated regulatory elements in S. japonica. Differential expression of genes was negatively correlated with DNA methylation with the highest methylation levels measured in both haploid gametophytes. Hypomethylated and highly expressed genes in diploid sporophytes included genes involved in morphogenesis and halogen metabolism. The data herein provide evidence that cytosine methylation, although occurring at a low level, is significantly contributing to the formation of different life-cycle stages, tissue differentiation and metabolism in brown algae.

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

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

Corrigendum: A Potential Role for Epigenetic Processes in the Acclimation Response to Elevated pCO2 in the Model Diatom Phaeodactylum tricornutum.

[This corrects the article DOI: 10.3389/fmicb.2018.03342.].

Comparative characterization of putative chitin deacetylases from Phaeodactylum tricornutum and Thalassiosira pseudonana highlights the potential for distinct chitin-based metabolic processes in diatoms.

Chitin is generally considered to be present in centric diatoms but not in pennate species. Many aspects of chitin biosynthetic pathways have not been explored in diatoms. We retrieved chitin metabolic genes from pennate (Phaeodactylum tricornutum) and centric (Thalassiosira pseudonana) diatom genomes. Chitin deacetylase (CDA) genes from each genome (PtCDA and TpCDA) were overexpressed in P. tricornutum. We performed comparative analysis of their sequence structure, phylogeny, transcriptional profiles, localization and enzymatic activities. The chitin relevant proteins show complex subcellular compartmentation. PtCDA was likely acquired by horizontal gene transfer from prokaryotes, whereas TpCDA has closer relationships with sequences in Opisthokonta. Using transgenic P. tricornutum lines expressing CDA-green fluorescent protein (GFP) fusion proteins, PtCDA predominantly localizes to Golgi apparatus whereas TpCDA localizes to endoplasmic reticulum/chloroplast endoplasmic reticulum membrane. CDA-GFP overexpression upregulated the transcription of chitin synthases and potentially enhanced the ability of chitin synthesis. Although both CDAs are active on GlcNAc5 , TpCDA is more active on the highly acetylated chitin polymer DA60. We have addressed the ambiguous characters of CDAs from P. tricornutum and T. pseudonana. Differences in localization, evolution, expression and activities provide explanations underlying the greater potential of centric diatoms for chitin biosynthesis. This study paves the way for in vitro applications of novel CDAs.