Identification of DNA Viruses in Ancient DNA from Herbarium Samples.

Herbaria encompass millions of plant specimens, mostly collected in the nineteenth and twentieth centuries that can represent a key resource for investigating the history and evolution of phytopathogens. In the last years, the application of high-throughput sequencing technologies for the analysis of ancient nucleic acids has revolutionized the study of ancient pathogens including viruses, allowing the reconstruction of historical genomic viral sequences, improving phylogenetic based molecular dating, and providing essential insight into plant virus ecology. In this chapter, we describe a protocol to reconstruct ancient plant and soil viral sequences starting from highly fragmented ancient DNA extracted from herbarium plants and their associated rhizospheric soil. Following Illumina high-throughput sequencing, sequence data are de novo assembled, and DNA viral sequences are selected, according to their similarity with known viruses.

Beyond the limits of the unassigned protist microbiome: inferring large-scale spatio-temporal patterns of Syndiniales marine parasites.

Marine protists are major components of the oceanic microbiome that remain largely unrepresented in culture collections and genomic reference databases. The exploration of this uncharted protist diversity in oceanic communities relies essentially on studying genetic markers from the environment as taxonomic barcodes. Here we report that across 6 large scale spatio-temporal planktonic surveys, half of the genetic barcodes remain taxonomically unassigned at the genus level, preventing a fine ecological understanding for numerous protist lineages. Among them, parasitic Syndiniales (Dinoflagellata) appear as the least described protist group. We have developed a computational workflow, integrating diverse 18S rDNA gene metabarcoding datasets, in order to infer large-scale ecological patterns at 100% similarity of the genetic marker, overcoming the limitation of taxonomic assignment. From a spatial perspective, we identified 2171 unassigned clusters, i.e., Syndiniales sequences with 100% similarity, exclusively shared between the Tropical/Subtropical Ocean and the Mediterranean Sea among all Syndiniales orders and 25 ubiquitous clusters shared within all the studied marine regions. From a temporal perspective, over 3 time-series, we highlighted 39 unassigned clusters that follow rhythmic patterns of recurrence and are the best indicators of parasite community's variation. These clusters withhold potential as ecosystem change indicators, mirroring their associated host community responses. Our results underline the importance of Syndiniales in structuring planktonic communities through space and time, raising questions regarding host-parasite association specificity and the trophic mode of persistent Syndiniales, while providing an innovative framework for prioritizing unassigned protist taxa for further description.

Genomic differentiation of three pico-phytoplankton species in the Mediterranean Sea.

For more than a decade, high-throughput sequencing has transformed the study of marine planktonic communities and has highlighted the extent of protist diversity in these ecosystems. Nevertheless, little is known relative to their genomic diversity at the species-scale as well as their major speciation mechanisms. An increasing number of data obtained from global scale sampling campaigns is becoming publicly available, and we postulate that metagenomic data could contribute to deciphering the processes shaping protist genomic differentiation in the marine realm. As a proof of concept, we developed a findable, accessible, interoperable and reusable (FAIR) pipeline and focused on the Mediterranean Sea to study three a priori abundant protist species: Bathycoccus prasinos, Pelagomonas calceolata and Phaeocystis cordata. We compared the genomic differentiation of each species in light of geographic, environmental and oceanographic distances. We highlighted that isolation-by-environment shapes the genomic differentiation of B. prasinos, whereas P. cordata is impacted by geographic distance (i.e. isolation-by-distance). At present time, the use of metagenomics to accurately estimate the genomic differentiation of protists remains challenging since coverages are lower compared to traditional population surveys. However, our approach sheds light on ecological and evolutionary processes occurring within natural marine populations and paves the way for future protist population metagenomic studies.

Towards omics-based predictions of planktonic functional composition from environmental data.

Marine microbes play a crucial role in climate regulation, biogeochemical cycles, and trophic networks. Unprecedented amounts of data on planktonic communities were recently collected, sparking a need for innovative data-driven methodologies to quantify and predict their ecosystemic functions. We reanalyze 885 marine metagenome-assembled genomes through a network-based approach and detect 233,756 protein functional clusters, from which 15% are functionally unannotated. We investigate all clusters' distributions across the global ocean through machine learning, identifying biogeographical provinces as the best predictors of protein functional clusters' abundance. The abundances of 14,585 clusters are predictable from the environmental context, including 1347 functionally unannotated clusters. We analyze the biogeography of these 14,585 clusters, identifying the Mediterranean Sea as an outlier in terms of protein functional clusters composition. Applicable to any set of sequences, our approach constitutes a step towards quantitative predictions of functional composition from the environmental context.

