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

From nets to networks: tools for deciphering phytoplankton metabolic interactions within communities and their global significance

Our oceans are populated with a wide diversity of planktonic organisms that form complex dynamic communities at the base of marine trophic networks. Within such communities are phytoplankton, unicellular photosynthetic taxa that provide an estimated half of global primary production and support biogeochemical cycles, along with other essential ecosystem services. One of the major challenges for microbial ecologists has been to try to make sense of this complexity. While phytoplankton distributions can be well explained by abiotic factors such as temperature and nutrient availability, there is increasing evidence that their ecological roles are tightly linked to their metabolic interactions with other plankton members through complex mechanisms (e.g. competition and symbiosis). Therefore, unravelling phytoplankton metabolic interactions is the key for inferring their dependency on, or antagonism with, other taxa and better integrating them into the context of carbon and nutrient fluxes in marine trophic networks. In this review, we attempt to summarize the current knowledge brought by ecophysiology, organismal imaging, in silico predictions and co-occurrence networks using 'omics data, highlighting successful combinations of approaches that may be helpful for future investigations of phytoplankton metabolic interactions within their complex communities.This article is part of the theme issue 'Connected interactions: enriching food web research by spatial and social interactions'.

A key gene, violaxanthin de-epoxidase-like 1, enhances fucoxanthin accumulation in Phaeodactylum tricornutum

BACKGROUND

Fucoxanthin has been widely investigated owing to its beneficial biological properties, and the model diatom Phaeodactylum tricornutum, possessing fucoxanthin (Fux) chlorophyll proteins as light-harvesting systems, is considered to have the potential to become a commercial cell factory for the pigment production.

RESULTS

Here, we compared the pigment contents in 10 different P. tricornutum strains from the globe, and found that strain CCMP631 (Pt6) exhibited the highest Fux content but with a low biomass. Comparison of mRNA levels revealed that higher Fux content in Pt6 was related with the higher expression of gene violaxanthin de-epoxidase-like (VDL) protein 1 (VDL1), which encodes the enzyme catalyzing the tautomerization of violaxanthin to neoxanthin in Fux biosynthesis pathway. Single nucleotide variants of VDL1 gene and allele-specific expression in strains Pt1 (the whole genome sequenced strain CCMP632) and Pt6 were analyzed, and overexpressing of each of the 4 VDL1 alleles, two from Pt1 and two from Pt6, in strain Pt1 leads to an increase in downstream product diadinoxanthin and channels the pigments towards Fux biosynthesis. All the 8 VDL1 overexpression (OE) lines showed significant increases by 8.2 to 41.7% in Fux content without compromising growth, and VDL1 Allele 2 OE lines even exhibited the higher cell density on day 8, with an increase by 24.2-28.7% in two Pt1VDL1-allele 2 OE lines and 7.1-11.1% in two Pt6VDL1-allele 2 OE lines, respectively.

CONCLUSIONS

The results reveal VDL1, localized in the plastid stroma, plays a key role in Fux over-accumulation in P. tricornutum. Overexpressing VDL1, especially allele 2, improved both the Fux content and growth rate, which provides a new strategy for the manipulation of Fux production in the future.

Allele-dependent expression and functionality of lipid enzyme phospholipid:diacylglycerol acyltransferase affect diatom carbon storage and growth

In the acyl-CoA-independent pathway of triacylglycerol (TAG) synthesis unique to plants, fungi, and algae, TAG formation is catalyzed by the enzyme phospholipid:diacylglycerol acyltransferase (PDAT). The unique PDAT gene of the model diatom Phaeodactylum tricornutum strain CCMP2561 boasts 47 single nucleotide variants within protein coding regions of the alleles. To deepen our understanding of TAG synthesis, we observed the allele-specific expression of PDAT by the analysis of 87 published RNA-sequencing (RNA-seq) data and experimental validation. The transcription of one of the two PDAT alleles, Allele 2, could be specifically induced by decreasing nitrogen concentrations. Overexpression of Allele 2 in P. tricornutum substantially enhanced the accumulation of TAG by 44% to 74% under nutrient stress; however, overexpression of Allele 1 resulted in little increase of TAG accumulation. Interestingly, a more serious growth inhibition was observed in the PDAT Allele 1 overexpression strains compared with Allele 2 counterparts. Heterologous expression in yeast (Saccharomyces cerevisiae) showed that enzymes encoded by PDAT Allele 2 but not Allele 1 had TAG biosynthetic activity, and 7 N-terminal and 3 C-terminal amino acid variants between the 2 allele-encoded proteins substantially affected enzymatic activity. P. tricornutum PDAT, localized in the innermost chloroplast membrane, used monogalactosyldiacylglycerol and phosphatidylcholine as acyl donors as demonstrated by the increase of the 2 lipids in PDAT knockout lines, which indicated a common origin in evolution with green algal PDATs. Our study reveals unequal roles among allele-encoded PDATs in mediating carbon storage and growth in response to nitrogen stress and suggests an unsuspected strategy toward lipid and biomass improvement for biotechnological purposes.

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.

