Genomes of the dinoflagellate Polarella glacialis encode tandemly repeated single-exon genes with adaptive functions

Stress responses in an Arctic microalga (Pelagophyceae) following sudden salinity change revealed by gene expression analysis

Marine microbes that have for eons been adapted to stable salinity regimes are confronted with sudden decreases in salinity in the Arctic Ocean. The episodic freshening is increasing due to climate change with melting multi-year sea-ice and glaciers, greater inflows from rivers, and increased precipitation. To investigate algal responses to lowered salinity, we analyzed the responses and acclimatation over 24 h in a non-model Arctic marine alga (pelagophyte CCMP2097) following transfer to realistic lower salinities. Using RNA-seq transcriptomics, here we show rapid differentially expressed genes related to stress oxidative responses, proteins involved in the photosystem and circadian clock, and those affecting lipids and inorganic ions. After 24 h the pelagophyte adjusted to the lower salinity seen in the overexpression of genes associated with freezing resistance, cold adaptation, and salt tolerance. Overall, a suite of ancient widespread pathways is recruited enabling the species to adjust to the stress of rapid salinity change.

Lessons from Extremophiles: Functional Adaptations and Genomic Innovations across the Eukaryotic Tree of Life

From hydrothermal vents, to glaciers, to deserts, research in extreme environments has reshaped our understanding of how and where life can persist. Contained within the genomes of extremophilic organisms are the blueprints for a toolkit to tackle the multitude of challenges of survival in inhospitable environments. As new sequencing technologies have rapidly developed, so too has our understanding of the molecular and genomic mechanisms that have facilitated the success of extremophiles. Although eukaryotic extremophiles remain relatively understudied compared to bacteria and archaea, an increasing number of studies have begun to leverage 'omics tools to shed light on eukaryotic life in harsh conditions. In this perspective paper, we highlight a diverse breadth of research on extremophilic lineages across the eukaryotic tree of life, from microbes to macrobes, that are collectively reshaping our understanding of molecular innovations at life's extremes. These studies are not only advancing our understanding of evolution and biological processes but are also offering a valuable roadmap on how emerging technologies can be applied to identify cellular mechanisms of adaptation to cope with life in stressful conditions, including high and low temperatures, limited water availability, and heavy metal habitats. We shed light on patterns of molecular and organismal adaptation across the eukaryotic tree of life and discuss a few promising research directions, including investigations into the role of horizontal gene transfer in eukaryotic extremophiles and the importance of increasing phylogenetic diversity of model systems.

Whole-genome duplication in an algal symbiont bolsters coral heat tolerance

The algal endosymbiont Durusdinium trenchii enhances the resilience of coral reefs under thermal stress. D. trenchii can live freely or in endosymbiosis, and the analysis of genetic markers suggests that this species has undergone whole-genome duplication (WGD). However, the evolutionary mechanisms that underpin the thermotolerance of this species are largely unknown. Here, we present genome assemblies for two D. trenchii isolates, confirm WGD in these taxa, and examine how selection has shaped the duplicated genome regions using gene expression data. We assess how the free-living versus endosymbiotic lifestyles have contributed to the retention and divergence of duplicated genes, and how these processes have enhanced the thermotolerance of D. trenchii. Our combined results suggest that lifestyle is the driver of post-WGD evolution in D. trenchii, with the free-living phase being the most important, followed by endosymbiosis. Adaptations to both lifestyles likely enabled D. trenchii to provide enhanced thermal stress protection to the host coral.

Assessing the global distribution and risk of harmful microalgae: A focus on three toxic Alexandrium dinoflagellates

Harmful dinoflagellates and their resulting blooms pose a threat to marine life and human health. However, to date, global maps of marine life often overlook harmful microorganisms. As harmful algal blooms (HABs) increase in frequency, severity, and extent, understanding the distribution of harmful dinoflagellates and their drivers is crucial for their management. We used MaxEnt, random forest, and ensemble models to map the habitats of the representative HABs species in the genus Alexandrium, including A. catenella, A. minutum, and A. pacificum. Since species occurrence records used in previous studies were solely morphology-based, potentially leading to misidentifications, we corrected these species' distribution records using molecular criteria. The results showed that the key environmental drivers included the distance to the coastline, bathymetry, sea surface temperature (SST), and dissolved oxygen. Alexandrium catenella thrives in temperate to cold zones and is driven by low SST and high oxygen levels. Alexandrium pacificum mainly inhabits the Temperate Northern Pacific and prefers warmer SST and lower oxygen levels. Alexandrium minutum thrives universally and adapts widely to SST and oxygen. By analyzing the habitat suitability of locations with recorded HAB occurrences, we found that high habitat suitability could serve as a reference indicator for bloom risk. Therefore, we have proposed a qualitative method to spatially assess the harmful algae risk according to the habitat suitability. On the global risk map, coastal temperate seas, such as the Mediterranean, Northwest Pacific, and Southern Australia, faced higher risks. Although HABs currently have restricted geographic distributions, our study found these harmful algae possess high environmental tolerance and can thrive across diverse habitats. HAB impacts could increase if climate changes or ocean conditions became more favorable. Marine transportation may also spread the harmful algae to new unaffected ecosystems. This study has pioneered the assessment of harmful algal risk based on habitat suitability.

Investigation of heterotrophs reveals new insights in dinoflagellate evolution

Dinoflagellates are diverse and ecologically important protists characterized by many morphological and molecular traits that set them apart from other eukaryotes. These features include, but are not limited to, massive genomes organized using bacterially-derived histone-like proteins (HLPs) and dinoflagellate viral nucleoproteins (DVNP) rather than histones, and a complex history of photobiology with many independent losses of photosynthesis, numerous cases of serial secondary and tertiary plastid gains, and the presence of horizontally acquired bacterial rhodopsins and type II RuBisCo. Elucidating how this all evolved depends on knowing the phylogenetic relationships between dinoflagellate lineages. Half of these species are heterotrophic, but existing molecular data is strongly biased toward the photosynthetic dinoflagellates due to their amenability to cultivation and prevalence in culture collections. These biases make it impossible to interpret the evolution of photosynthesis, but may also affect phylogenetic inferences that impact our understanding of character evolution. Here, we address this problem by isolating individual cells from the Salish Sea and using single cell, culture-free transcriptomics to expand molecular data for dinoflagellates to include 27 more heterotrophic taxa, resulting in a roughly balanced representation. Using these data, we performed a comprehensive search for proteins involved in chromatin packaging, plastid function, and photoactivity across all dinoflagellates. These searches reveal that 1) photosynthesis was lost at least 21 times, 2) two known types of HLP were horizontally acquired around the same time rather than sequentially as previously thought; 3) multiple rhodopsins are present across the dinoflagellates, acquired multiple times from different donors; 4) kleptoplastic species have nucleus-encoded genes for proteins targeted to their temporary plastids and they are derived from multiple lineages, and 5) warnowiids are the only heterotrophs that retain a whole photosystem, although some photosynthesis-related electron transport genes are widely retained in heterotrophs, likely as part of the iron-sulfur cluster pathway that persists in non-photosynthetic plastids.

Nuclear genome of dinoflagellates: Size variation and insights into evolutionary mechanisms

Recent progress in high-throughput sequencing technologies has dramatically increased availability of genome data for prokaryotes and eukaryotes. Dinoflagellates have distinct chromosomes and a huge genome size, which make their genomic analysis complicated. Here, we reviewed the nuclear genomes of core dinoflagellates, focusing on the genome and cell size. Till now, the genome sizes of several dinoflagellates (more than 25) have been measured by certain methods (e.g., flow cytometry), showing a range of 3-250 pg of genomic DNA per cell. In contrast to their relatively small cell size, their genomes are huge (about 1-80 times the human haploid genome). In the present study, we collected the genome and cell size data of dinoflagellates and compared their relationships. We found that dinoflagellate genome size exhibits a positive correlation with cell size. On the other hand, we recognized that the genome size is not correlated with phylogenetic relatedness. These may be caused by genome duplication, increased gene copy number, repetitive non-coding DNA, transposon expansion, horizontal gene transfer, organelle-to-nucleus gene transfer, and/or mRNA reintegration into the genome. Ultimate verification of these factors as potential causative mechanisms would require sequencing of more dinoflagellate genomes in the future.

