Assessing the viral content of uncultured picoeukaryotes in the global-ocean by single cell genomics.

Viruses are the most abundant biological entities on Earth and have fundamental ecological roles in controlling microbial communities. Yet, although their diversity is being increasingly explored, little is known about the extent of viral interactions with their protist hosts as most studies are limited to a few cultivated species. Here, we exploit the potential of single-cell genomics to unveil viral associations in 65 individual cells of 11 essentially uncultured stramenopiles lineages sampled during the Tara Oceans expedition. We identified viral signals in 57% of the cells, covering nearly every lineage and with narrow host specificity signal. Only seven out of the 64 detected viruses displayed homologies to known viral sequences. A search for our viral sequences in global ocean metagenomes showed that they were preferentially found at the DCM and within the 0.2-3 µm size fraction. Some of the viral signals were widely distributed, while others geographically constrained. Among the viral signals we detected an endogenous mavirus virophage potentially integrated within the nuclear genome of two distant uncultured stramenopiles. Virophages have been previously reported as a cell's defence mechanism against other viruses, and may therefore play an important ecological role in regulating protist populations. Our results point to single-cell genomics as a powerful tool to investigate viral associations in uncultured protists, suggesting a wide distribution of these relationships, and providing new insights into the global viral diversity.

Constant abundances of ubiquitous uncultured protists in the open sea assessed by automated microscopy.

Protists have fundamental ecological roles in marine environments and their diversity is being increasingly explored, yet little is known about the quantitative importance of specific taxa in these ecosystems. Here we optimized a newly developed automated system of image acquisition and image analysis to enumerate minute uncultured cells of different sizes targeted by fluorescence in situ hybridization. The automated counting routine was highly reproducible, well correlated with manual counts, and was then applied on surface and deep chlorophyll maximum samples from the Malaspina 2010 circumnavigation. The three targeted uncultured taxa (MAST-4, MAST-7 and MAST-1C) were found in virtually all samples from several ocean basins (Atlantic, Indian and Pacific) in fairly constant cell abundances, following typical lognormal distributions. Their global abundances averaged 49, 23 and 7 cells ml-1 , respectively, and altogether the three groups accounted for about 10%-20% of heterotrophic picoeukaryotes. Our innovative high-throughput cell counting routine allows for the first time a direct assessment of the biogeographic distribution of small protists (< 5 µm) and shows the ubiquity in sunlit oceans of three bacterivorous taxa, suggesting their key roles in marine ecosystems.

Single-cell genomics of multiple uncultured stramenopiles reveals underestimated functional diversity across oceans.

Single-celled eukaryotes (protists) are critical players in global biogeochemical cycling of nutrients and energy in the oceans. While their roles as primary producers and grazers are well appreciated, other aspects of their life histories remain obscure due to challenges in culturing and sequencing their natural diversity. Here, we exploit single-cell genomics and metagenomics data from the circumglobal Tara Oceans expedition to analyze the genome content and apparent oceanic distribution of seven prevalent lineages of uncultured heterotrophic stramenopiles. Based on the available data, each sequenced genome or genotype appears to have a specific oceanic distribution, principally correlated with water temperature and depth. The genome content provides hypotheses for specialization in terms of cell motility, food spectra, and trophic stages, including the potential impact on their lifestyles of horizontal gene transfer from prokaryotes. Our results support the idea that prominent heterotrophic marine protists perform diverse functions in ocean ecology.

X-Cells Are Globally Distributed, Genetically Divergent Fish Parasites Related to Perkinsids and Dinoflagellates.

"X-cells" have long been associated with tumor-like formations (xenomas) in marine fish, including many of commercial interest. The name was first used to refer to the large polygonal cells that were found in epidermal xenomas from flatfish from the Pacific Northwest [1]. Similar looking cells from pseudobranchial xenomas had previously been reported from cod in the Atlantic [2] and Pacific Oceans [3]. X-cell pathologies have been reported from five teleost orders: Pleuronectiformes (flatfish), Perciformes (perch-like fish), Gadiformes (cods), Siluriformes (catfish), and Salmoniformes (salmonids). Various explanations have been elicited for their etiology, including being adenomas or adenocarcinomas [4, 5], virally transformed fish cells [6-8], or products of coastal pollution [9, 10]. It was hypothesized that X-cells were protozoan parasites [1, 11-13], and although recent molecular analyses have confirmed this, they have failed to place them in any phylum [14-18], demonstrating weak phylogenetic associations with the haplosporidians [16] or the alveolates [15]. Here, we sequenced rRNA genes from European and Japanese fish that are known to develop X-cell xenomas. We also generated a metagenomic sequence library from X-cell xenomas of blue whiting and Atlantic cod and assembled 63 X-cell protein-coding genes for a eukaryote-wide phylogenomic analysis. We show that X-cells group in two highly divergent clades, robustly sister to the bivalve parasite Perkinsus. We formally describe these as Gadixcellia and Xcellia and provide a phylogenetic context to catalyze future research. We also screened Atlantic cod populations for xenomas and residual pathologies and show that X-cell infections are more prevalent and widespread than previously known.