Phaeodactylum tricornutum: An established model species for diatom molecular research and an emerging chassis for algal synthetic biology

Diatoms are prominent and highly diverse microalgae in aquatic environments. Compared with other diatom species, Phaeodactylum tricornutum is an "atypical diatom" displaying three different morphotypes and lacking the usual silica shell. Despite being of limited ecological relevance, its ease of growth in the laboratory and well-known physiology, alongside the steady increase in genome-enabled information coupled with effective tools for manipulating gene expression, have meant it has gained increased recognition as a powerful experimental model for molecular research on diatoms. We here present a brief overview of how over the last 25 years P. tricornutum has contributed to the unveiling of fundamental aspects of diatom biology, while also emerging as a new tool for algal process engineering and synthetic biology.

Comparison of two Phaeodactylum tricornutum ecotypes under nitrogen starvation and resupply reveals distinct lipid accumulation strategies but a common degradation process

Introduction

Phaeodactylum tricornutum is a model species frequently used to study lipid metabolism in diatoms. When exposed to a nutrient limitation or starvation, diatoms are known to accumulate neutral lipids in cytoplasmic lipid droplets (LDs). Those lipids are produced partly de novo and partly from the recycle of plastid membrane lipids. Under a nitrogen resupply, the accumulated lipids are catabolized, a phenomenon about which only a few data are available. Various strains of P. tricornutum have been isolated around the world that may differ in lipid accumulation patterns.

Methods

To get further information on this topic, two genetically distant ecotypes of P. tricornutum (Pt1 and Pt4) have been cultivated under nitrogen deprivation during 11 days followed by a resupply period of 3 days. The importance of cytoplasmic LDs relative to the plastid was assessed by a combination of confocal laser scanning microscopy and cell volume estimation using bright field microscopy pictures.

Results and discussion

We observed that in addition to a basal population of small LDs (0.005 μm3 to 0.7 μm3) present in both strains all along the experiment, Pt4 cells immediately produced two large LDs (up to 12 μm3 after 11 days) while Pt1 cells progressively produced a higher number of smaller LDs (up to 7 μm3 after 11 days). In this work we showed that, in addition to intracellular available space, lipid accumulation may be limited by the pre-starvation size of the plastid as a source of membrane lipids to be recycled. After resupplying nitrogen and for both ecotypes, a fragmentation of the largest LDs was observed as well as a possible migration of LDs to the vacuoles that would suggest an autophagic degradation. Altogether, our results deepen the understanding of LDs dynamics and open research avenues for a better knowledge of lipid degradation in diatoms.

Transporter Proteins as Ecological Assets and Features of Microbial Eukaryotic Pangenomes

Here we review two connected themes in evolutionary microbiology: (a) the nature of gene repertoire variation within species groups (pangenomes) and (b) the concept of metabolite transporters as accessory proteins capable of providing niche-defining "bolt-on" phenotypes. We discuss the need for improved sampling and understanding of pangenome variation in eukaryotic microbes. We then review the factors that shape the repertoire of accessory genes within pangenomes. As part of this discussion, we outline how gene duplication is a key factor in both eukaryotic pangenome variation and transporter gene family evolution. We go on to outline how, through functional characterization of transporter-encoding genes, in combination with analyses of how transporter genes are gained and lost from accessory genomes, we can reveal much about the niche range, the ecology, and the evolution of virulence of microbes. We advocate for the coordinated systematic study of eukaryotic pangenomes through genome sequencing and the functional analysis of genes found within the accessory gene repertoire.