Facultative lifestyle drives diversity of coral algal symbionts

The photosynthetic symbionts of corals sustain biodiverse reefs in nutrient-poor, tropical waters. Recent genomic data illuminate the evolution of coral symbionts under genome size constraints and suggest that retention of the facultative lifestyle, widespread among these algae, confers a selective advantage when compared with a strict symbiotic existence. We posit that the coral symbiosis is analogous to a 'bioreactor' that selects winner genotypes and allows them to rise to high numbers in a sheltered habitat prior to release by the coral host. Our observations lead to a novel hypothesis, the 'stepping-stone model', which predicts that local adaptation under both the symbiotic and free-living stages, in a stepwise fashion, accelerates coral alga diversity and the origin of endemic strains and species.

Population genomic analyses reveal that salinity and geographic isolation drive diversification in a free-living protist

Protists make up the vast diversity of eukaryotic life and play a critical role in biogeochemical cycling and in food webs. Because of their small size, cryptic life cycles, and large population sizes, our understanding of speciation in these organisms is very limited. We performed population genomic analyses on 153 strains isolated from eight populations of the recently radiated dinoflagellate genus Apocalathium, to explore the drivers and mechanisms of speciation processes. Species of this genus inhabit both freshwater and saline habitats, lakes and seas, and are found in cold temperate environments across the world. RAD sequencing analyses revealed that the populations were overall highly differentiated, but morphological similarity was not congruent with genetic similarity. While geographic isolation was to some extent coupled to genetic distance, this pattern was not consistent. Instead, we found evidence that the environment, specifically salinity, is a major factor in driving ecological speciation in Apocalathium. While saline populations were unique in loci coupled to genes involved in osmoregulation, freshwater populations appear to lack these. Our study highlights that adaptation to freshwater through loss of osmoregulatory genes may be an important speciation mechanism in free-living aquatic protists.

Flexible B12 ecophysiology of Phaeocystis antarctica due to a fusion B12-independent methionine synthase with widespread homologues

Coastal Antarctic marine ecosystems are significant in carbon cycling because of their intense seasonal phytoplankton blooms. Southern Ocean algae are primarily limited by light and iron (Fe) and can be co-limited by cobalamin (vitamin B12). Micronutrient limitation controls productivity and shapes the composition of blooms which are typically dominated by either diatoms or the haptophyte Phaeocystis antarctica. However, the vitamin requirements and ecophysiology of the keystone species P. antarctica remain poorly characterized. Using cultures, physiological analysis, and comparative omics, we examined the response of P. antarctica to a matrix of Fe-B12 conditions. We show that P. antarctica is not auxotrophic for B12, as previously suggested, and identify mechanisms underlying its B12 response in cultures of predominantly solitary and colonial cells. A combination of proteomics and proteogenomics reveals a B12-independent methionine synthase fusion protein (MetE-fusion) that is expressed under vitamin limitation and interreplaced with the B12-dependent isoform under replete conditions. Database searches return homologues of the MetE-fusion protein in multiple Phaeocystis species and in a wide range of marine microbes, including other photosynthetic eukaryotes with polymorphic life cycles as well as bacterioplankton. Furthermore, we find MetE-fusion homologues expressed in metaproteomic and metatranscriptomic field samples in polar and more geographically widespread regions. As climate change impacts micronutrient availability in the coastal Southern Ocean, our finding that P. antarctica has a flexible B12 metabolism has implications for its relative fitness compared to B12-auxotrophic diatoms and for the detection of B12-stress in a more diverse set of marine microbes.

Massive genome reduction predates the divergence of Symbiodiniaceae dinoflagellates

Dinoflagellates in the family Symbiodiniaceae are taxonomically diverse, predominantly symbiotic lineages that are well-known for their association with corals. The ancestor of these taxa is believed to have been free-living. The establishment of symbiosis (i.e. symbiogenesis) is hypothesized to have occurred multiple times during Symbiodiniaceae evolution, but its impact on genome evolution of these taxa is largely unknown. Among Symbiodiniaceae, the genus Effrenium is a free-living lineage that is phylogenetically positioned between two robustly supported groups of genera within which symbiotic taxa have emerged. The apparent lack of symbiogenesis in Effrenium suggests that the ancestral features of Symbiodiniaceae may have been retained in this lineage. Here, we present de novo assembled genomes (1.2-1.9 Gbp in size) and transcriptome data from three isolates of Effrenium voratum and conduct a comparative analysis that includes 16 Symbiodiniaceae taxa and the other dinoflagellates. Surprisingly, we find that genome reduction, which is often associated with a symbiotic lifestyle, predates the origin of Symbiodiniaceae. The free-living lifestyle distinguishes Effrenium from symbiotic Symbiodiniaceae vis-à-vis their longer introns, more-extensive mRNA editing, fewer (~30%) lineage-specific gene sets, and lower (~10%) level of pseudogenization. These results demonstrate how genome reduction and the adaptation to distinct lifestyles intersect to drive diversification and genome evolution of Symbiodiniaceae.

Contribution of the satellitome to the exceptionally large genome of dinoflagellates: The case of the harmful alga Alexandrium minutum

Dinoflagellates are known to possess an exceptionally large genome organized in permanently condensed chromosomes. Focusing on the contribution of satellite DNA (satDNA) to the whole DNA content of genomes and its potential role in the architecture of the chromosomes, we present the characterization of the satellitome of Alexandriun minutum strain VGO577. To achieve this, we analyzed Illumina reads using graph-based clustering and performed complementary bioinformatic analyses. In this way, we discovered 180 satDNAs occupying 17.38 % of the genome. The 12 most abundant satDNAs represent the half of the satellitome but no satDNA is overrepresented, with the most abundant contributing ∼1.56 % of the genome. The largest repeat unit is 517 bp long but more than the half of the satDNAs (101) have repeat units shorter than 20 bp. We used FISH to map a selected set of 26 satDNAs. Although some satDNAs generate discrete hybridization signals at specific chromosomal locations (hybridization sites, HS), our cytological analysis showed that most satDNAs are dispersed throughout the genome, probably forming short arrays. Two satDNAs co-localize with the 45S rDNA. With the exception of telomeric DNA, no other satDNA yields HS on all chromosomes. In addition, we analyzed nine satDNAs yielding HS in VGO577 in four other A. minutum strains. Polymorphism at the intraspecific level was found for the presence/absence and/or abundance of some satDNAs, suggesting the amplification/deletion of these satDNAs following geographic separation or during culture maintenance of the strains. We also discuss how these results contribute to the understanding of chromosome architecture and evolution of dinoflagellate genomes.

Multi-omics analysis reveals the molecular response to heat stress in a "red tide" dinoflagellate

BACKGROUND

"Red tides" are harmful algal blooms caused by dinoflagellate microalgae that accumulate toxins lethal to other organisms, including humans via consumption of contaminated seafood. These algal blooms are driven by a combination of environmental factors including nutrient enrichment, particularly in warm waters, and are increasingly frequent. The molecular, regulatory, and evolutionary mechanisms that underlie the heat stress response in these harmful bloom-forming algal species remain little understood, due in part to the limited genomic resources from dinoflagellates, complicated by the large sizes of genomes, exhibiting features atypical of eukaryotes.

RESULTS

We present the de novo assembled genome (~ 4.75 Gbp with 85,849 protein-coding genes), transcriptome, proteome, and metabolome from Prorocentrum cordatum, a globally abundant, bloom-forming dinoflagellate. Using axenic algal cultures, we study the molecular mechanisms that underpin the algal response to heat stress, which is relevant to current ocean warming trends. We present the first evidence of a complementary interplay between RNA editing and exon usage that regulates the expression and functional diversity of biomolecules, reflected by reduction in photosynthesis, central metabolism, and protein synthesis. These results reveal genomic signatures and post-transcriptional regulation for the first time in a pelagic dinoflagellate.

CONCLUSIONS

Our multi-omics analyses uncover the molecular response to heat stress in an important bloom-forming algal species, which is driven by complex gene structures in a large, high-G+C genome, combined with multi-level transcriptional regulation. The dynamics and interplay of molecular regulatory mechanisms may explain in part how dinoflagellates diversified to become some of the most ecologically successful organisms on Earth.

