Unlocking the Complex Cell Biology of Coral-Dinoflagellate Symbiosis: A Model Systems Approach.

Symbiotic interactions occur in all domains of life, providing organisms with resources to adapt to new habitats. A prime example is the endosymbiosis between corals and photosynthetic dinoflagellates. Eukaryotic dinoflagellate symbionts reside inside coral cells and transfer essential nutrients to their hosts, driving the productivity of the most biodiverse marine ecosystem. Recent advances in molecular and genomic characterization have revealed symbiosis-specific genes and mechanisms shared among symbiotic cnidarians. In this review, we focus on the cellular and molecular processes that underpin the interaction between symbiont and host. We discuss symbiont acquisition via phagocytosis, modulation of host innate immunity, symbiont integration into host cell metabolism, and nutrient exchange as a fundamental aspect of stable symbiotic associations. We emphasize the importance of using model systems to dissect the cellular complexity of endosymbiosis, which ultimately serves as the basis for understanding its ecology and capacity to adapt in the face of climate change.

New plastids, old proteins: repeated endosymbiotic acquisitions in kareniacean dinoflagellates.

Dinoflagellates are a diverse group of ecologically significant micro-eukaryotes that can serve as a model system for plastid symbiogenesis due to their susceptibility to plastid loss and replacement via serial endosymbiosis. Kareniaceae harbor fucoxanthin-pigmented plastids instead of the ancestral peridinin-pigmented ones and support them with a diverse range of nucleus-encoded plastid-targeted proteins originating from the haptophyte endosymbiont, dinoflagellate host, and/or lateral gene transfers (LGT). Here, we present predicted plastid proteomes from seven distantly related kareniaceans in three genera (Karenia, Karlodinium, and Takayama) and analyze their evolutionary patterns using automated tree building and sorting. We project a relatively limited ( ~ 10%) haptophyte signal pointing towards a shared origin in the family Chrysochromulinaceae. Our data establish significant variations in the functional distributions of these signals, emphasizing the importance of micro-evolutionary processes in shaping the chimeric proteomes. Analysis of plastid genome sequences recontextualizes these results by a striking finding the extant kareniacean plastids are in fact not all of the same origin, as two of the studied species (Karlodinium armiger, Takayama helix) possess plastids from different haptophyte orders than the rest.

Conspicuous chloroplast with light harvesting-photosystem I/II megacomplex in marine Prorocentrum cordatum.

Marine photosynthetic (micro)organisms drive multiple biogeochemical cycles and display a large diversity. Among them, the bloom-forming, free-living dinoflagellate Prorocentrum cordatum CCMP 1329 (formerly P. minimum) stands out with its distinct cell biological features. Here, we obtained insights into the structural properties of the chloroplast and the photosynthetic machinery of P. cordatum using microscopic and proteogenomic approaches. High-resolution FIB/SEM analysis revealed a single large chloroplast (∼40% of total cell volume) with a continuous barrel-like structure, completely lining the inner face of the cell envelope and enclosing a single reticular mitochondrium, the Golgi apparatus, as well as diverse storage inclusions. Enriched thylakoid membrane fractions of P. cordatum were comparatively analyzed with those of the well-studied model-species Arabidopsis (Arabidopsis thaliana) using 2D BN DIGE. Strikingly, P. cordatum possessed a large photosystem-light harvesting megacomplex (>1.5 MDa), which is dominated by photosystems I and II (PSI, PSII), chloroplast complex I, and chlorophyll a-b binding light harvesting complex proteins. This finding parallels the absence of grana in its chloroplast and distinguishes from the predominant separation of PSI and PSII complexes in A. thaliana, indicating a different mode of flux balancing. Except for the core elements of the ATP synthase and the cytb6f-complex, the composition of the other complexes (PSI, PSII, and pigment-binding proteins, PBPs) of P. cordatum differed markedly from those of A. thaliana. Furthermore, a high number of PBPs was detected, accounting for a large share of the total proteomic data (∼65%) and potentially providing P. cordatum with flexible adaptation to changing light regimes.