Temperature sensitivity of carbon concentrating mechanisms in the diatom Phaeodactylum tricornutum

Marine diatoms are key primary producers across diverse habitats in the global ocean. Diatoms rely on a biophysical carbon concentrating mechanism (CCM) to supply high concentrations of CO₂ around their carboxylating enzyme, RuBisCO. The necessity and energetic cost of the CCM are likely to be highly sensitive to temperature, as temperature impacts CO₂ concentration, diffusivity, and the kinetics of CCM components. Here, we used membrane inlet mass spectrometry (MIMS) and modeling to capture temperature regulation of the CCM in the diatom Phaeodactylum tricornutum (Pt). We found that enhanced carbon fixation rates by Pt at elevated temperatures were accompanied by increased CCM activity capable of maintaining RuBisCO close to CO₂ saturation but that the mechanism varied. At 10 and 18 °C, diffusion of CO₂ into the cell, driven by Pt's 'chloroplast pump' was the major inorganic carbon source. However, at 18 °C, upregulation of the chloroplast pump enhanced (while retaining the proportion of) both diffusive CO₂ and active HCO₃^- uptake into the cytosol, and significantly increased chloroplast HCO₃^- concentrations. In contrast, at 25 °C, compared to 18 °C, the chloroplast pump had only a slight increase in activity. While diffusive uptake of CO₂ into the cell remained constant, active HCO₃^- uptake across the cell membrane increased resulting in Pt depending equally on both CO₂ and HCO₃^- as inorganic carbon sources. Despite changes in the CCM, the overall rate of active carbon transport remained double that of carbon fixation across all temperatures tested. The implication of the energetic cost of the Pt CCM in response to increasing temperatures was discussed.

Whole-genome scanning reveals environmental selection mechanisms that shape diversity in populations of the epipelagic diatom Chaetoceros

Diatoms form a diverse and abundant group of photosynthetic protists that are essential players in marine ecosystems. However, the microevolutionary structure of their populations remains poorly understood, particularly in polar regions. Exploring how closely related diatoms adapt to different environments is essential given their short generation times, which may allow rapid adaptations, and their prevalence in marine regions dramatically impacted by climate change, such as the Arctic and Southern Oceans. Here, we address genetic diversity patterns in Chaetoceros, the most abundant diatom genus and one of the most diverse, using 11 metagenome-assembled genomes (MAGs) reconstructed from Tara Oceans metagenomes. Genome-resolved metagenomics on these MAGs confirmed a prevalent distribution of Chaetoceros in the Arctic Ocean with lower dispersal in the Pacific and Southern Oceans as well as in the Mediterranean Sea. Single-nucleotide variants identified within the different MAG populations allowed us to draw a landscape of Chaetoceros genetic diversity and revealed an elevated genetic structure in some Arctic Ocean populations. Gene flow patterns of closely related Chaetoceros populations seemed to correlate with distinct abiotic factors rather than with geographic distance. We found clear positive selection of genes involved in nutrient availability responses, in particular for iron (e.g., ISIP2a, flavodoxin), silicate, and phosphate (e.g., polyamine synthase), that were further supported by analysis of Chaetoceros transcriptomes. Altogether, these results highlight the importance of environmental selection in shaping diatom diversity patterns and provide new insights into their metapopulation genomics through the integration of metagenomic and environmental data.

Increased genetic diversity loss and genetic differentiation in a model marine diatom adapted to ocean warming compared to high CO2

Although high CO2 and warming could act interactively on marine phytoplankton, little is known about the molecular basis for this interaction on an evolutionary scale. Here we explored the adaptation to high CO2 in combination with warming in a model marine diatom Phaeodactylum tricornutum. Whole-genome re-sequencing identifies, in comparison to populations grown under control conditions, a larger genetic diversity loss and a higher genetic differentiation in the populations adapted for 2 years to warming than in those adapted to high CO2. However, this diversity loss was less under high CO2 combined with warming, suggesting that the evolution driven by warming was constrained by high CO2. By integrating genomics, transcriptomics, and physiological data, we found that the underlying molecular basis for this constraint is associated with the expression of genes involved in some key metabolic pathways or biological processes, such as the glyoxylate pathway, amino acid and fatty acid metabolism, and diel variability. Our results shed new light on the evolutionary responses of marine phytoplankton to multiple environmental changes in the context of global change and provide new insights into the molecular basis underpinning interactions among those multiple drivers.