Mixotrophic protists display contrasted biogeographies in the global ocean.

Mixotrophy, or the ability to acquire carbon from both auto- and heterotrophy, is a widespread ecological trait in marine protists. Using a metabarcoding dataset of marine plankton from the global ocean, 318,054 mixotrophic metabarcodes represented by 89,951,866 sequences and belonging to 133 taxonomic lineages were identified and classified into four mixotrophic functional types: constitutive mixotrophs (CM), generalist non-constitutive mixotrophs (GNCM), endo-symbiotic specialist non-constitutive mixotrophs (eSNCM), and plastidic specialist non-constitutive mixotrophs (pSNCM). Mixotrophy appeared ubiquitous, and the distributions of the four mixotypes were analyzed to identify the abiotic factors shaping their biogeographies. Kleptoplastidic mixotrophs (GNCM and pSNCM) were detected in new zones compared to previous morphological studies. Constitutive and non-constitutive mixotrophs had similar ranges of distributions. Most lineages were evenly found in the samples, yet some of them displayed strongly contrasted distributions, both across and within mixotypes. Particularly divergent biogeographies were found within endo-symbiotic mixotrophs, depending on the ability to form colonies or the mode of symbiosis. We showed how metabarcoding can be used in a complementary way with previous morphological observations to study the biogeography of mixotrophic protists and to identify key drivers of their biogeography.

Author Correction: Clade-specific diversification dynamics of marine diatoms since the Jurassic.

In the version of this Article originally published, the authors did not give credit to David G. Mann for the four microscopic images used in Fig. 1a. This has now been amended in all versions of the Article.

Clade-specific diversification dynamics of marine diatoms since the Jurassic.

Diatoms are one of the most abundant and diverse groups of phytoplankton and play a major role in marine ecosystems and the Earth's biogeochemical cycles. Here we combine DNA metabarcoding data from the Tara Oceans expedition with palaeoenvironmental data and phylogenetic models of diversification to analyse the diversity dynamics of marine diatoms. We reveal a primary effect of variation in carbon dioxide partial pressure (pCO2) on early diatom diversification, followed by a major burst of diversification in the late Eocene epoch, after which diversification is chiefly affected by sea level, an influx of silica availability and competition with other planktonic groups. Our results demonstrate a remarkable heterogeneity of diversification dynamics across diatoms and suggest that a changing climate will favour some clades at the expense of others.

A de novo approach to disentangle partner identity and function in holobiont systems.

BACKGROUND: Study of meta-transcriptomic datasets involving non-model organisms represents bioinformatic challenges. The production of chimeric sequences and our inability to distinguish the taxonomic origins of the sequences produced are inherent and recurrent difficulties in de novo assembly analyses. As the study of holobiont meta-transcriptomes is affected by challenges invoked above, we propose an innovative bioinformatic approach to tackle such difficulties and tested it on marine models as a proof of concept. RESULTS: We considered three holobiont models, of which two transcriptomes were previously published and a yet unpublished transcriptome, to analyze and sort their raw reads using Short Read Connector, a k-mer based similarity method. Before assembly, we thus defined four distinct categories for each holobiont meta-transcriptome: host reads, symbiont reads, shared reads, and unassigned reads. Afterwards, we observed that independent de novo assemblies for each category led to a diminution of the number of chimeras compared to classical assembly methods. Moreover, the separation of each partner's transcriptome offered the independent and comparative exploration of their functional diversity in the holobiont. Finally, our strategy allowed to propose new functional annotations for two well-studied holobionts (a Cnidaria-Dinophyta, a Porifera-Bacteria) and a first meta-transcriptome from a planktonic Radiolaria-Dinophyta system forming widespread symbiotic association for which our knowledge is considerably limited. CONCLUSIONS: In contrast to classical assembly approaches, our bioinformatic strategy generates less de novo assembled chimera and allows biologists to study separately host and symbiont data from a holobiont mixture. The pre-assembly separation of reads using an efficient tool as Short Read Connector is an effective way to tackle meta-transcriptomic challenges and offers bright perpectives to study holobiont systems composed of either well-studied or poorly characterized symbiotic lineages and ultimately expand our knowledge about these associations.

Analysis of the genomic basis of functional diversity in dinoflagellates using a transcriptome-based sequence similarity network.