Temperature regulation of marine heterotrophic prokaryotes increases latitudinally as a breach between bottom-up and top-down controls.

Planktonic heterotrophic prokaryotes make up the largest living biomass and process most organic matter in the ocean. Determining when and where the biomass and activity of heterotrophic prokaryotes are controlled by resource availability (bottom-up), predation and viral lysis (top-down) or temperature will help in future carbon cycling predictions. We conducted an extensive survey across subtropical and tropical waters of the Atlantic, Indian and Pacific Oceans during the Malaspina 2010 Global Circumnavigation Expedition and assessed indices for these three types of controls at 109 stations (mostly from the surface to 4,000 m depth). Temperature control was approached by the apparent activation energy in eV (ranging from 0.46 to 3.41), bottom-up control by the slope of the log-log relationship between biomass and production rate (ranging from -0.12 to 1.09) and top-down control by an index that considers the relative abundances of heterotrophic nanoflagellates and viruses (ranging from 0.82 to 4.83). We conclude that temperature becomes dominant (i.e. activation energy >1.5 eV) within a narrow window of intermediate values of bottom-up (0.3-0.6) and top-down 0.8-1.2) controls. A pervasive latitudinal pattern of decreasing temperature regulation towards the Equator, regardless of the oceanic basin, suggests that the impact of global warming on marine microbes and their biogeochemical function will be more intense at higher latitudes. Our analysis predicts that 1°C ocean warming will result in increased biomass of heterotrophic prokaryoplankton only in waters with <26°C of mean annual surface temperature.

Accessing the genomic information of unculturable oceanic picoeukaryotes by combining multiple single cells.

Pico-sized eukaryotes play key roles in the functioning of marine ecosystems, but we still have a limited knowledge on their ecology and evolution. The MAST-4 lineage is of particular interest, since it is widespread in surface oceans, presents ecotypic differentiation and has defied culturing efforts so far. Single cell genomics (SCG) are promising tools to retrieve genomic information from these uncultured organisms. However, SCG are based on whole genome amplification, which normally introduces amplification biases that limit the amount of genomic data retrieved from a single cell. Here, we increase the recovery of genomic information from two MAST-4 lineages by co-assembling short reads from multiple Single Amplified Genomes (SAGs) belonging to evolutionary closely related cells. We found that complementary genomic information is retrieved from different SAGs, generating co-assembly that features >74% of genome recovery, against about 20% when assembled individually. Even though this approach is not aimed at generating high-quality draft genomes, it allows accessing to the genomic information of microbes that would otherwise remain unreachable. Since most of the picoeukaryotes still remain uncultured, our work serves as a proof-of-concept that can be applied to other taxa in order to extract genomic data and address new ecological and evolutionary questions.

Rarity in aquatic microbes: placing protists on the map.

Most microbial richness at any given time tends to be represented by low-abundance (rare) taxa, which are collectively referred to as the "rare biosphere". Here we review works on the rare biosphere using high-throughput sequencing (HTS), with a particular focus on unicellular eukaryotes or protists. Evidence thus far indicates that the rare biosphere encompasses dormant as well as metabolically active microbes that could potentially play key roles in ecosystem functioning. Rare microbes appear to have biogeography, and sometimes the observed patterns can be similar to what is observed among abundant taxa, suggesting similar community-structuring mechanisms. There is limited evidence indicating that the rare biosphere contains taxa that are phylogenetically distantly related to abundant counterparts; therefore, the rare biosphere may act as a reservoir of deep-branching phylogenetic diversity. The potential role of the rare biosphere as a bank of redundant functions that can help to maintain continuous ecosystem function following oscillations in taxonomic abundances is hypothesized as its main ecological role. Future studies focusing on rare microbes are crucial for advancing our knowledge of microbial ecology and evolution and unveiling their links with ecosystem function.

Geographic distance and ecosystem size determine the distribution of smallest protists in lacustrine ecosystems.

Understanding the spatial distribution of aquatic microbial diversity and the underlying mechanisms causing differences in community composition is a challenging and central goal for ecologists. Recent insights into protistan diversity and ecology are increasing the debate over their spatial distribution. In this study, we investigate the importance of spatial and environmental factors in shaping the small protists community structure in lakes. We analyzed small protists community composition (beta-diversity) and richness (alpha-diversity) at regional scale by different molecular methods targeting the gene coding for 18S rRNA gene (T-RFLP and 454 pyrosequencing). Our results show a distance-decay pattern for rare and dominant taxa and the spatial distribution of the latter followed the prediction of the island biogeography theory. Furthermore, geographic distances between lakes seem to be the main force shaping the protists community composition in the lakes studied here. Finally, the spatial distribution of protists was discussed at the global scale (11 worldwide distributed lakes) by comparing these results with those present in the public database. UniFrac analysis showed 18S rRNA gene OTUs compositions significantly different among most of lakes, and this difference does not seem to be related to the trophic status.