What Was Old Is New Again: The Pennate Diatom Haslea ostrearia (Gaillon) Simonsen in the Multi-Omic Age

The marine pennate diatom Haslea ostrearia has long been known for its characteristic blue pigment marennine, which is responsible for the greening of invertebrate gills, a natural phenomenon of great importance for the oyster industry. For two centuries, this taxon was considered unique; however, the recent description of a new blue Haslea species revealed unsuspected biodiversity. Marennine-like pigments are natural blue dyes that display various biological activities—e.g., antibacterial, antioxidant and antiproliferative—with a great potential for applications in the food, feed, cosmetic and health industries. Regarding fundamental prospects, researchers use model organisms as standards to study cellular and physiological processes in other organisms, and there is a growing and crucial need for more, new and unconventional model organisms to better correspond to the diversity of the tree of life. The present work, thus, advocates for establishing H. ostrearia as a new model organism by presenting its pros and cons—i.e., the interesting aspects of this peculiar diatom (representative of benthic-epiphytic phytoplankton, with original behavior and chemodiversity, controlled sexual reproduction, fundamental and applied-oriented importance, reference genome, and transcriptome will soon be available); it will also present the difficulties encountered before this becomes a reality as it is for other diatom models (the genetics of the species in its infancy, the transformation feasibility to be explored, the routine methods needed to cryopreserve strains of interest).

Comparative Proteomic Analysis of the Diatom Phaeodactylum tricornutum Reveals New Insights Into Intra- and Extra-Cellular Protein Contents of Its Oval, Fusiform, and Triradiate Morphotypes

Phaeodactylum tricornutum is an atypical diatom since it can display three main morphotypes: fusiform, triradiate, and oval. Such pleomorphism is possible thanks to an original metabolism, which is tightly regulated in order to acclimate to environmental conditions. Currently, studies dedicated to the comparison of each morphotype issued from one specific strain are scarce and little information is available regarding the physiological significance of this morphogenesis. In this study, we performed a comparative proteomic analysis of the three morphotypes from P. tricornutum. Cultures highly enriched in one dominant morphotype (fusiform, triradiate, or oval) of P. tricornutum Pt3 strain were used. Pairwise comparisons highlighted biological processes, which are up- and down-regulated in the oval (e.g., purine and cellular amino acid metabolism) and triradiate morphotypes (e.g., oxido-reduction and glycolytic processes) compared to the fusiform one used as a reference. Intersection analysis allowed us to identify the specific features of the oval morphotype. Results from this study confirmed previous transcriptomic RNA sequencing observation showing that the oval cells present a distinct metabolism with specific protein enrichment compared to fusiform and triradiate cells. Finally, the analysis of the secretome of each morphotype was also performed.

Metabolite Quantification by Fourier Transform Infrared Spectroscopy in Diatoms: Proof of Concept on Phaeodactylum tricornutum

Diatoms are feedstock for the production of sustainable biocommodities, including biofuel. The biochemical characterization of newly isolated or genetically modified strains is seminal to identify the strains that display interesting features for both research and industrial applications. Biochemical quantification of organic macromolecules cellular quotas are time-consuming methodologies which often require large amount of biological sample. Vibrational spectroscopy is an essential tool applied in several fields of research. A Fourier transform infrared (FTIR) microscopy-based imaging protocol was developed for the simultaneous cellular quota quantification of lipids, carbohydrates, and proteins of the diatom Phaeodactylum tricornutum. The low amount of sample required for the quantification allows the high throughput quantification on small volume cultures. A proof of concept was performed (1) on nitrogen-starved experimental cultures and (2) on three different P. tricornutum wild-type strains. The results are supported by the observation in situ of lipid droplets by confocal and brightfield microscopy. The results show that major differences exist in the regulation of lipid metabolism between ecotypes of P. tricornutum.

