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

Alex Veglia, Kalia Bistolas, Christian Voolstra, Benjamin Hume, Hans-Joachim Ruscheweyh, Serge Planes, Denis Allemand, Emilie Boissin, Patrick Wincker, Julie Poulain, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

Alex Veglia, Kalia Bistolas, Christian Voolstra, Benjamin Hume, Hans-Joachim Ruscheweyh, et al. "Endogenous viral elements reveal associations between a non-retroviral RNA virus and symbiotic dinoflagellate genomes." Communications Biology, 2023, 6 (1), pp.566. 10.1038/s42003-023-04917-9. hal-04176296

HAL Id: hal-04176296

https://hal.science/hal-04176296

Submitted on 12 Aug 2023

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>
<figure><img src="image_3.png" /></figure>

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

Alex J. Veglia ID ¹,²³, Kalia S. I. Bistolas ID ²,²³✉, Christian R. Voolstra ID ³, Benjamin C. C. Hume ID ³, Hans-Joachim Ruscheweyh ID ⁴, Serge Planes⁵, Denis Allemand ID ⁶, Emilie Boissin ID ⁵, Patrick Wincker ID ⁷,⁸, Julie Poulain ⁷,⁸, Clémentine Moulin⁹, Guillaume Bourdin ID ¹⁰, Guillaume Iwankow⁵, Sarah Romac¹¹, Sylvain Agostini ID ¹², Bernard Banaigs⁵, Emmanuel Boss ID ¹⁰, Chris Bowler ID ¹³, Colomban de Vargas¹¹, Eric Douville ID ¹⁴, Michel Flores ID ¹⁵, Didier Forcioli ID ¹⁶,¹⁷, Paola Furla ID ¹⁶,¹⁷, Pierre E. Galand ID ¹⁸, Eric Gilson ID ¹⁶,¹⁹, Fabien Lombard ID ²⁰, Stéphane Pesant²¹, Stéphanie Reynaud ID ⁶, Shinichi Sunagawa ID ⁴, Olivier P. Thomas²², Romain Troublé ID ⁹, Didier Zoccola ID ⁶, Adrienne M. S. Correa ID ¹ & Rebecca L. Vega Thurber ID ²

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.
Endogenous viral elements, or “EVEs,” arise when whole or fragmented viral genomes are incorporated into host cell germlines. Once integrated, EVEs may propagate across successive host generations, potentially becoming fixed in a population through natural selection or drift¹,². Therefore, the presence and content of EVEs can provide clues into the evolutionary relationships among host species and shed light on ancient and modern virus-host interactions³. To date, most EVEs described in metazoan and plant genomes are retroviral, as this viral group must integrate their genome (as a provirus) into the genome of the host to replicate. Retroviruses thus possess and encode all of the molecular machinery (e.g. reverse transcriptases, integrases) required to integrate autonomously⁴. Remarkably, however, sequences from viruses that do not encode reverse transcriptases or exploit integration as a component of an obligate replication strategy—even viruses with no DNA stage—have also recently been detected as EVEs in diverse eukaryotic genomes⁵⁻¹¹. These non-retroviral RNA EVEs have been reported in hosts ranging from unicellular algae to chiropteran (bat) genomes¹²⁻¹⁸. Though the mechanisms behind non-retroviral integration continue to be explored, viral sequences may be introduced via nonhomologous recombination and repair, through interactions with host-provisioned integrases and reverse transcriptases supplied on mobile elements (e.g. retroelements), or by utilizing co-infecting viruses⁶,⁷.

Endogenization of any viral sequence (including non-retroviral EVEs) may have positive, neutral or negative effects on a host¹⁹⁻²¹. While many EVEs are functionally defective or deleterious and ultimately removed from a population via purifying selection, retained EVEs may remodel the genomic architecture of their hosts or introduce sources of genetic innovation later co-opted for host function (i.e. exaptation²²,²³). Such ‘domesticated’ EVEs can be co-opted by hosts and utilized as regulatory elements, transcription factors, or functional proteins with purposes ranging from organism development to synaptic plasticity in the mammalian brain²⁴⁻²⁸. In particular, non-retroviral EVEs potentially serve as antiviral prototypes that help hosts combat infection by exogenous viruses currently circulating in the population¹⁴,²⁹⁻³¹. Mechanisms underpinning EVE-derived immunity can include cell receptor interference, nucleic acid sequence recognition (e.g., RNAi), or even replication sabotage through production of faulty virus proteins from EVEs³². If expressed, EVEs may have a significant influence on the health, physiology and/or behavior of their hosts in natural and experimental systems³¹,³³,³⁴.

Investigating the distribution, sequence identity, and function of EVEs can yield insight into virus-host interactions across generations. EVEs catalogue a subset of the viruses that a host lineage has encountered and can link homologous extant viruses to contemporary hosts or known disease states³¹. Because integrated elements may accrue mutations at a slower rate than exogenous viral genomes⁶,³⁵, EVEs can fill gaps in virus-host networks and act as synapomorphies, indicating the minimum time that a virus may have interacted with a host. As ‘genomic fossils’, EVEs have helped paleovirologists date the minimum origin of *Circoviridae*³⁶, *Hepadnaviridae*³⁷, *Bornaviridae*³⁸, *Flaviviridae*³⁹, *Lentiviridae*⁴⁰,⁴¹, and *Spumaviridae*⁴² infections within metazoans²,²⁴,³⁵,⁴³,⁴⁴ (reviewed by Barreat and Katzourakis in 2022⁴⁵).

Coral holobionts – the cnidarian animal and its resident microbial assemblage, including dinoflagellates in the family Symbiodiniaceae, bacteria, archaea, fungi, and viruses – are an ecologically and economically valuable, multipartite non-model system⁴⁶,⁴⁷. Symbiodiniaceae are key obligate nutritional symbionts of corals and support their hosts in the construction of reef frameworks⁴⁸. However, environmental stress can break down

coral-Symbiodiniaceae partnerships, resulting in bleaching – the mass loss of Symbiodiniaceae cells⁴⁹. Some bleaching signs (paling of a coral colony) are hypothesized to also result from viral lysis of Symbiodiniaceae²¹,⁵⁰⁻⁵⁴, but direct evidence supporting this hypothesis remains limited. Overall, the role of viruses in coral colony health and disease requires further examination.

Non-retroviral +ssRNA dinoRNAV sequences were first reported in stony corals based on five metatranscriptomic sequences and corroborated by Symbiodiniaceae EST libraries⁵⁵. Subsequent studies indicated that similar +ssRNA viruses are commonly detected in coral RNA viromes and metatranscriptomes, as well as via targeted amplicon assays⁵⁴,⁵⁶⁻⁵⁸. These viruses exhibit synteny and significant homology to *Heterocapsa circularisquama* RNA virus (HcRNAV⁵⁷), the sole recognized representative of the genus *Dinornavirus* and a known pathogen of free-living dinoflagellates⁵⁹. Both HcRNAV and dinoRNAV sequences detected in coral holobiont tissues contain two ORFs – a Major Capsid Protein (MCP) and RNA dependent RNA polymerase (RdRp). Furthermore, icosahedral virus-like particle (VLP) arrays resembling HcRNAV (but with 40% smaller individual particle diameters) have been imaged in the Symbiodiniaceae-dense coral gastrodermis tissue and in Symbiodiniaceae themselves⁶⁰. Levin et al. (2017)⁵⁷ assembled the 5.2 kb genome of a putative dinoRNAV from a poly(A)-selected metatranscriptome generated from cultured *Symbiodinium*. The assembly contained a 5' dinoflagellate spliced leader ("dinoSL")⁶¹ – a component of >95% of Symbiodiniaceae mRNAs, speculated to illustrate molecular mimicry – and exhibited >1000-fold higher expression in a thermosensitive *Cladocopium* C1 population relative to a thermotolerant population of this Symbiodiniaceae strain at ambient temperatures (27 °C⁴⁸,⁵⁷). Together, the findings from these studies suggest that Symbiodiniaceae are target hosts of reef-associated dinoRNAVs.

This study (1) systematically searched for putative endogenized dinoRNAVs in metagenomes from in situ (symbiotic) coral colonies and seawater, as well as in available genomes of Symbiodiniaceae and aposymbiotic (symbiont-free) cnidarians, (2) investigated the evolutionary relationship of putative dinoRNAV EVEs to exogenous reef-associated dinoRNAV sequences, and (3) made preliminary inferences regarding the distribution and possible function of these dinoRNAV EVEs based on their detection, prevalence, and genomic context.

## Results and discussion

**Evidence of endogenized dinoRNAVs in coral holobiont metagenomes.** Putative dinoRNAV EVEs were detected in metagenomes generated from 42 cnidarian holobionts out of 269 sampled across the South Pacific Ocean (Supplementary Data 1). The majority of endogenized dinoRNAVs were identified in hydrocoral metagenomes (*Millepora* spp.; 70.5%, *n* = 105) which predominantly harbored *Symbiodinium* dinoflagellates but EVE-like sequences were also observed in scleractinian coral metagenomes (*Pocillopora* spp.; 29.5%, *n* = 15.) which predominantly harbored *Cladocopium* and *Durusdinium* dinoflagellates (Fig. 1a, c). No dinoRNAV-like sequences were detected among *Porites* spp. metagenomes (Figs. 1, 2). Hydrocoral metagenomes were sequenced at equivalent depths as scleractinian corals and had comparable levels of annotation (Supplementary Fig. 1, Supplementary Data 2); thus, higher dinoRNAV EVE prevalence in hydrocoral libraries was likely not a result of methodological bias. Of the 11 evaluated South Pacific islands, dinoRNAV EVEs were identified in samples from eight (Guam, Gambier, Moorea, Cook, Niue, Malpelo, Coïba, and Las Perlas), spanning 18 unique sites (Fig. 1b, d). Among *Pocillopora* spp. metagenomes, putative dinoRNAV EVEs were only
<figure><img src="image_4.png" /><figcaption>**Fig. 1 Islands and species (cnidarian and dinoflagellate) correlating with dinoRNAV EVE-like sequence detection among Tara Pacific metagenomes.**</figcaption></figure>