Genome-wide identification, comparative analysis and functional roles in flavonoid biosynthesis of cytochrome P450 superfamily in pear (Pyrus spp.)

BACKGROUND

The cytochrome P450 (CYP) superfamily is the largest enzyme metabolism family in plants identified to date, and it is involved in many biological processes, including secondary metabolite biosynthesis, hormone metabolism and stress resistance. However, the P450 gene superfamily has not been well studied in pear (Pyrus spp.).

RESULTS

Here, the comprehensive identification and a comparative analysis of P450 superfamily members were conducted in cultivated and wild pear genomes. In total, 338, 299 and 419 P450 genes were identified in Chinese white pear, European pear and the wild pear, respectively. Based on the phylogenetic analyses, pear P450 genes were divided into ten clans, comprising 48 families. The motif and gene structure analyses further supported this classification. The expansion of the pear P450 gene family was attributed to whole-genome and single-gene duplication events. Several P450 gene clusters were detected, which have resulted from tandem and proximal duplications. Purifying selection was the major force imposed on the long-term evolution of P450 genes. Gene dosage balance, subfunctionalization and neofunctionalization jointly drove the retention and functional diversification of P450 gene pairs. Based on the association analysis between transcriptome expression profiles and flavonoid content during fruit development, three candidate genes were identified as being closely associated with the flavonoid biosynthesis, and the expression of one gene was further verified using qRT-PCR and its function was validated through transient transformation in pear fruit.

CONCLUSIONS

The study results provide insights into the evolution and biological functions of P450 genes in pear.

Gene duplication is the primary driver of intraspecific genomic divergence in coral algal symbionts

Dinoflagellates in the order Suessiales include the family Symbiodiniaceae, which have essential roles as photosymbionts in corals, and their cold-adapted sister group, Polarella glacialis. These diverse taxa exhibit extensive genomic divergence, although their genomes are relatively small (haploid size < 3 Gbp) when compared with most other free-living dinoflagellates. Different strains of Symbiodiniaceae form symbiosis with distinct hosts and exhibit different regimes of gene expression, but intraspecific whole-genome divergence is poorly understood. Focusing on three Symbiodiniaceae species (the free-living Effrenium voratum and the symbiotic Symbiodinium microadriaticum and Durusdinium trenchii) and the free-living outgroup P. glacialis, for which whole-genome data from multiple isolates are available, we assessed intraspecific genomic divergence with respect to sequence and structure. Our analysis, based on alignment and alignment-free methods, revealed a greater extent of intraspecific sequence divergence in Symbiodiniaceae than in P. glacialis. Our results underscore the role of gene duplication in generating functional innovation, with a greater prevalence of tandemly duplicated single-exon genes observed in the genomes of free-living species than in symbionts. These results demonstrate the remarkable intraspecific genomic divergence in dinoflagellates under the constraint of reduced genome sizes, shaped by genetic duplications and symbiogenesis events during the diversification of Symbiodiniaceae.

Extreme environments offer an unprecedented opportunity to understand microbial eukaryotic ecology, evolution, and genome biology

Research in extreme environments has substantially expanded our understanding of the ecology and evolution of life on Earth, but a major group of organisms has been largely overlooked: microbial eukaryotes (i.e., protists). In this Perspective, we summarize data from over 80 studies of protists in extreme environments and identify focal lineages that are of significant interest for further study, including clades within Echinamoebida, Heterolobosea, Radiolaria, Haptophyta, Oomycota, and Cryptophyta. We argue that extreme environments are prime sampling targets to fill gaps in the eukaryotic tree of life and to increase our understanding of the ecology, metabolism, genome architecture, and evolution of eukaryotic life.

Genomic copy number variability at the genus, species and population levels impacts in situ ecological analyses of dinoflagellates and harmful algal blooms

The application of meta-barcoding, qPCR, and metagenomics to aquatic eukaryotic microbial communities requires knowledge of genomic copy number variability (CNV). CNV may be particularly relevant to functional genes, impacting dosage and expression, yet little is known of the scale and role of CNV in microbial eukaryotes. Here, we quantify CNV of rRNA and a gene involved in Paralytic Shellfish Toxin (PST) synthesis (sxtA4), in 51 strains of 4 Alexandrium (Dinophyceae) species. Genomes varied up to threefold within species and ~7-fold amongst species, with the largest (A. pacificum, 130 ± 1.3 pg cell-1 /~127 Gbp) in the largest size category of any eukaryote. Genomic copy numbers (GCN) of rRNA varied by 6 orders of magnitude amongst Alexandrium (102- 108 copies cell-1) and were significantly related to genome size. Within the population CNV of rRNA was 2 orders of magnitude (105 - 107 cell-1) in 15 isolates from one population, demonstrating that quantitative data based on rRNA genes needs considerable caution in interpretation, even if validated against locally isolated strains. Despite up to 30 years in laboratory culture, rRNA CNV and genome size variability were not correlated with time in culture. Cell volume was only weakly associated with rRNA GCN (20-22% variance explained across dinoflagellates, 4% in Gonyaulacales). GCN of sxtA4 varied from 0-102 copies cell-1, was significantly related to PSTs (ng cell-1), displaying a gene dosage effect modulating PST production. Our data indicate that in dinoflagellates, a major marine eukaryotic group, low-copy functional genes are more reliable and informative targets for quantification of ecological processes than unstable rRNA genes.

Endogenous viral elements reveal associations between a non-retroviral RNA virus and symbiotic dinoflagellate genomes

Endogenous viral elements (EVEs) offer insight into the evolutionary histories and hosts of contemporary viruses. This study leveraged DNA metagenomics and genomics to detect and infer the host of a non-retroviral dinoflagellate-infecting +ssRNA virus (dinoRNAV) common in coral reefs. As part of the Tara Pacific Expedition, this study surveyed 269 newly sequenced cnidarians and their resident symbiotic dinoflagellates (Symbiodiniaceae), associated metabarcodes, and publicly available metagenomes, revealing 178 dinoRNAV EVEs, predominantly among hydrocoral-dinoflagellate metagenomes. Putative associations between Symbiodiniaceae and dinoRNAV EVEs were corroborated by the characterization of dinoRNAV-like sequences in 17 of 18 scaffold-scale and one chromosome-scale dinoflagellate genome assembly, flanked by characteristically cellular sequences and in proximity to retroelements, suggesting potential mechanisms of integration. EVEs were not detected in dinoflagellate-free (aposymbiotic) cnidarian genome assemblies, including stony corals, hydrocorals, jellyfish, or seawater. The pervasive nature of dinoRNAV EVEs within dinoflagellate genomes (especially Symbiodinium), as well as their inconsistent within-genome distribution and fragmented nature, suggest ancestral or recurrent integration of this virus with variable conservation. Broadly, these findings illustrate how +ssRNA viruses may obscure their genomes as members of nested symbioses, with implications for host evolution, exaptation, and immunity in the context of reef health and disease.

Convergent evolution and horizontal gene transfer in Arctic Ocean microalgae

Microbial communities in the world ocean are affected strongly by oceanic circulation, creating characteristic marine biomes. The high connectivity of most of the ocean makes it difficult to disentangle selective retention of colonizing genotypes (with traits suited to biome specific conditions) from evolutionary selection, which would act on founder genotypes over time. The Arctic Ocean is exceptional with limited exchange with other oceans and ice covered since the last ice age. To test whether Arctic microalgal lineages evolved apart from algae in the global ocean, we sequenced four lineages of microalgae isolated from Arctic waters and sea ice. Here we show convergent evolution and highlight geographically limited HGT as an ecological adaptive force in the form of PFAM complements and horizontal acquisition of key adaptive genes. Notably, ice-binding proteins were acquired and horizontally transferred among Arctic strains. A comparison with Tara Oceans metagenomes and metatranscriptomes confirmed mostly Arctic distributions of these IBPs. The phylogeny of Arctic-specific genes indicated that these events were independent of bacterial-sourced HGTs in Antarctic Southern Ocean microalgae.

Mutational Bias and Natural Selection Driving the Synonymous Codon Usage of Single-Exon Genes in Rice (Oryza sativa L.)