Phylogenetic affiliation of SSU rRNA genes generated by massively parallel sequencing: new insights into the freshwater protist diversity.

Recent advances in next-generation sequencing (NGS) technologies spur progress in determining the microbial diversity in various ecosystems by highlighting, for example, the rare biosphere. Currently, high-throughput pyrotag sequencing of PCR-amplified SSU rRNA gene regions is mainly used to characterize bacterial and archaeal communities, and rarely to characterize protist communities. In addition, although taxonomic assessment through phylogeny is considered as the most robust approach, similarity and probabilistic approaches remain the most commonly used for taxonomic affiliation. In a first part of this work, a tree-based method was compared with different approaches of taxonomic affiliation (BLAST and RDP) of 18S rRNA gene sequences and was shown to be the most accurate for near full-length sequences and for 400 bp amplicons, with the exception of amplicons covering the V5-V6 region. Secondly, the applicability of this method was tested by running a full scale test using an original pyrosequencing dataset of 18S rRNA genes of small lacustrine protists (0.2-5 µm) from eight freshwater ecosystems. Our results revealed that i) fewer than 5% of the operational taxonomic units (OTUs) identified through clustering and phylogenetic affiliation had been previously detected in lakes, based on comparison to sequence in public databases; ii) the sequencing depth provided by the NGS coupled with a phylogenetic approach allowed to shed light on clades of freshwater protists rarely or never detected with classical molecular ecology approaches; and iii) phylogenetic methods are more robust in describing the structuring of under-studied or highly divergent populations. More precisely, new putative clades belonging to Mamiellophyceae, Foraminifera, Dictyochophyceae and Euglenida were detected. Beyond the study of protists, these results illustrate that the tree-based approach for NGS based diversity characterization allows an in-depth description of microbial communities including taxonomic profiling, community structuring and the description of clades of any microorganisms (protists, Bacteria and Archaea).

Short-term dynamics of diversity patterns: evidence of continual reassembly within lacustrine small eukaryotes.

The short-term variation in the community structure of freshwater small eukaryotes (0.2-5 µm) was investigated in a mesotrophic lake every 2-3 days over one summer by coupling three molecular methods: 454 amplicon pyrosequencing, qPCR and TSA-FISH. The pyrosequencing approach unveiled a much more extensive small-eukaryotic diversity (991 OTUs) than has been described previously. The vast majority of the diversity described was represented by rare OTUs (≤ 0.01% of reads) belonging primarily to Cryptomycota, Dikarya and photosynthetic organisms, which were never detected as abundant in any of the samples. The small eukaryote community was characterized by a continual and important reassembly. These rearrangements involved the 20 'core taxa' (≥ 1% of reads), and, were essentially due to a handful of OTUs that were detected in intermediate abundance (0.01-1% of reads) and sporadically in dominant taxa. Putative bacterivorous (Ciliophora and Cercozoa) as well as parasitic and saprotrophic taxa (Perkinsozoa and Cryptomycota) were involved in these changes of diversity. A putative infection of microalgae by a lacustrine perkinsozoan was also reported for the first time in this study. Open questions regarding both the patterns that govern the rapid small eukaryote reassemblies and the possible biogeography of these organisms arise from this study.

Diversity and dynamics of free-living and particle-associated Betaproteobacteria and Actinobacteria in relation to phytoplankton and zooplankton communities.

The diversity of attached and free-living Actinobacteria and Betaproteobacteria, based on 16S rRNA gene sequences, was investigated in a mesotrophic lake during two periods of contrasting phytoplankton dominance. Comparison analyses showed a phylogenetic difference between attached and free-living communities for the two bacterial groups. For Betaproteobacteria, the betaI clade was detected at all sampling dates in free-living and attached bacterial communities and was the dominant clade contributing to 57.8% of the total retrieved operational taxonomic units (OTUs). For Actinobacteria, the acIV cluster was found to be dominant, followed by acI contributing to 45% and 25% of the total retrieved OTUs, respectively. This study allows the determination of eight new putative clades among the Betaproteobacteria termed lbI-lbVIII and a new putative clade named acLBI belonging to the Actinobacteria. The seasonal dynamics of phytoplankton and zooplankton communities have been reflected as changes in distinct bacterial phylotypes for both attached and free-living communities. For attached communities, relationships were observed between Actinobacteria and Chrysophyceae, and between Betaproteobacteria and Dinophyceae and Chlorophyceae biomass. On the other hand, within free-living communities, few actinobacterial clades were found to be dependent on either nutrients or phytoplankton communities, whereas Betaproteobacteria were mainly associated with biological parameters (i.e. phytoplankton and copepod communities).