a Count of scaffolds with putative endogenized dinoRNAV-like sequences among Tara Pacific metagenomes, grouped by island and spaced longitudinally by location sampled. b Sampling sites of Tara Pacific metagenomes explored for endogenized dinoRNAV-like sequences in this study. Internal circles indicate dominant Symbiodiniaceae genera based on ITS2 type profiles, outer ring denotes coral host(s) sampled at each island. c Symbiodiniaceae ITS2 type profile metabarcoding as delineated via Symportal<sup>119</sup> within island and host. d Sample design of Tara Pacific libraries queried for dinoRNAV EVEs. [x] and black circles on map indicate island locations or species where no dinoRNAV-like sequences were detected. Icons derived from the Noun Project.

<figure><img src="image_5.png" /><figcaption>**Fig. 2 Quantity of putative endogenized dinoRNAV per metagenome or genome.** Total quantity of putative endogenized dinoRNAV EVEs identified, broadly organized by sample source (metagenome or genome), and number of libraries or assemblies queried (after source name, Supplementary Data 6). Opaque circles denote the sum total of dinoRNAV EVE-like sequences identified from each source, while transparent circles denote individual counts of putative dinoRNAV EVEs per library. Icons created with BioRender.</figcaption></figure>
identified on the Central American coast (CAMR, Coastal Pacific Longhurst Province) and were absent in Melanesia, Micronesia, and Polynesia; at these latter sites, dinoRNAVs were largely found in *Millepora* hydrocoral metagenomes. Importantly, endogenized dinoRNAV open reading frames (ORFs) appeared to be immediately adjacent to ORFs identified as dinoflagellate (typically Symbiodiniaceae) genes—they were not proximal to coral genes or those of other cellular organisms abundant in these metagenomes (Supplementary Data 3).

We examined the Symbiodiniaceae ITS2 profiles associated with each metagenome and found that putative dinoRNAV EVEs were primarily associated with *Symbiodinium*, *Cladocopium*, and *Durusdinium*, which exhibited variation on both host and regional scales (Fig. 1c, Supplementary Data 4). DinoRNAV EVEs were more common in *Symbiodinium*-dominated cnidarians ($F_{2,1044} = 25.8, p < 0.0001$, nested ANOVA; Supplementary Fig. 2, Supplementary Data 5) relative to cnidarians hosting other Symbiodiniaceae genera, regardless of host. This suggested that dinoRNAV integration may be particularly recurrent or conserved within the genus *Symbiodinium* (Fig. 1).

To determine if these putative viral integrations were specific to cnidarian holobiont metagenomes and ensure that they were not artifacts of shared sample processing and sequencing procedures of the Tara Pacific pipeline, we also analyzed seawater metagenomes and publicly available metagenomes from the stony coral-dinoflagellate holobiont, *Acropora* spp. (Supplementary Data 1B,⁶²⁻⁶⁹). Examination of 120 Tara Oceans pelagic seawater metagenomes⁷⁰ yielded no sequences sharing homology to dinoRNAVs. The concentration of Symbiodiniaceae cells within cnidarian tissues is considerably higher than that of the surrounding seawater⁷¹⁻⁷⁴. On average, only 1.46 ± 0.08% of assembled contigs in seawater metagenomes were annotated as Symbiodiniaceae. Thus, lack of detection of dinoRNAV-like sequences from seawater metagenomes is likely due to reduced genomic signal of Symbiodiniaceae in the water column, rather than a lack of EVEs associated with Symbiodiniaceae lineages in seawater. However, it also must be noted that these Tara Oceans seawater metagenomes were not collected concurrently with coral samples⁷⁵. Analysis of the 30 non-Tara *Acropora* holobiont metagenomes identified 29 more putative dinoRNAV EVEs (Fig. 2; Supplementary Data 6). These dinoRNAV EVEs were again neighboring dinoflagellate ORFs. While the Caribbean *Acropora* metagenomes analyzed contained too few reads to resolve the dominant Symbiodiniaceae present, earlier studies of the same coral colonies identified *Symbiodinium* spp. as the primary symbiont present⁷⁶.

The identification of endogenized dinoRNAV-like sequences in cnidarian holobiont metagenomes, combined with the proximity of dinoRNAV-like ORFs to dinoflagellate-like sequences across metagenomes harboring diverse dinoflagellate consortia, collectively indicate that dinoRNAV EVEs are widespread among Symbiodiniaceae genera (Fig. 2 cyan dots).

**Endogenized DinoRNAVs detected in Symbiodiniaceae genomes.** To further test the hypothesis that dinoRNAVs on reefs infect dinoflagellate symbionts and not cnidarians, we examined 18 scaffold-scale genome assemblies representing the dinoflagellate families Symbiodiniaceae and Suessiaceae as well as 25 cnidarian genomes spanning 10 genera (Supplementary Data 1B⁶²⁻⁶⁹; Fig. 2; Table 1). Alignments revealed no evidence of endogenized dinoRNAVs in any of the 151,782 aposymbiotic (dinoflagellate-free) cnidarian scaffolds. In contrast, the same approach uncovered 351 (of 593,433) dinoflagellate scaffolds with evidence of endogenized dinoRNAVs (Fig. 2; Table 1). The identified 351 dinoRNAV EVE-containing scaffolds were

observed across 17 of the 18 dinoflagellate genome assemblies (Table 1). DinoRNAV EVEs were also observed in two assemblies from the free-living dinoflagellate genus, *Polarella* (family Suesiaceae), which is closely related to the family Symbiodiniaceae, and served as an outgroup in this study⁷⁷,⁷⁸. Interestingly, assemblies belonging to *Symbiodinium*, the most ancestral Symbiodiniaceae genus⁴⁸, contained a higher number of scaffolds with putative dinoRNAV EVEs ($\bar{x}=28.11, stdev=10.7$) relative to assemblies of other Symbiodiniaceae genera ($\bar{x}=8.71, stdev=11$; Fig. 2 cyan dots; Table 1). This result may clarify why observations of dinoRNAV-like ORFs were more common in metagenomes dominated by *Symbiodinium* (Fig. 1c). The dinoflagellate genome assembly with no detected dinoRNAV EVEs belonged to a relatively incomplete assembly of *Cladocopium* C15, which had the second lowest N50 and lowest BUSCO completeness score of all genomes examined (completeness 11.6%, relative to the average 24.54%; Table 1, Supplementary Data 7). The lower coverage/completeness of the *Cladocopium* C15 assembly indicates a reduced window into this genome. It is therefore possible that when a more complete assembly is generated, dinoRNAV EVE-like sequences will be detectable from this dinoflagellate. However, a linear model suggested that there was no relationship between dinoRNAV EVE detection and assembly statistics (i.e. query length, N50, or completeness; see Supplementary Data 8 for linear model output). Instead, dinoflagellate genus was the strongest predictor of dinoRNAV detection in a genome (LM results: Genus $F = 5.74, p = 0.012$) and dinoRNAV detections were significantly higher in *Symbiodinium* than *Cladocopium* genomes (pairwise estimated difference = -27.77 ± 5.91, $p = 0.01$; Supplementary Data 9). Furthermore, since we were unable to detect dinoRNAV EVEs in *Porites* metagenomes—a coral species primarily harboring *Cladocopium* C15 symbionts – we hypothesize that dinoRNAV endogenization was either less common in this lineage of Symbiodiniaceae or integrations have been lost over evolutionary time⁷⁹,⁸⁰.

**Incomplete ORFs and possible duplications indicate endogenization of DinoRNAVs.** The repeated observation of putative dinoRNAV EVEs in dinoflagellate scaffolds and contigs from metagenomes and genomes suggests these sequences are either (1) conserved sequence artifacts of Symbiodiniaceae-dinoRNAV interactions, and/or (2) evidence of highly prevalent dinoflagellate viruses, commonly integrated and propagated via their single-celled hosts. If the observed dinoRNAV-like sequences represent active infections capable of generating virions during egress, we would, at minimum, expect essential ORFs associated with replication (RNA-dependent RNA polymerase, *RdRp*) and virion structure (Major Capsid Protein, *MCP*) to be endogenized on the same scaffold. We would additionally expect to observe overall conservation of ORF length/composition (with a lack of internal stop codons or substantial deletions) when aligning the dinoRNAV-like sequences detected here with known exogenous dinoRNAV sequences.