The genome of a nonphotosynthetic diatom provides insights into the metabolic shift to heterotrophy and constraints on the loss of photosynthesis

Although most of the tens of thousands of diatom species are photoautotrophs, a small number of heterotrophic species no longer photosynthesize. We sequenced the genome of a nonphotosynthetic diatom, Nitzschia Nitz4, to determine how carbon metabolism was altered in the wake of this trophic shift. Nitzschia Nitz4 has retained its plastid and plastid genome, but changes associated with the transition to heterotrophy were cellular-wide and included losses of photosynthesis-related genes from the nuclear and plastid genomes, elimination of isoprenoid biosynthesis in the plastid, and remodeling of mitochondrial glycolysis to maximize adenosine triphosphte (ATP) yield. The genome contains a β-ketoadipate pathway that may allow Nitzschia Nitz4 to metabolize lignin-derived compounds. Diatom plastids lack an oxidative pentose phosphate pathway (oPPP), leaving photosynthesis as the primary source of NADPH to support essential biosynthetic pathways in the plastid and, by extension, limiting available sources of NADPH in nonphotosynthetic plastids. The genome revealed similarities between nonphotosynthetic diatoms and apicomplexan parasites for provisioning NADPH in their plastids and highlighted the ancestral absence of a plastid oPPP as a potentially important constraint on loss of photosynthesis, a hypothesis supported by the higher frequency of transitions to parasitism or heterotrophy in lineages that have a plastid oPPP.

Potential for Biomass Production and Remediation by Cultivation of the Marine Model Diatom Phaeodactylum tricornutum in Oil Field Produced Wastewater Media

While oilfield produced water (PW) is one of the largest, unclaimed wastewater streams of the oil industry, it could potentially be used as a cultivation medium for microalgae. Microalgae could help with the remediation of this water while also delivering biomass that can be transformed into valuable byproducts such as biofuels. The coupling of these two purposes is expected to cut production costs of biofuels while aiding environmental protection. In this study, we compared the cultivation capacity of the marine model diatom Phaeodactylum tricornutum in media at varying salinities and in media composed of PW from two oilfields in the Central Valley of California that differed drastically in the concentration of inorganic and organic constituents. Specifically, we measured the carrying capacity of these media, the maximum growth rates of P. tricornutum, its cellular lipid accumulation capacity, and its capacity to remediate the most polluted PW source. Our study shows that P. tricornutum can successfully adjust to the tested cultivation media through processes of short-term acclimation and long-term adaptation. Furthermore, the cultivation of P. tricornutum in the most heavily polluted PW source led to significant increases in cell yield and improved photosynthetic capacity during the stationary phase, which could be attributed chiefly to the higher levels of nitrate present in this PW source. Chemical water analyses also demonstrated the capability of P. tricornutum to remediate major nutrient content and potentially harmful elements like fluorine and copper. Because P. tricornutum is amenable to advanced genetic engineering, which could be taken advantage of to improve its cultivation resilience and productivity in an economic setting, we propose this study as a step towards essential follow-up studies that will identify the genetic regulation behind its growth in oilfield PW media and its remediation of the PW constituents.

Evolution of Mutation Rate in Astronomically Large Phytoplankton Populations

Genetic diversity is expected to be proportional to population size, yet, there is a well-known, but unexplained lack of genetic diversity in large populations-the "Lewontin's paradox." Larger populations are expected to evolve lower mutation rates, which may help to explain this paradox. Here, we test this conjecture by measuring the spontaneous mutation rate in a ubiquitous unicellular marine phytoplankton species Emiliania huxleyi (Haptophyta) that has modest genetic diversity despite an astronomically large population size. Genome sequencing of E. huxleyi mutation accumulation lines revealed 455 mutations, with an unusual GC-biased mutation spectrum. This yielded an estimate of the per site mutation rate µ = 5.55×10-10 (CI 95%: 5.05×10-10 - 6.09×10-10), which corresponds to an effective population size Ne ∼ 2.7×106. Such a modest Ne is surprising for a ubiquitous and abundant species that accounts for up to 10% of global primary productivity in the oceans. Our results indicate that even exceptionally large populations do not evolve mutation rates lower than ∼10-10 per nucleotide per cell division. Consequently, the extreme disparity between modest genetic diversity and astronomically large population size in the plankton species cannot be explained by an unusually low mutation rate.