The relative abundance of single-exon genes (SEGs) in higher plants is perplexing. Uncovering the synonymous codon usage pattern of SEGs will benefit for further understanding their underlying evolutionary mechanism in plants. Using internal correspondence analysis (ICA), we reveal a significant difference in synonymous codon usage between SEGs and multiple-exon genes (MEGs) in rice. But the effect is weak, accounting for only 2.61% of the total codon usage variability. SEGs and MEGs contain remarkably different base compositions, and are under clearly differential selective constraints, with the former having higher GC content, and evolving relatively faster during evolution. In the group of SEGs, the variability in synonymous codon usage among genes is partially due to the variations in GC content, gene function, and gene expression level, which accounts for 22.03%, 5.99%, and 3.32% of the total codon usage variability, respectively. Therefore, mutational bias and natural selection should work on affecting the synonymous codon usage of SEGs in rice. These findings may deepen our knowledge for the mechanisms of origination, differentiation and regulation of SEGs in plants.

Genetic and Structural Diversity of Prokaryotic Ice-Binding Proteins from the Central Arctic Ocean

Ice-binding proteins (IBPs) are a group of ecologically and biotechnologically relevant enzymes produced by psychrophilic organisms. Although putative IBPs containing the domain of unknown function (DUF) 3494 have been identified in many taxa of polar microbes, our knowledge of their genetic and structural diversity in natural microbial communities is limited. Here, we used samples from sea ice and sea water collected in the central Arctic Ocean as part of the MOSAiC expedition for metagenome sequencing and the subsequent analyses of metagenome-assembled genomes (MAGs). By linking structurally diverse IBPs to particular environments and potential functions, we reveal that IBP sequences are enriched in interior ice, have diverse genomic contexts and cluster taxonomically. Their diverse protein structures may be a consequence of domain shuffling, leading to variable combinations of protein domains in IBPs and probably reflecting the functional versatility required to thrive in the extreme and variable environment of the central Arctic Ocean.

Expanding the taxonomic and environmental extent of an underexplored carbon metabolism-oxalotrophy

Oxalate serves various functions in the biological processes of plants, fungi, bacteria, and animals. It occurs naturally in the minerals weddellite and whewellite (calcium oxalates) or as oxalic acid. The environmental accumulation of oxalate is disproportionately low compared to the prevalence of highly productive oxalogens, namely plants. It is hypothesized that oxalotrophic microbes limit oxalate accumulation by degrading oxalate minerals to carbonates via an under-explored biogeochemical cycle known as the oxalate-carbonate pathway (OCP). Neither the diversity nor the ecology of oxalotrophic bacteria is fully understood. This research investigated the phylogenetic relationships of the bacterial genes oxc, frc, oxdC, and oxlT, which encode key enzymes for oxalotrophy, using bioinformatic approaches and publicly available omics datasets. Phylogenetic trees of oxc and oxdC genes demonstrated grouping by both source environment and taxonomy. All four trees included genes from metagenome-assembled genomes (MAGs) that contained novel lineages and environments for oxalotrophs. In particular, sequences of each gene were recovered from marine environments. These results were supported with marine transcriptome sequences and description of key amino acid residue conservation. Additionally, we investigated the theoretical energy yield from oxalotrophy across marine-relevant pressure and temperature conditions and found similar standard state Gibbs free energy to "low energy" marine sediment metabolisms, such as anaerobic oxidation of methane coupled to sulfate reduction. These findings suggest further need to understand the role of bacterial oxalotrophy in the OCP, particularly in marine environments, and its contribution to global carbon cycling.

Strong population genomic structure of the toxic dinoflagellate Alexandrium minutum inferred from meta-transcriptome samples

Despite theoretical expectations, marine microeukaryote population are often highly structured and the mechanisms behind such patterns remain to be elucidated. These organisms display huge census population sizes, yet genotyping usually requires clonal strains originating from single cells, hindering proper population sampling. Estimating allelic frequency directly from population wide samples, without any isolation step, offers an interesting alternative. Here, we validate the use of meta-transcriptome environmental samples to determine the population genetic structure of the dinoflagellate Alexandrium minutum. Strain and meta-transcriptome based results both indicated a strong genetic structure for A. minutum in Western Europe, to the level expected between cryptic species. The presence of numerous private alleles, and even fixed polymorphism, would indicate ancient divergence and absence of gene flow between populations. Single nucleotide polymorphisms (SNPs) displaying strong allele frequency differences were distributed throughout the genome, which might indicate pervasive selection from standing genetic variation (soft selective sweeps). However, a few genomic regions displayed extremely low diversity that could result from the fixation of adaptive de novo mutations (hard selective sweeps) within the populations.

Genetic association of toxin production in the dinoflagellate Alexandrium minutum

Dinoflagellates of the genus Alexandrium are responsible for harmful algal blooms and produce paralytic shellfish toxins (PSTs). Their very large and complex genomes make it challenging to identify the genes responsible for toxin synthesis. A family-based genomic association study was developed to determine the inheritance of toxin production in Alexandrium minutum and identify genomic regions linked to this production. We show that the ability to produce toxins is inheritable in a Mendelian way, while the heritability of the toxin profile is more complex. We developed the first dinoflagellate genetic linkage map. Using this map, several major results were obtained: 1. A genomic region related to the ability to produce toxins was identified. 2. This region does not contain any polymorphic sxt genes, known to be involved in toxin production in cyanobacteria. 3. The sxt genes, known to be present in a single cluster in cyanobacteria, are scattered on different linkage groups in A. minutum. 4. The expression of two sxt genes not assigned to any linkage group, sxtI and sxtG, may be regulated by the genomic region related to the ability to produce toxins. Our results provide new insights into the organization of toxicity-related genes in A. minutum, suggesting a dissociated genetic mechanism for the production of the different analogues and the ability to produce toxins. However, most of the newly identified genes remain unannotated. This study therefore proposes new candidate genes to be further explored to understand how dinoflagellates synthesize their toxins.

Multiple waves of viral invasions in Symbiodiniaceae algal genomes

Dinoflagellates from the family Symbiodiniaceae are phototrophic marine protists that engage in symbiosis with diverse hosts. Their large and distinct genomes are characterized by pervasive gene duplication and large-scale retroposition events. However, little is known about the role and scale of horizontal gene transfer (HGT) in the evolution of this algal family. In other dinoflagellates, high levels of HGTs have been observed, linked to major genomic transitions, such as the appearance of a viral-acquired nucleoprotein that originated via HGT from a large DNA algal virus. Previous work showed that Symbiodiniaceae from different hosts are actively infected by viral groups, such as giant DNA viruses and ssRNA viruses, that may play an important role in coral health. Latent viral infections may also occur, whereby viruses could persist in the cytoplasm or integrate into the host genome as a provirus. This hypothesis received experimental support; however, the cellular localization of putative latent viruses and their taxonomic affiliation are still unknown. In addition, despite the finding of viral sequences in some genomes of Symbiodiniaceae, viral origin, taxonomic breadth, and metabolic potential have not been explored. To address these questions, we searched for putative viral-derived proteins in thirteen Symbiodiniaceae genomes. We found fifty-nine candidate viral-derived HGTs that gave rise to twelve phylogenies across ten genomes. We also describe the taxonomic affiliation of these virus-related sequences, their structure, and their genomic context. These results lead us to propose a model to explain the origin and fate of Symbiodiniaceae viral acquisitions.

The enigmatic clock of dinoflagellates, is it unique?

Dinoflagellate clocks are unique as they show no resemblance to any known model eukaryotic or prokaryotic clock architecture. Dinoflagellates are unicellular, photosynthetic, primarily marine eukaryotes are known for their unique biology and rhythmic physiology. Their physiological rhythms are driven by an internal oscillator whose molecular underpinnings are yet unknown. One of the primary reasons that slowed the progression of their molecular studies is their extremely large and repetitive genomes. Dinoflagellates are primary contributors to the global carbon cycle and oxygen levels, therefore, comprehending their internal clock architecture and its interaction with their physiology becomes a subject of utmost importance. The advent of high throughput Omics technology provided the momentum to understand the molecular architecture and functioning of the dinoflagellate clocks. We use these extensive databases to perform meta-analysis to reveal the status of clock components in dinoflagellates. In this article, we will delve deep into the various “Omics” studies that catered to various breakthroughs in the field of circadian biology in these organisms that were not possible earlier. The overall inference from these omics studies points toward an uncommon eukaryotic clock model, which can provide promising leads to understand the evolution of molecular clocks.