However, both DIAMOND and gene prediction analyses generally depicted dinoRNAV-like ORFs in isolation on separate scaffolds. While 28 *MCP* and 73 *RdRp* dinoRNAV ORFs were annotated, both ORFs were present on a Symbiodiniaceae scaffold – potentially representing whole dinoRNAV genome integrations – in only 14 instances. Thirteen of these 14 were from *Symbiodinium* genomes, whereas one scaffold was from *Breviolum minutum*, a member of the second most ancestral dinoflagellate genus (Table 1)⁴⁸. To assess the conservation of putative dinoRNAV EVE sequence length/composition, we aligned the genomic and single ORF EVEs to reference exogenous dinoRNAV sequences. The reference genome for reef-associated dinoRNAVs is ~5 Kbp
<table><caption>Table 1 DinoRNAV EVE-like open reading frames from representative Symbiodiniaceae and Suesiaceae dinoflagellate scaffold-level genome assemblies.</caption><thead><tr><th rowspan="2">Dinoflagellate Species (strain)</th><th rowspan="2">Total # Scaffolds</th><th rowspan="2">Host</th><th rowspan="2">Location</th><th rowspan="2">BUSCO score</th><th colspan="3">dinoRNAV EVE ORFs on scaffolds</th></tr><tr><th>RdRp</th><th>MCP</th><th>Both</th></tr></thead><tbody><tr><td>Symbiodiniaceae</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Symbiodinium linuccheae (CCMP2456)<sup>83</sup></td><td>37,772</td><td>Plexaura hornamalla</td><td>Bermuda</td><td>21.8%</td><td>39</td><td>0</td><td>1</td></tr><tr><td>Symbiodinium microadriaticum (04-503SC1.03)<sup>83</sup></td><td>57,558</td><td>Orbicella faveolata</td><td>Florida, USA</td><td>41.6%</td><td>30</td><td>1</td><td>3</td></tr><tr><td>Symbiodinium microadriaticum (CassKB8)<sup>83</sup></td><td>67,937</td><td>Cassiopea sp.</td><td>Hawaii, USA</td><td>73.3%</td><td>29</td><td>1</td><td>3</td></tr><tr><td>Symbiodinium microadriaticum (CCMP2467)<sup>113</sup></td><td>9688</td><td>Stylophora pistillata</td><td>Red Sea</td><td>15.6%</td><td>29</td><td>1</td><td>3</td></tr><tr><td>Symbiodinium natans (CCMP2548)<sup>83**</sup></td><td>2855</td><td>N/A (isolated from seawater)</td><td>Hawaii, USA</td><td>15.5%</td><td>14</td><td>1</td><td>3</td></tr><tr><td>Symbiodinium necroappetens (CCMP2469)<sup>83*</sup></td><td>104,583</td><td>Condylactis gigantea</td><td>Jamaica</td><td>22.8%</td><td>37</td><td>4</td><td>2</td></tr><tr><td>Symbiodinium pilosum (CCMP2461)<sup>83**</sup></td><td>48,302</td><td>Zoanthus sociatus</td><td>Jamaica</td><td>19.8%</td><td>15</td><td>0</td><td>0</td></tr><tr><td>Symbiodinium sp. A5 (formerly S. tridacnidorum)<sup>83</sup></td><td>6245</td><td>Heliofungia actiniformis</td><td>Australia</td><td>21.1%</td><td>17</td><td>1</td><td>0</td></tr><tr><td>Symbiodinium tridacnidorum (sh18 A3 Y106)<sup>114</sup></td><td>16,176</td><td>Tridacna crocea</td><td>Japan</td><td>19.8%</td><td>20</td><td>1</td><td>0</td></tr><tr><td>Brevolium minutum (Mf1.05b)<sup>115</sup></td><td>21,899</td><td>Orbicella faveolata</td><td>Florida, USA</td><td>14.2%</td><td>21</td><td>3</td><td>1</td></tr><tr><td>Cladocopium C15<sup>116</sup></td><td>34,589</td><td>Porites lutea</td><td>Australia</td><td>11.6%</td><td>0</td><td>0</td><td>0</td></tr><tr><td>Cladocopium sp. Clacro (formerly C. goreaui)<sup>89</sup></td><td>41,289</td><td>Acropora tenuis</td><td>Australia</td><td>27.7%</td><td>4</td><td>0</td><td>0</td></tr><tr><td>Cladocopium sp C92 (Y103)<sup>114</sup></td><td>6686</td><td>Frugum sp.</td><td>Japan</td><td>19.5%</td><td>2</td><td>0</td><td>0</td></tr><tr><td>Durusdinium trenchii<sup>117</sup></td><td>19,593</td><td>Favia speciosa</td><td>Japan</td><td>28.7%</td><td>10</td><td>1</td><td>0</td></tr><tr><td>Fugacium kawagutti (CSI56 CCMP2468)<sup>89</sup></td><td>16,959</td><td>N/A (Free-living)</td><td>Hawaii, USA</td><td>8.3%</td><td>8</td><td>1</td><td>0</td></tr><tr><td>Fugacium kawagutti (CCMP2468)<sup>118</sup></td><td>30,040</td><td>N/A (Free-living)</td><td>Hawaii, USA</td><td>17.9%</td><td>9</td><td>1</td><td>0</td></tr><tr><td>n = 314 scaffolds with DinoRNAV EVE-like sequences</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Suesiaceae</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Polarella glacialis (CCMP1383)<sup>78**</sup></td><td>33,494</td><td>N/A (Free-living, isolated from seawater)</td><td>Antarctica</td><td>20.8%</td><td>20</td><td>0</td><td>0</td></tr><tr><td>Polarella glacialis (CCMP2088)<sup>78**</sup></td><td>37,768</td><td>N/A (Free-living, isolated from seawater)</td><td>Arctic</td><td>21.8%</td><td>18</td><td>0</td><td>0</td></tr><tr><td>n = 38 total scaffolds with DinoRNAV EVE-like sequences</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

Total counts of dinoflagellate scaffolds in genomes queried with individual endogenized dinoRNAV ORFs (RdRp, MCP) or both nearby each other, potentially indicating full viral genome integration. Table supplies associated dinoflagellate host species and location of isolation. RdRp = RNA-dependent RNA polymerase; MCP = major capsid protein. Assembly coverage and completeness are measured via BUSCO score (% completeness, or %C)<sup>102</sup>. Indicates species with documented opportunistic life history. ** Indicates species with documented free-living life history per principal species description. Gonzalez-Pech et al.<sup>83</sup>, Aranda et al.<sup>115</sup>, Shoguchi et al.<sup>114</sup>, Shoguchi et al.<sup>115</sup>, Liu et al.<sup>89</sup>, Shoguchi et al.<sup>117</sup>, Lin et al.<sup>118</sup>, Stephens et al.<sup>78</sup>. Further genome citations (including accession numbers) and BUSCO completion metrics can be found in Supplementary Data 7 and in superscript.
<figure><img src="image_6.png" /><figcaption>**Fig. 3** Amino acid alignment of putative endogenous and exogenous dinoRNAV-like +ssRNA virus Major Capsid Protein (MCP) sequences against transcriptome reference. Putatively endogenous dinoRNAV Major Capsid Protein (MCP) amino acid sequences were aligned against exogenous references, including: (1) Symbiodiniaceae +ssRNA virus MCP ORFs recovered from a *Cladocopium* sp. transcriptome (Levin et al, 2017), and (2) dinoRNAV MCP amplicons from fractionated coral tissue (Montalvo-Proaño et al. 2017). Conserved regions were observed between exogenous and putatively endogenous viral sequences (labeled Regions 1-3).</figcaption></figure>

long and contains a 1,071 bp noncoding region between ORFs, with a 124-nucleotide internal ribosomal binding site⁵⁷. In this study, for 13 of the scaffolds in which dinoRNAV ORFs were detected, the putative noncoding region between the *MCP* and *RdRp* EVEs ranged from ~200-800 bp (except for a scaffold belonging to *S. linuchea* CCMP2456, which contained a ~79 kbp noncoding region, and was excluded in further alignments). No internal ribosomal binding sites were detected within the putative dinoRNAV EVEs identified in dinoflagellate genomes. A nucleotide-based alignment to Levin et al.'s (2017)⁵⁷ reference dinoRNAV genome indicated that the putative dinoRNAV EVEs presented here contained substantial insertions and/or deletions (Supplementary Fig. 3). Translated exogenous dinoRNAV MCP ORFs are reported to be ~358 aa in length⁵⁷; Fig. 3 top sequences), but dinoRNAV-like MCP sequences recovered in this study ranged from 116-605aa in length. Furthermore, comparisons of these endogenous MCPs to exogenous reference sequences revealed internal stop codons and overall low similarity (Fig. 3). Amino acid-based alignment of endogenous dinoRNAV MCPs to metatranscriptome- and amplicon-generated exogenous reference sequences⁵⁷,⁵⁸ revealed indels and regions of low similarity between three conserved regions across both endogenous and exogenous MCP sequences (red boxes in Fig. 3).

Interestingly, multiple whole dinoRNAV integrations were sometimes observed in a single dinoflagellate genome. For example, genome assemblies of four different *S. microadriaticum* strains contained two or three whole dinoRNAV EVEs each (Table 1; Fig. 2). Pairwise alignments measuring shared nucleotide identity of whole dinoRNAV EVEs across Symbiodiniaceae scaffolds revealed that the *S. microadriaticum* genomes and the *S. necroappetens* genome share two whole genome dinoRNAV EVEs (provisionally dinoRNAV-A and dinoRNAV-

B; Supplementary Fig. 3; Clustal-Omega)⁸¹. *S. microadriaticum* dinoRNAV-B was identical in all strains and shared 97% identity with the *S. necroappetens* dinoRNAV-B, yet proximal genes varied (Supplementary Data 10, 11). Importantly, the inconsistent composition and fragmented nature of both the genomic and single ORF dinoRNAV EVEs reported here supports the hypothesis that these sequences are not capable of generating replicative virions and are best interpreted as multiple integrations of dinoRNAVs into a host genome.