A widespread coral-infecting apicomplexan with chlorophyll biosynthesis genes.

Apicomplexa is a group of obligate intracellular parasites that includes the causative agents of human diseases such as malaria and toxoplasmosis. Apicomplexans evolved from free-living phototrophic ancestors, but how this transition to parasitism occurred remains unknown. One potential clue lies in coral reefs, of which environmental DNA surveys have uncovered several lineages of uncharacterized basally branching apicomplexans1,2. Reef-building corals have a well-studied symbiotic relationship with photosynthetic Symbiodiniaceae dinoflagellates (for example, Symbiodinium3), but the identification of other key microbial symbionts of corals has proven to be challenging4,5. Here we use community surveys, genomics and microscopy analyses to identify an apicomplexan lineage-which we informally name 'corallicolids'-that was found at a high prevalence (over 80% of samples, 70% of genera) across all major groups of corals. Corallicolids were the second most abundant coral-associated microeukaryotes after the Symbiodiniaceae, and are therefore core members of the coral microbiome. In situ fluorescence and electron microscopy confirmed that corallicolids live intracellularly within the tissues of the coral gastric cavity, and that they possess apicomplexan ultrastructural features. We sequenced the genome of the corallicolid plastid, which lacked all genes for photosystem proteins; this indicates that corallicolids probably contain a non-photosynthetic plastid (an apicoplast6). However, the corallicolid plastid differs from all other known apicoplasts because it retains the four ancestral genes that are involved in chlorophyll biosynthesis. Corallicolids thus share characteristics with both their parasitic and their free-living relatives, which suggests that they are evolutionary intermediates and implies the existence of a unique biochemistry during the transition from phototrophy to parasitism.

Orchestrated translation specializes dinoflagellate metabolism three times per day.

Many cells specialize for different metabolic tasks at different times over their normal ZT cycle by changes in gene expression. However, in most cases, circadian gene expression has been assessed at the mRNA accumulation level, which may not faithfully reflect protein synthesis rates. Here, we use ribosome profiling in the dinoflagellate Lingulodinium polyedra to identify thousands of transcripts showing coordinated translation. All of the components in carbon fixation are concurrently regulated at ZT0, predicting the known rhythm of carbon fixation, and many enzymes involved in DNA replication are concurrently regulated at ZT12, also predicting the known rhythm in this process. Most of the enzymes in glycolysis and the TCA cycle are also regulated together, suggesting rhythms in these processes as well. Surprisingly, a third cluster of transcripts show peak translation at approximately ZT16, and these transcripts encode enzymes involved in transcription, translation, and amino acid biosynthesis. The latter has physiological consequences, as measured free amino acid levels increase at night and thus represent a previously undocumented rhythm in this model. Our results suggest that ribosome profiling may be a more accurate predictor of changed metabolic state than transcriptomics.

Ribosome profiling in the Symbiodiniacean dinoflagellate Fugacium kawagutii shows coordinated protein synthesis of enzymes in different pathways at different times of day.

Dinoflagellates respond to daily changes in light and dark by changes in cellular metabolism, yet the mechanisms used are still unclear. For example, Fugacium (previously Symbiodinium) kawagutii shows little difference in the transcriptome between day and night suggesting little transcriptional control over gene expression. Here, we have performed ribosome profiling at 2 h intervals over a daily light-dark cycle to assess the degree to which protein synthesis rates might change over the daily cycle. The number of F. kawagutii coding sequences with significant differences in the number of ribosome-protected fragments (RPF) over the 24-h cycle was 2923 using JTK_Cycle and 3655 using ECHO. The majority of the regulated transcripts showed peak translation at the onset of the dark period. The regulated sequences were assigned to different KEGG pathways and transcripts that were translated at roughly the same time were termed concurrently regulated. Both analyses revealed concurrent regulation of many transcripts whose gene products were involved in spliceosome or lysosome biogenesis with peak translation rates around the onset of the dark period, while others, involved in nitrate metabolism and ribosomal proteins, were preferentially translated around the onset of the day phase or the end of the night phase, respectively. In addition, some sequences involved in DNA synthesis were preferentially translated at the end of the day. We conclude that light-dark cycles seem able to synchronize translation of some transcripts encoding proteins involved in a range of different cellular processes, and propose that these changes may help the cells adapt and alter their metabolism as a function of the time of day.