Active meiosis during dinoflagellate blooms: A 'sex for proliferation' hypothesis

In dinoflagellates, sexual reproduction is best known to be induced by adverse environmental conditions and culminate in encystment for survival ('sex for encystment'). Although increasing laboratory observations indicate that sex can lead to production of vegetative cells bypassing encystment, the occurrence of this alternative pathway in natural populations and its ecological roles remain poorly understood. Here we report evidence that sex in dinoflagellates can potentially be an instrument for bloom proliferation or extension. By bloom metatranscriptome profiling, we documented elevated expression of meiosis genes in two evolutionarily distinct species (Prorocentrum shikokuense and Karenia mikimotoi) during bloom, a timing unexpected of the 'sex for encystment' scenario. To link these genes to meiosis, we induced encystment and cyst germination in the cyst-forming species Scrippsiella acuminata, and found that five of these genes were upregulated during cyst germination, when meiosis occurs. Integrating data from all three species revealed that SPO11, MND1, and DMC1 were likely common between cyst-forming and non-encysting sex in dinoflagellates. Furthermore, flow cytometric analyses revealed consecutive rounds of DNA halving during blooms of P. shikokuense and K. mikimotoi, evidencing meiosis. These data provided novel evidence that sexual reproduction in dinoflagellates might serve to promote cell proliferation, and along with the consequent enhancement of genetic diversity facilitating resistance against pathogens and environmental stress, to boost or extend a bloom ('sex for proliferation'). The putative meiosis-specific genes and insights reported here will prove to be helpful for rigorously testing the hypothesis and addressing whether the two modes of sex are genetically predisposed (i.e. species-specific) or environmentally induced (switchable within species), and if the latter what triggers the switch.

Improved Cladocopium goreaui Genome Assembly Reveals Features of a Facultative Coral Symbiont and the Complex Evolutionary History of Dinoflagellate Genes

Dinoflagellates of the family Symbiodiniaceae are crucial photosymbionts in corals and other marine organisms. Of these, Cladocopium goreaui is one of the most dominant symbiont species in the Indo-Pacific. Here, we present an improved genome assembly of C. goreaui combining new long-read sequence data with previously generated short-read data. Incorporating new full-length transcripts to guide gene prediction, the C. goreaui genome (1.2 Gb) exhibits a high extent of completeness (82.4% based on BUSCO protein recovery) and better resolution of repetitive sequence regions; 45,322 gene models were predicted, and 327 putative, topologically associated domains of the chromosomes were identified. Comparison with other Symbiodiniaceae genomes revealed a prevalence of repeats and duplicated genes in C. goreaui, and lineage-specific genes indicating functional innovation. Incorporating 2,841,408 protein sequences from 96 taxonomically diverse eukaryotes and representative prokaryotes in a phylogenomic approach, we assessed the evolutionary history of C. goreaui genes. Of the 5246 phylogenetic trees inferred from homologous protein sets containing two or more phyla, 35–36% have putatively originated via horizontal gene transfer (HGT), predominantly (19–23%) via an ancestral Archaeplastida lineage implicated in the endosymbiotic origin of plastids: 10–11% are of green algal origin, including genes encoding photosynthetic functions. Our results demonstrate the utility of long-read sequence data in resolving structural features of a dinoflagellate genome, and highlight how genetic transfer has shaped genome evolution of a facultative symbiont, and more broadly of dinoflagellates.

An overview of transcription in dinoflagellates

Dinoflagellates are a vital diverse family of unicellular algae widespread in various aquatic environments. Typically large genomes and permanently condensed chromosomes without histones make these organisms unique among eukaryotes in terms of chromatin structure and gene expression. Genomic and transcriptomic sequencing projects have provided new insight into the genetic foundation of dinoflagellate behaviors. Genes in tandem arrays, trans-splicing of mRNAs and lower levels of transcriptional regulation compared to other eukaryotes all contribute to the differences seen. Here we present a general overview of transcription in dinoflagellates based on previously described work.

The role of zinc in the adaptive evolution of polar phytoplankton

Zinc is an essential trace metal for oceanic primary producers with the highest concentrations in polar oceans. However, its role in the biological functioning and adaptive evolution of polar phytoplankton remains enigmatic. Here, we have applied a combination of evolutionary genomics, quantitative proteomics, co-expression analyses and cellular physiology to suggest that model polar phytoplankton species have a higher demand for zinc because of elevated cellular levels of zinc-binding proteins. We propose that adaptive expansion of regulatory zinc-finger protein families, co-expanded and co-expressed zinc-binding proteins families involved in photosynthesis and growth in these microalgal species and their natural communities were identified to be responsible for the higher zinc demand. The expression of their encoding genes in eukaryotic phytoplankton metatranscriptomes from pole-to-pole was identified to correlate not only with dissolved zinc concentrations in the upper ocean but also with temperature, suggesting that environmental conditions of polar oceans are responsible for an increased demand of zinc. These results suggest that zinc plays an important role in supporting photosynthetic growth in eukaryotic polar phytoplankton and that this has been critical for algal colonization of low-temperature polar oceans.

Hypothesis: Trans-splicing Generates Evolutionary Novelty in the Photosynthetic Amoeba Paulinella

Plastid primary endosymbiosis has occurred twice, once in the Archaeplastida ancestor and once in the Paulinella (Rhizaria) lineage. Both events precipitated massive evolutionary changes, including the recruitment and activation of genes that are horizontally acquired (HGT) and the redeployment of existing genes and pathways in novel contexts. Here we address the latter aspect in Paulinella micropora KR01 (hereafter, KR01) that has independently evolved spliced leader (SL) trans-splicing (SLTS) of nuclear-derived transcripts. We investigated the role of this process in gene regulation, novel gene origination, and endosymbiont integration. Our analysis shows that 20% of KR01 genes give rise to transcripts with at least one (but in some cases, multiple) sites of SL addition. This process, which often occurs at canonical cis-splicing acceptor sites (internal introns), results in shorter transcripts that may produce 5'-truncated proteins with novel functions. SL-truncated transcripts fall into four categories that may show: (i) altered protein localization, (ii) altered protein function, structure, or regulation, (iii) loss of valid alternative start codons, preventing translation, or (iv) multiple SL addition sites at the 5'-terminus. The SL RNA genes required for SLTS are putatively absent in the heterotrophic sister lineage of photosynthetic Paulinella species. Moreover, a high proportion of transcripts derived from genes of endosymbiotic gene transfer (EGT) and HGT origin contain SL sequences. We hypothesize that truncation of transcripts by SL addition may facilitate the generation and expression of novel gene variants and that SLTS may have enhanced the activation and fixation of foreign genes in the host genome of the photosynthetic lineages, playing a key role in primary endosymbiont integration.

Recent developments in horizontal gene transfer with the adaptive innovation of fermented foods

Horizontal gene transfer (HGT) has contributed significantly to the adaptability of bacteria, yeast and mold in fermented foods, whose evidence has been found in several fermented foods. Although not every HGT has biological significance, it plays an important role in improving the quality of fermented foods. In this review, how HGT facilitated microbial domestication and adaptive evolution in fermented foods was discussed. HGT can assist in the industrial innovation of fermented foods, and this adaptive evolution strategy can improve the quality of fermented foods. Additionally, the mechanism underlying HGT in fermented foods were analyzed. Furthermore, the critical bottlenecks involved in optimizing HGT during the production of fermented foods and strategies for optimizing HGT were proposed. Finally, the prospect of HGT for promoting the industrial innovation of fermented foods was highlighted. The comprehensive report on HGT in fermented foods provides a new trend for domesticating preferable starters for food fermentation, thus optimizing the quality and improving the industrial production of fermented foods.

Salinity tolerance mechanisms of an Arctic Pelagophyte using comparative transcriptomic and gene expression analysis

Little is known at the transcriptional level about microbial eukaryotic adaptations to short-term salinity change. Arctic microalgae are exposed to low salinity due to sea-ice melt and higher salinity with brine channel formation during freeze-up. Here, we investigate the transcriptional response of an ice-associated microalgae over salinities from 45 to 8. Our results show a bracketed response of differential gene expression when the cultures were exposed to progressively decreasing salinity. Key genes associated with salinity changes were involved in specific metabolic pathways, transcription factors and regulators, protein kinases, carbohydrate active enzymes, and inorganic ion transporters. The pelagophyte seemed to use a strategy involving overexpression of Na+-H+ antiporters and Na+ -Pi symporters as salinity decreases, but the K+ channel complex at higher salinities. Specific adaptation to cold saline arctic conditions was seen with differential expression of several antifreeze proteins, an ice-binding protein and an acyl-esterase involved in cold adaptation.

