Single T. nucleata cells were collected from Station ALOHA (A Long-Term Oligotrophic Habitat Assessment; 2° 45′N, 158° 00′W) located 100 km north of Oahu, Hawaii aboard the R/V Kilo Moana. Plankton cells were gently collected by conducting drift tows (ship not underway) using a 1 mm mesh plankton net (Sea-Gear Corporation, Melbourne, FL, United States) from the side of the vessel for a duration of 10 min. Plankton cells were sorted manually from the plankton net tow using a dissecting microscope (Leica Wild M10 Stereo Dissecting Microscope). Solitary radiolarians were individually transferred into a multi-welled culture dish filled with 0.22 um-filtered seawater, incubated for 12 h, during which time they were moved twice to newly filtered seawater. This procedure allowed hosts to clear themselves of any prey and other organisms that might cause contamination in downstream analysis. All specimens subsequently used for sequencing recovered fully in the lab and exhibited large numbers of photosymbionts present in their pseudopodial networks (Figure 1).

After undergoing a thorough wash with filtered seawater, four individual T. nucleata cells were transferred into a single sterile 2.0 ml cyrovial tube. Approximately 1.5 ml of RLT buffer (Qiagen, #79216, add 1% β-Mercaptoethanol) was added to the cyrovial, along with 0.2 μl of sterile beads (BioSpec, #1107915) and flash frozen in liquid nitrogen. Cryovials were stored at −80°C before processing. Further homogenization was facilitated through mechanical bead beating for 5 min and 20 shakes/min (TissueLyser II, Qiagen, #85300). Immediately following cell lysis, the lysate was transferred to a Qiagen RNeasy Micro Column and total RNA was isolated from the sample using Qiagen RNeasy Micro Kit following the manufacturer’s protocol. Total RNA integrity was assessed using a Fragment Analyzer (Advanced Analytical Technology). Total RNA had a RNA Quality Number (RQN) of 6.3 and quantified with Qubit 3.0 RNA High Sensitivity Assay. Total RNA was stored at −80°C until further processing.

cDNA synthesis and amplification were performed following the manufacturer’s protocol for the SMART-Seq V4 Ultra Low Input RNA kit (Clontech) with 18 cycles in the primary PCR. Most of the rRNA was removed during the process via poly-A tail selection. The resulting cDNA library was quality checked with Agilent 2100 Bioanalyzer using a High Sensitivity DNA chip and quantified with Qubit 3.0 DNA High Sensitivity Assay. The cDNA library was multiplexed and prepared for Illumina sequencing using Nextera XT Index Kit (Ilumina, # FC-131-1001) followed by a bead clean up (AMPure, Beckman Coulter #A63881). The sequencing was performed on a single lane of an Illumina HiSeq-2500 sequencer (UPC Genome Core, USC, Los Angeles, CA, United States). A total of 58.3 million paired-end reads were generated. The original data can be accessed through the Sequence Read Archive using the accession number SRP154215.

All reads were quality trimmed and filtered; Illumina adapter sequences were removed using Trimmomatic v. 0.32 (Bolger et al., 2014). Reads that could be mapped onto the PhiX genome were also removed (artifact from sequencing). All remaining reads were assembled using Trinity v. 2.1.1 (Grabherr et al., 2011) with minimum transcript length set as 250 bp. Reads were mapped back to the assembly and read counts were calculated using RSEM v. 1.2.29 (Li and Dewey, 2011) with bowtie2 (Langmead and Salzberg, 2012). Transcripts with less than 5 aligned reads were discarded. Ribosomal RNA transcripts were identified by BLASTN against the SILVA database release 123. Chloroplast and mitochondrial transcripts were identified by BLASTX and BLASTN against Heterocapsa triquetra, Breviolum minutum (formerly known as Symbiodinium minutum) and Scrippsiella spp. chloroplast and mitochondrial DNA and protein sequences found in Genbank (see Supplementary Material for sequences). These sequences were set aside and not analyzed (see Supplementary Table 1 for summary statistics of different types of transcripts).

The remaining sequences, presumably mRNA transcripts, were searched against a custom database containing all Marine Microbial Eukaryote Transcriptome Sequencing Project (MMETSP) transcriptomes and additional aquatic protistan genomes and transcriptomes (details described in Hu et al., 2018) using DIAMOND v. 0.8.28 (Buchfink et al., 2015). For each transcript, all hits with a bitscore >90% of the highest bitscore were kept. If all those hits were from dinoflagellate sources, the transcript was considered a dinoflagellate transcript. The same procedure was performed to obtain radiolarian transcripts. The GC content of the resulting transcripts were plotted to show a clear difference between dinoflagellate and radiolarian transcripts.

A Gaussian Mixture Model with two components was fitted to the distribution of GC contents of all mRNA transcripts. All transcripts were then predicted to belong to one of the two components based on the model. Only predictions with >99% confidence were accepted, the rest were deemed uncertain and were not analyzed further.

A chimera identification and removal procedure was carried out for all dinoflagellate transcripts using methods and scripts described in Yang and Smith (2013) and using 30 Peridiniales transcriptomes found in MMETSP as reference. This procedure removed 468 transcripts.

Dinoflagellate transcripts were then aligned to Brandtodinium nutricula RCC 3387 (MMETSP1462) transcriptome assembly (Johnson et al., 2018) with BLASTN using a 90% identity cutoff and 50% coverage cutoff for the shorter sequence. Transcripts aligned to the B. nutricula transcriptome were considered Brandtodinium sequences. These cutoffs were determined using test transcriptomes of related dinoflagellates such as Scrippsiella spp., Heterocapsa spp., and Azadinium spinosum to ensure that the false positive rate of being identified as Brandtodinium was low (around 1%) while maximizing the number of transcripts identified as Brandtodinium (see Supplementary Table 2). Sequences not aligned to B. nutricula most likely belonged to a parasitic dinoflagellate, as described below, but could also contain some Brandtodinium transcripts since recruitment of Brandtodinium sequences using B. nutricula transcriptome was unlikely to be complete.

Coding and protein sequences were predicted from all mRNA transcripts using TransDecoder v. 2.0 (Haas et al., 2013). Functional annotation of the predicted proteins were obtained by HMM searches against Pfam and TIGRFAM databases, and by searches against KEGG database using KAAS server (Moriya et al., 2007).

To validate the identity of the host and symbionts, nearly full length 18S rRNA sequences of the three organisms were PCR amplified from the extracted holobiont cDNA using primers 5′-GTACAAAGGACAGGGACGCA-3′ and 5′-TCCTGCCAGTAGTCATACGC-3′ for T. nucleata, 5′-ACACGGCAAAACTGCGAATG-3′ and 5′-GCGGCAGCTTTCAGGAACT-3′ for Brandtodinium, and 5′-ATACGGCGAAACTGCGAATG-3′ and 5′-TCCGCAGAAAAACTGGGTAA-3′ for the putative parasite. Primers were designed based on transcriptome assemblies. Amplicons were sequenced using Sanger sequencing. Sequences were aligned with other relevant dinoflagellate and radiolarian sequences that were also close to full length using MUSCLE (Edgar, 2004). The alignment was manually curated to remove poorly aligned positions defined as positions where more than half of the sequences were gaps. A phylogenetic tree was generated using the maximum likelihood method with default settings using MEGA7 (Kumar et al., 2016).

Original reads of the B. nutricula RCC3387 (MMETSP1462) transcriptome were downloaded from iMicrobe¹. Its transcriptome assembly were downloaded from Johnson et al. (2018). The same pipeline of read alignment, read count generation, rRNA/chloroplast/mitochondria sequence identification, protein prediction and annotation described above was also applied to the MMETSP1462 transcriptome. To compare gene expression of the two different Brandtodinium species/strains, homologous gene clusters were first generated by using the MCL algorithm (Enright et al., 2002). Both transcriptome assemblies were combined, an all-vs.-all BLASTN was conducted with 90% identity cutoff and 50% coverage cutoff of the shorter sequence. Alignment identity multiplied by coverage of the shorter sequence was used as edge weights and an inflation parameter of 1.5 was used (Fischer et al., 2011). Only clusters with transcripts from both organisms were extracted and analyzed. Read counts of transcripts from the same organism and cluster were summed and used as input for differential gene expression analysis using edgeR v. 3.6 (Robinson et al., 2010) with a BCV of 0.4 (Chen et al., 2014). Gene clusters that were significantly (FDR < 0.05) differentially expressed between the two organisms were extracted.

Three 18S rRNA sequences with significant read counts were found in the de novo transcriptome assembly. In addition to the two expected sequences, which were those of T. nucleata and of Brandtodinium sp., the 18S rRNA sequence of a second dinoflagellate species was also found. The sequence of this second dinoflagellate is almost identical to those of the non-photosynthetic, potentially parasitic dinoflagellate previously found in T. nucleata (Gast, 2006). Additionally, the closest (98–99% similar) known relatives of this putative parasite were those found in two nassellarian species, Androcyclas gamphonyca and Ceratospyris hyperborea (Dolven et al., 2007). Phylogenetic analysis revealed that these putative parasites were not related to Syndiniales, which are widespread and well-known protistan parasites (Guillou et al., 2008), but are more closely related to other Peridiniales (Figure 2A). The 18S rRNA sequence of the photosynthetic symbiont was almost identical to those of B. nutricula (formerly known as Scrippsiella nutricula) found in radiolarians (Gast and Caron, 1996; Figure 2A). Sequences of many Brandtodinium species from the Probert et al. (2014) study were shorter and therefore not included in our phylogenetic tree, but they were also more than 99% similar to our photosymbiont sequence (Figure 2A). The 18S rRNA sequences of the two dinoflagellates within T. nucleata were ∼97% similar to each other.

The 18S rRNA sequence of T. nucleata was >99% similar to that of T. nucleata described in Zettler et al. (1998). We compared these two sequences with reference sequences of the three main groups of polycystine radiolarians according to Biard et al. (2015). T. nucleata did not appeared to be closely related to any of the reference sequences or any of the Thalassicolla sequences generated by Biard et al. (2015). Out of the three main groups, T. nucleata was more similar to Collosphaeridae than the other two (Figure 2B).

De novo assembly of the holobiont transcriptome generated 435,985 transcripts, totaling 310,563,069 bp, with a N50 of 996 bp. After filtering out transcripts with low read support, and separating transcripts of rRNA, chloroplast and mitochondrial origin, 262,445 transcripts remained. The GC content of filtered transcripts clearly exhibited a bimodial distribution (Figure 3A). The likely taxonomy of the origin of these transcripts was obtained by database searches. Transcripts that were most similar to dinoflagellates in the database had much higher GC% (66.0%) than those most similar to radiolarians (41.8%), and their combined GC% distributions were similar to that of the whole transcriptome (Figure 3B). Therefore, transcripts were statistically separated into three bins: the high GC% dinoflagellate bin, which consisted of 99,580 sequences, and the low GC% T. nucleata bin, which consisted of 140,628 sequences, and a third bin whose source could not be confidently assigned (Table 1).

Brandtodinium transcripts were further separated from the dinoflagellate bin by aligning sequences to the transcriptome assembly of B. nutricula RCC 3387 in its free-living state. Sequences (31,079) that were highly similar (>90%) to those of B. nutricula RCC 3387 were considered Brandtodinium transcripts. The recruitment of sequences to the Brandtodinium reference is thought to be an underestimate as the reference transcriptome originates from a separate strain cultivated in a different living condition. Therefore, the sequences that were not recruited were considered a mixture of transcripts from Brandtodinium and the putative parasite. A number of other bioinformatic approaches including expression level-based and tetranucleotide frequency-based methods were attempted to separate the transcripts of the two dinoflagellates, but were unsuccessful because of the nature of transcriptome assembly (shorter sequences) and the similarity between the two organisms. A summary of the statistics of the different transcript bins is listed in Table 1.

Relative gene expression levels of Brandtodinium sp. were compared with those of B. nutricula in the free-living state. Significant differences were observed in several pathways and functions. In all, 19,987 homologous gene clusters were compared between the two transcriptomes, and 8,598 of them had significantly different expression levels between the two (FDR < 0.05). Among them, 5,091 gene clusters had lower expression levels in hospite, while 3,507 had higher expression levels in hospite.

Almost all genes encoding eukaryotic ribosomal proteins were expressed at a lower level in hospite compared to the free-living state (Figure 4A). Other pathways related to RNA and protein production such as spliceosome, and aminoacyl-tRNA biosynthesis were also mostly down-regulated in hospite (Figure 4B,C). In contrast to eukaryotic ribosomal genes, genes encoding chloroplast and mitochondrial ribosomes had similar expression levels between the two lifestyles (Figure 4D). Among photosynthesis genes, more gene clusters were up-regulated (13) than were down-regulated (6) in hospite (Figure 4E). The Calvin cycle, which is responsible for photosynthetic carbon fixation, shares a lot of genes with other carbon metabolism pathways such as glycolysis and pentose phosphate pathway. Most Calvin cycle genes were expressed at lower levels in hospite, however, the most important and most highly expressed gene in the pathway, rbcL encoding RuBisCO, was unchanged between the two living states (labeled in Figure 4F). Most of the genes involved in carbon metabolism were expressed at a lower level in hospite. Specifically, most glycolysis/gluconeogenesis genes, including the fructose-1,6-bisphosphatase (fbp) which is the specific gene for gluconeogenesis, were down-regulated (Figure 4G). The same was true for most TCA cycle genes including citrate synthase (CS, Figure 4H). The majority of fatty acid biosynthesis genes exhibited the same pattern (Figure 4I).

Genes involved in amino acid biosynthesis were mostly down-regulated in hospite (Figure 4J), but a small number of them displayed increased expression levels, including the three genes responsible for serine biosynthesis, serABC, and glyA; the latter is responsible for glycine biosynthesis (Figure 4J, 5). Genes involved in nitrogen metabolism had varied expression patterns which appeared to be related to the different substrate with which they interact. Nitrate transporter, nitrate reductase, and nitrite reductase (NR) were all down-regulated in hospite. In contrast, genes responsible for ammonium uptake, namely glutatmine synthetase and glutamate synthase (GS and GOGAT), were up-regulated (Figure 4K). Among genes involved in glycan biosynthesis, for example, N-glycan biosynthesis (Figure 4L), more genes were up-regulated than were down-regulated. A similar trend was observed for O-glycan and glycan precursor biosynthesis (data not shown).

Thalassicolla nucleata transcripts accounted for more than half of the filtered assembled transcriptome (Table 1), yet more than half of the transcripts did not have any functional annotation. Genes encoding structural proteins such as tubulin were among the most highly expressed genes. 56 transcripts encoding c-type lectins were found in the T. nucleata transcriptome and one of them was among the 100 most highly expressed transcripts. Several nitrogen metabolism genes including two glutamate dehydrogenases, two ammonium transporters, a glutamine synthetase, and an amino acid transporter were among the 500 most highly expressed transcripts (see Supplementary Table 3). Nitrate reduction genes were not found in the T. nucleata transcriptome.

The presence of photosymbionts and parasites in T. nucleata has been documented multiple times. Indeed, our identification of the photosymbiont, a Brandtodinium sp., was consistent with previous studies (Gast and Caron, 1996, 2001; Probert et al., 2014), as the rRNA gene sequence of the photosymbiont from the host in this study was nearly identical to sequences obtained in those studies. However, most observations of parasites in T. nucleata and other radiolarians have been most closely related to Syndiniales species, with only one exception. Gast (2006) reported a colorless dinoflagellate being released from disintegrated central capsules of T. nucleata. The putative parasite was not related to Syndiniales, but rather was closer to Peridiniales. That observation was confirmed in our study. The 18S rRNA sequence of this putative parasite was most closely related to two other dinoflagellates found in two different polycystine radiolarians (Dolven et al., 2007). This suggests the existence of a genus within Peridiniales that, similar to some Syndiniales species, specializes as parasites of radiolarians (Figure 2). It is not clear how widespread these parasites are, however, the study of Gast and Caron (1996) was conducted on specimens collected from the North Atlantic while specimens in the present study were collected in the North Pacific, implying a very broad geographic distribution. In an 18S rRNA survey of 30 different polycystine radiolarians, no parasite 18S sequences were reported (Probert et al., 2014), even though the primers used in the study should have amplified 18S rRNA of the parasites. One explanation for this is that the parasites were not detected simply because their particular hosts were not sampled, suggesting the relationships between the parasites and hosts are more specific than those between the photosymbionts and hosts. Nevertheless, our results showed that a Peridiniales parasite is a third player in the radiolarian-dinoflagellate association, at least in some species, and a potentially important interaction that should be considered in future studies.

Studying the transcriptome of the radiolarian holobiont is a challenge because these associations are difficult to access with their exclusively open ocean distributions. Additionally, specimens are often a complex mixture of organisms (e.g., captured prey in their pseudopodial networks) (Balzano et al., 2015). As a result, we presently lack reference genomes/transcriptomes, especially of the radiolarians. In our study, there was a very large difference of GC% (>24%) between the radiolarian and dinoflagellate sequences. That difference allowed us to separate radiolarian from dinoflagellate sequences with confidence. Depending on the GC% of the host, this approach may be suitable for studying the transcriptomes of other radiolarian holobionts; while the GC% of Brandtodinium spp. in various polycystine hosts is expected to be similar to the 66% observed here, the GC% of radiolarians varies widely from ∼45 to ∼55% in different species (Balzano et al., 2015). Recently, a k-mer based similarity method showed promise in separating transcripts of host and photosymbiont in Collozoum sp. holobiont (Meng et al., 2018). Such approaches may be good alternatives when there is not enough separation in GC%.

Our tranascriptomic analysis was complicated by the presence of a putative parasite that was closely related to Brandtodinium sp. Fortunately, the B. nutricula transcriptome from MMETSP allowed us to recruit Brandtodinium sp. transcripts and compare them to those derived from the transcriptome obtained from a free-living culture. However, this comparison between the two Brandtodinium species had its own challenges. This comparison was in no way an ideal or even typical comparative transcriptome study. The strains compared, the media used, the sample collection and library preparation protocols, and the sequencing platforms in the two studies were all different. These differences could all contribute to differences in gene expression in ways that have nothing to do with symbiosis. For example, the two different strains could have different gene pools and different baseline gene expression; B. nutricula was grown in L1 media, which had far higher nitrate and phosphate concentration than surface ocean waters. In an effort to take these differences into account in our analysis, we used a very large biological coefficient of variation (BCV) when comparing gene expression. A larger BCV is suitable for comparing different individuals of the same species for statistical analyses, allowing us to account for the larger systemic noises and biases caused by those differences. Additionally, we largely drew our inferences based on patterns observed in pathways and groups of genes, not individual genes. When discussing the reasons behind an observed gene expression pattern, we always try to consider the above-mentioned difference first (see below). Nevertheless, we recognize that conclusions of gene expression patterns in this study will need to be verified in better controlled experiments in the future.

The most remarkable differences in gene expression of Brandtodinium sp. between the free-living and symbiotic states was the dramatic decrease of genes involved in RNA and protein synthesis in hospite. An overwhelming majority of ribosome genes were down-regulated in hospite (48 out of 52 genes, Figure 4A). The expression patterns of spliceosome genes and aminoacyl-tRNA biosynthesis genes were similar to that of ribosome genes in that they were also largely down-regulated (Figure 4B,C). Taken together, these results imply that the growth of Brandtodinium sp. is suppressed inside the host, a finding that is consistent with reports from other host-photosymbiont systems. In cnidaria-Symbiodiniaceae symbiosis, Symbiodiniaceae grow much slower in hospite than in culture, and it is hypothesized that a reduced growth rate might be imposed and regulated by the host (reviewed by Davy et al., 2012). However, it is important to recognize that the laboratory conditions under which the free-living B. nutricula were grown might have been much more favorable for growth and reproduction compared to conditions in the host’s cytoplasm. The growth rates of the symbionts in radiolarian hosts have not been measured, and it remains an interesting question whether intracellular existence in the host is beneficial or detrimental with respect to symbiont growth rates.

In contrast to the eukaryotic ribosomal genes, prokaryotic ribosomal genes (i.e., chloroplast and mitochondrial ribosomal genes) did not exhibit large changes in expression between the free-living and in hospite conditions (Figure 4D). This result suggested that the photosymbiont maintained similar activities in its chloroplast and mitochondrion relative to its free-living state. Consistent with this interpretation, the expression levels of photosynthesis genes were similar between the two conditions, if not slightly higher in hospite (Figure 4E). Expression of carbon fixation genes was more difficult to infer because the Calvin cycle largely overlaps with other carbon metabolism pathways. Yet the expression of the gene encoding RuBisCO, arguably the most important gene expressed in these pathways, remained unchanged between the two conditions (Figure 4F). These data suggested that photosynthetic activities were likely similar between the two conditions.

The significant discordance between the expression patterns of photosynthesis genes and cell growth-related genes suggested a different fate for fixed carbon in hospite. We examined various metabolic pathways that could act as a carbon sink for the photosymbiont. Most of the differentially expressed genes in glycolysis/gluconeogenesis pathway were down-regulated in hospite (Figure 4G), including the key gene in making sugar, fructose-1,6-bisphosphatase, and the key gene in making glycerol, triose phosphate isomerase (Figure 5). These findings imply that neither sugar nor glycerol were the likely destination of carbon fixed by the photosymbiont. Similar patterns (largely down-regulation in hospite) were observed in TCA cycle genes (Figure 4H) and fatty acid biosynthesis genes (Figure 4I), suggesting that respiration and fatty acid/lipid synthesis were also unlikely routes of increased organic carbon flux. While the expression levels of amino acid biosynthesis were lower in hospite overall (Figure 4J), genes involved in the synthesis of specific amino acids including serine and glycine were higher (Figure 5). Collectively, these data seemed to imply that certain amino acids might be produced more than others in the photosymbiont. Though, we caution this deduction was based on the observation of only a few genes, and would require more evidence in future studies.

The nature of organic carbon and nutrient exchange between hosts and their photosymbionts is one of the most central questions in understanding these associations. It has been concluded in other symbiotic interactions involving dinoflagellates that these photosymbionts released a variety of compounds including glucose, glycerol, small organic acids, amino acids, and lipids to the hosts (Wang and Douglas, 1999; Whitehead and Douglas, 2003; Davy et al., 2012). Similar processes (or some subset of them) are presumably occurring in radiolarian holobionts because these associations have exceptionally high rates of photosynthetic activity, with symbionts contributing up to ∼80% of the overall carbon budgets (Caron et al., 1995). However, as noted above, patterns of gene expression in the present study did not support an important role for sugar, glycerol or fatty acid biosynthesis as major pathways for production and translocation of organic carbon to the host in hospite. There was only limited support for increased production of specific amino acids. We therefore hypothesize that amino acids might be at least one of the predominant forms of translocated compounds to the host in T. nucleata (Figure 5).

If amino acids constitute the major form of translocated carbon from photosymbiont to host in the radiolarian, then a mechanism for transporting nitrogen back into the photosymbiont should also exist. In agreement with this expectation, we observed dramatic changes in the expression of nitrogen metabolism genes of the photosymbiont in hospite. Genes related to nitrate reduction were all down-regulated, but genes involved in ammonium uptake were all up-regulated (Figure 4K, 5). This was at least partially due to the fact that B. nutricula was grown in L1 media with nitrate as the only nitrogen source. Nevertheless, our data indicate that ammonium was an important part of the nitrogen source for Brandtodinium sp. during symbiosis. Ammonium is the most likely form of nitrogen produced from prey digestion by the host. High expression levels of several nitrogen metabolism genes in the host seemed to corroborate this speculation. Oligotrophic surface waters are often nitrogen-limited, and nitrate reduction is expensive. Therefore, ammonium provided by the host is likely a vital source of nitrogen for the intracellular photosymbionts, and a mechanism by which the host modulates the metabolism of its photosymbionts.

Another likely important factor controlling intracellular gene expression and metabolism of the photosymbiont is the high CO₂/HCO₃⁻ concentration that may be present in the cytoplasm, promoting more efficient photosynthesis by the symbiont. We observed decreases in the expression of carbonic anhydrase genes in the symbiont in hospite (Figure 5). This finding implies a lessened requirement for a carbon-concentrating mechanism (CCM) during symbiosis, presumably due to the close proximity of CO₂ produced by host respiration. In some other symbiotic systems such as the sea anemone-Symbiodiniaceae symbiosis, there is evidence that the host also possesses CCMs (Weis and Reynolds, 1999; Davy et al., 2012). We found more than 30 carbonic anhydrase transcripts in the host transcriptome, although none of them were highly expressed. It is therefore unclear if T. nucleata employs CCMs to influence photosymbiont photosynthesis.

How hosts and symbionts recognize each other, or from a different perspective, how symbionts evade host digestion, is another fundamental question for heterotroph-alga symbiotic associations. Previous studies have concluded that interactions between c-type lectins and glycans play an important role in host-symbiont recognition in cnidarian-Symbiodiniaceae symbiosis (Wood-Charlson et al., 2006; Davy et al., 2012). Many c-type lectin genes have been identified in metazoan (Wood-Charlson and Weis, 2009) and radiolarian hosts (Balzano et al., 2015). We identified 56 c-type lectin transcripts in T. nucleata, with at least one of them being expressed at very high level. We also observed genes involved in glycan biosynthesis in the symbiont expressed at higher levels in hospite (Figure 4L). These findings are consistent with the possibility that recognition between the radiolarian host, T. nucleata, and its photosymbiont, Brandtodinium sp., might also involve lectin/glycan interactions.