Symbiosis modulates gene expression of symbionts, but not coral hosts, under thermal challenge.

Increasing ocean temperatures are causing dysbiosis between coral hosts and their symbionts. Previous work suggests that coral host gene expression responds more strongly to environmental stress compared to their intracellular symbionts; however, the causes and consequences of this phenomenon remain untested. We hypothesized that symbionts are less responsive because hosts modulate symbiont environments to buffer stress. To test this hypothesis, we leveraged the facultative symbiosis between the scleractinian coral Oculina arbuscula and its symbiont Breviolum psygmophilum to characterize gene expression responses of both symbiotic partners in and ex hospite under thermal challenges. To characterize host and in hospite symbiont responses, symbiotic and aposymbiotic O. arbuscula were exposed to three treatments: (1) control (18°C), (2) heat (32°C), and (3) cold (6°C). This experiment was replicated with B. psygmophilum cultured from O. arbuscula to characterize ex hospite symbiont responses. Both thermal challenges elicited classic environmental stress responses (ESRs) in O. arbuscula regardless of symbiotic state, with hosts responding more strongly to cold challenge. Hosts also exhibited stronger responses than in hospite symbionts. In and ex hospite B. psygmophilum both down-regulated gene ontology pathways associated with photosynthesis under thermal challenge; however, ex hospite symbionts exhibited greater gene expression plasticity and differential expression of genes associated with ESRs. Taken together, these findings suggest that O. arbuscula hosts may buffer environments of B. psygmophilum symbionts; however, we outline the future work needed to confirm this hypothesis.

Symbiodiniaceae diversity varies by host and environment across thermally distinct reefs.

Endosymbiotic dinoflagellates (Symbiodiniaceae) influence coral thermal tolerance at both local and regional scales. In isolation, the effects of host genetics, environment, and thermal disturbances on symbiont communities are well understood, yet their combined effects remain poorly resolved. Here, we investigate Symbiodiniaceae across 1300 km in Australia's Coral Sea Marine Park to disentangle these interactive effects. We identified Symbiodiniaceae to species-level resolution for three coral species (Acropora cf humilis, Pocillopora verrucosa, and Pocillopora meandrina) by sequencing two genetic markers of the symbiont (ITS2 and psbAncr), paired with genotype-by-sequencing of the coral host (DArT-seq). Our samples predominantly returned sequences from the genus Cladocopium, where Acropora cf humilis affiliated with C3k, Pocillopora verrucosa with C. pacificum, and Pocillopora meandrina with C. latusorum. Multivariate analyses revealed that Acropora symbionts were driven strongly by local environment and thermal disturbances. In contrast, Pocillopora symbiont communities were both partitioned 2.5-fold more by host genetic structure than by environmental structure. Among the two Pocillopora species, the effects of environment and host genetics explained four times more variation in symbionts for P. meandrina than P. verrucosa. The concurrent bleaching event in 2020 had variable impacts on symbiont communities, consistent with patterns in P. verrucosa and A. cf humilis, but not P. meandrina. Our findings demonstrate how symbiont macroscale community structure responses to environmental gradients depend on host species and their respective population structure. Integrating host, symbiont, and environmental data will help forecast the adaptive potential of corals and their symbionts amidst a rapidly changing environment.

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.