Dinoflagellates are one of the most abundant and functionally diverse groups of eukaryotes. Despite an overall scarcity of genomic information for dinoflagellates, constantly emerging high-throughput sequencing resources can be used to characterize and compare these organisms. We assembled de novo and processed 46 dinoflagellate transcriptomes and used a sequence similarity network (SSN) to compare the underlying genomic basis of functional features within the group. This approach constitutes the most comprehensive picture to date of the genomic potential of dinoflagellates. A core-predicted proteome composed of 252 connected components (CCs) of putative conserved protein domains (pCDs) was identified. Of these, 206 were novel and 16 lacked any functional annotation in public databases. Integration of functional information in our network analyses allowed investigation of pCDs specifically associated with functional traits. With respect to toxicity, sequences homologous to those of proteins found in species with toxicity potential (e.g., sxtA4 and sxtG) were not specific to known toxin-producing species. Although not fully specific to symbiosis, the most represented functions associated with proteins involved in the symbiotic trait were related to membrane processes and ion transport. Overall, our SSN approach led to identification of 45,207 and 90,794 specific and constitutive pCDs of, respectively, the toxic and symbiotic species represented in our analyses. Of these, 56% and 57%, respectively (i.e., 25,393 and 52,193 pCDs), completely lacked annotation in public databases. This stresses the extent of our lack of knowledge, while emphasizing the potential of SSNs to identify candidate pCDs for further functional genomic characterization.

The epibiotic life of the cosmopolitan diatom Fragilariopsis doliolus on heterotrophic ciliates in the open ocean.

Diatoms are a diverse and ecologically important group of phytoplankton. Although most species are considered free living, several are known to interact with other organisms within the plankton. Detailed imaging and molecular characterization of any such partnership is, however, limited, and an appraisal of the large-scale distribution and ecology of such consortia was never attempted. Here, observation of Tara Oceans samples from the Benguela Current led to the detection of an epibiotic association between a pennate diatom and a tintinnid ciliate. We identified the diatom as Fragilariopsis doliolus that possesses a unique feature to form barrel-shaped chains, associated with seven different genera of tintinnids including five previously undescribed associations. The organisms were commonly found together in the Atlantic and Pacific Ocean basins, and live observations of the interaction have been recorded for the first time. By combining confocal and scanning electron microscopy of individual consortia with the sequencing of high-resolution molecular markers, we analyzed their distribution in the global ocean, revealing morpho-genetically distinct tintinnid haplotypes and biogeographically structured diatom haplotypes. The diatom was among the most abundant in the global ocean. We show that the consortia were particularly prevalent in nutrient-replete conditions, rich in potential predators. These observations support the hypothesis of a mutualistic symbiosis, wherein diatoms acquire increased motility and tintinnids benefit from silicification through increased protection, and highlight that such associations may be more prevalent than currently appreciated.

The evolution of diatoms and their biogeochemical functions.

In contemporary oceans diatoms are an important group of eukaryotic phytoplankton that typically dominate in upwelling regions and at high latitudes. They also make significant contributions to sporadic blooms that often occur in springtime. Recent surveys have revealed global information about their abundance and diversity, as well as their contributions to biogeochemical cycles, both as primary producers of organic material and as conduits facilitating the export of carbon and silicon to the ocean interior. Sequencing of diatom genomes is revealing the evolutionary underpinnings of their ecological success by examination of their gene repertoires and the mechanisms they use to adapt to environmental changes. The rise of the diatoms over the last hundred million years is similarly being explored through analysis of microfossils and biomarkers that can be traced through geological time, as well as their contributions to seafloor sediments and fossil fuel reserves. The current review aims to synthesize current information about the evolution and biogeochemical functions of diatoms as they rose to prominence in the global ocean.This article is part of the themed issue 'The peculiar carbon metabolism in diatoms'.

Acclimation of a low iron adapted Ostreococcus strain to iron limitation through cell biomass lowering.

Iron is an essential micronutrient involved in many biological processes and is often limiting for primary production in large regions of the World Ocean. Metagenomic and physiological studies have identified clades or ecotypes of marine phytoplankton that are specialized in iron depleted ecological niches. Although less studied, eukaryotic picophytoplankton does contribute significantly to primary production and carbon transfer to higher trophic levels. In particular, metagenomic studies of the green picoalga Ostreococcus have revealed the occurrence of two main clades distributed along coast-offshore gradients, suggesting niche partitioning in different nutrient regimes. Here, we present a study of the response to iron limitation of four Ostreococcus strains isolated from contrasted environments. Whereas the strains isolated in nutrient-rich waters showed high iron requirements, the oceanic strains could cope with lower iron concentrations. The RCC802 strain, in particular, was able to maintain high growth rate at low iron levels. Together physiological and transcriptomic data indicate that the competitiveness of RCC802 under iron limitation is related to a lowering of iron needs though a reduction of the photosynthetic machinery and of protein content, rather than to cell size reduction. Our results overall suggest that iron is one of the factors driving the differentiation of physiologically specialized Ostreococcus strains in the ocean.