Metagenome-assembled genomes of phytoplankton microbiomes from the Arctic and Atlantic Oceans

BACKGROUND

Phytoplankton communities significantly contribute to global biogeochemical cycles of elements and underpin marine food webs. Although their uncultured genomic diversity has been estimated by planetary-scale metagenome sequencing and subsequent reconstruction of metagenome-assembled genomes (MAGs), this approach has yet to be applied for complex phytoplankton microbiomes from polar and non-polar oceans consisting of microbial eukaryotes and their associated prokaryotes.

RESULTS

Here, we have assembled MAGs from chlorophyll a maximum layers in the surface of the Arctic and Atlantic Oceans enriched for species associations (microbiomes) with a focus on pico- and nanophytoplankton and their associated heterotrophic prokaryotes. From 679 Gbp and estimated 50 million genes in total, we recovered 143 MAGs of medium to high quality. Although there was a strict demarcation between Arctic and Atlantic MAGs, adjacent sampling stations in each ocean had 51-88% MAGs in common with most species associations between Prasinophytes and Proteobacteria. Phylogenetic placement revealed eukaryotic MAGs to be more diverse in the Arctic whereas prokaryotic MAGs were more diverse in the Atlantic Ocean. Approximately 70% of protein families were shared between Arctic and Atlantic MAGs for both prokaryotes and eukaryotes. However, eukaryotic MAGs had more protein families unique to the Arctic whereas prokaryotic MAGs had more families unique to the Atlantic.

CONCLUSION

Our study provides a genomic context to complex phytoplankton microbiomes to reveal that their community structure was likely driven by significant differences in environmental conditions between the polar Arctic and warm surface waters of the tropical and subtropical Atlantic Ocean. Video Abstract.

Alignment-Free Analysis of Whole-Genome Sequences From Symbiodiniaceae Reveals Different Phylogenetic Signals in Distinct Regions

Dinoflagellates of the family Symbiodiniaceae are predominantly essential symbionts of corals and other marine organisms. Recent research reveals extensive genome sequence divergence among Symbiodiniaceae taxa and high phylogenetic diversity hidden behind subtly different cell morphologies. Using an alignment-free phylogenetic approach based on sub-sequences of fixed length k (i.e. k-mers), we assessed the phylogenetic signal among whole-genome sequences from 16 Symbiodiniaceae taxa (including the genera of Symbiodinium, Breviolum, Cladocopium, Durusdinium and Fugacium) and two strains of Polarella glacialis as outgroup. Based on phylogenetic trees inferred from k-mers in distinct genomic regions (i.e. repeat-masked genome sequences, protein-coding sequences, introns and repeats) and in protein sequences, the phylogenetic signal associated with protein-coding DNA and the encoded amino acids is largely consistent with the Symbiodiniaceae phylogeny based on established markers, such as large subunit rRNA. The other genome sequences (introns and repeats) exhibit distinct phylogenetic signals, supporting the expected differential evolutionary pressure acting on these regions. Our analysis of conserved core k-mers revealed the prevalence of conserved k-mers (>95% core 23-mers among all 18 genomes) in annotated repeats and non-genic regions of the genomes. We observed 180 distinct repeat types that are significantly enriched in genomes of the symbiotic versus free-living Symbiodinium taxa, suggesting an enhanced activity of transposable elements linked to the symbiotic lifestyle. We provide evidence that representation of alignment-free phylogenies as dynamic networks enhances the ability to generate new hypotheses about genome evolution in Symbiodiniaceae. These results demonstrate the potential of alignment-free phylogenetic methods as a scalable approach for inferring comprehensive, unbiased whole-genome phylogenies of dinoflagellates and more broadly of microbial eukaryotes.

First record of the spatial organization of the nucleosome-less chromatin of dinoflagellates: The nonrandom distribution of microsatellites and bipolar arrangement of telomeres in the nucleus of Gambierdiscus australes (Dinophyceae)

Dinoflagellates are a group of protists whose exceptionally large genome is organized in permanently condensed nucleosome-less chromosomes. In this study, we examined the potential role of repetitive DNAs in both the structure of dinoflagellate chromosomes and the architecture of the dinoflagellate nucleus. Non-denaturing fluorescent in situ hybridization (ND-FSH) was used to determine the abundance and physical distribution of telomeric DNA and 16 microsatellites (1- to 4-bp repeats) in the nucleus of Gambierdiscus australes. The results showed an increased relative abundance of the different microsatellite motifs with increasing GC content. Two ND-FISH probes, (A)₂₀ and (AAT)₅ , did not yield signals whereas the remainder revealed a dispersed but nonrandom distribution of the microsatellites, mostly in clusters. The bean-shaped interphase nucleus of G. australes contained a region with a high density of trinucleotides. This nuclear compartment was located between the nucleolar organizer region (NOR), located on the concave side of the nucleus, and the convex side. Telomeric DNA was grouped in multiple foci and distributed in two polarized compartments: one associated with the NOR and the other peripherally located along the convex side of the nucleus. Changes in the position of the telomeres during cell division evidenced their dynamic distribution and thus that of the chromosomes during dinomitosis. These insights into the spatial organization of microsatellites and telomeres and thus into the nuclear architecture of G. australes will open up new lines of research into the structure and function of the nucleosome-less chromatin of dinoflagellates.

Diatoms and Their Microbiomes in Complex and Changing Polar Oceans

Diatoms, a key group of polar marine microbes, support highly productive ocean ecosystems. Like all life on earth, diatoms do not live in isolation, and they are therefore under constant biotic and abiotic pressures which directly influence their evolution through natural selection. Despite their importance in polar ecosystems, polar diatoms are understudied compared to temperate species. The observed rapid change in the polar climate, especially warming, has created increased research interest to discover the underlying causes and potential consequences on single species to entire ecosystems. Next-Generation Sequencing (NGS) technologies have greatly expanded our knowledge by revealing the molecular underpinnings of physiological adaptations to polar environmental conditions. Their genomes, transcriptomes, and proteomes together with the first eukaryotic meta-omics data of surface ocean polar microbiomes reflect the environmental pressures through adaptive responses such as the expansion of protein families over time as a consequence of selection. Polar regions and their microbiomes are inherently connected to climate cycles and their feedback loops. An integrated understanding built on "omics" resources centered around diatoms as key primary producers will enable us to reveal unifying concepts of microbial co-evolution and adaptation in polar oceans. This knowledge, which aims to relate past environmental changes to specific adaptations, will be required to improve climate prediction models for polar ecosystems because it provides a unifying framework of how interacting and co-evolving biological communities might respond to future environmental change.

Dinoflagellate Phosphopantetheinyl Transferase (PPTase) and Thiolation Domain Interactions Characterized Using a Modified Indigoidine Synthesizing Reporter

Photosynthetic dinoflagellates synthesize many toxic but also potential therapeutic compounds therapeutics via polyketide/non-ribosomal peptide synthesis, a common means of producing natural products in bacteria and fungi. Although canonical genes are identifiable in dinoflagellate transcriptomes, the biosynthetic pathways are obfuscated by high copy numbers and fractured synteny. This study focuses on the carrier domains that scaffold natural product synthesis (thiolation domains) and the phosphopantetheinyl transferases (PPTases) that thiolate these carriers. We replaced the thiolation domain of the indigoidine producing BpsA gene from Streptomyces lavendulae with those of three multidomain dinoflagellate transcripts and coexpressed these constructs with each of three dinoflagellate PPTases looking for specific pairings that would identify distinct pathways. Surprisingly, all three PPTases were able to activate all the thiolation domains from one transcript, although with differing levels of indigoidine produced, demonstrating an unusual lack of specificity. Unfortunately, constructs with the remaining thiolation domains produced almost no indigoidine and the thiolation domain for lipid synthesis could not be expressed in E. coli. These results combined with inconsistent protein expression for different PPTase/thiolation domain pairings present technical hurdles for future work. Despite these challenges, expression of catalytically active dinoflagellate proteins in E. coli is a novel and useful tool going forward.