Complex parental effects impact variation in larval thermal tolerance in a vertically transmitting coral.

Coral populations must be able to adapt to changing environmental conditions for coral reefs to persist under climate change. The adaptive potential of these organisms is difficult to forecast due to complex interactions between the host animal, dinoflagellate symbionts and the environment. Here we created 26 larval families from six Montipora capitata colonies from a single reef, showing significant, heritable variation in thermal tolerance. Our results indicate that 9.1% of larvae are expected to exhibit four times the thermal tolerance of the general population. Differences in larval thermotolerance were driven mainly by maternal contributions, but we found no evidence that these effects were driven by symbiont identity despite vertical transmission from the dam. We also document no evidence of reproductive incompatibility attributable to symbiont identity. These data demonstrate significant genetic variation within this population which provides the raw material upon which natural selection can act.

The second most abundant dinophyte in the ponds of a botanical garden is a species new to science.

In the microscopy realm, a large body of dark biodiversity still awaits to be uncovered. Unarmoured dinophytes are particularly neglected here, as they only present inconspicuous traits. In a remote German locality, we collected cells, from which a monoclonal strain was established, to study morphology using light and electron microscopy and to gain DNA sequences from the rRNA operon. In parallel, we detected unicellular eukaryotes in ponds of the Botanical Garden Munich-Nymphenburg by DNA-metabarcoding (V4 region of the 18S rRNA gene), weekly sampled over the course of a year. Strain GeoK*077 turned out to be a new species of Borghiella with a distinct position in molecular phylogenetics and characteristic coccoid cells of ovoid shape as the most important diagnostic trait. Borghiella ovum, sp. nov., was also present in artificial ponds of the Botanical Garden and was the second most abundant dinophyte detected in the samples. More specifically, Borghiella ovum, sp. nov., shows a clear seasonality, with high frequency during winter months and complete absence during summer months. The study underlines the necessity to assess the biodiversity, particularly of the microscopy realm more ambitiously, if even common species such as formerly Borghiella ovum are yet unknown to science.

DsRNA sequencing revealed a previously missed terminal sequence of a +ssRNA virus that infects dinoflagellate Heterocapsa circularisquama.

Heterocapsa circularisquama RNA virus (HcRNAV) is the only dinoflagellate-infecting RNA virus cultured. However, only two strains of HcRNAV have been registered with complete genome sequences (strains 34 and 109 for UA and CY types, respectively). To extend the genomic information of HcRNAV, we performed full-genome sequencing of an unsequenced strain of HcRNAV (strain A8) using the fragmented and primer-ligated double-stranded RNA (dsRNA) sequencing (FLDS) method. The complete genome of HcRNAV A8 with 4457 nucleotides (nt) was successfully determined, and sequence alignment of the major capsid protein gene suggested that A8 was a UA-type strain, consistent with its intraspecific host specificity. The complete sequence was found to be 80 nt longer at the 5' terminus than the registered sequences of HcRNAV strains (34 and 109), suggesting that FLDS is more reliable for determining the terminal sequence than conventional methods (5' Rapid Amplification of cDNA End). Our study contributes to a better understanding of dinoflagellate-infecting viruses with limited sequence data.

Fatty acid production and associated gene pathways are altered by increased salinity and dimethyl sulfoxide treatments during cryopreservation of Symbiodinium pilosum (Symbiodiniaceae).