**A Potential Mechanism for dinoRNAV Endogenization: Host-Provisioned Retroelements.** To assess if general genomic “neighborhoods” are conserved across dinoRNAV integrations (e.g. site location and synteny) and to better understand the genes proximal to EVEs on Symbiodiniaceae genomes, a chromosome-scale *Symbiodinium microadriaticum* genome assembly was evaluated (Fig. 4). The highest quality dinoflagellate genome assembly currently available revealed dinoRNAV-like ORFs on 18 of 94 chromosomes, with at least one *RdRp* on each, and some with multiple (two with *n* = 2 *RdRps*, three with *n* = 3 *RdRps*). On three of the chromosomes (# 30, 35, and 74), there were predicted ORFs annotated as dinoRNAV *MCPs* in close proximity to a *RdRp* ORF (separated by noncoding regions 319-656nt), indicative of a potential full-length dinoRNAV genome integration. These results corroborate detections of multiple genomic dinoRNAV EVEs in scaffold-scale assemblies of *Symbiodinium microadriaticum* genomes (Supplementary Fig. 3). The higher-resolution *S. microadriaticum* chromosome-level assembly facilitated the identification of an additional dinoRNAV genomic EVE (*n* = 4 for chromosome-level vs. *n* = 3 for scaffold-level, Supplementary Fig. 3), two of which were identified on Chromosome 74 and were separated by 2501 nucleotides. Of note,
<figure><img src="image_7.png" /><figcaption>**Fig. 4** Representative scaffolds and chromosome fragments containing putative dinoRNAVEs. Scaffolds annotation: MCP ORFs are indicated by light blue, *RdRp* indicated by navy blue ORFs (with complete description in Supplementary Data 10). Open reading frame (ORF) color broadly indicates cellular versus putative +ssRNA viral homology; yellow and some green (e.g., integrases, polyproteins) ORFs may be exploited mechanisms for viral integration. (+/-) base pair values represent sequence lengths between ORFs.</figcaption></figure>

Nand et al. (2021)⁸² reported a decreasing abundance and expression of genes towards the center of chromosomes (past ~2Mpb of a telomere), where there was an increase in repetitive elements; this is where 26 of 29 putative dinoRNAVEs were identified in the chromosome-level assembly. Furthermore, ORFs neighboring integrations often varied widely, both in proximity and predicted function, from collagen and RNA binding protein to reverse transcriptase and non-LTR retrotransposable elements. These ORFs potentially contributed to the endogenization of dinoRNAV via mechanisms such as retrotransposition (Fig. 4, Supplementary Data 10).

Retroposition through host-provisioned retroelements is one proposed mechanism of non-retroviral RNA virus integration into eukaryotic genomes⁶,⁷. An indicator of this form of integration is the nearby presence of a relict dinoflagellate spliced

leader (“dinoSL”), a 22nt sequence located at the 5’ end of mRNAs⁸³⁻⁸⁶. Such a sequence flanks the *RdRp* gene on some extant dinoRNAVs⁵⁷. We detected dinoSLs within 500 bp of 23.1% (six of 26) endogenized *RdRp* ORFs on *S. micoradriaticum* chromosomes, providing support for retroposition of these viral elements into Symbiodiniaceae genomes (Supplementary Data 11, 12). DinoRNAV gene integration may be facilitated by any of three major orders of retroelements associated with Symbiodiniaceae, including long terminal repeat (LTR) retrotransposons, short interspersed nuclear elements (SINEs), and long interspersed nuclear elements (LINEs⁸³,⁸⁷,⁸⁸). Evidence suggests that these LINEs are common and non-active remnants of an ancient proliferation of LINEs that preceded the diversification of *Suessiales*⁷⁸,⁸³,⁸⁹. *Symbiodinium* contains more LINEs relative to other Symbiodiniaceae genera, comprising 74.10-171.31 Mbp of
*Symbiodinium* genomes, relative to an average of 7.48 Mbp of the genomes of in other genera, indicating the loss of these retroelements across speciation events⁸²,⁸³. The loss of LINEs in more recently derived Symbiodiniaceae genera coincides with a decrease in dinoRNAV EVE detection in these genomes (Table 1). Conversely, the genomes of *Polarella*, the psychrophilic and free-living outgroup from which Symbiodiniaceae diversified ~160 million years ago, are LINE-rich and generally have comparable numbers of dinoRNAV EVEs to *Symbiodinium* (Table 1⁴⁸,⁷⁷,⁷⁸,⁸³). Together, this suggests that LINE activity during speciation may have facilitated dinoRNAV integration and the resulting EVEs may constitute dinornavirus “fossils.” This may explain their degree of sequence fragmentation and relatively low sequence similarity to modern extant dinoRNAVs (Fig. 3).

LINE-mediated retroposition is further supported by the observation of a LINE reverse transcriptase homolog ~17 kbp upstream of a *RdRp* EVE with a relict dinoSL on chromosome 45 (Supplementary Data 12) and a LINE retroelement 95 bp downstream of an EVE recovered from a *Pocillopora* metagenome (Fig. 4). Additionally, ~40% of annotated ORFs (35 of 88 annotated proteins) proximal to dinoRNAV ORFs on *S. microadriaticum* chromosomes were similar to non-LTR elements seen in other eukaryotic genomes sometimes <300 bp 5' upstream (Supplementary Data 12, Supplementary Fig. 4). Collectively, these findings implicate host provisioned retroelements, such as LINEs, as facilitators of dinoRNAV gene integration.

**DinoRNAV EVEs show homology to extant exogenous viruses.**
Modern, exogenous dinoRNAVs (Order: *Sobelivirales*) are highly divergent and hypothesized to form chronic infections within dinoflagellate hosts⁵⁴,⁵⁵,⁵⁷,⁵⁸. This chronic infection strategy likely provides opportunities for retroelement-driven endogenization into host genomes. Because many EVEs evolve at the rate of the host genome, rather than at the much faster rate of exogenous +ssRNA viral genomes, EVEs can serve as a snapshot of viral ancestry⁹⁰. We compared translated dinoRNAV EVEs to exogenous dinoRNAVs and other *Dinornavirus* taxa to assess the conservation of EVEs, the potential for host utilization of these elements, and their relatedness to contemporary dinoRNAVs. We found that amino acid translations of endogenous dinoRNAV *MCP* sequences contained conserved motifs observed in the exogenous *MCP* sequences (e.g. Regions 1–3 in Fig. 3), yet the associated phylogeny was highly polyphyletic along inferred ancestral nodes (Fig. 5a). Endogenous *MCP* ORFs also appear to be evolving under neutral selection (dN/dS=0.958).

Endogenized dinoRNAV *MCP* form their own clades within the *MCP* tree, each closely related to specific clades consisting of extant dinoRNAVs or environmental (i.e. unclassified) sobeliviruses with similar conserved motifs. The majority of dinoRNAV *MCP* EVEs shared similarity to extant *MCPs* identified from unfractionated stony coral holobionts via amplicon sequencing⁵⁸; these sequences formed an independent, disorganized clade (Fig. 5a clade containing yellow and blue sequences), relative to those recovered from dinoflagellate transcriptomes or those of other invertebrate hosts. Likewise, dinoRNAV *RdRp* EVEs identified via metagenomics appear most similar to HcRNAV, the defining member of family Alvernaviridae and a protist pathogen, further supporting the affiliation of this EVE with a dinoflagellate host. *MCP* and *RdRp* ORFs putatively derived from the same dinoflagellate genomes often shared clades (clades containing multiple blue or green sequences in Fig. 5a, b), perhaps indicative of duplications within genomes or multiple integration events of particular dinoRNAV lineages within host genera. The detection of putative dinoRNAV *RdRp* ORFs within *Polarella* genomes is therefore indicative of either the antiquity of dinoRNAV-dinoflagellate interactions and/or a propensity for recent dinoRNAV integration across Dinophyceae families. However, the

exclusion of the *P. glacialis* dinoRNAV-*RdRp* from *RdRps* of other dinoflagellate clades (pink, Fig. 5B) further illustrates the congruence between EVEs and their host genomes. Overall, the evident homology to contemporary *Dinornaviruses* support these integrations as Alvernaviridae within order *Sobelivirales*.

The expression and functional potential of endogenized dinoRNAV elements (if any) remains unclear. With no isolated Symbiodiniaceae-infecting dinoRNAV strains available, investigation into EVE functionality is limited to in silico approaches. Sequence data mining efforts identified RNA sequences either sharing sequence similarity with dinoRNAVs, or containing whole dinoRNAV-like ORFs that also annotated as dinoflagellate transcripts (i.e. with cellular ORFs or sequence similarity) in seven out of nine publicly accessible dinoflagellate transcriptomes (Supplementary Data 13). Additionally, two transcripts from an exogenous dinoRNAV infection identified in *Cladocopium* transcriptomes carried *MCP* ORFs of +ssRNA viral sequences ('TR74740_c13-g1 Il' and 'TR74740_c13-g1_i2'⁵⁷, red text in Fig. 5a) and form a clade with putative *Symbiodinium* dinoRNAV EVEs (Fig. 5). Likewise, the *RdRp* ORF of 'TR74740_c13-g1 Il' and the *RdRp* of 'GAKY01194223.1'— a transcript derived from a cultured *Symbiodinium microadriaticum* A1 transcriptome—shared some areas of similarity to putative endogenous dinoRNAVs (Fig. 5b⁵⁷,⁹¹). Importantly, both RNA transcripts also shared features characteristic of dinoflagellates, such as a 5' dinoSL⁶¹ or dinoflagellate sequence space flanking the dinoRNAV itself⁹¹. Furthermore, 'TR74740_c13-g1 Il' appeared to be in the top 0.03% of expressed transcripts at under certain thermal conditions, and GAKY01194223.1 appeared to exhibit moderately differential expression at the extremes of temperature and ionic stress in a cultured host⁵⁷,⁹¹.

While viral *RdRps* have been leveraged by eukaryotes in multiple pathways⁹², the apparent fragmentation of the putative dinoRNAV EVEs in silico may indicate a role in triggering antiviral mechanisms within their hosts³¹,⁹³. Given that the *Symbiodinium* genome contains all core RNAi protein machinery, including Argonaute and Dicer, and that GAKY01194223.1 folds into several hairpins (ΔG = -142.5 kcal/mol; Supplementary Fig. 5 examples), Symbiodiniaceae may use the putative EVE ncRNA identified here to develop host immunity against extant, exogenous dinoRNAVs. Furthermore, Symbiodiniaceae harboring dinoRNAV EVEs also contained numerous non-retroviral EVEs of other viral families (Supplementary Data 11, Fig. 7) in close proximity, such as *Herpesviridae*, *Baculoviridae*, *Poxviridae*, *Iridoviridae*, *Phycodnaviridae*, *Pandoraviridae* and *Pithoviridae*, ssDNA viruses of the family *Shotokuvirae*, -ssRNA viruses from the family *Rhabdoviridae* and +ssRNA viruses from the family *Coronaviridae* (Supplementary Fig. 6). Metagenomes corroborate findings of similar *RdRps* from these viral families (Supplementary Fig. 6). This provides support for host-mediated integration (e.g. retroposition) as a means of defense for single celled organisms, though further research is needed⁹⁴.

**Conclusions.** Over recent decades, endogenous viral elements (EVEs) have enabled investigators to better understand the evolutionary history of viruses (“paleovirology”) in diverse terrestrial systems, uncovering ancient and modern virus-host interactions. Our study further demonstrates how in silico identification of EVEs can provide ecological context for enigmatic viral genomes in non-model, multipartite systems such as coral holobionts, impacting how we study coral reefs and their viral consortia. Here, we detected heritable integrations of multiple putative dinoRNAV genes in Symbiodiniaceae scaffolds from cnidarian metagenomes, as well as in diverse genomes of cultured Symbiodiniaceae; no integrations were detected from seawater metagenomes nor diverse aposymbiotic cnidarian genomes. The
a

<figure><img src="image_8.png" /></figure>

b

<figure><img src="image_9.png" /></figure>

apparent pervasive nature of dinoRNAV-like sequences among dinoflagellate genomes (especially the genus *Symbiodinium*) suggests widespread and recurrent/ancestral integration of these EVEs. We propose that host-provisioned mechanisms drive dinoRNAV integration into single-celled dinoflagellate genomes as EVEs. The findings presented in this study further validate the dinoRNAV-Symbiodiniaceae virus-host pair, enhancing our

understanding of ecologically and economically important cnidarian holobionts and opening the door to examining the role of EVEs in reef health.

**Methods**

**Identification and computational validation of dinoRNAV EVEs leveraging meta'omics.** The Tara Pacific Expedition (2016-2018) sampled coral reefs to
**Fig. 5** Phylogenies of dinoRNAV major capsid protein (MCP) and RNA-dependent RNA polymerase (RdRp) ORFs. **a** dinoRNAV Major Capsid Protein (MCP) ORF phylogeny. Maximum-likelihood tree of MCP amino acid sequences generated with a LG + F + G4 substitution model and 50,000 parametric bootstraps, illustrating the similarity of putative dinoRNAV EVEs (this study) to extant dinoRNAVs from stony coral colonies. **b** RNA-dependent RNA polymerase (RdRp) ORF phylogeny. Maximum-likelihood tree of RdRp amino acid sequences generated with a Blosum62 + G4 substitution model and 50,000 parametric bootstraps, demonstrating the similarity of metagenomic dinoRNAV EVE RdRps to RdRps of the sole recognized *Dinornavirus*, *Heterocapsa circularisquama* RNA virus (HcRNAV), as well as alignment to each other. ORFs were recovered from host metagenomes, transcriptomes, genomes, and extant +ssRNA reference viruses from amplicon libraries (a only). Both trees include *Dinornavirus* reference sequences and visualized in iTOL.

investigate reef health and ecology using multiple methods, including amplicon sequencing and metagenomics (see Pesant et al. 2020⁹⁵ and https://doi.org/10.5281/zenodo.4068293 for coral reef sampling and processing methods). In this study, we explored metagenomes generated from hydrocorals (n = 60 *Millepora*), stony corals (n = 108 *Porites*, n = 101 *Pocillopora*) sampled from 11 islands (three replicate sites per island) across the South Pacific Ocean during the Tara Pacific Expedition for dinoRNAV EVES (Fig. 1, Supplementary Data 1A, 1B⁹⁵). Amplicon libraries of the dinoflagellate *Internal Transcribed Spacer 2* (ITS2) gene fragment were sequenced in tandem with the metagenomes, to characterize the dominant Symbiodiniaceae harbored by hydrozoan and stony coral colonies⁹⁵.

To confirm that these dinoRNAV EVE sequences were affiliated with coral holobionts and reduce the possibility that they are technical artifacts, publicly available metagenome libraries were analyzed (Supplementary Data 1B). These additional libraries included 120 assembled pelagic water samples presumed to include pelagic dinoflagellate sequences from the Tara Oceans dataset (2009-2013⁷⁰) and 30 MiSeq metagenomes from unfractionated samples of the stony coral genus *Acropora*, which were processed and sequenced via a different pipeline (Supplementary Data 1B, Supplementary Fig. 7). Publicly accessible transcriptomes from nine Symbiodiniaceae assemblies (Supplementary Data 1B) were also queried to determine if dinoRNAV-like sequences were present in poly(A)-selected dinoflagellate transcriptomes and resembled EVEs in terms of proximal gene composition and presence of a characteristic pre-mRNA spliced leader (dinoSL) sequence (as in Levin et al, 2017⁵⁷). Details regarding the collection of samples, generation of metagenomes and associated Symbiodiniaceae amplicon libraries, and associated bioinformatic analyses are provided in Supplementary Fig. 7).

Metagenomic and transcriptomic scaffolds were annotated against a curated database of dinoRNAV-like sequences (Supplementary Data 14) via BLASTx (e-value < 1 × 10⁻⁵; see Supplementary Fig. 7 for workflow⁹⁶). Alignments to the custom database with a bit score <50 and percent shared amino acid identity <30% were excluded from further analysis. A length penalty was not imposed during this step due to the limited length of assembled scaffolds (average N50 = 3341 ± 127 nt across all queried libraries). Open reading frames (ORFs) from selected scaffolds were called via Prodigal (v2.6.3⁹⁷) and annotated against the NCBI-nr database (e-value < 0.001; DIAMOND v2.0.6⁹⁸) to confirm homology to dinoRNAVs and to identify adjacent dinoflagellate sequences (e-value < 1 × 10⁻⁵, bit≥50). In the absence of complete ORFs (potentially due to the limited size of scaffolds, partial integrations, etc.), homology was confirmed through comparison of the initial alignments to the curated database and 300nt of upstream/downstream flanking sequences (bedtools v2.30.⁹⁹) against the NCBI-nr database (e-value < 0.001; DIAMOND v2.0.6⁹⁸). This served as further curation and verification, as EVEs can exist in fragmented or degraded states. Non-normalized quality-controlled reads were mapped via bbmap (v.38.84¹⁰⁰), and putative EVEs were assessed for uniform read coverage across scaffolds, reducing the probability of chimeric assembly. RNA secondary structure was predicted via mfold (v.3.5¹⁰¹).

**dinoRNAV EVEs in dinoflagellate and aposymbiotic cnidarian genomes.**

Publicly available dinoflagellate and aposymbiotic (dinoflagellate-free) cnidarian genome assemblies were queried to resolve the putative host(s) of dinoRNAVs, to assess homology among detected dinoRNAVs within coral holobionts, and to compare genes proximal to dinoRNAV EVEs in different host species/strains. A chromosome-scale dinoflagellate genome assembly generated from a *Symbiodinium microadriaticum* culture (Accession: GSE152150)⁸², and scaffold-scale genome assemblies were examined for dinoRNAV EVES (Supplementary Data 1B, Supplementary Fig. 7). Scaffold-scale genome assemblies were from the closely related families Symbiodiniaceae and Suessiaceae, and included representatives from the genera *Symbiodinium* (n = 9), *Breviolum* (n = 1), *Cladocopium* (n = 3), *Durusdinium* (n = 1), *Fugacium* (n = 2), and *Polarella* (n = 2), as well as 25 aposymbiotic cnidarian genome assemblies, including the stony coral genera *Acropora* (n = 13), *Astreopora* (n = 1), *Galaxea* (n = 1), *Montastraea* (n = 1), *Montipora* (n = 3), *Orbicella* (n = 1), *Pocillopora* (n = 2), *Porites* (n = 1), and *Stylophora* (n = 1), and the jellyfish *Clytia* (n = 1; Fig. 2, Supplementary Data 1B). All publicly available genome assemblies had undergone a form of microbial decontamination, trimming, and quality control prior to assembly, minimizing risk of microbial contamination. Genome completeness and quality further were assessed via BUSCO (v3)¹⁰² with the Eukaryota dataset and QUAST (v5.0.2¹⁰³), respectively. Scaffolds/chromosomes containing putative dinoRNAV EVEs were

identified by aligning sequences to the protein version of the Reference Viral DataBase (RVDB v.19¹⁰⁴) using DIAMOND BLASTx (v0.9.30)⁹⁸. The same exclusion criteria were maintained for alignments of metagenomic scaffolds, also omitting alignments <100 amino acids. Regions of dinoflagellate genomes exhibiting similarity to the *MCP* or *RdRp* of reef-associated dinoRNAV reference genomes⁵⁷ or other closely related +ssRNA viruses (Supplementary Data 14) were extracted and re-aligned to the NCBI-nr database to further confirm viral homology.

We tested the relationship between the number of identified dinoRNAV EVE-containing scaffolds, dinoflagellate genera, and genome quality metrics using a linear model. Model selection was performed with an F-test (package car, v.3.0-12) and assumptions were visually checked. Pairwise comparisons between genera were conducted using the package emmeans (v.1.7.2). Putative whole dinoRNAV-like genomes within scaffolds were identified based on the presence of *MCP* and *RdRp*-like sequences on the same scaffold no further than 1.5 Kbp apart (Table 1; Supplementary Fig. 3). IRESPred¹⁰⁵ was utilized to identify internal ribosomal entry sites (IRES) with default parameters on putative dinoRNAV EVE with whole sequence integrations.

ORFs were predicted and annotated from dinoRNAV EVE-containing scaffolds and all dinoflagellate chromosomes using Prodigal⁹⁷ and MAKER2 annotation pipeline¹⁰⁶ with the AUGUSTUS gene prediction software¹⁰⁷. Translated ORFs were then aligned to a hybrid database containing the UniProt/Swiss-Prot database and protein version of RVDB (v.19; DIAMOND-BLASTp). ORFs on putative dinoRNAV EVE-containing scaffolds and chromosomes were further annotated using InterProScan (v5.4.8-83.0, Pfam analysis with default parameters) to identify sequences proximal to putative dinoRNAV integrations. The presence of dinoflagellate spliced leaders ("dinoSLs") were examined within 500nt of dinoRNAV EVEs using BLASTn with default parameters (except word size=9, excluding two ambiguous positions as specified in Gonzalez-Pech et al. 2021⁸³).

**Phylogenetic analysis of dinoRNAV EVES.** Amino acid-based phylogenetic trees were generated with dinoRNAV EVE ORFs (MCP and RdRp) from scaffold-scale genomic assemblies, metagenomes, transcriptomes, and sequences from exogenous and closely related +ssRNA reference viruses (Supplementary Data 1A, B, Supplementary Data 14). Sequences were aligned using the best fit algorithm determined by MAFFT (v7.464)¹⁰⁸ and reviewed and trimmed manually in MEGA (v7)¹⁰⁹. Maximum-likelihood trees were generated with IQTREE2¹¹⁰ using the model determined by ModelFinder¹¹¹ and 50,000 parametric bootstraps¹¹² with nearest neighbor interchange optimization. ORFs from the chromosome-level assembly for S. microadriaticum culture CCMP2467 were not included in the phylogeny in order to avoid redundancy with those from the analogous scaffold-level assembly. To calculate dN/dS, ORFs were aligned in Clustal Omega (v.1.2.4), refined in MUSCLE (v.3.6), before using pal2nal (v.14) for codon-based nucleic acid alignment. Evolutionary trajectory was then assessed via CODEML (PAML package, v.4.10.5).

**Statistics and reproducibility.** As indicated throughout the article, metagenomes (n = 269) and genomes (n = 18) served as technical replicates, and ORFs or full EVE sequences served as comparative evoevolutionary units (replicates described in Table 1) when available. Negative controls (seawater metagenomes, coral host, etc) were also evaluated. All statistical packages are reported in methods or Supplementary Fig. 7.

**Reporting summary.** Further information on research design and collection permits are available in the Nature Research Reporting Summary linked to this article.

**Data availability**

Metadata are accessible in zenodo: https://zenodo.org/record/6299409#Y.-ClwuzMKml. Metagenomes are available via https://doi.org/10.5281/zenodo.7839794. Seawater metagenomes are available through the European bioinformatics institute (Tara Oceans; ERP001736) and NCBI (PRJEB1787). NCBI accession numbers for individual holobiont species metagenomes, genome assemblies and reference sequences can be found in Supplementary Data 1B, 3 and 14, respectively.
Received: 20 April 2022; Accepted: 24 April 2023;
Published online: 01 June 2023

References

1. Johnson, W. E. Endogenous retroviruses in the genomics era. *Annu. Rev. Virol.* **2**, 135–159 (2015).

2. Johnson, W. E. Origins and evolutionary consequences of ancient endogenous retroviruses. *Nat. Rev. Microbiol.* **17**, 355–370 (2019).

3. Johnson, W. E. Endless forms most viral. *PLoS Genet.* **6**, e1001210 (2010).

4. Stoye, J. P. Studies of endogenous retroviruses reveal a continuing evolutionary saga. *Nat. Rev. Microbiol.* **10**, 395–406 (2012).

5. Gallot-Lavallée, L. & Blanc, G. A glimpse of nucleo-cytoplasmic large DNA virus biodiversity through the eukaryotic genomics window. *Viruses* **9**, 17 (2017).

6. Flynn, P. J. & Moreau, C. S. Assessing the diversity of endogenous viruses throughout ant genomes. *Front Microbiol* **10**, 1139 (2019).

7. Horie, M. et al. Endogenous non-retroviral RNA virus elements in mammalian genomes. *Nature* **463**, 84–87 (2010).

8. Katzourakis, A. & Gifford, R. J. Endogenous viral elements in animal genomes. *PLOS Genet.* **6**, e1001191 (2010).

9. Chiba, S. et al. Widespread endogenization of genome sequences of non-retroviral RNA viruses into plant genomes. *PLOS Pathog.* **7**, e1002146 (2011).

10. Chu, H., Jo, Y. & Cho, W. K. Evolution of endogenous non-retroviral genes integrated into plant genomes. *Curr. Plant Biol.* **1**, 55–59 (2014).

11. Kojima, S. et al Virus-like insertions with sequence signatures similar to those of endogenous nonretroviral RNA viruses in the human genome. *Proc Natl Acad Sci.* **118**, https://doi.org/10.1073/pnas.2010758118 (2021).

12. Ballinger, M. J., Bruenn, J. A. & Taylor, D. J. Phylogeny, integration and expression of sigma virus-like genes in. *Drosoph. Mol. Phylogenet Evol.* **65**, 251–258 (2012).

13. Tromas, N., Zwart, M. P., Forment, J. & Elena, S. F. Shrinkage of genome size in a plant RNA virus upon transfer of an essential viral gene into the host genome. *GBE* **6**, 538–550 (2014).

14. Palatini, U. et al. Comparative genomics shows that viral integrations are abundant and express piRNAs in the arboviral vectors *Aedes aegypti* and *Aedes albopictus*. *BMC Genomics* **18**, 512 (2017).

15. Wang, L. et al. Endogenous viral elements in algal genomes. *Acta Ocean. Sin.* **33**, 102–107 (2014).

16. Jebb, D. et al. Six reference-quality genomes reveal evolution of bat adaptations. *Nature* **583**, 578–584 (2020).

17. Moniruzzaman, M., Weinheimer, A. R., Martinez-Gutierrez, C. A. & Aylward, F. O. Widespread endogenization of giant viruses shapes genomes of green algae. *Nature* **588**, 141–145 (2020).

18. Skirmuntt, E. C., Escalera-Zamudio, M., Teeling, E. C., Smith, A. & Katzourakis, A. The potential role of Endogenous Viral Elements in the evolution of bats as reservoirs for zoonotic viruses. *Annu. Rev. Virol.* **7**, 103–119 (2020). 092818–015613.

19. Roossinck, M. J. The good viruses: viral mutualistic symbioses. *Nat. Rev. Microbiol.* **9**, 99–108 (2011).

20. Harrison, E. & Brockhurst, M. A. Ecological and evolutionary benefits of temperate phage: what does or doesn't kill you makes you stronger. *BioEssays* **39**, 1700112 (2017).

21. Correa, A. M. S. et al. Revisiting the rules of life for viruses of microorganisms. *Nat. Rev. Microbiol.* **19**, 501–513 (2021).

22. Jern, P. & Coffin, J. M. Effects of retroviruses on host genome function. *Annu Rev. Genet.* **42**, 709–732 (2008).

23. Oliveira, N. M., Satija, H., Kouwenhoven, I. A. & Eiden, M. V. Changes in viral protein function that accompany retroviral endogenization. *Proc. Natl Acad. Sci.* **104**, 17506–17511 (2007).

24. Feschotte, C. & Gilbert, C. Endogenous viruses: insights into viral evolution and impact on host biology. *Nat. Rev. Genet.* **13**, 283–296 (2012).

25. Frank, J. A. & Feschotte, C. Co-option of endogenous viral sequences for host cell function. *COVIRO* **25**, 81–89 (2017).

26. Mortelmans, K., Wang-Johanning, F. & Johannning, G. L. The role of human endogenous retroviruses in brain development and function. *APMIS* **124**, 105–115 (2016).

27. Sofuku, K., & Honda, T. Influence of endogenous viral sequences on gene expression. *IntechOpen*. https://doi.org/10.5772/intechopen.71864 (2018)

28. Takahashi, H., Fukuhara, T., Kitazawa, H., & Kormelink, R. Virus latency and the impact on plants. *Front. Microbiol.* **10** https://doi.org/10.3389/fmicb.2019.02764 (2019).

29. Whitfield, Z. J. et al. The diversity, structure, and function of heritable adaptive immunity sequences in the *Aedes aegypti* genome. *Curr. Biol.* **27**, 3511–3519.e7 (2017).

30. ter Horst, A. M., Nigg, J. C., Dekker, F. M. & Falk, B. W. Endogenous viral elements are widespread in arthropod genomes and commonly give rise to PIWI-interacting RNAs. *J. Virol.* **93**, e02124–18 (2019).

31. Suzuki, Y. et al. Non-retroviral Endogenous Viral Element limits cognate virus replication in *Aedes aegypti* ovaries. *Curr. Biol.* **30**, 3495–3506.e6 (2020).

32. Aswad, S. & Katzourakis, A. Paleovirology and virally derived immunity. *Trends Ecol. Evol.* **27**, 627–36 (2012).

33. Parker, B. J. & Brisson, J. A. A laterally transferred viral gene modifies aphid wing plasticity. *Curr. Biol.* **29**, 2098–2103.e5 (2019).

34. Wilson, W., Francis, I., Ryan, K. & Davy, S. Temperature induction of viruses in symbiotic dinoflagellates. *Aquat. Micro. Ecol.* **25**, 99–102 (2001).

35. Aiewsakun, P. & Katzourakis, A. Endogenous viruses: connecting recent and ancient viral evolution. *Virology* **479–480**, 26–37 (2015).

36. Belyi, V. A., Levine, A. J. & Skalka, A. M. Sequences from ancestral single-stranded DNA viruses in vertebrate genomes: the parvoviridae and circoviridae are more than 40 to 50 million years old. *J. Virol.* **84**, 12458–62 (2010).

37. Gilbert, C. & Feschotte, C. Genomic fossils calibrate the long-term evolution of hepadnaviruses. *PLoS Biol.* **8**, e1000495 (2010).

38. Kawasaki, J., Kojima, S., Mukai, Y., Tomonaga, K. & Horie, M. 100-My history of bornavirus infections hidden in vertebrate genomes. *Proc. Natl Acad. Sci.* **118**, e2026235118 (2021).

39. Li, Y. Q. et al. Discovery of Flaviviridae-derived endogenous viral elements in shrew genomes provide novel insights into Pestivirus ancient history. *bioRxiv.* 02.11.480044. https://doi.org/10.1101/2022.02.11.480044 (2022)

40. Cui, J. & Holmes, E. C. Endogenous lentiviruses in the ferret genome. *J. Virol.* **86**, 3383–5 (2012).

41. Keckesova, Z., Ylinen, L. M., Towers, G. J., Gifford, R. J. & Katzourakis, A. Identification of a RELIK orthologue in the European hare (Lepus europaeus) reveals a minimum age of 12 million years for the lagomorph lentiviruses. *Virology* **5**, 7–11 (2009). 384.

42. Katzourakis, A., Gifford, R. J., Tristem, M., Gilbert, M. T. & Pybus, O. G. Macroevolution of complex retroviruses. *Science* **18**, 1512 (2009). 325.

43. Katzourakis, A. Paleovirology: inferring viral evolution from host genome sequence data. *Philos. Trans. R. Soc.* **368**, 20120493 (2013).

44. Patel, M. R., Emerman, M. & Malik, H. S. Paleovirology—ghosts and gifts of viruses past. *COVIRO* **1**, 304–309 (2011).

45. Barreat, J. G. N. & Katzourakis, A. Paleovirology of the DNA viruses of eukaryotes. *Trends Microbiol* **30**, 281–292 (2022).

46. Knowlton, N. & Rohwer, F. Multispecies microbial mutualisms on coral reefs: the host as a habitat. *Am. Nat.* **162**, S51–S62 (2003).

47. Matthews, J. L. et al. Symbiodiniaceae-bacteria interactions: rethinking metabolite exchange in reef-building corals as multi-partner metabolic networks. *Environ. Microbiol.* **22**, 1675–1687 (2020).

48. Lajeunesse, T. C. et al. Systematic revision of symbiodiniaceae highlights the antiquity and diversity of coral endosymbionts. *Curr. Biol.* **28**, 2570–2580.e6 (2018).

49. Glynn, P. W. Coral reef bleaching: facts, hypotheses and implications. *Glob. Change Biol.* **2**, 495–509 (1996).

50. van Oppen, M. J. H., Leong, J.-A. & Gates, R. D. Coral-virus interactions: a double-edged sword? *Symbiosis* **47**, 1–8 (2009).

51. Correa, A. M. S. et al. Viral outbreak in corals associated with an in situ bleaching event: atypical herpes- like viruses and a new Megavirus infecting. *Symbiodinium*. *Front Microbiol* **7**, 127 (2016).

52. Vega Thurber, R., Payet, J. P., Thurber, A. R. & Correa, A. M. S. Virus-host interactions and their roles in coral reef health and disease. *Nat. Rev. Microbiol* **15**, 205–216 (2017).

53. Messyasz, A. et al. Coral bleaching phenotypes associated with differential abundances of Nucleocytoplasmic Large DNA Viruses. *Front Mar. Sci.* **7**, 555474 (2020).

54. Grupstra, C. G. B. et al. Thermal stress triggers productive viral infection of a key coral reef symbiont. *ISME J.* **16**, 1430–1441 (2022).

55. Correa, A. M. S., Welsh, R. M. & Vega Thurber, R. L. Unique nucleocytoplasmic dsDNA and +ssRNA viruses are associated with the dinoflagellate endosymbionts of corals. *ISME J.* **7**, 13–27 (2013).

56. Weynberg, K. D., Wood-Charlson, E. M., Suttle, C. A. & van Oppen, M. J. H. Generating viral metagenomes from the coral holobiont. *Front Microbiol* **5**, 206 (2014).

57. Levin, R. A., Voolstra, C. R., Weynberg, K. D. & van Oppen, M. J. H. Evidence for a role of viruses in the thermal sensitivity of coral photosymbionts. *ISME J.* **11**, 808–812 (2017).

58. Montalvo-Proaño, J., Buerger, P., Weynberg, K. D. & van Oppen, M. J. H. A PCR-Based Assay targeting the major capsid protein gene of a Dinorna-Like ssRNA virus that infects coral photosymbionts. *Front. Microbiol* **8**, 1665 (2017).

59. Nagasaki, K. et al. Comparison of genome sequences of single-stranded RNA viruses infecting the bivalve-killing dinoflagellate *Heterocapsa circularisquama*. *Appl. Environ. Microbiol* **71**, 8888–8894 (2005).
60. Lawrence, S. A., Davy, J. E., Aeby, G. S., Wilson, W. H. & Davy, S. K. Quantification of virus-like particles suggests viral infection in corals affected by *Porites* tissue loss. *Coral Reefs* **33**, 687–691 (2014).

61. Zhang, H., Zhuang, Y., Gill, J. & Lin, S. Proof that dinoflagellate Spliced Leader (DinoSL) is a useful hook for fishing dinoflagellate transcripts from mixed microbial samples: *Symbiodinium kawagutii* as a case study. *Protist* **164**, 510–527 (2013).

62. Voolstra, C. R. et al. Comparative analysis of the genomes of *Stylophora pistillata* and *Acropora digitifera* provides evidence for extensive differences between species of corals. *Sci. Rep.* **7**, 17583 (2017).

63. Buitrigo-López, C., Mariappan, K. G., Cárdenas, A., Gegner, H. M. & Voolstra, C. R. The Genome of the Cauliflower Coral *Pocillopora verrucosa*. *Genome Biol. Evol.* **1**, 12, 1911–1917 (2020).

64. Cunningham, R., Bay, R. A., Gillette, P., Baker, A. C. & Traylor-Knowles, N. Comparative analysis of the *Pocillopora damicornis* genome highlights role of immune system in coral evolution. *Sci. Rep.* **8**, 16134 (2018).

65. Helmkampf, M., Bellinger, M. R., Geib, S. M., Sim, S. B. & Takabayashi, M. Draft Genome of the Rice Coral Montipora capitata Obtained from Linked-Read Sequencing. *Genome Biol. Evol.* **1**, 11, 2045–2054 (2019).

66. Kitchen, S. A. et al. Genomic variants among threatened acropora corals. *G3 (Bethesda)* **9**, 1633–1646 (2019).

67. RefuGe 2020 Consortium. The RefuGe 2020 Consortium—using “omics” approaches to explore the adaptability and resilience of coral holobionts to environmental change. *Front. Mar. Sci.* **2**. https://doi.org/10.3389/fmars.2015.00068 (2015).

68. Shinzato, C. et al. Eighteen coral genomes reveal the evolutionary origin of acropora strategies to accommodate environmental changes. *Mol. Biol. Evol.* **4**, 16–30 (2021). 38.

69. Ying, H. et al. Comparative genomics reveals the distinct evolutionary trajectories of the robust and complex coral lineages. *Genome Biol.* **2**, 175 (2018). 19.

70. Tara Oceans Consortium Coordinators. et al. Open science resources for the discovery and analysis of Tara Oceans data. *Sci. Data* **2**, 150023 (2015).

71. Littman, R. A., van Oppen, M. J. H. & Willis, B. L. Methods for sampling free-living *Symbiodinium* (zooxanthellae) and their distribution and abundance at Lizard Island (Great Barrier Reef). *J. Exp. Mar. Biol. Ecol.* **364**, 48–53 (2008).

72. Scheufen, T., Iglesias-Prieto, R. & Enriquez, S. Changes in the number of symbionts and *Symbiodinium* cell pigmentation modulate differentially coral light absorption and photosynthetic performance. *Front. Mar. Sci.* **4**, 309 (2017).

73. Fujise, L. et al. Unlocking the phylogenetic diversity, primary habitats, and abundances of free-living Symbiodiniaceae on a coral reef. *Mol. Ecol.* **30**, 343–360 (2021).

74. Grupstra, C. G. B., Rabbitt, K. M., Howe-Kerr, L. I. & Correa, A. M. S. Fish predation on corals promotes the dispersal of coral symbionts. *Anim. Microbiome* **3**, 25 (2021).

75. Sunagawa, S. et al. Ocean plankton. Structure and function of the global ocean microbiome. *Science* **348**, 1261359 (2015).

76. Muller, E. M., Bartels, E. & Baums, I. B. Bleaching causes loss of disease resistance within the threatened coral species *Acropora cervicornis*. *eLife* **7**, e35066 (2018).

77. Janouškovec, J. et al. Major transitions in dinoflagellate evolution unveiled by phylotranscriptomics. *Proc. Natl Acad. Sci.* **114**, E171–E180 (2017).

78. Stephens, T. G. et al. Genomes of the dinoflagellate *Polarella glacialis* encode tandemly repeated single-exon genes with adaptive functions. *BMC Biol.* **18**, 56 (2020).

79. Tan, Y. T. R. et al. Endosymbiont diversity and community structure in *Porites lutea* from Southeast Asia are driven by a suite of environmental variables. *Symbiosis* **80**, 269–277 (2020).

80. Qin, Z. et al. Diversity of Symbiodiniaceae in 15 coral species from the Southern South China Sea: potential relationship with coral thermal adaptability. *Front Microbiol* **10**, 2343 (2019).

81. Sievers, F. et al. Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. *Mol. Syst. Biol.* **7**, 539 (2011).

82. Nand, A. et al. Genetic and spatial organization of the unusual chromosomes of the dinoflagellate *Symbiodinium microadriaticum*. *Nat. Genet.* **53**, 618–629 (2021).

83. González-Pech, R. A. et al. Comparison of 15 dinoflagellate genomes reveals extensive sequence and structural divergence in family Symbiodiniaceae and genus *Symbiodinium*. *BMC Biol.* **19**, 73 (2021).

84. Zhang, H., Hou, Y., Miranda, L. & Lin, S. Spliced leader RNA trans-splicing in dinoflagellates. *Proc. Natl Acad. Sci.* **104**, 4618–4623 (2007).

85. Lidie, K. B. & van Dolah, F. M. Spliced leader RNA-mediated trans-splicing in a dinoflagellate, *Karenia brevis*. *J. Eukaryot. Microbiol.* **54**, 427–35 (2007).

86. Song, B., Chen, S. & Chen, W. Dinoflagellates, a unique lineage for retrogene research. *Front. Microbiol.* **9**, 1556 (2018).

87. Elbarbay, R. A., Lucas, B. A., Maquat, L. E. Retrotransposons as regulators of gene expression. *Science.* **12**, 351(6274):aac7247 (2016)

88. Mita, P. & Boeke, J. D. How retrotransposons shape genome regulation. *Curr. Opin. Genet. Dev.* **37**, 90–100 (2016).

89. Liu, H. et al. *Symbiodinium* genomes reveal adaptive evolution of functions related to coral-dinoflagellate symbiosis. *Commun. Biol.* **1**, 95 (2018).

90. Holmes, E. C. The evolution of endogenous viral elements. *Cell Host Microbe* **10**, 368–377 (2011).

91. Baumgarten, S. et al. Integrating microRNA and mRNA expression profiling in *Symbiodinium microadriaticum*, a dinoflagellate symbiont of reef-building corals. *BMC Genomics* **14**, 704 (2013).

92. Lipardi, C. & Paterson, B. M. Identification of an RNA-dependent RNA polymerase in *Drosophila* involved in RNAi and transposon suppression. *Proc. Natl Acad. Sci.* **106**, 15645–15650 (2009).

93. Blair, C. D., Olson, K. E., & Bonizzoni, M. The Widespread occurrence and potential biological roles of endogenous viral elements in insect genomes. *Curr. Issues Mol. Biol.* **34**, 13–30 (2020).

94. Yan, N. & Chen, Z. Intrinsic antiviral immunity. *Nat. Immunol.* **13**, 214–222 (2012).

95. Pesant, S. et al. (2020). Tara Pacific samples provenance and environmental context—version 2. Zenodo. https://doi.org/10.5281/zenodo.4068293 (2020).

96. Altschul, S. F., Gish, W., Miller, W., Myers, E. W. & Lipman, D. J. Basic local alignment search tool. *J. Mol. Biol.* **215**, 403–410 (1990).

97. Hyatt, D. et al. Prodigal: prokaryotic gene recognition and translation initiation site identification. *BMC Bioinform.* **11**, 119 (2010).

98. Buchfink, B., Xie, C. & Huson, D. H. Fast and sensitive protein alignment using DIAMOND. *Nat. Methods* **12**, 59–60 (2015).

99. Quinlan, A. R. & Hall, I. M. BEDTools: a flexible suite of utilities for comparing genomic features. *Bioinformatics* **26**, 841–842 (2010).

100. Bushnell, B., Rood, J. & Singer, E. BBMerge—accurate paired shotgun read merging via overlap. *PLOS ONE* **12**, e0185056 (2017).

101. Zuker, M. Mfold web server for nucleic acid folding and hybridization prediction. *Nucleic Acids Res.* **31**, 3406–3415 (2003).

102. Simão, F. A., Waterhouse, R. M., Ioannidis, P., Kriventseva, E. V. & Zdobnov, E. M. BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. *Bioinformatics* **31**, 3210–3212 (2015).

103. Gurevich, A., Saveliev, V., Vyahhi, N. & Tesler, G. QUAST: quality assessment tool for genome assemblies. *Bioinformatics* **29**, 1072–1075 (2013).

104. Bigot, T., Temmam, S., Pérot, P. & Eloit, M. RVDB-prot, a reference viral protein database and its HMM profiles. *F1000 Res.* **8**, 530 (2020).

105. Kolekar, P. et al. IRESPred: web server for prediction of cellular and viral internal ribosome entry site. *IRES. Sci. Rep.* **6**, 27436 (2016).

106. Holt, C. & Yandell, M. MAKER2: an annotation pipeline and genome-database management tool for second-generation genome projects. *BMC Bioinform.* **12**, 491 (2011).

107. Stanke, M. et al. AUGUSTUS: ab initio prediction of alternative transcripts. *Nucleic Acids Res.* **34**, W435–W439 (2006). (Web Server).

108. Katoh, K. & Standley, D. M. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. *Mol. Biol. Evol.* **30**, 772–780 (2013).

109. Kumar, S., Stecher, G. & Tamura, K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. *Mol. Biol. Evol.* **33**, 1870–1874 (2016).

110. Minh, B. Q. et al. IQ-TREE 2: new models and efficient methods for phylogenetic inference in the genomic era. *Mol. Biol. Evol.* **37**, 1530–1534 (2020).

111. Kalyaanamoorthy, S., Minh, B. Q., Wong, T. K. F., von Haeseler, A. & Jermiin, L. S. ModelFinder: fast model selection for accurate phylogenetic estimates. *Nat. Methods* **14**, 587–589 (2017).

112. Hoang, D. T., Chernomor, O., von Haeseler, A., Minh, B. Q. & Vinh, L. S. UFBoot2: improving the ultrafast bootstrap approximation. *Mol. Biol. Evol.* **35**, 518–522 (2018).

113. Aranda, M. et al. Genomes of coral dinoflagellate symbionts highlight evolutionary adaptations conducive to a symbiotic lifestyle. *Sci. Rep.* **6**, 39734 (2016).

114. Shoguchi, E. et al. Two divergent *Symbiodinium* genomes reveal conservation of a gene cluster for sunscreen biosynthesis and recently lost genes. *BMC Genomics* **19**, 458 (2018).

115. Shoguchi, E. et al. Draft assembly of the *Symbiodinium minutum* nuclear genome reveals dinoflagellate gene structure. *Curr. Biol.* **23**, 1399–1408 (2013).

116. Robbins, S. J. et al. A genomic view of the reef-building coral *Porites lutea* and its microbial symbionts. *Nat. Microbiol.* **4**, 2090–2100 (2019).

117. Shoguchi, E. et al. A new dinoflagellate genome illuminates a conserved gene cluster involved in sunscreen biosynthesis. *GBE* **13**, evaa235 (2021).

118. Lin, S. et al. The *Symbiodinium kawagutii* genome illuminates dinoflagellate gene expression and coral symbiosis. *Science* **350**, 691–694 (2015).

119. Hume, B. C. C. et al. SymPortal: A novel analytical framework and platform for coral algal symbiont next-generation sequencing *ITS2* profiling. *Mol. Ecol. Resour.* **19**, 1063–1080 (2019).
## Acknowledgements

Special thanks to the Tara Ocean Foundation, the R/V Tara crew and the Tara Pacific Expedition Participants (https://doi.org/10.5281/zenodo.3777760). Thank you also to Carsten G.B. Grupstra and Clark Hamor for their input regarding analyses. We are keen to thank the commitment of the following institutions for their financial and scientific support that made this unique Tara Pacific Expedition possible: CNRS, PSL, CSM, EPHE, Genoscope, CEA, Inserm, Université Côte d'Azur, ANR, agnès b., UNESCO-IOC, the Veolia Foundation, the Prince Albert II de Monaco Foundation, Région Bretagne, Billerudkorsnas, AmerisourceBergen Company, Lorient Agglomération, Oceans by Disney, L'Oréal, Biotherm, France Collectivités, Fonds Français pour l'Environnement Mondial (FFEM), Etienne Bourgois, FRANCE GENOMIQUE (#ANR-10-INBS-09 to P.W.), and the Tara Ocean Foundation teams. Tara Pacific would not exist without the continuous support of the participating institutes. This research is further supported by NSF OCE #2145472 to AMSC, NSF DOB Grant 2025457 to RLVT, and with additional support from NSF PRFB #1907184 to KSIB.

## Author contributions

A.J.V. and K.S.I.B. contributed equally to this study. K.S.I.B. identified endogenization in metagenomes with support from RVT; A.J.V. confirmed the presence in genomes, with support from AMSC. All other authors (C.R.V., B.C.C.H., H.J.R., S.P., D.A., E.B., C.M., G.B., P.W., J.P., G.I., S.R., S.A., B.B., E.B., C.B., C.d.V., E.D., M.F., D.F., P.F., P.G., E.G., F.L., S.P., S.R., S.S., O.T., R.T., D.Z.) are contributing members of the Tara Pacific Expedition, collecting biosamples for metagenomics used in this study or guiding the expedition, with essential bioinformatic support in metagenomic assembly from HJR/SS, sequencing support from PW and JP, and ITS2 analysis from CRV/BCCH. A.V., K.S.I.B., A.M.S.C., and R.V.T. wrote the manuscript. This is publication number #26 of the Tara Pacific Consortium.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s42003-023-04917-9.

**Correspondence** and requests for materials should be addressed to Kalia S. I. Bistolas.

**Peer review information** *Communications Biology* thanks Guan-Zhu Han, Raúl González-Pech and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Primary Handling Editor: Zhijuan Qiu and Gene Chong.

**Reprints and permission information** is available at http://www.nature.com/reprints

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_10.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2023

<footnote n="1">¹BioSciences Department, Rice University, Houston, TX, USA. ²Microbiology Department, Oregon State University, Corvallis, OR, USA.</footnote>

<footnote n="3">³Department of Biology, University of Konstanz, Konstanz, Germany. ⁴Department of Biology, Institute of Microbiology and Swiss Institute of Bioinformatics, Vladimir-Prelog-Weg 4, ETH Zürich, CH-8093 Zürich, Switzerland. ⁵PSL Research University: EPHE-UPVD-CNRS, USR 3278 CRIOBE, Laboratoire d'Excellence CORAIL, Université de Perpignan, 52 Avenue Paul Alduy, 66860 Perpignan, Cedex, France. ⁶Centre Scientifique de Monaco, 8 Quai Antoine Ier, Monaco MC-98000, Principality of Monaco. ⁷Génomique Métabolique, Genoscope, Institut François Jacob, CEA, CNRS, Univ Evry, Université Paris-Saclay, Evry, France. ⁸Research Federation for the study of Global Ocean Systems Ecology and Evolution, FR2022/ Tara Oceans-GOSEE, 3 rue Michel-Ange, 75016, Paris, France. ⁹Fondation Tara Océan, Base Tara, 8 rue de Prague, 75012 Paris, France.</footnote>

<footnote n="10">¹⁰School of Marine Sciences, University of Maine, Orono, ME, USA. ¹¹Sorbonne Université, CNRS, Station Biologique de Roscoff, AD2M, UMR 7144 ECOMAP, Roscoff, France. ¹²Shimoda Marine Research Center, University of Tsukuba, 5-10-1 Shimoda, Shizuoka, Japan. ¹³Institut de Biologie de l'Ecole Normale Supérieure (IBENS), Ecole normale supérieure, CNRS, INSERM, Université PSL, 75005 Paris, France. ¹⁴Laboratoire des Sciences du Climat et de l'Environnement, LSCE/IPSL, CEA-CNRS-UVSQ, Université Paris-Saclay, F-91191 Gif-sur-Yvette, France. ¹⁵Weizmann Institute of Science, Department of Earth and Planetary Sciences, 76100 Rehovot, Israel. ¹⁶Université Côte d'Azur, CNRS, INSERM, IRCAN, Medical School, Nice, France. ¹⁷Laboratoire International Associé Université Côte d'Azur-Centre Scientifique de Monaco, LIA ROPSE, Monaco, France. ¹⁸Sorbonne Université, CNRS, Laboratoire d'Ecogéochimie des Environnements Benthiques (LECOB), Observatoire Océanologique de Banyuls, 66650 Banyuls sur Mer, France. ¹⁹Department of Medical Genetics, CHU of Nice, Nice, France. ²⁰Sorbonne Université, Institut de la Mer de Villefranche sur Mer, Laboratoire d'Océanographie de Villefranche, F-06230 Villefranche-sur-Mer, France.</footnote>

<footnote n="21">²¹European Molecular Biology Laboratory, European Bioinformatics Institute, Wellcome Genome Campus, Hinxton, Cambridge CB10 1SD, UK.</footnote>

<footnote n="22">²²School of Biological and Chemical Sciences, Ryan Institute, University of Galway, University Road H91 TK33, Galway, Ireland. ²³These authors contributed equally: Alex J. Veglia, Kalia S.I. Bistolas. ²⁵email: ksb97@cornell.edu</footnote>