Mitotic recombination between homologous chromosomes drives genomic diversity in diatoms

Diatoms, an evolutionarily successful group of microalgae, display high levels of intraspecific genetic variability in natural populations. However, the contribution of various mechanisms generating such diversity is unknown. Here we estimated the genetic micro-diversity within a natural diatom population and mapped the genomic changes arising within clonally propagated diatom cell cultures. Through quantification of haplotype diversity by next-generation sequencing and amplicon re-sequencing of selected loci, we documented a rapid accumulation of multiple haplotypes accompanied by the appearance of novel protein variants in cell cultures initiated from a single founder cell. Comparison of the genomic changes between mother and daughter cells revealed copy number variation and copy-neutral loss of heterozygosity leading to the fixation of alleles within individual daughter cells. The loss of heterozygosity can be accomplished by recombination between homologous chromosomes. To test this hypothesis, we established an endogenous readout system and estimated that the frequency of interhomolog mitotic recombination was under standard growth conditions 4.2 events per 100 cell divisions. This frequency is increased under environmental stress conditions, including treatment with hydrogen peroxide and cadmium. These data demonstrate that copy number variation and mitotic recombination between homologous chromosomes underlie clonal variability in diatom populations. We discuss the potential adaptive evolutionary benefits of the plastic response in the interhomolog mitotic recombination rate, and we propose that this may have contributed to the ecological success of diatoms.

Acclimation and adaptation to elevated pCO2 increase arsenic resilience in marine diatoms

Arsenic pollution is a widespread threat to marine life, but the ongoing rise pCO₂ levels is predicted to decrease bio-toxicity of arsenic. However, the effects of arsenic toxicity on marine primary producers under elevated pCO₂ are not well characterized. Here, we studied the effects of arsenic toxicity in three globally distributed diatom species (Phaeodactylum tricornutum, Thalassiosira pseudonana, and Chaetoceros mulleri) after short-term acclimation (ST, 30 days), medium-term exposure (MT, 750 days), and long-term (LT, 1460 days) selection under ambient (400 µatm) and elevated (1000 and 2000 µatm) pCO₂. We found that elevated pCO₂ alleviated arsenic toxicity even after short acclimation times but the magnitude of the response decreased after mid and long-term adaptation. When fed with these elevated pCO₂ selected diatoms, the scallop Patinopecten yessoensis had significantly lower arsenic content (3.26-52.83%). Transcriptomic and biochemical analysis indicated that the diatoms rapidly developed arsenic detoxification strategies, which included upregulation of transporters associated with shuttling harmful compounds out of the cell to reduce arsenic accumulation, and upregulation of proteins involved in synthesizing glutathione (GSH) to chelate intracellular arsenic to reduce arsenic toxicity. Thus, our results will expand our knowledge to fully understand the ecological risk of trace metal pollution under increasing human activity induced ocean acidification.

Re-examination of two diatom reference genomes using long-read sequencing

Background

The marine diatoms Thalassiosira pseudonana and Phaeodactylum tricornutum are valuable model organisms for exploring the evolution, diversity and ecology of this important algal group. Their reference genomes, published in 2004 and 2008, respectively, were the product of traditional Sanger sequencing. In the case of T. pseudonana, optical restriction site mapping was employed to further clarify and contextualize chromosome-level scaffolds. While both genomes are considered highly accurate and reasonably contiguous, they still contain many unresolved regions and unordered/unlinked scaffolds.

Results

We have used Oxford Nanopore Technologies long-read sequencing to update and validate the quality and contiguity of the T. pseudonana and P. tricornutum genomes. Fine-scale assessment of our long-read derived genome assemblies allowed us to resolve previously uncertain genomic regions, further characterize complex structural variation, and re-evaluate the repetitive DNA content of both genomes. We also identified 1862 previously undescribed genes in T. pseudonana. In P. tricornutum, we used transposable element detection software to identify 33 novel copia-type LTR-RT insertions, indicating ongoing activity and rapid expansion of this superfamily as the organism continues to be maintained in culture. Finally, Bionano optical mapping of P. tricornutum chromosomes was combined with long-read sequence data to explore the potential of long-read sequencing and optical mapping for resolving haplotypes.