The Symbiodinium genus is ancestral among other Symbiodiniaceae lineages with species that are both symbiotic and free living. Changes in marine ecosystems threaten their existence and crucial ecological roles. Cryopreservation offers an avenue for their long-term storage for future habitat restoration after coral bleaching. In our previous study we demonstrated that high salinity treatments of Symbiodiniaceae isolates led to changes in their fatty acid (FA) profiles and higher cell viabilities after cryopreservation. In this study, we investigated the role of increased salinity on FA production and the genes involved in FA biosynthesis and degradation pathways during the cryopreservation of Symbiodinium pilosum. Overall, there was a twofold increase in mass of FAs produced by S. pilosum after being cultured in medium with increased salinity (54 parts per thousand; ppt). Dimethyl sulfoxide (Me2SO) led to a ninefold increase of FAs in standard salinity (SS) treatment, compared to a fivefold increase in increased salinity (IS) treatments. The mass of the FA classes returned to baseline during recovery. Transcriptomic analyses showed an acyl carrier protein gene was significantly upregulated after Me2SO treatment in the SS cultures. Cytochrome P450 reductase genes were significantly down regulated after Me2SO addition in SS treatment preventing FA degradation. These changes in the expression of FA biosynthesis and degradation genes contributed to more FAs in SS treated isolates. Understanding how increased salinity changes FA production and the roles of specific genes in regulating FA pathways will help improve current freezing protocols for Symbiodiniaceae and other marine microalgae.

Cross resistance to brevetoxin-3 by kdr and super-kdr mutations in house flies.

The dinoflagellate Karenia brevis is a causative agent of red tides in the Gulf of Mexico and generates a potent family of structurally related brevetoxins that act via the voltage-sensitive Na+ channel. This project was undertaken to better understand the neurotoxicology and kdr cross-resistance to brevetoxins in house flies by comparing the susceptible aabys strain to ALkdr (kdr) and JPskdr (super-kdr). When injected directly into the hemocoel, larvae exhibited rigid, non-convulsive paralysis consistent with prolongation of sodium channel currents, the known mechanism of action of brevetoxins. In neurophysiological studies, the firing frequency of susceptible larval house fly central nervous system preparations showed a > 200% increase 10 min after treatment with 1 nM brevetoxin-3. This neuroexcitation is consistent with the spastic paralytic response seen after hemocoel injections. Target site mutations in the voltage-sensitive sodium channel of house flies, known to confer knockdown resistance (kdr and super-kdr) against pyrethroids, attenuated the effect of brevetoxin-3 in baseline firing frequency and toxicity assays. The rank order of sensitivity to brevetoxin-3 in both assays was aabys > ALkdr > JPskdr. At the LD50 level, resistance ratios for the knockdown resistance strains were 6.9 for the double mutant (super-kdr) and 2.3 for the single mutant (kdr). The data suggest that knockdown resistance mutations may be one mechanism by which flies survive brevetoxin-3 exposure during red tide events.

DNA and RNA Stability of Marine Microalgae in Cold-Stored Sediments and Its Implications in Metabarcoding Analyses.

The ever-increasing applications of metabarcoding analyses for environmental samples demand a well-designed assessment of the stability of DNA and RNA contained in cells that are deposited or buried in marine sediments. We thus conducted a qPCR quantification of the DNA and RNA in the vegetative cells of three microalgae entrapped in facsimile marine sediments and found that >90% of DNA and up to 99% of RNA for all microalgal species were degraded within 60 days at 4 °C. A further examination of the potential interference of the relic DNA of the vegetative cells with resting cyst detection in sediments was performed via a metabarcoding analysis in artificial marine sediments spiked with the vegetative cells of two Kareniaceae dinoflagellates and the resting cysts of another three dinoflagellates. The results demonstrated a dramatic decrease in the relative abundances of the two Kareniaceae dinoflagellates in 120 days, while those of the three resting cysts increased dramatically. Together, our results suggest that a positive detection of microalgae via metabarcoding analysis in DNA or RNA extracted from marine sediments strongly indicates the presence of intact or viable cysts or spores due to the rapid decay of relic DNA/RNA. This study provides a solid basis for the data interpretation of metabarcoding surveys, particularly in resting cyst detection.

Tightly Constrained Genome Reduction and Relaxation of Purifying Selection during Secondary Plastid Endosymbiosis.