Gene clusters for biosynthesis of mycosporine-like amino acids in dinoflagellate nuclear genomes: Possible recent horizontal gene transfer between species of Symbiodiniaceae (Dinophyceae)

Global warming increases the temperature of the ocean surface, which can disrupt dinoflagellate-coral symbioses and result in coral bleaching. Photosynthetic dinoflagellates of the family Symbiodiniaceae include bleaching-tolerant and bleaching-sensitive coral symbionts. Therefore, understanding the molecular mechanisms for changing symbiont diversity is potentially useful to assist recovery of coral holobionts (corals and their associated microbes, including multiple species of Symbiodiniaceae), although sexual reproduction has not been observed in the Symbiodiniaceae. Recent molecular phylogenetic analyses estimate that the Symbiodiniaceae appeared 160 million years ago and diversified into 15 groups, five genera of which now have available draft genomes (i.e., Symbiodinium, Durusdinium, Breviolum, Fugacium, and Cladocopium). Comparative genomic analyses have suggested that crown groups have fewer gene families than early-diverging groups, although many genes that were probably acquired via gene duplications and horizontal gene transfers (HGTs) have been found in each decoded genome. Because UV stress is likely a contributor to coral bleaching, and because the highly conserved gene cluster for mycosporine-like amino acid (MAA) biosynthesis has been found in thermal-tolerant symbiont genomes, I reviewed genomic features of the Symbiodiniaceae, focusing on possible acquisition of a biosynthetic gene cluster for MAAs, which absorb UV radiation. On the basis of highly conserved noncoding sequences, I hypothesized that HGTs have occurred among members of the Symbiodiniaceae and have contributed to the diversification of Symbiodiniaceae-host relationships. Finally, I proposed that bleaching tolerance may be strengthened by multiple MAAs from both symbiotic dinoflagellates and corals.

Codon Usage Bias in Phytoplankton

Non-random usage of synonymous codons, known as “codon bias”, has been described in many organisms, from bacteria to Drosophila, but little is known about it in phytoplankton. This phenomenon is thought to be driven by selection for translational efficiency. As the efficacy of selection is proportional to the effective population size, species with large population sizes, such as phytoplankton, are expected to have strong codon bias. To test this, we measured codon bias in 215 strains from Haptophyta, Chlorophyta, Ochrophyta (except diatoms that were studied previously), Dinophyta, Cryptophyta, Ciliophora, unicellular Rhodophyta and Chlorarachniophyta. Codon bias is modest in most groups, despite the astronomically large population sizes of marine phytoplankton. The strength of the codon bias, measured with the effective number of codons, is the strongest in Haptophyta and the weakest in Chlorarachniophyta. The optimal codons are GC-ending in most cases, but several shifts to AT-ending codons were observed (mainly in Ochrophyta and Ciliophora). As it takes a long time to reach a new equilibrium after such shifts, species having AT-ending codons show a lower frequency of optimal codons compared to other species. Genetic diversity, calculated for species with more than three strains sequenced, is modest, indicating that the effective population sizes are many orders of magnitude lower than the astronomically large census population sizes, which helps to explain the modest codon bias in marine phytoplankton. This study represents the first comparative analysis of codon bias across multiple major phytoplankton groups.

Spatial organization of dinoflagellate genomes: Novel insights and remaining critical questions

As is true for many other aspects, genome architecture, evolution, and function in dinoflagellates are enigmatic and, in the meantime, continuous inspiration for scientific quests. Recent third-generation sequencing and Hi-C linkage analyses brought new insights into the spatial organization of symbiodiniacean genomes, revealing the topologically associated domains, discrete gene clusters and their cis and trans orientations, and relationships with transcription. Where do these new findings bring us in dinoflagellate genomics? Here, we aim to place these new results in the backdrop of the long history of research on this topic and in the context of what critical questions remain to be pursued in the future. The new data suggest, pending verification of other complete chromosome assemblies, a potential evolutionary trend in chromosome number decrease and length increase within the Symbiodiniaceae. While questions remain about the mechanics of the three-dimensional chromosome structure and cell cycle-related DNA replication, the mechanisms of gene transcription and genome size evolution, these latest findings set new starting points for further inquiries.

Bioluminescence and Photoreception in Unicellular Organisms: Light-Signalling in a Bio-Communication Perspective

Bioluminescence, the emission of light catalysed by luciferases, has evolved in many taxa from bacteria to vertebrates and is predominant in the marine environment. It is now well established that in animals possessing a nervous system capable of integrating light stimuli, bioluminescence triggers various behavioural responses and plays a role in intra- or interspecific visual communication. The function of light emission in unicellular organisms is less clear and it is currently thought that it has evolved in an ecological framework, to be perceived by visual animals. For example, while it is thought that bioluminescence allows bacteria to be ingested by zooplankton or fish, providing them with favourable conditions for growth and dispersal, the luminous flashes emitted by dinoflagellates may have evolved as an anti-predation system against copepods. In this short review, we re-examine this paradigm in light of recent findings in microorganism photoreception, signal integration and complex behaviours. Numerous studies show that on the one hand, bacteria and protists, whether autotrophs or heterotrophs, possess a variety of photoreceptors capable of perceiving and integrating light stimuli of different wavelengths. Single-cell light-perception produces responses ranging from phototaxis to more complex behaviours. On the other hand, there is growing evidence that unicellular prokaryotes and eukaryotes can perform complex tasks ranging from habituation and decision-making to associative learning, despite lacking a nervous system. Here, we focus our analysis on two taxa, bacteria and dinoflagellates, whose bioluminescence is well studied. We propose the hypothesis that similar to visual animals, the interplay between light-emission and reception could play multiple roles in intra- and interspecific communication and participate in complex behaviour in the unicellular world.

Insights into Alexandrium minutum Nutrient Acquisition, Metabolism and Saxitoxin Biosynthesis through Comprehensive Transcriptome Survey

Simple Summary

Alexandrium minutum is one of the causing organisms for the occurrence of harmful algae bloom (HABs) in marine ecosystems. This species produces saxitoxin, one of the deadliest neurotoxins which can cause human mortality. However, molecular information such as genes and proteins catalog on this species is still lacking. Therefore, this study has successfully characterized several new molecular mechanisms regarding A. minutum environmental adaptation and saxitoxin biosynthesis. Ultimately, this study provides a valuable resource for facilitating future dinoflagellates’ molecular response to environmental changes.

Abstract

The toxin-producing dinoflagellate Alexandrium minutum is responsible for the outbreaks of harmful algae bloom (HABs). It is a widely distributed species and is responsible for producing paralytic shellfish poisoning toxins. However, the information associated with the environmental adaptation pathway and toxin biosynthesis in this species is still lacking. Therefore, this study focuses on the functional characterization of A. minutum unigenes obtained from transcriptome sequencing using the Illumina Hiseq 4000 sequencing platform. A total of 58,802 (47.05%) unigenes were successfully annotated using public databases such as NCBI-Nr, UniprotKB, EggNOG, KEGG, InterPRO and Gene Ontology (GO). This study has successfully identified key features that enable A. minutum to adapt to the marine environment, including several carbon metabolic pathways, assimilation of various sources of nitrogen and phosphorus. A. minutum was found to encode homologues for several proteins involved in saxitoxin biosynthesis, including the first three proteins in the pathway of saxitoxin biosynthesis, namely sxtA, sxtG and sxtB. The comprehensive transcriptome analysis presented in this study represents a valuable resource for understanding the dinoflagellates molecular metabolic model regarding nutrient acquisition and biosynthesis of saxitoxin.

A Global Approach to Estimating the Abundance and Duplication of Polyketide Synthase Domains in Dinoflagellates