Conclusion

Despite its potential to yield highly contiguous scaffolds, long-read sequencing is not a panacea. Even for relatively small nuclear genomes such as those investigated herein, repetitive DNA sequences cause problems for current genome assembly algorithms. Determining whether a long-read derived genomic assembly is ‘better’ than one produced using traditional sequence data is not straightforward. Our revised reference genomes for P. tricornutum and T. pseudonana nevertheless provide additional insight into the structure and evolution of both genomes, thereby providing a more robust foundation for future diatom research.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07666-3.

Comprehensive and Functional Analysis of Horizontal Gene Transfer Events in Diatoms

Diatoms are a diverse group of mainly photosynthetic algae, responsible for 20% of worldwide oxygen production, which can rapidly respond to favorable conditions and often outcompete other phytoplankton. We investigated the contribution of horizontal gene transfer (HGT) to its ecological success. A large-scale phylogeny-based prokaryotic HGT detection procedure across nine sequenced diatoms showed that 3-5% of their proteome has a horizontal origin and a large influx occurred at the ancestor of diatoms. More than 90% of HGT genes are expressed, and species-specific HGT genes in Phaeodactylum tricornutum undergo strong purifying selection. Genes derived from HGT are implicated in several processes including environmental sensing and expand the metabolic toolbox. Cobalamin (vitamin B12) is an essential cofactor for roughly half of the diatoms and is only produced by bacteria. Five consecutive genes involved in the final synthesis of the cobalamin biosynthetic pathway, which could function as scavenging and repair genes, were detected as HGT. The full suite of these genes was detected in the cold-adapted diatom Fragilariopsis cylindrus. This might give diatoms originating from the Southern Ocean, a region typically depleted in cobalamin, a competitive advantage. Overall, we show that HGT is a prevalent mechanism that is actively used in diatoms to expand its adaptive capabilities.

Integrative omics identification, evolutionary and structural analysis of low affinity nitrate transporters in diatoms, diNPFs

Diatoms are one of the major and most diverse groups of phytoplankton, with chimeric genomes harbouring a combination of genes of bacterial, animal and plant origin. They have developed sophisticated mechanisms to face environmental variations. In marine environments, nutrients concentration shows significant temporal and spatial variability, influencing phytoplankton growth. Among nutrients, nitrogen, present at micromolar levels, is often a limiting resource. Here, we report a comprehensive characterization of the Nitrate Transporter 1/Peptide Transporter Family (NPF) in diatoms, diNPFs. NPFs are well characterized in many organisms where they recognize a broad range of substrates, ranging from short-chained di- and tri-peptides in bacteria, fungi and mammals to a wide variety of molecules including nitrate in higher plants. Scarce information is available for diNPFs. We integrated-omics, phylogenetic, structural and expression analyses, to infer information on their role in diatoms. diNPF genes diverged to produce two distinct clades with strong sequence and structural homology with either bacterial or plant NPFs, with different predicted sub-cellular localization, suggesting that the divergence resulted in functional diversification. Moreover, transcription analysis of diNPF genes under different laboratory and environmental growth conditions suggests that diNPF diversification led to genetic adaptations that might contribute to diatoms ability to flourish in diverse environmental conditions.

New paradigm in diatom omics and genetic manipulation

Diatoms are one of the most heterogeneous eukaryotic plankton known for regulating earth's biogeochemical cycles and maintaining the marine ecosystems ever since the late Eocene epoch. The advent of multidisciplinary omics approach has both epitomized and revolutionized the nature of their chimeric genetic toolkit, ecophysiology, and metabolic adaptability as well as their interaction with other communities. In addition, advanced functional annotation of transcriptomic and proteomic data using cutting edge bioinformatics tools together with high-resolution genome-scale mathematical modeling has effectively proven as the catapult in solving genetic bottlenecks in microbial as well as diatom exploration. In this review, a corroborative summation of the robust work done in manipulating, engineering, and sequencing of the diatom genomes besides underpinning the holistic application of omics in transcription and translation has been discussed in order to shrewd their multifarious novel potential in the field of biotechnology and provide an insight into their dynamic evolutionary relevance.

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.