Endosymbiosis, the establishment of a former free-living prokaryotic or eukaryotic cell as an organelle inside a host cell, can dramatically alter the genomic architecture of the endosymbiont. Plastids or chloroplasts, the light-harvesting organelle of photosynthetic eukaryotes, are excellent models to study this phenomenon because plastid origin has occurred multiple times in evolution. Here, we investigate the genomic signature of molecular processes acting through secondary plastid endosymbiosis-the origination of a new plastid from a free-living eukaryotic alga. We used phylogenetic comparative methods to study gene loss and changes in selective regimes on plastid genomes, focusing on green algae that have given rise to three independent lineages with secondary plastids (euglenophytes, chlorarachniophytes, and Lepidodinium). Our results show an overall increase in gene loss associated with secondary endosymbiosis, but this loss is tightly constrained by the retention of genes essential for plastid function. The data show that secondary plastids have experienced temporary relaxation of purifying selection during secondary endosymbiosis. However, this process is tightly constrained, with selection relaxed only relative to the background in primary plastids. Purifying selection remains strong in absolute terms even during the endosymbiosis events. Selection intensity rebounds to pre-endosymbiosis levels following endosymbiosis events, demonstrating the changes in selection efficiency during different origin phases of secondary plastids. Independent endosymbiosis events in the euglenophytes, chlorarachniophytes, and Lepidodinium differ in their degree of relaxation of selection, highlighting the different evolutionary contexts of these events. This study reveals the selection-drift interplay during secondary endosymbiosis and evolutionary parallels during organellogenesis.

Mesoscale oceanographic meanders influence protist community function and structure in the southern Indian Ocean.

The interface between the nutrient-rich Southern Ocean and oligotrophic Indian Ocean creates unique environmental conditions that can strongly influence biological processes. We investigated protist communities across a mesoscale meander of the Subtropical Front within the Southern Indian Ocean. 18S V9 rDNA metabarcoding suggests a diverse protist community in which the dinoflagellates and parasitic Syndiniales were abundant. Diversity was highest in frontal waters of the mesoscale meander, with differences in community structure inside and outside the meander. While the overall community was dominated by mixotrophic taxa, the frontal boundary of the meander had increased abundances of heterotrophic taxa, with potential implications for net atmospheric CO2 drawdown. Pulse amplitude modulated (PAM) fluorimetry revealed significant differences in the photophysiology of phytoplankton communities inside and outside the meander. By using single-cell PAM microscopy, we identified physiological differences between dinoflagellate and coccolithophore taxa, which may have contributed to changes in photophysiology observed at community level. Overall, our results demonstrate that frontal areas have a strong impact on the composition of protist communities in the Southern Ocean with important implications for understanding biological processes in this region.

Assessing the contribution of bacteria to the heat tolerance of experimentally evolved coral photosymbionts.

Coral reefs are extremely vulnerable to ocean warming, which triggers coral bleaching-the loss of endosymbiotic microalgae (Symbiodiniaceae) from coral tissues, often leading to death. To enhance coral climate resilience, the symbiont, Cladocopium proliferum was experimentally evolved for >10 years under elevated temperatures resulting in increased heat tolerance. Bacterial 16S rRNA gene metabarcoding showed the composition of intra- and extracellular bacterial communities of heat-evolved strains was significantly different from that of wild-type strains, suggesting bacteria responded to elevated temperatures, and may even play a role in C. proliferum thermal tolerance. To assess whether microbiome transplantation could enhance heat tolerance of the sensitive wild-type C. proliferum, we transplanted bacterial communities from heat-evolved to the wild-type strain and subjected it to acute heat stress. Microbiome transplantation resulted in the incorporation of only 30 low-abundance strains into the microbiome of wild-type cultures, while the relative abundance of 14 pre-existing strains doubled in inoculated versus uninoculated samples. Inoculation with either wild-type or heat-evolved bacterial communities boosted C. proliferum growth, although no difference in heat tolerance was observed between the two inoculation treatments. This study provides evidence that Symbiodiniaceae-associated bacterial communities respond to heat selection and may contribute to coral adaptation to climate change.