Benthic protists: the under-charted majority.

Marine protist diversity inventories have largely focused on planktonic environments, while benthic protists have received relatively little attention. We therefore hypothesize that current diversity surveys have only skimmed the surface of protist diversity in marine sediments, which may harbor greater diversity than planktonic environments. We tested this by analyzing sequences of the hypervariable V4 18S rRNA from benthic and planktonic protist communities sampled in European coastal regions. Despite a similar number of OTUs in both realms, richness estimations indicated that we recovered at least 70% of the diversity in planktonic protist communities, but only 33% in benthic communities. There was also little overlap of OTUs between planktonic and benthic communities, as well as between separate benthic communities. We argue that these patterns reflect the heterogeneity and diversity of benthic habitats. A comparison of all OTUs against the Protist Ribosomal Reference database showed that a higher proportion of benthic than planktonic protist diversity is missing from public databases; similar results were obtained by comparing all OTUs against environmental references from NCBI's Short Read Archive. We suggest that the benthic realm may therefore be the world's largest reservoir of marine protist diversity, with most taxa at present undescribed.

Insights into global diatom distribution and diversity in the world's ocean.

Diatoms (Bacillariophyta) constitute one of the most diverse and ecologically important groups of phytoplankton. They are considered to be particularly important in nutrient-rich coastal ecosystems and at high latitudes, but considerably less so in the oligotrophic open ocean. The Tara Oceans circumnavigation collected samples from a wide range of oceanic regions using a standardized sampling procedure. Here, a total of ∼12 million diatom V9-18S ribosomal DNA (rDNA) ribotypes, derived from 293 size-fractionated plankton communities collected at 46 sampling sites across the global ocean euphotic zone, have been analyzed to explore diatom global diversity and community composition. We provide a new estimate of diversity of marine planktonic diatoms at 4,748 operational taxonomic units (OTUs). Based on the total assigned ribotypes, Chaetoceros was the most abundant and diverse genus, followed by Fragilariopsis, Thalassiosira, and Corethron We found only a few cosmopolitan ribotypes displaying an even distribution across stations and high abundance, many of which could not be assigned with confidence to any known genus. Three distinct communities from South Pacific, Mediterranean, and Southern Ocean waters were identified that share a substantial percentage of ribotypes within them. Sudden drops in diversity were observed at Cape Agulhas, which separates the Indian and Atlantic Oceans, and across the Drake Passage between the Atlantic and Southern Oceans, indicating the importance of these ocean circulation choke points in constraining diatom distribution and diversity. We also observed high diatom diversity in the open ocean, suggesting that diatoms may be more relevant in these oceanic systems than generally considered.

Plankton networks driving carbon export in the oligotrophic ocean.

The biological carbon pump is the process by which CO2 is transformed to organic carbon via photosynthesis, exported through sinking particles, and finally sequestered in the deep ocean. While the intensity of the pump correlates with plankton community composition, the underlying ecosystem structure driving the process remains largely uncharacterized. Here we use environmental and metagenomic data gathered during the Tara Oceans expedition to improve our understanding of carbon export in the oligotrophic ocean. We show that specific plankton communities, from the surface and deep chlorophyll maximum, correlate with carbon export at 150 m and highlight unexpected taxa such as Radiolaria and alveolate parasites, as well as Synechococcus and their phages, as lineages most strongly associated with carbon export in the subtropical, nutrient-depleted, oligotrophic ocean. Additionally, we show that the relative abundance of a few bacterial and viral genes can predict a significant fraction of the variability in carbon export in these regions.

Marine protist diversity in European coastal waters and sediments as revealed by high-throughput sequencing.

Although protists are critical components of marine ecosystems, they are still poorly characterized. Here we analysed the taxonomic diversity of planktonic and benthic protist communities collected in six distant European coastal sites. Environmental deoxyribonucleic acid (DNA) and ribonucleic acid (RNA) from three size fractions (pico-, nano- and micro/mesoplankton), as well as from dissolved DNA and surface sediments were used as templates for tag pyrosequencing of the V4 region of the 18S ribosomal DNA. Beta-diversity analyses split the protist community structure into three main clusters: picoplankton-nanoplankton-dissolved DNA, micro/mesoplankton and sediments. Within each cluster, protist communities from the same site and time clustered together, while communities from the same site but different seasons were unrelated. Both DNA and RNA-based surveys provided similar relative abundances for most class-level taxonomic groups. Yet, particular groups were overrepresented in one of the two templates, such as marine alveolates (MALV)-I and MALV-II that were much more abundant in DNA surveys. Overall, the groups displaying the highest relative contribution were Dinophyceae, Diatomea, Ciliophora and Acantharia. Also, well represented were Mamiellophyceae, Cryptomonadales, marine alveolates and marine stramenopiles in the picoplankton, and Monadofilosa and basal Fungi in sediments. Our extensive and systematic sequencing of geographically separated sites provides the most comprehensive molecular description of coastal marine protist diversity to date.