Mating type specific transcriptomic response to sex inducing pheromone in the pennate diatom Seminavis robusta

Sexual reproduction is a fundamental phase in the life cycle of most diatoms. Despite its role as a source of genetic variation, it is rarely reported in natural circumstances and its molecular foundations remain largely unknown. Here, we integrate independent transcriptomic datasets to prioritize genes responding to sex inducing pheromones (SIPs) in the pennate diatom Seminavis robusta. We observe marked gene expression changes associated with SIP treatment in both mating types, including an inhibition of S phase progression, chloroplast division, mitosis, and cell wall formation. Meanwhile, meiotic genes are upregulated in response to SIP, including a sexually induced diatom specific cyclin. Our data further suggest an important role for reactive oxygen species, energy metabolism, and cGMP signaling during the early stages of sexual reproduction. In addition, we identify several genes with a mating type specific response to SIP, and link their expression pattern with physiological specialization, such as the production of the attraction pheromone diproline in mating type − (MT−) and mate-searching behavior in mating type + (MT+). Combined, our results provide a model for early sexual reproduction in pennate diatoms and significantly expand the suite of target genes to detect sexual reproduction events in natural diatom populations.

Sexual reproduction potential implied by functional analysis of SPO11 in Phaeodactylum tricornutum

Phaeodactylum tricornutum is a model microalgae that is widely used to study diatom physiology and ecology. Since the meiotic process and sexual cycle have never been observed directly, P. tricornutum has been considered to be an asexual species. However, phylogenetic analysis of the P. tricornutum genome has revealed a series of meiosis-specific gene homologues in this species. We identified two copies of differently transcribed SPO11 homologs that contain the conserved motifs of Winged-helix and Toprim domains. The homolog PtSPO11-3 interacts with TopoVIB in yeast two-hybrid analysis, whereas the homolog PtSPO11-2 could rescue the sporulation defect of a Spo11 yeast mutant strain. PtSPO11-2 was also found to be significantly up-regulated at low temperatures in P. tricornutum and its key catalytic residue was important to the homolog's function in sporulation. The results herein provide positive clue that meiosis and sexual reproduction could exist in this diatom.

The Seminavis robusta genome provides insights into the evolutionary adaptations of benthic diatoms

Benthic diatoms are the main primary producers in shallow freshwater and coastal environments, fulfilling important ecological functions such as nutrient cycling and sediment stabilization. However, little is known about their evolutionary adaptations to these highly structured but heterogeneous environments. Here, we report a reference genome for the marine biofilm-forming diatom Seminavis robusta, showing that gene family expansions are responsible for a quarter of all 36,254 protein-coding genes. Tandem duplications play a key role in extending the repertoire of specific gene functions, including light and oxygen sensing, which are probably central for its adaptation to benthic habitats. Genes differentially expressed during interactions with bacteria are strongly conserved in other benthic diatoms while many species-specific genes are strongly upregulated during sexual reproduction. Combined with re-sequencing data from 48 strains, our results offer insights into the genetic diversity and gene functions in benthic diatoms.

Diatom Molecular Research Comes of Age: Model Species for Studying Phytoplankton Biology and Diversity

Diatoms are the world's most diverse group of algae, comprising at least 100,000 species. Contributing ∼20% of annual global carbon fixation, they underpin major aquatic food webs and drive global biogeochemical cycles. Over the past two decades, Thalassiosira pseudonana and Phaeodactylum tricornutum have become the most important model systems for diatom molecular research, ranging from cell biology to ecophysiology, due to their rapid growth rates, small genomes, and the cumulative wealth of associated genetic resources. To explore the evolutionary divergence of diatoms, additional model species are emerging, such as Fragilariopsis cylindrus and Pseudo-nitzschia multistriata Here, we describe how functional genomics and reverse genetics have contributed to our understanding of this important class of microalgae in the context of evolution, cell biology, and metabolic adaptations. Our review will also highlight promising areas of investigation into the diversity of these photosynthetic organisms, including the discovery of new molecular pathways governing the life of secondary plastid-bearing organisms in aquatic environments.