Vertical structure of small eukaryotes in three lakes that differ by their trophic status: a quantitative approach.

In lakes, the diversity of eukaryotic picoplankton has been recently studied by the analysis of 18S ribosomal RNA gene sequences; however, quantitative data are rare. In this study, the vertical structure and abundance of the small eukaryotic size fraction (0.2-5 µm) were investigated in three lakes by tyramide signal amplification-fluorescent in situ hybridization targeting six phylogenetic groups: Chlorophyta, Haptophyta, Cercozoa, LKM11, Perkinsozoa and fungi. The groups targeted in this study are found in all lakes; however, both the abundance and structure of small eukaryotes are dependent on the system's productivity and depth. These data highlighted the presence of Chlorophyta contributing on an average to 19.3%, 14.7% and 41.2% of total small eukaryotes in lakes Bourget, Aydat and Pavin, respectively. This study also revealed the unexpected importance of Haptophyta, reaching 62.8% of eukaryotes in the euphotic zone of Lake Bourget. The high proportions of these pigmented cells highlight the underestimation of these groups by PCR-based methods. The presence of pigmented Chlorophyta in the deepest zones of the lakes suggests a mixotrophic behaviour of these taxa. We also confirmed the presence of putative parasites such as Perkinsozoa (5.1% of small eukaryotes in Lake Pavin and Bourget) and, with lower abundances, fungi (targeted by the MY1574 probe). Cells targeted by LKM11 probes represented the second group of abundance within heterotrophs. Open questions regarding the functional roles of the targeted groups arise from this study, especially regarding parasitism and mixotrophy, which are interactions poorly taken into account in planktonic food web models.

Community structure and dynamics of small eukaryotes targeted by new oligonucleotide probes: new insight into the lacustrine microbial food web.

The seasonal dynamics of the small eukaryotic fraction (cell diameter, 0.2 to 5 microm) was investigated in a mesotrophic lake by tyramide signal amplification-fluorescence in situ hybridization targeting seven different phylogenetic groups: Chlorophyceae, Chrysophyceae, Cryptophyceae, Cercozoa, LKM11, Perkinsozoa (two clades), and Fungi. The abundance of small eukaryotes ranged from 1,692 to 10,782 cells ml(-1). The dominant groups were the Chrysophyceae and the Chlorophyceae, which represented 19.6% and 17.9% of small eukaryotes, respectively. The results also confirmed the quantitative importance of putative parasites, Fungi and Perkinsozoa, in the small heterotrophic eukaryotic assemblage. The relative abundances recorded for the Perkinsozoa group reached as much as 31.6% of total targeted eukaryotes during the summer. The dynamics of Perkinsozoa clade 1 coincided with abundance variations in Peridinium and Ceratium spp. (Dinoflagellates), while the dynamics of Perkinsozoa clade 2 was linked to the presence of Dinobryon spp. (Chrysophyceae). Fungi, represented by chytrids, reached maximal abundance in December (569 cells ml(-1)) and were mainly correlated with the dynamics of diatoms, especially Melosira varians. A further new finding of this study is the recurrent presence of Cercozoa (6.2%) and LKM11 (4.5%) cells. This quantitative approach based on newly designed probes offers a promising means of in-depth analysis of microbial food webs in lakes, especially by revealing the phylogenetic composition of the small heterotrophic flagellate assemblage, for which an important fraction of cells are generally unidentified by classical microscopy (on average, 96.8% of the small heterotrophic flagellates were identified by the specific probes we used in this study).

Isolation and characterization of mesotrione-degrading Bacillus sp. from soil.

Dissipation kinetics of mesotrione, a new triketone herbicide, sprayed on soil from Limagne (Puy-de-Dôme, France) showed that the soil microflora were able to biotransform it. Bacteria from this soil were cultured in mineral salt solution supplemented with mesotrione as sole source of carbon for the isolation of mesotrione-degrading bacteria. The bacterial community structure of the enrichment cultures was analyzed by temporal temperature gradient gel electrophoresis (TTGE). The TTGE fingerprints revealed that mesotrione had an impact on bacterial community structure only at its highest concentrations and showed mesotrione-sensitive and mesotrione-adapted strains. Two adapted strains, identified as Bacillus sp. and Arthrobacter sp., were isolated by colony hybridization methods. Biodegradation assays showed that only the Bacillus sp. strain was able to completely and rapidly biotransform mesotrione. Among several metabolites formed, 2-amino-4-methylsulfonylbenzoic acid (AMBA) accumulated in the medium. Although sulcotrione has a chemical structure closely resembling that of mesotrione, the isolates were unable to degrade it.