Original Research ARTICLE
Bacteria in Ostreococcus tauri cultures – friends, foes or hitchhikers?
- 1Institut Pasteur, Microbial Evolutionary Genomics, Paris, France
- 2CNRS, UMR 3525, Paris, France
- 3CNRS, UMR 7232, Biologie Intégrative des Organismes Marins, Observatoire Océanologique, Banyuls-sur-Mer, France
- 4Sorbonne Universités, UPMC Université Paris 06, UMR 7232, BIOM, Observatoire Océanologique, Banyuls-sur-Mer, France
Marine phytoplankton produce half of the oxygen we breathe and their astounding diversity is just starting to be unraveled. Many microbial phytoplankton are thought to be phototrophic, depending solely on inorganic sources of carbon and minerals for growth rather than preying on other planktonic cells. However, there is increasing evidence that symbiotic associations, to a large extent with bacteria, are required for vitamin or nutrient uptake for many eukaryotic microalgae. Here, we use in silico approaches to look for putative symbiotic interactions by analysing the gene content of microbial communities associated with 13 different Ostreococcus tauri (Chlorophyta, Mamilleophyceae) cultures sampled from the Mediterranean Sea. While we find evidence for bacteria in all cultures, there is no ubiquitous bacterial group, and the most prevalent group, Flavobacteria, is present in 10 out of 13 cultures. Among seven of the microbiomes, we detected genes predicted to encode type 3 secretion systems (T3SS, in 6/7 microbiomes) and/or putative type 6 secretion systems (T6SS, in 4/7 microbiomes). Phylogenetic analyses show that the corresponding genes are closely related to genes of systems identified in bacterial-plant interactions, suggesting that these T3SS might be involved in cell-to-cell interactions with O. tauri.
Eukaryotes acquired photosynthesis on multiple occasions from endosymbiosis (Archibald, 2012; De Clerck et al., 2012), resulting in an astounding phylogenetic diversity of phytoplanktonic eukaryotes (Not et al., 2012 for a review). The coexistence of so many species competing for the same resources does not fit theoretical prediction that in a stable environment, the best competitor wins, a puzzle coined as the “paradox of the plankton” (Hutchinson, 1961). However spatial and temporal environmental heterogeneities affect the environmental stability hypothesis (e.g., Levins and Culver, 1971), and interactions among competitors (e.g., Gross, 2008), as well as non-competitive interactions with other species (e.g., mutualism, commensalism) may increase their probability of coexistence. The “phycosphere” – the region immediately surrounding and influenced by phytoplankton cells – is an important bacterial habitat (Bell and Mitchell, 1972; Blackburn et al., 1998). Heterotrophic bacterial communities are sustained by phytoplankton exudates and play an important role in remineralization of nitrogen (N) and phosphate (P). Many protist algae are mixotrophic, gaining nutrients either by photosynthesis or by heterotrophy depending on environmental conditions (Flynn et al., 2013 for a review), adding further complexity to plankton community assemblages. Facilitative interactions have often been suggested by phycologists whose algal laboratory cultures were often most successful when they did not eliminate all bacteria from the cultures (Cole, 1982 for a review), and the generality of bacterial-seaweed associations is now well known (Goecke et al., 2010; Egan et al., 2013; Hollants et al., 2013 for reviews). For example, some protists are known to fix nitrogen or carbon via their cyanobacterial endosymbionts (reviewed in Nowack and Melkonian, 2010; Thompson and Zehr, 2013) and some chemical pathways mediating algal-bacterial interactions were identified (e.g., Seyedsayamdost et al., 2011; Patzelt et al., 2013; Syrpas et al., 2014). Whether some phytoplanktonic eukaryotes evolved specific interactions with bacteria beyond ecological facilitation is still a matter of debate, and little information is available on close algal-bacterial interactions for unicellular eukaryotes. The study of interactions between planktonic microbes has been long hampered by our lack of knowledge of these unicellular organisms, especially for the smallest sized planktonic eukaryotes, the picoeukaryotes (cell diameter size <2 μm). These microorganisms often lack morphologically informative characters and are difficult to isolate and maintain in culture. For example, species of the genus Ostreococcus are hard to discriminate, and it is thus difficult to study species-specific interactions between partners one cannot identify (Subirana et al., 2013).
Early observations of physical attachment between some diatom species and bacteria date back to the first cytological observations (e.g., the diatom Skeletonema costatum in Droop and Elson, 1966). A pioneering study of bacteria-phytoplankton interactions screened the bacterial content of microalgal cultures by standard microbiological techniques (Berland and Maestrini, 1969). With the development of molecular biology tools, ribosomal RNA genes sequencing and barcoding approaches allowed establishing links between phytoplankton and bacterial community dynamics in natural communities (Fukami et al., 1992; Rooney-Varga et al., 2005) and culture collections of diatoms and dinoflagellates (Schäfer et al., 2002; Jasti et al., 2005; Sapp et al., 2007). Recently, integrative approaches associating sequence-based identification of cells, cytometry, cell sorting and tracer experiments with 15N and 13C, demonstrated a mutualistic interaction between N2-fixing cyanobacteria and a phytoplanktonic prymnesiophyte (Thompson et al., 2012). Cell sorting and single cell sequencing of natural isolates now enables to identify candidate bacteria–protist partnerships without cultivation, and these techniques promise to uncover overlooked interactions (Martinez-Garcia et al., 2012).
The development of next generation sequencing now enables to sequence microalgae and associated microbiomes and to investigate the molecular toolkit of putative interactions from gene content analyses. For example, symbiosis involving nitrogen fixation requires the presence of the nitrogenase operon in the bacterial partner, and complementation for vitamin biosynthesis requires the presence of the genes responsible for the vitamin B pathway in the bacteria along with suitable algal uptake mechanisms (Table 1). In some plant-bacterium interactions, physical cell-to cell-interactions occur via specialized protein secretion systems such as the type 3 secretion systems (T3SS; Kosarewicz et al., 2012), or the T4SS (Smillie et al., 2010). These systems are sophisticated molecular needles that enable the translocation of bacterial effectors from the bacterial cytoplasm to the eukaryotic cell (Figure 1). These systems are involved in both antagonistic and beneficial interactions between bacteria and eukaryotes.
TABLE 1. The Molecular toolkit of bacteria – eukaryote interactions: examples of bacterial target genes for screening genomes and metagenomes.
FIGURE 1. Assembled structure of the type 3 secretion system (modified from Abby and Rocha, 2012). This system is made of an extracellular part for close contact with host cell, and a membrane-associated part containing the core secretion apparatus. The components with bright colors correspond to the most conserved components, and are thus the one searched by sequence similarity to assess the system’s presence in the assembled sequences.
While picoalgae of the genus Ostreococcus are distributed worldwide (Piganeau and Moreau, 2007; Demir-Hilton et al., 2011), Ostreococcus tauri, first isolated from a French Mediterranean coastal lagoon and described as the smallest eukaryotic species known (Courties et al., 1994; Chrétiennot-Dinet et al., 1995) has been found mainly in coastal regions and lagoons (Subirana et al., 2013). Within a 3-year period (2006–2008) of sustained effort, we could isolate 17 new wild-type clonal lines of O. tauri and characterize them with a few genetic markers (Grimsley et al., 2010). Interestingly, none of the O. tauri cultures we retained were completely axenic despite initial treatment with antibiotics. Here, we postulate that the bacterial microbiome consistently present in successful O. tauri cultures could be required for the health of these cultures. We present an analysis of the bacterial microbiome associated with these cultures and apply recently developed tools to mine the O. tauri microbiome for protein secretion systems involved in bacteria-eukaryotes interaction.
Material and Methods
We analyzed data from 13 O. tauri strains that were sampled from surface water in five locations by the North-West Mediterranean sea previously described (Blanc-Mathieu et al., 2013). Cultures were isolated by serial filtrations, addition of Keller’s salts as a supplement (Keller et al., 1987), and growth in a culture chamber in the laboratory. These strains were established from clonal culture by plating out in gelled medium and re-inoculating cells picked from a single colony in liquid medium to obtain cell densities above 107 ml-1. One strain was the control O. tauri laboratory strain RCC4221, cloned from the RCC745 culture (Courties et al., 1994), and the other 12 were isolated more recently as described in Grimsley et al. (2010). Despite the treatment by antibiotics as described in Grimsley et al. (2010); kanamycin 20 μg/ml, penicillin 25 μg/ml, and neomycin 20 μg/ml final concentrations, none of the strains were found to be completely axenic. This is the case not only for the 13 strains analyzed here, but also for 100s of other isolations made by plating out for single algal cells, including O. mediterraneus, Micromonas sp. and Bathycoccus prasinos (Nigel Grimsley, unpublished observations). Three micrograms of total DNA was extracted from each culture as previously described (Derelle et al., 2006). Genomic DNA of the strains was randomly sheared into ∼250-bp fragments. The libraries created from these fragments were sequenced on an Illumina GAIIx and Hiseq system at the Joint Genome Institute1 Community Sequencing Program (CSP-129). Sequence data for strains RCC1108, RCC1114, RCC1115, RCC1116, RCC1558, RCC1559, and RCC4221 were 76 bp paired-end reads (depth of coverage for the algal genome 160-fold to 340-fold) and sequence data of strains RCC1110, RCC1112, RCC1117, RCC1118, RCC1123, and RCC1561 were 101 bp paired-end reads (depth of coverage for the algal genome 780-fold to 1130-fold).
Extraction of Bacterial Sequences and Taxonomic Affiliation
Paired-end data of each strain were mapped to the reference nuclear genome of O. tauri (GenBank accession numbers: CAID01000001-CAID01000020, (Blanc-Mathieu et al., under review) with the Burrow-Wheeler Aligner (BWA) with parameters n= 6 l= 35 k= 3 e= 3 (Li and Durbin, 2009). Read pairs with no read mappings to the genome sequence were identified based on SAM flags and extracted in fastq files with the seqtk package2. Paired-end reads were assembled with Velvet (Zerbino and Birney, 2008) using default parameters and various k-mer sizes. The assembly with the highest median contig length (N50) was retained.
To remove remaining Ostreococcus sequences that might have been too divergent to be mapped on to the reference sequence, contigs with more than 80% amino-acids identity over 90 bp with available Mamiellales coding sequences were discarded using the PRASINOID interface3 (Vaulot et al., 2012). The remaining contigs were analyzed via MG-RAST (Meyer et al., 2008) on the Refseq protein database. Taxonomic affiliations for each microbiome were downloaded from the MG-RAST server and analyzed using in-house scripts to retrieve all contigs with an alignment against Refseq longer than 30 amino-acids and with more than 80% sequence identity.
Estimation of Ubiquity and Abundance of Bacteria
The abundance of each bacterial group was measured as the sum of reads affiliated to one group, divided by the total number of affiliated reads for each microbiome. Ubiquity was defined as the number of occurrences of a genus in the 13 metagenomes.
To estimate the bacteria to Ostreococcus-cells ratio, we assumed that the number of reads affiliated to Ostreococcus, rO is equal to the product of the number of Ostreococcus cells in the sample, CO, by the genome size, GO, by a constant αO (representing the product of the DNA extraction efficiency and the sequencing efficiency):
The number of reads affiliated to bacteria, RB, is equal to the sum of the product of CBi the number of bacterial cells from each bacterial strain i, by GBi their average genome size, by a constant αBi. Assuming an equal α and genome size between bacteria (marine bacteria have a genome in the 2–7 Mb size range), the complete number of bacterial reads can be expressed as the product of CB, GB, and αB.
From these two equations we can estimate the ratio of bacterial to Ostreococcus cells as a function of two parameters: the relative genome size differences, the relative number of reads and the relative α parameter between bacteria and Ostreococcus DNA.
Statistical analyses and ubiquity abundance plots were done with R4.
Search for Nitrogenases
To screen for the presence of nitrogenases, we used the amino-acid sequences of the different types of nitrogenases described in Raymond et al., (2004). We processed the output of the blastx of this dataset against the assemblies to retain all hits with amino-acid identity greater or equal to 60% and total alignment length higher than 100 amino-acids. Ten cultures contained hits with these criteria. Further analysis of theses hits did not confirm that the encoding genes were nitrogenases, but genes belonging to related gene families, like hydrogenases, so that we did not proceed to further analysis.
Identification and Phylogenetic Analysis of Protein Secretion Systems
Tools were recently developed to identify type 3, type 4, and type 6 secretion systems from similarity search of essential components and gene content/gene organization criteria (Abby and Rocha, 2012; Gama et al., 2012; Guglielmini et al., 2014). They were used with the MacSyFinder framework (Abby et al., 2014) to detect these systems in the assembled contigs of the microbiomes. Phylogenetically relevant components of T3SS detected in the microbiome were analyzed (i.e., SctJ, SctN, SctQ, SctR, SctS, SctT, SctU, and SctV; see Figure 1). Their sequences were introduced into the appropriate pre-existing gene families (dataset of Abby and Rocha, 2012). We aligned their sequences with Muscle (default parameters) and selected informative sites with BMGE (BLOSUM30 similarity matrix, gap rate cut-off = 0.20, sliding window size = 3, entropy score cut-off = 0.5; Edgar, 2004; Criscuolo and Gribaldo, 2010). Then, these alignments were concatenated, and a phylogenetic tree including microbiome sequences was built with RAxML (Le and Gascuel matrix + 4-categories-discretized Gamma distribution for rate variation among sites + empirical frequencies of amino-acids) with 100 rapid bootstraps (Stamatakis, 2006; Le and Gascuel, 2008).
Which Bacteria Live in O. tauri Cultures?
MG-RAST taxonomic affiliation of contigs based on sequence identity to the Refseq protein database listed 149 distinct bacterial taxonomic affiliations at the genus level. A few contigs were assigned to higher taxonomic ranks like alphaproteobacteria or gammaproteobacteria. We estimated the average abundance of each group as the proportion of reads affiliated to a bacterial group to the total number of reads assigned to bacteria. Based on these estimates, we found no bacterial group to be both ubiquitous and abundant throughout the 13 microbiomes (Figure 2). The most ubiquitous bacterial group is Flavobacterium and is found in 10 out of 13 cultures, with an average abundance of 18% of reads in the 10 microbiomes. The most abundant group of bacteria is Limnobacter as this genus is represented by 60% of the reads in the three microbiomes where it was detected. Other abundant groups are Pseudomonas, Roseovarius, and Oceanocaulis (Figure 2), each group represent more than 10% of the reads. The most abundant bacterial groups for each microbiome are reported in Table 2. They belong to seven different genera for the 13 microbiomes; Flavobacterium, Pseudomonas, and Limnobacter dominate in the six microbiomes with more than one bacterium for 10 O. tauri cells, while Sphingomonas, Robiginitalea, Oceanocaulis, and Roseovarius are the most abundant bacterial groups in the seven microbiomes with less than one bacterium for 10 O. tauri cells (see below). These genera belong to the Bacteroidetes and Proteobacteria phyla, with three orders present for the latter: alphaproteobacteria, betaproteobacteria, and gammaproteobacteria.
FIGURE 2. Ubiquity versus abundance (average % of reads across microbiomes where the group is present) of bacterial groups after MGRAST analysis of the 13 microbiomes.
TABLE 2. Number of reads assigned to Ostreococcus tauri and bacteria for each cultured strain and estimation of the number of bacterial cells for 10 O. tauri cells from Eqn (2) for the most abundant bacteria.
How Many Bacteria Per O. tauri Cell?
We used the number of reads affiliated to O. tauri and to the most abundant bacterial groups to estimate the number of bacterial to microalgal cells using Eqn (2). This estimation relies on the assumption that the DNA extraction protocol and the sequencing are not biased towards O. tauri or bacteria. In a 45 marine bacterial genomes dataset, the average genome size was estimated to be ∼4 Mb (Moran and Armbrust, 2007). We used this value as a proxy for genome size of bacteria associated to O. tauri strains. Since marine bacterial genome size varies between 2 and 7 Mb, we do not expect more than a twofold difference in our estimate. The number of bacteria for 10 O. tauri cells varies by two orders of magnitude between cultures: from 0.02 to 4 (Table 2).
Looking for Bacterial Factors Involved in Interactions with Eukaryotes
Type 3 secretion systems have evolved and diversified into recognizable sub-types to interact with different kinds of eukaryotic cells (animal vs. plant), and participate in different types of interaction with eukaryotes (antagonistic vs. beneficial; Troisfontaines and Cornelis, 2005; Abby and Rocha, 2012). The human pathogen Salmonella uses two types of T3SS (SPI-1 and SPI-2) at different stages of host infection (Valdez et al., 2009), while nitrogen-fixing Rhizobiales use a particular T3SS (Rhizobiales type) to establish symbiotic interactions with host plants (Dai et al., 2008; Kambara et al., 2009). Therefore, the presence of T3SS components in a bacterial dataset is a hallmark of bacterium-eukaryote interactions, and the phylogenetic typing of these components can help identifying the kind of interaction, and the type of targeted eukaryotic cell. Among T4SS – classically dedicated to conjugation – some are involved in plant pathogenesis (e.g., tumor formation, Pitzschke and Hirt, 2010) but also in plant symbiosis (Hubber et al., 2004). T6SS allow the translocation of effectors from bacteria to eukaryotic cells in antagonistic relationships, but were also proved to target bacteria for bacterial competition (Pukatzki et al., 2007; Hood et al., 2010; Kapitein and Mogk, 2013). T6SS are virulence factors for several phytopathogens, and were observed in plant symbionts genomes (Amadou et al., 2008; Wu et al., 2008). We looked for signs of these putative factors of bacterium-eukaryote interaction in the microbiota associated to O. tauri. Using both gene content and close linkage distance between genes as criteria for the detection of T4SS and T6SS (Gama et al., 2012; Guglielmini et al., 2014), there was no evidence of T4SS seemingly involved in protein secretion, but we could detect six occurrences of putative T6SS in four different microbiomes (Table S1).
Using the same kind of detection methods for type 3 secretion systems (Abby and Rocha, 2012), we could detect seven occurrences of putative T3SS in contigs from six different microbiomes (Figure 3). We took advantage of a previous study to sub-type T3SS: we included the T3SS’ components detected in O. tauri’s microbiome in a reference phylogeny annotated with T3SS sub-types and corresponding host types (Figure 4 of Abby and Rocha, 2012). Interestingly, detected T3SS consistently placed within, or as sister-groups of plant-associated T3SS sub-types with high rapid bootstrap supports (Figure 4). A first sub-type placed within the “Hrp1” T3SS family, which includes many systems of plant pathogens and some involved in plant symbiosis (i.e., Pseudomonas syringae, Dickeya dadantii). This system, whose contig was assigned to the Pseudomonas genus, was found in a single microbiome and showed sequence and genetic architecture highly similar to the T3SS observed in the genome of P. brassicacearum, a root-associated plant symbiont (Ortet et al., 2011). Three T3SS occurrences seemed to correspond to a single system that placed within the clade of plant-symbionts, Rhizobiales T3SS. This system showed high similarity with a Mesorhizobium system (Rhizobiales species, 50–70% identity in a Blast analysis). The three corresponding contigs were attributed to three different alphaproteobacterial species: 2 Rhizobiales and 1 Rhodobacterales. Intriguingly, a third type of T3SS found in two microbiomes fell outside of the defined T3SS families, but constituted a sister-group of the Hrp1 family (Figure 4). The closest system found in the phylogeny was that of Herbaspirillum seropedicae, a plant symbiont, and both taxonomic attribution of contigs and similarity searches of the system using Blast and MG-RAST pointed at the species Limnobacter sp., a Burkholderiales.
FIGURE 3. Genetic architecture of microbiome T3SS. The genetic architecture of T3SS detected in microbiomes’ contigs and of the closest systems found in a complete genome is displayed along the corresponding phylogeny. The color of boxes corresponds to the colors in Figure 1.
Overall, one of the detected T6SS gene cluster was found on a contig assigned to the Pseudomonas genus, in the same microbiome where a Pseudomonas T3SS was inferred (RCC1108, see Table S1; Figure 4). Five of the six detected T6SS were found in three microbiomes showing also evidence of T3SS gene clusters.
FIGURE 4. The phylogeny of T3SS reflects host-cell type. The T3SS detected in O. tauri microbiomes were integrated in a reference phylogeny of the T3SS (Abby and Rocha, 2012). T3SS phylogeny shows a diversification of T3SS into host-cell adapted sub-types: animal associated, and plant associated (“ecology” panel, and column). While most of these sub-types are mostly involved in pathogenic interactions with plants or animals, some Hrp1 T3SS and all Rhizobiales T3SS are involved in symbiotic associations with plants.
Despite the efforts of many laboratories over the last century to define the media and growth conditions required for different marine algae, most have so far remained recalcitrant to growth in a completely defined medium, and require seawater to grow over essential extra nutrients (e.g., phosphate and nitrate). Seawater is a complex solution of chemicals and organisms that can vary considerably in its composition between geographically distant regions, complicating the development of appropriate culture media. Collections are thus often located in marine biology laboratories close to coastlines. The majority of algal cultures cannot be maintained axenically, rendering physiological analyses of their nutritional requirements more difficult. Indeed, many unicellular algae are mixotrophic, and can satisfy part of their nutritional requirements by ingesting bacteria.
Despite antibiotics treatments and the isolation of single cells from colonies in soft agarose, none of the O. tauri strains were axenic. This was the case not only for the 13 strains analyzed here, but also for 100s of other isolations made by plating out for single algal cells (unpublished observations). In addition, bacterial cultures issued from media of such cultures as well as observations of algal strains by flow cytometry almost always confirmed the presence of common seawater bacteria (unpublished data). These observations strongly suggest that either O. tauri adhere to bacterial cells during the cultivation process, or that O. tauri require some unidentified substance from bacterial cells for growth. Recent experimental evidence and genome analysis suggests that O. tauri is vitamin B12-dependent (Helliwell et al., 2011). We provided three vitamins (thiamine [B1]), biotin [H], and cyanocobalamin [B12]) as described for standard Keller’s medium (Keller et al., 1987), making unlikely the selection of bacteria for vitamin B12 production. In order to identify putative supplements required by isolated O. tauri cells for growth, future work could focus in the use of sterile artificial seawater, and step-by-step introduction of candidate substances. But considering the great deal of effort already put in attempts to define suitable culture media for this kind of algae (Keller et al., 1987), it might be more efficient to first isolate associated bacteria and investigate their influence on the physiology of the microalgae (Le Chevanton et al., 2013).
The bioinformatic analyses performed here confirmed the presence of a diverse collection of common marine bacteria in O. tauri cultures. Flavobacteria was largely found in our samples. This important class of Bacteroidetes often constitutes a significant portion of marine microbial communities and has been reported in microalgal cultures (Berland and Maestrini, 1969; Mann et al., 2013). Similarly, Bacteroidetes have also been reported in the surface waters of NW Mediterranean sea (Lami et al., 2009). Flavobacteria are found both free-living and attached to organic aggregates and are considered as major mineralizers of organic matter (Kolton et al., 2013). Interestingly, the type 3 and type 6 secretion systems have been detected in the microbiomes with higher bacterial prevalence. Microbiomes for strains that have a higher ratio of number of bacteria to O. tauri cells (Table 2), are more likely to contain a detectable T3SS (Wilcoxon signed rank test p < 0.01). The taxonomic affiliation of the contigs containing a predicted T3SS correspond to the most abundant bacterial group in the microbiome except for RCC1110 and RCC1561, whose T3SS contigs were assigned to alphaproteobacteria while the most abundant bacterial group is a Flavobacterium. The contigs containing the secretion systems were attributed to genera with evidence of species interacting with eukaryotes via protein secretion systems. Finally, no single strain of bacterium was found in all of the O. tauri cultures. Since all these bacteria were isolated from the same host species, it could seem unlikely at first that very specific physical interactions or nutritional requirements exist between O. tauri and its microbiome. However, several biases could partially explain the heterogeneity observed between the microbiomes both in terms of taxonomic diversity, and therefore gene content. Firstly as the effective detection of the secretions systems rely on sequence similarity search and the genetic organization of their components, it heavily depends on the sequencing and assembly quality. In the context of NGS approaches for metagenomics, whose short reads are difficult to assemble, it is likely that we missed occurrences of systems due to contig assembly errors and biases. Secondly, the filtration steps in O. tauri isolation is likely to sort out aggregates of bacterial cells, and as a consequence, some bacteria of interest to understand the growth of O. tauri. Finally, antibiotics treatment had also an impact on the bacterial populations we analysed in this study, and this might also explain some of the heterogeneity in terms of bacterial diversity and gene content. To specifically find preferred associated bacteria we should repeat this work without using antibiotics, even if these conditions, it may be difficult to isolate Ostreococcus strains.
In conclusion, we provide evidence of pervasive bacterial presence in O. tauri cultures, despite initial antibiotic treatment. We provide evidence for putative plant-associated T3SS in six microbiomes, and several cases of T6SS in four microbiomes (three displayed both systems). For now there are no studies showing a clear association between T6SS sub-types and their function, thus it is hard to define from genomic analyses whether the systems we detected are targeting bacteria or eukaryotes. But in both cases, they might be parts of interactions between bacteria and eukaryotes, even indirectly: via bacterial competition, T6SS were found to serve as colonization factors in the plant pathogen Agrobacterium tumefaciens, (Ma et al., 2014) and to provide plants a protection against pathogens in P. fluorescens (Decoin et al., 2014). On the other hand, the analysis of detected T3SS gave a clearer picture as it clearly shows they are involved in interaction with plant cells. The three sub-types they belong to – or the groups they are closer to, all contain systems typical of plant symbiosis, and pathogenicity in the case of the Hrp1 group. Further experimental work is required to determine the impact of these secretion systems in O. tauri growth, while keeping in mind that interactions are dynamic, and that the same bacteria may change between “friend,” “foe,” or “hitch-hiker” over time or environmental conditions (Andrade-Domínguez et al., 2014).
Conflict of Interest Statement
The Review Editor, Fabrice Not, declares that, despite having collaborated with author, Nigel Grimsley, the review process was handled objectively and no conflict of interest exists. The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We would like to thank the GenoToul bioinformatic platform (http://bioinfo.genotoul.fr/) for access to computing facilities and Romain Blanc-Mathieu, Sophie Sanchez-Ferandin, Marc Garcia-Garcerà and the Genomics of Phytoplankton group in Banyuls sur mer for stimulating discussions. We would like to acknowledge the Joint Genome Institute (http://www.jgi.doe.gov/) for sequencing. This work was funded by ANR-13-JSV6-0005 to GP, the Institut Pasteur, the French Centre National de la Recherche Scientifique and the European Research Council (grant EVOMOBILOME, number 281605).
The Supplementary Material for this article can be found online at: http://www.frontiersin.org/journal/10.3389/fmicb.2014.00505/abstract
Abby, S. S., Néron, B., Ménager, H., Touchon, M., and Rocha, E. P. C. (2014). MacSyFinder: a program to mine genomes for molecular systems with an application to CRISPR-Cas systems. PLoS ONE 9:e110726. doi: 10.1371/journal.pone.0110726
Abby, S. S., and Rocha, E. P. C. (2012). The non-flagellar type III secretion system evolved from the bacterial flagellum and diversified into host-cell adapted systems. PLoS Genet. 8:e1002983. doi: 10.1371/journal.pgen.1002983
Amadou, C., Pascal, G., Mangenot, S., Glew, M., Bontemps, C., Capela, D.,et al. (2008). Genome sequence of the beta-rhizobium Cupriavidus taiwanensis and comparative genomics of rhizobia. Genome Res. 18, 1472–1483. doi: 10.1101/gr.076448.108
Archibald, J. M. (2012). “Chapter three – the evolution of algae by secondary and tertiary endosymbiosis,” in Advances in Botanical Research, Vol. 64, ed. G. Piganeau (Waltham: Academic Press), 87–118.
Blanc-Mathieu, R., Sanchez-Ferandin, S., Eyre-Walker, A., and Piganeau, G. (2013). Organellar inheritance in the green lineage: insights from Ostreococcus tauri. Genome Biol. Evol. 5, 1503–1511. doi: 10.1093/gbe/evt106
Chrétiennot-Dinet, M. J., Courties, C., Vaquer, A., Neveux, J., Claustre, H., Lautier, J.,et al. (1995). A new marine picoeucaryote: Ostreococcus tauri gen. et sp. nov. (Chlorophyta, Prasinophyceae). Phycologia 34, 285–292. doi: 10.2216/i0031-8884-34-4-285.1
Criscuolo, A., and Gribaldo, S. (2010). BMGE (Block Mapping and Gathering with Entropy): a new software for selection of phylogenetic informative regions from multiple sequence alignments. BMC Evol. Biol. 10:210. doi: 10.1186/1471-2148-10-210
Croft, M. T., Lawrence, A. D., Raux-Deery, E., Warren, M. J., and Smith, A. G. (2005). Algae acquire vitamin B12 through a symbiotic relationship with bacteria. Nature 438, 90–93. doi: 10.1038/nature04056
Dai, W. J., Zeng, Y., Xie, Z. P., and Staehelin, C. (2008). Symbiosis-promoting and deleterious effects of NopT, a novel type 3 effector of Rhizobium sp. strain NGR234. J. Bacteriol. 190, 5101–5110. doi: 10.1128/JB.00306-08
Decoin, V., Barbey, C., Bergeau, D., Latour, X., Feuilloley, M. G., Orange, N.,et al. (2014). A type VI secretion system is involved in pseudomonas fluorescens bacterial competition. PLoS ONE 9:e89411. doi: 10.1371/journal.pone.0089411
De Clerck, O., Bogaert, K. A., and Leliaert, F. (2012). “Chapter two – diversity and evolution of Algae: primary endosymbiosis,” in Advances in Botanical Research, ed. G Piganeau (Waltham: Academic Press), 55–86.
Demir-Hilton, E., Sudek, S., Cuvelier, M. L., Gentemann, C. L., Zehr, J. P., and Worden, A. Z. (2011). Global distribution patterns of distinct clades of the photosynthetic picoeukaryote Ostreococcus. ISME J. 5, 1095–1107. doi: 10.1038/ismej.2010.209
Derelle, E., Ferraz, C., Rombauts, S., Rouzé, P., Worden, A. Z., Robbens, S.,et al. (2006). Genome analysis of the smallest free-living eukaryote Ostreococcus tauri unveils many unique features. Proc. Natl. Acad. Sci. U.S.A. 103, 11647–11652. doi: 10.1073/pnas.0604795103
Egan, S., Harder, T., Burke, C., Steinberg, P., Kjelleberg, S., and Thomas, T. (2013). The seaweed holobiont: understanding seaweed–bacteria interactions. FEMS Microbiol. Rev. 37, 462–476. doi: 10.1111/1574-6976.12011
Flynn, K. J., Stoecker, D. K., Mitra, A., Raven, J. A., Glibert, P. M., Hansen, P. J.,et al. (2013). Misuse of the phytoplankton–zooplankton dichotomy: the need to assign organisms as mixotrophs within plankton functional types. J. Plankton Res. 35, 3–11. doi: 10.1093/plankt/fbs062
Gama, J. A., Abby, S. S., Vieira-Silva, S., Dionisio, F., and Rocha, E. P. (2012). Immune subversion and quorum-sensing shape the variation in infectious dose among bacterial pathogens. PLoS Pathog. 8:e1002503. doi: 10.1371/journal.ppat.1002503
Guglielmini, J., Neron, B., Abby, S. S., Garcillan-Barcia, M. P., de la Cruz, F., and Rocha, E. P. (2014). Key components of the eight classes of type IV secretion systems involved in bacterial conjugation or protein secretion. Nucleic Acids Res. 42, 5715–5727. doi: 10.1093/nar/gku194
Helliwell, K. E., Wheeler, G. L., Leptos, K. C., Goldstein, R. E., and Smith, A. G. (2011). Insights into the evolution of vitamin B12 auxotrophy from sequenced algal genomes. Mol. Biol. Evol. 28, 2921–2933. doi: 10.1093/molbev/msr124
Hollants, J., Leliaert, F., De Clerck, O., and Willems, A. (2013). What we can learn from sushi: a review on seaweed–bacterial associations. FEMS Microbiol. Ecol. 83, 1–16. doi: 10.1111/j.1574-6941.2012.01446.x
Hood, R. D., Singh, P., Hsu, F., Güvener, T., Carl, M. A., Trinidad, R. R.,et al. (2010). A type VI secretion system of Pseudomonas aeruginosa targets a toxin to bacteria. Cell Host Microbe. 7, 25–37. doi: 10.1016/j.chom.2009.12.007
Hubber, A., Vergunst, A. C., Sullivan, J. T., Hooykaas, P. J., and Ronson, C. W. (2004). Symbiotic phenotypes and translocated effector proteins of the Mesorhizobium loti strain R7A VirB/D4 type IV secretion system. Mol. Microbiol. 54, 561–574. doi: 10.1111/j.1365-2958.2004.04292.x
Jasti, S., Sieracki, M. E., Poulton, N. J., Giewat, M. W., and Rooney-Varga, J. N. (2005). Phylogenetic diversity and specificity of bacteria closely associated with Alexandrium spp. and other phytoplankton. Appl. Environ. Microbiol. 71, 3483–3494. doi: 10.1128/AEM.71.7.3483-3494.2005
Kambara, K., Ardissone, S., Kobayashi, H., Saad, M. M., Schumpp, O., Broughton, W. J.,et al. (2009). Rhizobia utilize pathogen-like effector proteins during symbiosis. Mol. Microbiol. 71, 92–106. doi: 10.1111/j.1365-2958.2008.06507.x
Kapitein, N., and Mogk, A. (2013). Deadly syringes: type VI secretion system activities in pathogenicity and interbacterial competition. Curr. Opin. Microbiol. 16, 52–58. doi: 10.1016/j.mib.2012.11.009
Kolton, M., Sela, N., Elad, Y., and Cytryn, E. (2013). Comparative genomic analysis indicates that niche adaptation of terrestrial flavobacteria is strongly linked to plant glycan metabolism. PLoS ONE 8:e76704. doi: 10.1371/journal.pone.0076704
Lami, R., Ghiglione, J., Desdevises, Y., West, N. J., and Lebaron, P. (2009). Annual patterns of presence and activity of marine bacteria monitored by 16S rDNA–16S rRNA fingerprints in the coastal NW Mediterranean Sea. Aquat. Microb. Ecol. 54, 199–210. doi: 10.3354/ame01264
Le Chevanton, M., Garnier, M., Bougaran, G., Schreiber, N., Lukomska, E., Bérard, J.,et al. (2013). Screening and selection of growth-promoting bacteria for Dunaliella cultures. Algal Res. 2, 212–222. doi: 10.1016/j.algal.2013.05.003
Ma, L. S., Hachani, A., Lin, J. S., Filloux, A., and Lai, E. M. (2014). Agrobacterium tumefaciens deploys a superfamily of type VI secretion DNase effectors as weapons for interbacterial competition in planta. Cell Host Microbe 16, 94–104. doi: 10.1016/j.chom.2014.06.002
Mann, A. J., Hahnke, R. L., Huang, S., Werner, J., Xing, P., Barbeyron, T.,et al. (2013). The Genome of the alga-associated marine flavobacterium Formosa agariphila KMM 3901(T) reveals a broad potential for degradation of algal polysaccharides. Appl. Environ. Microbiol. 79, 6813–6822. doi: 10.1128/AEM.01937-13
Martinez-Garcia, M., Brazel, D., Poulton, N. J., Swan, B. K., Gomez, M. L., Masland, D.,et al. (2012). Unveiling in situ interactions between marine protists and bacteria through single cell sequencing. ISME J. 6, 703–707. doi: 10.1038/ismej.2011.126
Meyer, F., Paarmann, D., D’Souza, M., Olson, R., Glass, E. M., Kubal, M.,et al. (2008). The Metagenomics RAST server – A public resource for the automatic phylogenetic and functional analysis of metagenomes. BMC Bioinformatics 9:386. doi: 10.1186/1471-2105-9-386
Not, F., Siano, R., Kooistra, W. H. C. F., Simon, N., Vaulot, D., and Probert, I. (2012). “Chapter one – diversity and ecology of eukaryotic marine phytoplankton,” in Advances in Botanical Research, Vol. 64, ed. G. Piganeau (Waltham: Academic Press), 1–53.
Ortet, P., Barakat, M., Lalaouna, D., Fochesato, S., Barbe, V., Vacherie, B.,et al. (2011). Complete genome sequence of a beneficial plant root-associated bacterium, Pseudomonas brassicacearum. J. Bacteriol. 193, 3146. doi: 10.1128/JB.00411-11
Patzelt, D., Wang, H., Buchholz, I., Rohde, M., Groebe, L., Pradella, S.,et al. (2013). You are what you talk: quorum sensing induces individual morphologies and cell division modes in Dinoroseobacter shibae. ISME J. 7, 2274–2286. doi: 10.1038/ismej.2013.107
Piganeau, G., and Moreau, H. (2007). Screening the Sargasso Sea metagenome for data to investigate genome evolution in Ostreococcus (Prasinophyceae, Chlorophyta). Gene 406, 184–190. doi: 10.1016/j.gene.2007.09.015
Pukatzki, S., Ma, A. T., Revel, A. T., Sturtevant, D., and Mekalanos, J. J. (2007). Type VI secretion system translocates a phage tail spike-like protein into target cells where it cross-links actin. Proc. Natl. Acad. Sci. U.S.A. 104, 15508–15513. doi: 10.1073/pnas.0706532104
Pukatzki, S., Ma, A. T., Sturtevant, D., Krastins, B., Sarracino, D., Nelson, W. C.,et al. (2006). Identification of a conserved bacterial protein secretion system in Vibrio cholerae using the Dictyostelium host model system. Proc. Natl. Acad. Sci. U.S.A. 103, 1528–1533. doi: 10.1073/pnas.0510322103
Rooney-Varga, J. N., Giewat, M. W., Savin, M. C., Sood, S., LeGresley, M., and Martin, J. L. (2005). Links between phytoplankton and bacterial community dynamics in a coastal marine environment. Microb. Ecol. 49, 163–175. doi: 10.1007/s00248-003-1057-0
Sapp, M., Schwaderer, A. S., Wiltshire, K. H., Hoppe, H. G., Gerdts, G., and Wichels, A. (2007). Species-specific bacterial communities in the phycosphere of microalgae? Microb. Ecol. 53, 683–699. doi: 10.1007/s00248-006-9162-5
Seyedsayamdost, M. R., Carr, G., Kolter, R., and Clardy, J. (2011). Roseobacticides: small molecule modulators of an algal-bacterial symbiosis. J. Am. Chem. Soc. 133, 18343–18349. doi: 10.1021/ja207172s
Subirana, L., Péquin, B., Michely, S., Escande, M. L., Meilland, J., Derelle, E.,et al. (2013). Morphology, Genome Plasticity, and Phylogeny in the Genus Ostreococcus Reveal a Cryptic Species, O. mediterraneus sp. nov. (Mamiellales, Mamiellophyceae). Protist 164, 643–659. doi: 10.1016/j.protis.2013.06.002
Syrpas, M., Ruysbergh, E., Blommaert, L., Vanelslander, B., Sabbe, K., Vyverman, W.,et al. (2014). Haloperoxidase mediated quorum quenching by Nitzschia cf pellucida: study of the metabolization of N-Acyl homoserine lactones by a benthic diatom. Mar. Drugs 12, 352–367. doi: 10.3390/md12010352
Thompson, A. W., Foster, R. A., Krupke, A., Carter, B. J., Musat, N., Vaulot, D.,et al. (2012). Unicellular cyanobacterium symbiotic with a single-celled eukaryotic alga. Science 337, 1546–1550. doi: 10.1126/science.1222700
Vaulot, D., Lepère, C., Toulza, E., De la Iglesia, R., Poulain, J., Gaboyer, F.,et al. (2012). Metagenomes of the Picoalga Bathycoccus from the Chile Coastal Upwelling. PLoS ONE 7:e39648. doi: 10.1371/journal.pone.0039648
Wagner-Döbler, I., Ballhausen, B., Berger, M., Brinkhoff, T., Buchholz, I., Bunk, B.,et al. (2009). The complete genome sequence of the algal symbiont Dinoroseobacter shibae: a hitchhiker’s guide to life in the sea. ISME J. 4, 61–77. doi: 10.1038/ismej.2009.94
Wu, H. Y., Chung, P. C., Shih, H. W., Wen, S. R., and Lai, E. M. (2008). Secretome analysis uncovers an Hcp-family protein secreted via a type VI secretion system in Agrobacterium tumefaciens. J. Bacteriol. 190, 2841–2850. doi: 10.1128/JB.01775-07
Keywords: phytoplankton, bacterial symbiosis, secretion system, illumina sequencing, bacterial diversity, microbiome, phycosphere, Ostreococcus tauri
Citation: Abby SS, Touchon M, De Jode A, Grimsley N and Piganeau G (2014) Bacteria in Ostreococcus tauri cultures – friends, foes or hitchhikers? Front. Microbiol. 5:505. doi: 10.3389/fmicb.2014.00505
Received: 18 July 2014; Accepted: 08 September 2014;
Published online: 07 November 2014.
Edited by:Monica Medina, Pennsylvania State University, USA
Reviewed by:Fabrice Not, Centre National de la Recherche Scientifique, France
Scott Clingenpeel, United States Department of Energy Joint Genome Institute, USA
Copyright © 2014 Abby, Touchon, De Jode, Grimsley and Piganeau. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Gwenael Piganeau, CNRS, UMR 7232, Biologie Intégrative des Organismes Marins, Observatoire Océanologique, Sorbonne Universités – University Pierre and Marie Curie, 66650 Banyuls-sur-mer, France e-mail: firstname.lastname@example.org