Ocean plankton. Environmental characteristics of Agulhas rings affect interocean plankton transport.

Agulhas rings provide the principal route for ocean waters to circulate from the Indo-Pacific to the Atlantic basin. Their influence on global ocean circulation is well known, but their role in plankton transport is largely unexplored. We show that, although the coarse taxonomic structure of plankton communities is continuous across the Agulhas choke point, South Atlantic plankton diversity is altered compared with Indian Ocean source populations. Modeling and in situ sampling of a young Agulhas ring indicate that strong vertical mixing drives complex nitrogen cycling, shaping community metabolism and biogeochemical signatures as the ring and associated plankton transit westward. The peculiar local environment inside Agulhas rings may provide a selective mechanism contributing to the limited dispersal of Indian Ocean plankton populations into the Atlantic.

Ocean plankton. Eukaryotic plankton diversity in the sunlit ocean.

Marine plankton support global biological and geochemical processes. Surveys of their biodiversity have hitherto been geographically restricted and have not accounted for the full range of plankton size. We assessed eukaryotic diversity from 334 size-fractionated photic-zone plankton communities collected across tropical and temperate oceans during the circumglobal Tara Oceans expedition. We analyzed 18S ribosomal DNA sequences across the intermediate plankton-size spectrum from the smallest unicellular eukaryotes (protists, >0.8 micrometers) to small animals of a few millimeters. Eukaryotic ribosomal diversity saturated at ~150,000 operational taxonomic units, about one-third of which could not be assigned to known eukaryotic groups. Diversity emerged at all taxonomic levels, both within the groups comprising the ~11,200 cataloged morphospecies of eukaryotic plankton and among twice as many other deep-branching lineages of unappreciated importance in plankton ecology studies. Most eukaryotic plankton biodiversity belonged to heterotrophic protistan groups, particularly those known to be parasites or symbiotic hosts.

Ocean plankton. Determinants of community structure in the global plankton interactome.

Species interaction networks are shaped by abiotic and biotic factors. Here, as part of the Tara Oceans project, we studied the photic zone interactome using environmental factors and organismal abundance profiles and found that environmental factors are incomplete predictors of community structure. We found associations across plankton functional types and phylogenetic groups to be nonrandomly distributed on the network and driven by both local and global patterns. We identified interactions among grazers, primary producers, viruses, and (mainly parasitic) symbionts and validated network-generated hypotheses using microscopy to confirm symbiotic relationships. We have thus provided a resource to support further research on ocean food webs and integrating biological components into ocean models.

Seasonal diversity and dynamics of haptophytes in the Skagerrak, Norway, explored by high-throughput sequencing.

Microalgae in the division Haptophyta play key roles in the marine ecosystem and in global biogeochemical processes. Despite their ecological importance, knowledge on seasonal dynamics, community composition and abundance at the species level is limited due to their small cell size and few morphological features visible under the light microscope. Here, we present unique data on haptophyte seasonal diversity and dynamics from two annual cycles, with the taxonomic resolution and sampling depth obtained with high-throughput sequencing. From outer Oslofjorden, S Norway, nano- and picoplanktonic samples were collected monthly for 2 years, and the haptophytes targeted by amplification of RNA/cDNA with Haptophyta-specific 18S rDNA V4 primers. We obtained 156 operational taxonomic units (OTUs), from c. 400.000 454 pyrosequencing reads, after rigorous bioinformatic filtering and clustering at 99.5%. Most OTUs represented uncultured and/or not yet 18S rDNA-sequenced species. Haptophyte OTU richness and community composition exhibited high temporal variation and significant yearly periodicity. Richness was highest in September-October (autumn) and lowest in April-May (spring). Some taxa were detected all year, such as Chrysochromulina simplex, Emiliania huxleyi and Phaeocystis cordata, whereas most calcifying coccolithophores only appeared from summer to early winter. We also revealed the seasonal dynamics of OTUs representing putative novel classes (clades HAP-3-5) or orders (clades D, E, F). Season, light and temperature accounted for 29% of the variation in OTU composition. Residual variation may be related to biotic factors, such as competition and viral infection. This study provides new, in-depth knowledge on seasonal diversity and dynamics of haptophytes in North Atlantic coastal waters.