Interannual dynamics of putative parasites (Syndiniales Group II) in a coastal ecosystem.

Temporal dynamics of Syndiniales Group II were investigated combining 18S rDNA amplicon sequencing and direct microscopy counts (fluorescence in situ hybridization-tyramide signal amplification [FISH-TSA]) during 5 years. The study was undertaken in meso-eutrophic coastal ecosystem, dominated by diatoms, the haptophyte Phaeocystis globosa and exhibiting relatively low dinoflagellate abundance (max. 18.6 × 103 cells L-1 ). Consistent temporal patterns of Syndiniales Group II were observed over consecutive years highlighting the existence of local populations. According to sequencing data, Syndiniales Group II showed increasing abundance and richness in summer and autumn. Dinospores counted by microscopy, were present at low abundances and were punctuated by transient peaks. In summer dinospore highest abundance (559 × 103 L-1 ) and prevalence (38.5%) coincided with the peak abundance of the dinoflagellate Prorocentrum minimum (13 × 103 L-1 ) while in autumn Syndiniales Group II likely had more diversified hosts. Although, several peaks of dinospore and read abundances coincided, there was no consistent relation between them. Ecological assembly processes at a seasonal scale revealed that stochastic processes were the main drivers (80%) of the Group II community assembly, though deterministic processes were noticeable (20%) in June and July. This latter observation may reflect the specific Syndiniales-dinoflagellate interactions in summer.

Diverse RNA Viruses Associated with Diatom, Eustigmatophyte, Dinoflagellate, and Rhodophyte Microalgae Cultures.

Unicellular microalgae are of immense ecological importance with growing commercial potential in industries such as renewable energy, food, and pharmacology. Viral infections can have a profound impact on the growth and evolution of their hosts. However, very little is known of the diversity within, and the effect of, unicellular microalgal RNA viruses. In addition, identifying RNA viruses in these organisms that could have originated more than a billion years ago constitutes a robust data set to dissect molecular events and address fundamental questions in virus evolution. We assessed the diversity of RNA viruses in eight microalgal cultures, including representatives from the diatom, eustigmatophyte, dinoflagellate, red algae, and euglenid groups. Using metatranscriptomic sequencing combined with bioinformatic approaches optimized to detect highly divergent RNA viruses, we identified 10 RNA virus sequences, with nine constituting new viral species. Most of the newly identified RNA viruses belonged to the double-stranded Totiviridae, Endornaviridae, and Partitiviridae, greatly expanding the reported host range for these families. Two new species belonging to the single-stranded RNA viral clade Marnaviridae, commonly associated with microalgal hosts, were also identified. This study highlights that a substantial diversity of RNA viruses likely exists undetected within the unicellular microalgae. It also highlights the necessity for RNA viral characterization and for investigation of the effects of viral infections on microalgal physiology, biology, and growth, considering their environmental and industrial roles. IMPORTANCE Our knowledge of the diversity of RNA viruses infecting microbial algae-the microalgae-is minimal. However, describing the RNA viruses infecting these organisms is of primary importance at both the ecological and economic scales because of the fundamental roles these organisms play in aquatic environments and their growing value across a range of industrial fields. Using metatranscriptomic sequencing, we aimed to reveal the RNA viruses present in cultures of eight microalgae species belonging to the diatom, dinoflagellate, eustigmatophyte, rhodophyte, and euglena major clades of algae. Accordingly, we identified 10 new divergent RNA virus species belonging to RNA virus families as diverse as the double-stranded Totiviridae, Endornaviridae, and Partitiviridae and the single-stranded Marnaviridae. By expanding the known diversity of RNA viruses infecting unicellular eukaryotes, this study contributes to a better understanding of the early evolution of the virosphere and will inform the use of microalgae in industrial applications.