Many dinoflagellate species make toxins in a myriad of different molecular configurations but the underlying chemistry in all cases is presumably via modular synthases, primarily polyketide synthases. In many organisms modular synthases occur as discrete synthetic genes or domains within a gene that act in coordination thus forming a module that produces a particular fragment of a natural product. The modules usually occur in tandem as gene clusters with a syntenic arrangement that is often predictive of the resultant structure. Dinoflagellate genomes however are notoriously complex with individual genes present in many tandem repeats and very few synthetic modules occurring as gene clusters, unlike what has been seen in bacteria and fungi. However, modular synthesis in all organisms requires a free thiol group that acts as a carrier for sequential synthesis called a thiolation domain. We scanned 47 dinoflagellate transcriptomes for 23 modular synthase domain models and compared their abundance among 10 orders of dinoflagellates as well as their co-occurrence with thiolation domains. The total count of domain types was quite large with over thirty-thousand identified, 29 000 of which were in the core dinoflagellates. Although there were no specific trends in domain abundance associated with types of toxins, there were readily observable lineage specific differences. The Gymnodiniales, makers of long polyketide toxins such as brevetoxin and karlotoxin had a high relative abundance of thiolation domains as well as multiple thiolation domains within a single transcript. Orders such as the Gonyaulacales, makers of small polyketides such as spirolides, had fewer thiolation domains but a relative increase in the number of acyl transferases. Unique to the core dinoflagellates, however, were thiolation domains occurring alongside tetratricopeptide repeats that facilitate protein-protein interactions, especially hexa and hepta-repeats, that may explain the scaffolding required for synthetic complexes capable of making large toxins. Clustering analysis for each type of domain was also used to discern possible origins of duplication for the multitude of single domain transcripts. Single domain transcripts frequently clustered with synonymous domains from multi-domain transcripts such as the BurA and ZmaK like genes as well as the multi-ketosynthase genes, sometimes with a large degree of apparent gene duplication, while fatty acid synthesis genes formed distinct clusters. Surprisingly the acyl-transferases and ketoreductases involved in fatty acid synthesis (FabD and FabG, respectively) were found in very large clusters indicating an unprecedented degree of gene duplication for these genes. These results demonstrate a complex evolutionary history of core dinoflagellate modular synthases with domain specific duplications throughout the lineage as well as clues to how large protein complexes can be assembled to synthesize the largest natural products known.

Comparison of 15 dinoflagellate genomes reveals extensive sequence and structural divergence in family Symbiodiniaceae and genus Symbiodinium

Background

Dinoflagellates in the family Symbiodiniaceae are important photosynthetic symbionts in cnidarians (such as corals) and other coral reef organisms. Breakdown of the coral-dinoflagellate symbiosis due to environmental stress (i.e. coral bleaching) can lead to coral death and the potential collapse of reef ecosystems. However, evolution of Symbiodiniaceae genomes, and its implications for the coral, is little understood. Genome sequences of Symbiodiniaceae remain scarce due in part to their large genome sizes (1–5 Gbp) and idiosyncratic genome features.

Results

Here, we present de novo genome assemblies of seven members of the genus Symbiodinium, of which two are free-living, one is an opportunistic symbiont, and the remainder are mutualistic symbionts. Integrating other available data, we compare 15 dinoflagellate genomes revealing high sequence and structural divergence. Divergence among some Symbiodinium isolates is comparable to that among distinct genera of Symbiodiniaceae. We also recovered hundreds of gene families specific to each lineage, many of which encode unknown functions. An in-depth comparison between the genomes of the symbiotic Symbiodinium tridacnidorum (isolated from a coral) and the free-living Symbiodinium natans reveals a greater prevalence of transposable elements, genetic duplication, structural rearrangements, and pseudogenisation in the symbiotic species.

Conclusions

Our results underscore the potential impact of lifestyle on lineage-specific gene-function innovation, genome divergence, and the diversification of Symbiodinium and Symbiodiniaceae. The divergent features we report, and their putative causes, may also apply to other microbial eukaryotes that have undergone symbiotic phases in their evolutionary history.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-00994-6.

Rapid protein evolution, organellar reductions, and invasive intronic elements in the marine aerobic parasite dinoflagellate Amoebophrya spp

BACKGROUND

Dinoflagellates are aquatic protists particularly widespread in the oceans worldwide. Some are responsible for toxic blooms while others live in symbiotic relationships, either as mutualistic symbionts in corals or as parasites infecting other protists and animals. Dinoflagellates harbor atypically large genomes (~ 3 to 250 Gb), with gene organization and gene expression patterns very different from closely related apicomplexan parasites. Here we sequenced and analyzed the genomes of two early-diverging and co-occurring parasitic dinoflagellate Amoebophrya strains, to shed light on the emergence of such atypical genomic features, dinoflagellate evolution, and host specialization.

RESULTS

We sequenced, assembled, and annotated high-quality genomes for two Amoebophrya strains (A25 and A120), using a combination of Illumina paired-end short-read and Oxford Nanopore Technology (ONT) MinION long-read sequencing approaches. We found a small number of transposable elements, along with short introns and intergenic regions, and a limited number of gene families, together contribute to the compactness of the Amoebophrya genomes, a feature potentially linked with parasitism. While the majority of Amoebophrya proteins (63.7% of A25 and 59.3% of A120) had no functional assignment, we found many orthologs shared with Dinophyceae. Our analyses revealed a strong tendency for genes encoded by unidirectional clusters and high levels of synteny conservation between the two genomes despite low interspecific protein sequence similarity, suggesting rapid protein evolution. Most strikingly, we identified a large portion of non-canonical introns, including repeated introns, displaying a broad variability of associated splicing motifs never observed among eukaryotes. Those introner elements appear to have the capacity to spread over their respective genomes in a manner similar to transposable elements. Finally, we confirmed the reduction of organelles observed in Amoebophrya spp., i.e., loss of the plastid, potential loss of a mitochondrial genome and functions.

CONCLUSION

These results expand the range of atypical genome features found in basal dinoflagellates and raise questions regarding speciation and the evolutionary mechanisms at play while parastitism was selected for in this particular unicellular lineage.

Integrated omics unveil the secondary metabolic landscape of a basal dinoflagellate

BACKGROUND

Some dinoflagellates cause harmful algal blooms, releasing toxic secondary metabolites, to the detriment of marine ecosystems and human health. Our understanding of dinoflagellate toxin biosynthesis has been hampered by their unusually large genomes. To overcome this challenge, for the first time, we sequenced the genome, microRNAs, and mRNA isoforms of a basal dinoflagellate, Amphidinium gibbosum, and employed an integrated omics approach to understand its secondary metabolite biosynthesis.

RESULTS

We assembled the ~ 6.4-Gb A. gibbosum genome, and by probing decoded dinoflagellate genomes and transcriptomes, we identified the non-ribosomal peptide synthetase adenylation domain as essential for generation of specialized metabolites. Upon starving the cells of phosphate and nitrogen, we observed pronounced shifts in metabolite biosynthesis, suggestive of post-transcriptional regulation by microRNAs. Using Iso-Seq and RNA-seq data, we found that alternative splicing and polycistronic expression generate different transcripts for secondary metabolism.

CONCLUSIONS

Our genomic findings suggest intricate integration of various metabolic enzymes that function iteratively to synthesize metabolites, providing mechanistic insights into how dinoflagellates synthesize secondary metabolites, depending upon nutrient availability. This study provides insights into toxin production associated with dinoflagellate blooms. The genome of this basal dinoflagellate provides important clues about dinoflagellate evolution and overcomes the large genome size, which has been a challenge previously.

Sex in Symbiodiniaceae dinoflagellates: genomic evidence for independent loss of the canonical synaptonemal complex

Dinoflagellates of the Symbiodiniaceae family encompass diverse symbionts that are critical to corals and other species living in coral reefs. It is well known that sexual reproduction enhances adaptive evolution in changing environments. Although genes related to meiotic functions were reported in Symbiodiniaceae, cytological evidence of meiosis and fertilisation are however yet to be observed in these taxa. Using transcriptome and genome data from 21 Symbiodiniaceae isolates, we studied genes that encode proteins associated with distinct stages of meiosis and syngamy. We report the absence of genes that encode main components of the synaptonemal complex (SC), a protein structure that mediates homologous chromosomal pairing and class I crossovers. This result suggests an independent loss of canonical SCs in the alveolates, that also includes the SC-lacking ciliates. We hypothesise that this loss was due in part to permanently condensed chromosomes and repeat-rich sequences in Symbiodiniaceae (and other dinoflagellates) which favoured the SC-independent class II crossover pathway. Our results reveal novel insights into evolution of the meiotic molecular machinery in the ecologically important Symbiodiniaceae and in other eukaryotes.