A holistic approach to marine eco-systems biology

Global diversity and distribution of aerobic anoxygenic phototrophs in the tropical and subtropical oceans

The aerobic anoxygenic phototrophic (AAP) bacteria are common in most marine environments but their global diversity and biogeography remain poorly characterized. Here, we analyzed AAP communities across 113 globally-distributed surface ocean stations sampled during the Malaspina Expedition in the tropical and subtropical ocean. By means of amplicon sequencing of the pufM gene, a genetic marker for this functional group, we show that AAP communities along the surface ocean were mainly composed of members of the Halieaceae (Gammaproteobacteria), which were adapted to a large range of environmental conditions, and of different clades of the Alphaproteobacteria, which seemed to dominate under particular circumstances, such as in the oligotrophic gyres. AAP taxa were spatially structured within each of the studied oceans, with communities from adjacent stations sharing more taxonomic similarities. AAP communities were composed of a large pool of rare members and several habitat specialists. When compared to the surface ocean prokaryotic and picoeukaryotic communities, it appears that AAP communities display an idiosyncratic global biogeographical pattern, dominated by selection processes and less influenced by dispersal limitation. Our study contributes to the understanding of how AAP communities are distributed in the horizontal dimension and the mechanisms underlying their distribution across the global surface ocean.

CANT-HYD: A Curated Database of Phylogeny-Derived Hidden Markov Models for Annotation of Marker Genes Involved in Hydrocarbon Degradation

Many pathways for hydrocarbon degradation have been discovered, yet there are no dedicated tools to identify and predict the hydrocarbon degradation potential of microbial genomes and metagenomes. Here we present the Calgary approach to ANnoTating HYDrocarbon degradation genes (CANT-HYD), a database of 37 HMMs of marker genes involved in anaerobic and aerobic degradation pathways of aliphatic and aromatic hydrocarbons. Using this database, we identify understudied or overlooked hydrocarbon degradation potential in many phyla. We also demonstrate its application in analyzing high-throughput sequence data by predicting hydrocarbon utilization in large metagenomic datasets from diverse environments. CANT-HYD is available at https://github.com/dgittins/CANT-HYD-HydrocarbonBiodegradation.

OUP accepted manuscript

Oxygen (O₂) is the ultimate oxidant on Earth and its respiration confers such an energetic advantage that microorganisms have evolved the capacity to scavenge O₂ down to nanomolar concentrations. The respiration of O₂ at extremely low levels is proving to be common to diverse microbial taxa, including organisms formerly considered strict anaerobes. Motivated by recent advances in O₂ sensing and DNA/RNA sequencing technologies, we performed a systematic review of environmental metatranscriptomes revealing that microbial respiration of O₂ at nanomolar concentrations is ubiquitous and drives microbial activity in seemingly anoxic aquatic habitats. These habitats were key to the early evolution of life and are projected to become more prevalent in the near future due to anthropogenic-driven environmental change. Here, we summarize our current understanding of aerobic microbial respiration under apparent anoxia, including novel processes, their underlying biochemical pathways, the involved microorganisms, and their environmental importance and evolutionary origin.

A novel, divergent alkane monooxygenase (alkB) clade involved in crude oil biodegradation

Alkanes are ubiquitous in marine ecosystems and originate from diverse sources ranging from natural oil seeps to anthropogenic inputs and biogenic production by cyanobacteria. Enzymes that degrade cyanobacterial alkanes (typically C15-C17 compounds) such as the alkane monooxygenase (AlkB) are widespread, but it remains unclear whether or not AlkB variants exist that specialize in degradation of crude oil from natural or accidental spills, a much more complex mixture of long-chain hydrocarbons. In the present study, large-scale analysis of available metagenomic and genomic data from the Gulf of Mexico (GoM) oil spill revealed a novel, divergent AlkB clade recovered from genomes with no cultured representatives that was dramatically increased in abundance in crude-oil impacted ecosystems. In contrast, the AlkB clades associated with biotransformation of cyanobacterial alkanes belonged to 'canonical' or hydrocarbonoclastic clades, and based on metatranscriptomics data and compared to the novel clade, were much more weakly expressed during crude oil biodegradation in laboratory mesocosms. The absence of this divergent AlkB clade in metagenomes of uncontaminated samples from the global ocean survey but not from the GoM as well as its frequent horizontal gene transfer indicated a priming effect of the Gulf for crude oil biodegradation likely driven by natural oil seeps.

High variability in SSU rDNA gene copy number among planktonic foraminifera revealed by single-cell qPCR

Metabarcoding has become the workhorse of community ecology. Sequencing a taxonomically informative DNA fragment from environmental samples gives fast access to community composition across taxonomic groups, but it relies on the assumption that the number of sequences for each taxon correlates with its abundance in the sampled community. However, gene copy number varies among and within taxa, and the extent of this variability must therefore be considered when interpreting community composition data derived from environmental sequencing. Here we measured with single-cell qPCR the SSU rDNA gene copy number of 139 specimens of five species of planktonic foraminifera. We found that the average gene copy number varied between of ~4000 to ~50,000 gene copies between species, and individuals of the same species can carry between ~300 to more than 350,000 gene copies. This variability cannot be explained by differences in cell size and considering all plausible sources of bias, we conclude that this variability likely reflects dynamic genomic processes acting during the life cycle. We used the observed variability to model its impact on metabarcoding and found that the application of a correcting factor at species level may correct the derived relative abundances, provided sufficiently large populations have been sampled.

BLight: efficient exact associative structure for k-mers

MOTIVATION

A plethora of methods and applications share the fundamental need to associate information to words for high-throughput sequence analysis. Doing so for billions of k-mers is commonly a scalability problem, as exact associative indexes can be memory expensive. Recent works take advantage of overlaps between k-mers to leverage this challenge. Yet, existing data structures are either unable to associate information to k-mers or are not lightweight enough.

RESULTS

We present BLight, a static and exact data structure able to associate unique identifiers to k-mers and determine their membership in a set without false positive that scales to huge k-mer sets with a low memory cost. This index combines an extremely compact representation along with very fast queries. Besides, its construction is efficient and needs no additional memory. Our implementation achieves to index the k-mers from the human genome using 8 GB of RAM (23 bits per k-mer) within 10 min and the k-mers from the large axolotl genome using 63 GB of memory (27 bits per k-mer) within 76 min. Furthermore, while being memory efficient, the index provides a very high throughput: 1.4 million queries per second on a single CPU or 16.1 million using 12 cores. Finally, we also present how BLight can practically represent metagenomic and transcriptomic sequencing data to highlight its wide applicative range.

AVAILABILITY AND IMPLEMENTATION

We wrote the BLight index as an open source C++ library under the AGPL3 license available at github.com/Malfoy/BLight. It is designed as a user-friendly library and comes along with code usage samples.

Microbiome Studies from Saudi Arabia over the Last 10 Years: Achievements, Gaps, and Future Directions

In the past ten years, microbiome studies have shown tremendous potentiality for implementation of understanding microbiome structures and functions of various biomes and application of this knowledge for human betterment. Saudi Arabia is full of geographical, ecological, ethnical, and industrial diversities and scientific capacities. Therefore, there is a great potential in Saudi Arabia to conduct and implement microbiome-based research and applications. However, there is no review available on where Saudi Arabia stands with respect to global microbiome research trends. This review highlights the metagenome-assisted microbiome research from Saudi Arabia compared to the global focuses on microbiome research. Further, it also highlights the gaps and areas that should be focused on by Saudi microbiome researchers and the possible initiatives to be taken by Saudi government and universities. This literature review shows that the global trends of microbiome research cover a broad spectrum of human and animal health conditions and diseases, environmental and antimicrobial resistance surveillance, surveillance of food and food processing, production of novel industrial enzymes and bioactive pharmaceutical products, and space applications. However, Saudi microbiome studies are mostly confined to very few aspects of health (human and animal) and environment/ecology in last ten years, without much application. Therefore, Saudi Arabia should focus more on applied microbiome research through government, academic, and industry initiatives and global cooperation to match the global trends.

Prophage Tracer: precisely tracing prophages in prokaryotic genomes using overlapping split-read alignment

The life cycle of temperate phages includes a lysogenic cycle stage when the phage integrates into the host genome and becomes a prophage. However, the identification of prophages that are highly divergent from known phages remains challenging. In this study, by taking advantage of the lysis-lysogeny switch of temperate phages, we designed Prophage Tracer, a tool for recognizing active prophages in prokaryotic genomes using short-read sequencing data, independent of phage gene similarity searching. Prophage Tracer uses the criterion of overlapping split-read alignment to recognize discriminative reads that contain bacterial (attB) and phage (attP) att sites representing prophage excision signals. Performance testing showed that Prophage Tracer could predict known prophages with precise boundaries, as well as novel prophages. Two novel prophages, dsDNA and ssDNA, encoding highly divergent major capsid proteins, were identified in coral-associated bacteria. Prophage Tracer is a reliable data mining tool for the identification of novel temperate phages and mobile genetic elements. The code for the Prophage Tracer is publicly available at https://github.com/WangLab-SCSIO/Prophage_Tracer.

Improved high throughput protocol for targeting eukaryotic symbionts in metazoan and eDNA samples

Eukaryote symbionts of animals are major drivers of ecosystems not only because of their diversity and host interactions from variable pathogenicity but also through different key roles such as commensalism and to different types of interdependence. However, molecular investigations of metazoan eukaryomes require minimising coamplification of homologous host genes. In this study we (1) identified a previously published “antimetazoan” reverse primer to theoretically enable amplification of a wider range of microeukaryotic symbionts, including more evolutionarily divergent sequence types, (2) evaluated in silico several antimetazoan primer combinations, and (3) optimised the application of the best performing primer pair for high throughput sequencing (HTS) by comparing one‐step and two‐step PCR amplification approaches, testing different annealing temperatures and evaluating the taxonomic profiles produced by HTS and data analysis. The primer combination 574*F – UNonMet_DB tested in silico showed the largest diversity of nonmetazoan sequence types in the SILVA database and was also the shortest available primer combination for broadly‐targeting antimetazoan amplification across the 18S rRNA gene V4 region. We demonstrate that the one‐step PCR approach used for library preparation produces significantly lower proportions of metazoan reads, and a more comprehensive coverage of host‐associated microeukaryote reads than the two‐step approach. Using higher PCR annealing temperatures further increased the proportion of nonmetazoan reads in all sample types tested. The resulting V4 region amplicons were taxonomically informative even when only the forward read is analysed. This region also revealed a diversity of known and putatively parasitic lineages and a wider diversity of host‐associated eukaryotes.

Human impacts on deep-sea sponge grounds: Applying environmental omics to monitoring

Sponges (Phylum Porifera) are the oldest extant Metazoans. In the deep sea, sponges can occur at high densities forming habitats known as sponge grounds. Sponge grounds can extend over large areas of up to hundreds of km² and are biodiversity hotspots. However, as human activities, including deep-water hydrocarbon extraction, continue to expand into areas harbouring sponge grounds, understanding how anthropogenic impacts affect sponges and the ecosystem services they provide at multiple biological scales (community, individual and (sub)cellular levels) is key for achieving sustainable management. This chapter (1) provides an update to the chapter of Advances in Marine Biology Volume 79 entitled "Potential Impacts of Offshore Oil and Gas Activities on Deep-Sea Sponges and the Habitats They Form" and (2) discusses the use of omics as a future tool for deep-sea ecosystem monitoring. While metagenomics and (meta)transcriptomics studies have contributed to improve our understanding of sponge biology in recent years, metabolomics analysis has mostly been used to identify natural products. The sponge metabolome, therefore, remains vastly unknown despite the fact that the metabolome is a key link between the genotype and phenotype, giving us a unique new insight to how key components of an ecosystem are functioning. As the fraction of the metabolome released into the seawater, the sponge exometabolome has only just started to be characterised in comparative environmental metabolomic studies. Yet, the sponge exometabolome constitute a unique opportunity for the identification of biomarkers of sponge health as compounds can be measured in seawater, bypassing the need for physical samples which can still be difficult to collect in the deep sea. Within sponge grounds, the characterisation of a shared sponge exometabolome could lead to the identification of biomarkers of ecosystem functioning and overall health. Challenges remain in establishing omics approaches in environmental monitoring but constant technological advances and reduction in costs means these techniques will become widely available in the future.

Biotechnological applications of marine bacteria in bioremediation of environments polluted with hydrocarbons and plastics

Marine ecosystems are some of the most adverse environments on Earth and contain a considerable portion of the global bacterial population, and some of these bacterial species play pivotal roles in several biogeochemical cycles. Marine bacteria have developed different molecular mechanisms to address fluctuating environmental conditions, such as changes in nutrient availability, salinity, temperature, pH, and pressure, making them attractive for use in diverse biotechnology applications. Although more than 99% of marine bacteria cannot be cultivated with traditional microbiological techniques, several species have been successfully isolated and grown in the laboratory, facilitating investigations of their biotechnological potential. Some of these applications may contribute to addressing some current global problems, such as environmental contamination by hydrocarbons and synthetic plastics. In this review, we first summarize and analyze recently published information about marine bacterial diversity. Then, we discuss new literature regarding the isolation and characterization of marine bacterial strains able to degrade hydrocarbons and petroleum-based plastics, and species able to produce biosurfactants. We also describe some current limitations for the implementation of these biotechnological tools, but also we suggest some strategies that may contribute to overcoming them. KEY POINTS: • Marine bacteria have a great metabolic capacity to degrade hydrocarbons in harsh conditions. • Marine environments are an important source of new bacterial plastic-degrading enzymes. • Secondary metabolites from marine bacteria have diverse potential applications in biotechnology.

Environmental vulnerability of the global ocean epipelagic plankton community interactome

Marine plankton form complex communities of interacting organisms at the base of the food web, which sustain oceanic biogeochemical cycles and help regulate climate. Although global surveys are starting to reveal ecological drivers underlying planktonic community structure and predicted climate change responses, it is unclear how community-scale species interactions will be affected by climate change. Here, we leveraged Tara Oceans sampling to infer a global ocean cross-domain plankton co-occurrence network-the community interactome-and used niche modeling to assess its vulnerabilities to environmental change. Globally, this revealed a plankton interactome self-organized latitudinally into marine biomes (Trades, Westerlies, Polar) and more connected poleward. Integrated niche modeling revealed biome-specific community interactome responses to environmental change and forecasted the most affected lineages for each community. These results provide baseline approaches to assess community structure and organismal interactions under climate scenarios while identifying plausible plankton bioindicators for ocean monitoring of climate change.

Ecosystems monitoring powered by environmental genomics: A review of current strategies with an implementation roadmap

A decade after environmental scientists integrated high-throughput sequencing technologies in their toolbox, the genomics-based monitoring of anthropogenic impacts on the biodiversity and functioning of ecosystems is yet to be implemented by regulatory frameworks. Despite the broadly acknowledged potential of environmental genomics to this end, technical limitations and conceptual issues still stand in the way of its broad application by end-users. In addition, the multiplicity of potential implementation strategies may contribute to a perception that the routine application of this methodology is premature or "in development", hence restraining regulators from binding these tools into legal frameworks. Here, we review recent implementations of environmental genomics-based methods, applied to the biomonitoring of ecosystems. By taking a general overview, without narrowing our perspective to particular habitats or groups of organisms, this paper aims to compare, review and discuss the strengths and limitations of four general implementation strategies of environmental genomics for monitoring: (a) Taxonomy-based analyses focused on identification of known bioindicators or described taxa; (b) De novo bioindicator analyses; (c) Structural community metrics including inferred ecological networks; and (d) Functional community metrics (metagenomics or metatranscriptomics). We emphasise the utility of the three latter strategies to integrate meiofauna and microorganisms that are not traditionally utilised in biomonitoring because of difficult taxonomic identification. Finally, we propose a roadmap for the implementation of environmental genomics into routine monitoring programmes that leverage recent analytical advancements, while pointing out current limitations and future research needs.

Phage-encoded ten-eleven translocation dioxygenase (TET) is active in C5-cytosine hypermodification in DNA

Chemical tailoring of canonical bases expands the functionality of DNA in the same manner that posttranscriptional and -translational modifications enhance functional diversity in RNA and proteins. We describe the activities of ten-eleven translocation dioxygenase (TET)-like iron(II)- and 2-oxo-glutarate–dependent 5mC dioxygenases that are encoded by several bacteriophages to enable hypermodification of C5-methyl cytosine bases in their DNA. Phage TETs act on methylation marks deposited within GpC sequences by functionally-associated cytosine 5-methyltransferases. The hydroxymethyl groups installed are further elaborated by tailoring enzymes, thereby decorating the phage DNA with diverse, complex modifications. These modifications are predicted to have protective roles against host defenses during viral infection.

TET/JBP (ten-eleven translocation/base J binding protein) enzymes are iron(II)- and 2-oxo-glutarate–dependent dioxygenases that are found in all kingdoms of life and oxidize 5-methylpyrimidines on the polynucleotide level. Despite their prevalence, few examples have been biochemically characterized. Among those studied are the metazoan TET enzymes that oxidize 5-methylcytosine in DNA to hydroxy, formyl, and carboxy forms and the euglenozoa JBP dioxygenases that oxidize thymine in the first step of base J biosynthesis. Both enzymes have roles in epigenetic regulation. It has been hypothesized that all TET/JBPs have their ancestral origins in bacteriophages, but only eukaryotic orthologs have been described. Here we demonstrate the 5mC-dioxygenase activity of several phage TETs encoded within viral metagenomes. The clustering of these TETs in a phylogenetic tree correlates with the sequence specificity of their genomically cooccurring cytosine C5-methyltransferases, which install the methyl groups upon which TETs operate. The phage TETs favor Gp5mC dinucleotides over the 5mCpG sites targeted by the eukaryotic TETs and are found within gene clusters specifying complex cytosine modifications that may be important for DNA packaging and evasion of host restriction.

Physiological determinants of biogeography: The importance of metabolic depression to heat tolerance

A quantitative understanding of physiological thermal responses is vital for forecasting species distributional shifts in response to climate change. Many studies have focused on metabolic rate as a global metric for analyzing the sublethal effects of changing environments on physiology. Thermal performance curves (TPCs) have been suggested as a viable analytical framework, but standard TPCs may not fully capture physiological responses, due in part to failure to consider the process of metabolic depression. We derived a model based on the nonlinear regression of biological temperature-dependent rate processes and built a heart rate data set for 26 species of intertidal molluscs distributed from 33°S to ~40°N. We then calculated physiological thermal performance limits with continuous heating using T 1 / 2 H , the temperature at which heart rate is decreased to 50% of the maximal rate, as a more realistic measure of upper thermal limits. Results indicate that heat-induced metabolic depression of cardiac performance is a common adaptive response that allows tolerance of harsh environments. Furthermore, our model accounted for the high inter-individual variability in the shape of cardiac TPCs. We then used these TPCs to calculate physiological thermal safety margins (pTSM), the difference between the maximal operative temperature (95th percentile of field temperatures) and T 1 / 2 H of each individual. Using pTSMs, we developed a physiological species distribution model (pSDM) to forecast future geographic distributions. pSDM results indicate that climate-induced species range shifts are potentially less severe than predicted by a simple correlative SDM. Species with metabolic depression below the optimum temperature will be more thermal resistant at their warm trailing edges. High intraspecific variability further suggests that models based on species-level vulnerability to environmental change may be problematic. This multi-scale, mechanistic understanding that incorporates metabolic depression and inter-individual variability in thermal response enables better predictions about the relationship between thermal stress and species distributions.

Density Peak clustering of protein sequences associated to a Pfam clan reveals clear similarities and interesting differences with respect to manual family annotation

BACKGROUND

The identification of protein families is of outstanding practical importance for in silico protein annotation and is at the basis of several bioinformatic resources. Pfam is possibly the most well known protein family database, built in many years of work by domain experts with extensive use of manual curation. This approach is generally very accurate, but it is quite time consuming and it may suffer from a bias generated from the hand-curation itself, which is often guided by the available experimental evidence.

RESULTS

We introduce a procedure that aims to identify automatically putative protein families. The procedure is based on Density Peak Clustering and uses as input only local pairwise alignments between protein sequences. In the experiment we present here, we ran the algorithm on about 4000 full-length proteins with at least one domain classified by Pfam as belonging to the Pseudouridine synthase and Archaeosine transglycosylase (PUA) clan. We obtained 71 automatically-generated sequence clusters with at least 100 members. While our clusters were largely consistent with the Pfam classification, showing good overlap with either single or multi-domain Pfam family architectures, we also observed some inconsistencies. The latter were inspected using structural and sequence based evidence, which suggested that the automatic classification captured evolutionary signals reflecting non-trivial features of protein family architectures. Based on this analysis we identified a putative novel pre-PUA domain as well as alternative boundaries for a few PUA or PUA-associated families. As a first indication that our approach was unlikely to be clan-specific, we performed the same analysis on the P53 clan, obtaining comparable results.

CONCLUSIONS

The clustering procedure described in this work takes advantage of the information contained in a large set of pairwise alignments and successfully identifies a set of putative families and family architectures in an unsupervised manner. Comparison with the Pfam classification highlights significant overlap and points to interesting differences, suggesting that our new algorithm could have potential in applications related to automatic protein classification. Testing this hypothesis, however, will require further experiments on large and diverse sequence datasets.

Connecting high-throughput biodiversity inventories: Opportunities for a site-based genomic framework for global integration and synthesis

High-throughput sequencing (HTS) is increasingly being used for the characterization and monitoring of biodiversity. If applied in a structured way, across broad geographical scales, it offers the potential for a much deeper understanding of global biodiversity through the integration of massive quantities of molecular inventory data generated independently at local, regional and global scales. The universality, reliability and efficiency of HTS data can potentially facilitate the seamless linking of data among species assemblages from different sites, at different hierarchical levels of diversity, for any taxonomic group and regardless of prior taxonomic knowledge. However, collective international efforts are required to optimally exploit the potential of site-based HTS data for global integration and synthesis, efforts that at present are limited to the microbial domain. To contribute to the development of an analogous strategy for the nonmicrobial terrestrial domain, an international symposium entitled "Next Generation Biodiversity Monitoring" was held in November 2019 in Nicosia (Cyprus). The symposium brought together evolutionary geneticists, ecologists and biodiversity scientists involved in diverse regional and global initiatives using HTS as a core tool for biodiversity assessment. In this review, we summarize the consensus that emerged from the 3-day symposium. We converged on the opinion that an effective terrestrial Genomic Observatories network for global biodiversity integration and synthesis should be spatially led and strategically united under the umbrella of the metabarcoding approach. Subsequently, we outline an HTS-based strategy to collectively build an integrative framework for site-based biodiversity data generation.

Description of Candidatus Mesopelagibacter carboxydoxydans and Candidatus Anoxipelagibacter denitrificans: Nitrate-reducing SAR11 genera that dominate mesopelagic and anoxic marine zones

The diverse and ubiquitous members of the SAR11 lineage (Alphaproteobacteria) represent up to 30-40% of the surface and mesopelagic oceanic microbial communities. However, the molecular and ecological mechanisms that differentiate closely related, yet distinct, SAR11 members that often co-occur under similar environmental conditions remain speculative. Recently, two mesopelagic and oxygen minimum zone (OMZ)-associated subclades of SAR11 (Ic and IIa.A) were described using single-cell amplified genomes (SAGs) linked to nitrate reduction in OMZs. In this current study, the collection of genomes belonging to these two subclades was expanded with thirteen new metagenome-assembled genomes (MAGs), thus providing a more detailed phylogenetic and functional characterization of these subclades. Gene content-based predictions of metabolic functions revealed similarities in central carbon metabolism between subclades Ic and IIa.A and surface SAR11 clades, with small variations in central pathways. These variations included more versatile sulfur assimilation pathways, as well as a previously predicted capacity for nitrate reduction that conferred unique versatility on mesopelagic-adapted clades compared to their surface counterparts. Finally, consistent with previously reported abundances of carbon monoxide (CO) in surface and mesopelagic waters, subclades Ia (surface) and Ic (mesopelagic) have the genetic potential to oxidize carbon monoxide (CO), presumably taking advantage of this abundant compound as an electron donor. Based on genomic analyses, environmental distribution and metabolic reconstruction, we propose two new SAR11 genera, Ca. Mesopelagibacter carboxydoxydans (subclade Ic) and Ca. Anoxipelagibacter denitrificans (subclade IIa.A), which represent members of the mesopelagic and OMZ-adapted SAR11 clades.

The Ocean barcode atlas: A web service to explore the biodiversity and biogeography of marine organisms

The Ocean Barcode Atlas (OBA) is a user friendly web service designed for biologists who wish to explore the biodiversity and biogeography of marine organisms locked in otherwise difficult to mine planetary scale DNA metabarcode data sets. Using just a web browser, a comprehensive picture of the diversity of a taxon or a barcode sequence is visualized graphically on world maps and interactive charts. Interactive results panels allow dynamic threshold adjustments and the display of diversity results in their environmental context measured at the time of sampling (temperature, oxygen, latitude, etc). Ecological analyses such as alpha and beta-diversity plots are produced via publication quality vector graphics representations. Currently, the Ocean Barcode Altas is deployed online with the (i) Tara Oceans eukaryotic 18S-V9 rDNA metabarcodes; (ii) Tara Oceans 16S/18S rRNA _mi Tags; and (iii) 16S-V4 V5 metabarcodes collected during the Malaspina-2010 expedition. Additional prokaryotic or eukaryotic plankton barcode data sets will be added upon availability, given they provide the required complement of barcodes (including raw reads to compute barcode abundance) associated with their contextual environmental variables. Ocean Barcode Atlas is a freely-available web service at: http://oba.mio.osupytheas.fr/ocean-atlas/.

Metagenomics for taxonomy profiling: tools and approaches

The study of metagenomics is an emerging field that identifies the total genetic materials in an organism along with the set of all genetic materials like deoxyribonucleic acid and ribose nucleic acid, which play a key role with the maintenance of cellular functions. The best part of this technology is that it gives more flexibility to environmental microbiologists to instantly pioneer the immense genetic variability of microbial communities. However, it is intensively complex to identify the suitable sequencing measures of any specific gene that can exclusively indicate the involvement of microbial metagenomes and be able to advance valuable results about these communities. This review provides an overview of the metagenomic advancement that has been advantageous for aggregation of more knowledge about specific genes, microbial communities and its metabolic pathways. More specific drawbacks of metagenomes technology mainly depend on sequence-based analysis. Therefore, this 'targeted based metagenomics' approach will give comprehensive knowledge about the ecological, evolutionary and functional sequence of significantly important genes that naturally exist in living beings either human, animal and microorganisms from distinctive ecosystems.

Comprehensive screening of genomic and metagenomic data reveals a large diversity of tetracycline resistance genes

Tetracyclines are broad-spectrum antibiotics used to prevent or treat a variety of bacterial infections. Resistance is often mediated through mobile resistance genes, which encode one of the three main mechanisms: active efflux, ribosomal target protection or enzymatic degradation. In the last few decades, a large number of new tetracycline-resistance genes have been discovered in clinical settings. These genes are hypothesized to originate from environmental and commensal bacteria, but the diversity of tetracycline-resistance determinants that have not yet been mobilized into pathogens is unknown. In this study, we aimed to characterize the potential tetracycline resistome by screening genomic and metagenomic data for novel resistance genes. By using probabilistic models, we predicted 1254 unique putative tetracycline resistance genes, representing 195 gene families (<70 % amino acid sequence identity), whereof 164 families had not been described previously. Out of 17 predicted genes selected for experimental verification, 7 induced a resistance phenotype in an Escherichia coli host. Several of the predicted genes were located on mobile genetic elements or in regions that indicated mobility, suggesting that they easily can be shared between bacteria. Furthermore, phylogenetic analysis indicated several events of horizontal gene transfer between bacterial phyla. Our results also suggested that acquired efflux pumps originate from proteobacterial species, while ribosomal protection genes have been mobilized from Firmicutes and Actinobacteria. This study significantly expands the knowledge of known and putatively novel tetracycline resistance genes, their mobility and evolutionary history. The study also provides insights into the unknown resistome and genes that may be encountered in clinical settings in the future.

Sequencing effort dictates gene discovery in marine microbial metagenomes

Massive metagenomic sequencing combined with gene prediction methods were previously used to compile the gene catalogue of the ocean and host-associated microbes. Global expeditions conducted over the past 15 years have sampled the ocean to build a catalogue of genes from pelagic microbes. Here we undertook a large sequencing effort of a perturbed Red Sea plankton community to uncover that the rate of gene discovery increases continuously with sequencing effort, with no indication that the retrieved 2.83 million non-redundant (complete) genes predicted from the experiment represented a nearly complete inventory of the genes present in the sampled community (i.e., no evidence of saturation). The underlying reason is the Pareto-like distribution of the abundance of genes in the plankton community, resulting in a very long tail of millions of genes present at remarkably low abundances, which can only be retrieved through massive sequencing. Microbial metagenomic projects retrieve a variable number of unique genes per Tera base-pair (Tbp), with a median value of 14.7 million unique genes per Tbp sequenced across projects. The increase in the rate of gene discovery in microbial metagenomes with sequencing effort implies that there is ample room for new gene discovery in further ocean and holobiont sequencing studies.

Investigating population-scale allelic differential expression in wild populations of Oithona similis (Cyclopoida, Claus, 1866)

Acclimation allowed by variation in gene or allele expression in natural populations is increasingly understood as a decisive mechanism, as much as adaptation, for species evolution. However, for small eukaryotic organisms, as species from zooplankton, classical methods face numerous challenges. Here, we propose the concept of allelic differential expression at the population-scale (psADE) to investigate the variation in allele expression in natural populations. We developed a novel approach to detect psADE based on metagenomic and metatranscriptomic data from environmental samples. This approach was applied on the widespread marine copepod, Oithona similis, by combining samples collected during the Tara Oceans expedition (2009-2013) and de novo transcriptome assemblies. Among a total of 25,768 single nucleotide variants (SNVs) of O. similis, 572 (2.2%) were affected by psADE in at least one population (FDR < 0.05). The distribution of SNVs under psADE in different populations is significantly shaped by population genomic differentiation (Pearson r = 0.87, p = 5.6 × 10-30), supporting a partial genetic control of psADE. Moreover, a significant amount of SNVs (0.6%) were under both selection and psADE (p < .05), supporting the hypothesis that natural selection and psADE tends to impact common loci. Population-scale allelic differential expression offers new insights into the gene regulation control in populations and its link with natural selection.

Diversity and distribution of marine heterotrophic bacteria from a large culture collection

Background

Isolation of marine microorganisms is fundamental to gather information about their physiology, ecology and genomic content. To date, most of the bacterial isolation efforts have focused on the photic ocean leaving the deep ocean less explored. We have created a marine culture collection of heterotrophic bacteria (MARINHET) using a standard marine medium comprising a total of 1561 bacterial strains, and covering a variety of oceanographic regions from different seasons and years, from 2009 to 2015. Specifically, our marine collection contains isolates from both photic (817) and aphotic layers (744), including the mesopelagic (362) and the bathypelagic (382), from the North Western Mediterranean Sea, the North and South Atlantic Ocean, the Indian, the Pacific, and the Arctic Oceans. We described the taxonomy, the phylogenetic diversity and the biogeography of a fraction of the marine culturable microorganisms to enhance our knowledge about which heterotrophic marine isolates are recurrently retrieved across oceans and along different depths.

Results

The partial sequencing of the 16S rRNA gene of all isolates revealed that they mainly affiliate with the classes Alphaproteobacteria (35.9%), Gammaproteobacteria (38.6%), and phylum Bacteroidetes (16.5%). In addition, Alteromonas and Erythrobacter genera were found the most common heterotrophic bacteria in the ocean growing in solid agar medium. When comparing all photic, mesopelagic, and bathypelagic isolates sequences retrieved from different stations, 37% of them were 100% identical. This percentage increased up to 59% when mesopelagic and bathypelagic strains were grouped as the aphotic dataset and compared to the photic dataset of isolates, indicating the ubiquity of some bacterial isolates along different ocean depths. Finally, we isolated three strains that represent a new species, and the genome comparison and phenotypic characterization of two of these strains (ISS653 and ISS1889) concluded that they belong to a new species within the genus Mesonia.

Conclusions

Overall, this study highlights the relevance of culture-dependent studies, with focus on marine isolated bacteria from different oceanographic regions and depths, to provide a more comprehensive view of the culturable marine bacteria as part of the total marine microbial diversity.

Alteromonas Myovirus V22 Represents a New Genus of Marine Bacteriophages Requiring a Tail Fiber Chaperone for Host Recognition

Marine phages play a variety of critical roles in regulating the microbial composition of our oceans. Despite constituting the majority of genetic diversity within these environments, there are relatively few isolates with complete genome sequences or in-depth analyses of their host interaction mechanisms, such as characterization of their receptor binding proteins (RBPs). Here, we present the 92,760-bp genome of the Alteromonas-targeting phage V22. Genomic and morphological analyses identify V22 as a myovirus; however, due to a lack of sequence similarity to any other known myoviruses, we propose that V22 be classified as the type phage of a new Myoalterovirus genus within the Myoviridae family. V22 shows gene homology and synteny with two different subfamilies of phages infecting enterobacteria, specifically within the structural region of its genome. To improve our understanding of the V22 adsorption process, we identified putative RBPs (gp23, gp24, and gp26) and tested their ability to decorate the V22 propagation strain, Alteromonas mediterranea PT11, as recombinant green fluorescent protein (GFP)-tagged constructs. Only GFP-gp26 was capable of bacterial recognition and identified as the V22 RBP. Interestingly, production of functional GFP-gp26 required coexpression with the downstream protein gp27. GFP-gp26 could be expressed alone but was incapable of host recognition. By combining size-exclusion chromatography with fluorescence microscopy, we reveal how gp27 is not a component of the final RBP complex but instead is identified as a new type of phage-encoded intermolecular chaperone that is essential for maturation of the gp26 RBP.IMPORTANCE Host recognition by phage-encoded receptor binding proteins (RBPs) constitutes the first step in all phage infections and the most critical determinant of host specificity. By characterizing new types of RBPs and identifying their essential chaperones, we hope to expand the repertoire of known phage-host recognition machineries. Due to their genetic plasticity, studying RBPs and their associated chaperones can shed new light onto viral evolution affecting phage-host interactions, which is essential for fields such as phage therapy or biotechnology. In addition, since marine phages constitute one of the most important reservoirs of noncharacterized genetic diversity on the planet, their genomic and functional characterization may be of paramount importance for the discovery of novel genes with potential applications.

Diversity and evolution of bacterial bioluminescence genes in the global ocean

Although bioluminescent bacteria are the most abundant and widely distributed of all light-emitting organisms, the biological role and evolutionary history of bacterial luminescence are still shrouded in mystery. Bioluminescence has so far been observed in the genomes of three families of Gammaproteobacteria in the form of canonical lux operons that adopt the CDAB(F)E(G) gene order. LuxA and luxB encode the two subunits of bacterial luciferase responsible for light-emission. Our deep exploration of public marine environmental databases considerably expands this view by providing a catalog of new lux homolog sequences, including 401 previously unknown luciferase-related genes. It also reveals a broader diversity of the lux operon organization, which we observed in previously undescribed configurations such as CEDA, CAED and AxxCE. This expanded operon diversity provides clues for deciphering lux operon evolution and propagation within the bacterial domain. Leveraging quantitative tracking of marine bacterial genes afforded by planetary scale metagenomic sampling, our study also reveals that the novel lux genes and operons described herein are more abundant in the global ocean than the canonical CDAB(F)E(G) operon.

Global Trends in Marine Plankton Diversity across Kingdoms of Life

The ocean is home to myriad small planktonic organisms that underpin the functioning of marine ecosystems. However, their spatial patterns of diversity and the underlying drivers remain poorly known, precluding projections of their responses to global changes. Here we investigate the latitudinal gradients and global predictors of plankton diversity across archaea, bacteria, eukaryotes, and major virus clades using both molecular and imaging data from Tara Oceans. We show a decline of diversity for most planktonic groups toward the poles, mainly driven by decreasing ocean temperatures. Projections into the future suggest that severe warming of the surface ocean by the end of the 21^st century could lead to tropicalization of the diversity of most planktonic groups in temperate and polar regions. These changes may have multiple consequences for marine ecosystem functioning and services and are expected to be particularly significant in key areas for carbon sequestration, fisheries, and marine conservation.

Video Abstract

•

Most epipelagic planktonic groups exhibit a poleward decline of diversity

•

No latitudinal diversity gradient was observed below the photic zone

•

Temperature emerges as the best predictor of epipelagic plankton diversity

•

Global warming may increase plankton diversity, particularly at high latitudes

Applications of weighted association networks applied to compositional data in biology

Next-generation sequencing technologies have generated, and continue to produce, an increasingly large corpus of biological data. The data generated are inherently compositional as they convey only relative information dependent upon the capacity of the instrument, experimental design and technical bias. There is considerable information to be gained through network analysis by studying the interactions between components within a system. Network theory methods using compositional data are powerful approaches for quantifying relationships between biological components and their relevance to phenotype, environmental conditions or other external variables. However, many of the statistical assumptions used for network analysis are not designed for compositional data and can bias downstream results. In this mini-review, we illustrate the utility of network theory in biological systems and investigate modern techniques while introducing researchers to frameworks for implementation. We overview (1) compositional data analysis, (2) data transformations and (3) network theory along with insight on a battery of network types including static-, temporal-, sample-specific- and differential-networks. The intention of this mini-review is not to provide a comprehensive overview of network methods, rather to introduce microbiology researchers to (semi)-unsupervised data-driven approaches for inferring latent structures that may give insight into biological phenomena or abstract mechanics of complex systems.

MGnify: the microbiome analysis resource in 2020

MGnify (http://www.ebi.ac.uk/metagenomics) provides a free to use platform for the assembly, analysis and archiving of microbiome data derived from sequencing microbial populations that are present in particular environments. Over the past 2 years, MGnify (formerly EBI Metagenomics) has more than doubled the number of publicly available analysed datasets held within the resource. Recently, an updated approach to data analysis has been unveiled (version 5.0), replacing the previous single pipeline with multiple analysis pipelines that are tailored according to the input data, and that are formally described using the Common Workflow Language, enabling greater provenance, reusability, and reproducibility. MGnify's new analysis pipelines offer additional approaches for taxonomic assertions based on ribosomal internal transcribed spacer regions (ITS1/2) and expanded protein functional annotations. Biochemical pathways and systems predictions have also been added for assembled contigs. MGnify's growing focus on the assembly of metagenomic data has also seen the number of datasets it has assembled and analysed increase six-fold. The non-redundant protein database constructed from the proteins encoded by these assemblies now exceeds 1 billion sequences. Meanwhile, a newly developed contig viewer provides fine-grained visualisation of the assembled contigs and their enriched annotations.

Tara Oceans: towards global ocean ecosystems biology

A planetary-scale understanding of the ocean ecosystem, particularly in light of climate change, is crucial. Here, we review the work of Tara Oceans, an international, multidisciplinary project to assess the complexity of ocean life across comprehensive taxonomic and spatial scales. Using a modified sailing boat, the team sampled plankton at 210 globally distributed sites at depths down to 1,000 m. We describe publicly available resources of molecular, morphological and environmental data, and discuss how an ecosystems biology approach has expanded our understanding of plankton diversity and ecology in the ocean as a planetary, interconnected ecosystem. These efforts illustrate how global-scale concepts and data can help to integrate biological complexity into models and serve as a baseline for assessing ecosystem changes and the future habitability of our planet in the Anthropocene epoch.

Mixotrophy in marine picocyanobacteria: use of organic compounds by Prochlorococcus and Synechococcus

Marine picocyanobacteria of the Prochlorococcus and Synechococcus genera have been longtime considered as autotrophic organisms. However, compelling evidence published over the last 15 years shows that these organisms can use different organic compounds containing key elements to survive in oligotrophic oceans, such as N (amino acids, amino sugars), S (dimethylsulfoniopropionate, DMSP), or P (ATP). Furthermore, marine picocyanobacteria can also take up glucose and use it as a source of carbon and energy, despite the fact that this compound is devoid of limiting elements and can also be synthesized by using standard metabolic pathways. This review will outline the main findings suggesting mixotrophy in the marine picocyanobacteria Prochlorococcus and Synechococcus, and its ecological relevance for these important primary producers.

Disentangling the mechanisms shaping the surface ocean microbiota

BACKGROUND

The ocean microbiota modulates global biogeochemical cycles and changes in its configuration may have large-scale consequences. Yet, the underlying ecological mechanisms structuring it are unclear. Here, we investigate how fundamental ecological mechanisms (selection, dispersal and ecological drift) shape the smallest members of the tropical and subtropical surface-ocean microbiota: prokaryotes and minute eukaryotes (picoeukaryotes). Furthermore, we investigate the agents exerting abiotic selection on this assemblage as well as the spatial patterns emerging from the action of ecological mechanisms. To explore this, we analysed the composition of surface-ocean prokaryotic and picoeukaryotic communities using DNA-sequence data (16S- and 18S-rRNA genes) collected during the circumglobal expeditions Malaspina-2010 and TARA-Oceans.

RESULTS

We found that the two main components of the tropical and subtropical surface-ocean microbiota, prokaryotes and picoeukaryotes, appear to be structured by different ecological mechanisms. Picoeukaryotic communities were predominantly structured by dispersal-limitation, while prokaryotic counterparts appeared to be shaped by the combined action of dispersal-limitation, selection and drift. Temperature-driven selection appeared as a major factor, out of a few selected factors, influencing species co-occurrence networks in prokaryotes but not in picoeukaryotes, indicating that association patterns may contribute to understand ocean microbiota structure and response to selection. Other measured abiotic variables seemed to have limited selective effects on community structure in the tropical and subtropical ocean. Picoeukaryotes displayed a higher spatial differentiation between communities and a higher distance decay when compared to prokaryotes, consistent with a scenario of higher dispersal limitation in the former after considering environmental heterogeneity. Lastly, random dynamics or drift seemed to have a more important role in structuring prokaryotic communities than picoeukaryotic counterparts.

CONCLUSIONS

The differential action of ecological mechanisms seems to cause contrasting biogeography, in the tropical and subtropical ocean, among the smallest surface plankton, prokaryotes and picoeukaryotes. This suggests that the idiosyncrasy of the main constituents of the ocean microbiota should be considered in order to understand its current and future configuration, which is especially relevant in a context of global change, where the reaction of surface ocean plankton to temperature increase is still unclear. Video Abstract.

Enzymes to unravel bioproducts architecture

Enzymes are essential and ubiquitous biocatalysts involved in various metabolic pathways and used in many industrial processes. Here, we reframe enzymes not just as biocatalysts transforming bioproducts but also as sensitive probes for exploring the structure and composition of complex bioproducts, like meat tissue, dairy products and plant materials, in both food and non-food bioprocesses. This review details the global strategy and presents the most recent investigations to prepare and use enzymes as relevant probes, with a focus on glycoside-hydrolases involved in plant deconstruction and proteases and lipases involved in food digestion. First, to expand the enzyme repertoire to fit bioproduct complexity, novel enzymes are mined from biodiversity and can be artificially engineered. Enzymes are further characterized by exploring sequence/structure/dynamics/function relationships together with the environmental factors influencing enzyme interactions with their substrates. Then, the most advanced experimental and theoretical approaches developed for exploring bioproducts at various scales (from nanometer to millimeter) using active and inactive enzymes as probes are illustrated. Overall, combining multimodal and multiscale approaches brings a better understanding of native-form or transformed bioproduct architecture and composition, and paves the way to mainstream the use of enzymes as probes.

DNA barcoding of macrofauna act as a tool for assessing marine ecosystem

Nowadays, marine ecosystem monitoring and assessment are increasingly depending on variety of molecular tools. With these background, DNA barcoding play a key role in species identification with increasing speed and accuracy, and although the suitability for developing genetic tools like genomic AMBI (gAMBI). Presently we have submitted 13 benthic polychaete species using mtCOI to GenBank. Of these, nine species were newly submitted, and hence they act as a benchmark and reference organism for identifying respective polychaete species worldwide in the near future. Based on that, our study results tend to be helpful for motivating among the researcher in order to implementing the genomic AMBI (gAMBI).

Clades of huge phages from across Earth's ecosystems

Bacteriophages typically have small genomes1 and depend on their bacterial hosts for replication2. Here we sequenced DNA from diverse ecosystems and found hundreds of phage genomes with lengths of more than 200 kilobases (kb), including a genome of 735 kb, which is-to our knowledge-the largest phage genome to be described to date. Thirty-five genomes were manually curated to completion (circular and no gaps). Expanded genetic repertoires include diverse and previously undescribed CRISPR-Cas systems, transfer RNAs (tRNAs), tRNA synthetases, tRNA-modification enzymes, translation-initiation and elongation factors, and ribosomal proteins. The CRISPR-Cas systems of phages have the capacity to silence host transcription factors and translational genes, potentially as part of a larger interaction network that intercepts translation to redirect biosynthesis to phage-encoded functions. In addition, some phages may repurpose bacterial CRISPR-Cas systems to eliminate competing phages. We phylogenetically define the major clades of huge phages from human and other animal microbiomes, as well as from oceans, lakes, sediments, soils and the built environment. We conclude that the large gene inventories of huge phages reflect a conserved biological strategy, and that the phages are distributed across a broad bacterial host range and across Earth's ecosystems.

Holistic pelagic biodiversity monitoring of the Black Sea via eDNA metabarcoding approach: From bacteria to marine mammals

As the largest semi-closed marine ecosystem in the world, the Black Sea has been heavily affected by human activities for a long time. Describing the biodiversity of multi-trophic biota in pelagic zone of the Black Sea and identifying the dominant environmental factors are prerequisites for protecting the sustainability of ecosystems. However, up to now, the taxonomic and distributional information about the Black Sea biota is not clear. Here, we employed a Tree-of-Life metabarcoding to analyze the biodiversity of eight communities in the Black Sea, investigated their biogeographical distribution, and further analyzed the influence of biological and abiotic factors on biota on large scales. We found that, (1) Over 8900 OTUs were detected in the Black Sea, of which 630 species were identified, covering the holistic biota from single-celled (bacteria 5620 OTUs 141 species; algae 1096 OTUs 185 species; protozoa 546 OTUs 146 species) to multicellular organisms (invertebrate metazoans 150 OTUs 34 species; fishes 1369 OTUs 76 species; large marine mammals 39 OTUs 5 species). (2) Higher trophic organisms (fishes and large mammals) distributed more evenly in space than the lower (microorganisms, protozoa and invertebrates). For lower trophic organisms, the vertical stratification was more obvious than the horizontal stratification (vertical p < 0.02, horizontal p < 0.05). (3) The bottom trophic organisms (bacteria and algae) of the food web significantly affected the distribution and composition of the others through biological interactions (Mantel p < 0.05). (4) At the level of abiotic factors, the effect of local species sorting on the composition of communities was 15% higher than that of mass dispersal effect. For the first time, this study monitored and profiled the holistic biodiversity in the pelagic zone of the Black Sea, and provided technological advances and preliminary knowledge for the ongoing Black Sea ecosystem protection efforts.

Diatoms Are Selective Segregators in Global Ocean Planktonic Communities

Diatoms are key phytoplankton in the modern ocean that are involved in numerous biotic interactions, ranging from symbiosis to predation and viral infection, which have considerable effects on global biogeochemical cycles. However, despite recent large-scale studies of plankton, we are still lacking a comprehensive picture of the diversity of diatom biotic interactions in the marine microbial community. Through the ecological interpretation of both inferred microbial association networks and available knowledge on diatom interactions compiled in an open-access database, we propose an ecosystems approach for exploring diatom interactions in the ocean.

Diatoms are a major component of phytoplankton, believed to be responsible for around 20% of the annual primary production on Earth. As abundant and ubiquitous organisms, they are known to establish biotic interactions with many other members of plankton. Through analyses of cooccurrence networks derived from the Tara Oceans expedition that take into account both biotic and abiotic factors in shaping the spatial distributions of species, we show that only 13% of diatom pairwise associations are driven by environmental conditions; the vast majority are independent of abiotic factors. In contrast to most other plankton groups, on a global scale, diatoms display a much higher proportion of negative correlations with other organisms, particularly toward potential predators and parasites, suggesting that their biogeography is constrained by top-down pressure. Genus-level analyses indicate that abundant diatoms are not necessarily the most connected and that species-specific abundance distribution patterns lead to negative associations with other organisms. In order to move forward in the biological interpretation of cooccurrence networks, an open-access extensive literature survey of diatom biotic interactions was compiled, of which 18.5% were recovered in the computed network. This result reveals the extent of what likely remains to be discovered in the field of planktonic biotic interactions, even for one of the best-known organismal groups.

IMPORTANCE Diatoms are key phytoplankton in the modern ocean that are involved in numerous biotic interactions, ranging from symbiosis to predation and viral infection, which have considerable effects on global biogeochemical cycles. However, despite recent large-scale studies of plankton, we are still lacking a comprehensive picture of the diversity of diatom biotic interactions in the marine microbial community. Through the ecological interpretation of both inferred microbial association networks and available knowledge on diatom interactions compiled in an open-access database, we propose an ecosystems approach for exploring diatom interactions in the ocean.

Genome Resolved Biogeography of Mamiellales

Among marine phytoplankton, Mamiellales encompass several species from the genera Micromonas, Ostreococcus and Bathycoccus, which are important contributors to primary production. Previous studies based on single gene markers described their wide geographical distribution but led to discussion because of the uneven taxonomic resolution of the method. Here, we leverage genome sequences for six Mamiellales species, two from each genus Micromonas, Ostreococcus and Bathycoccus, to investigate their distribution across 133 stations sampled during the Tara Oceans expedition. Our study confirms the cosmopolitan distribution of Mamiellales and further suggests non-random distribution of species, with two triplets of co-occurring genomes associated with different temperatures: Ostreococcuslucimarinus, Bathycoccusprasinos and Micromonaspusilla were found in colder waters, whereas Ostreococcus spp. RCC809, Bathycoccus spp. TOSAG39-1 and Micromonascommoda were more abundant in warmer conditions. We also report the distribution of the two candidate mating-types of Ostreococcus for which the frequency of sexual reproduction was previously assumed to be very low. Indeed, both mating types were systematically detected together in agreement with either frequent sexual reproduction or the high prevalence of a diploid stage. Altogether, these analyses provide novel insights into Mamiellales' biogeography and raise novel testable hypotheses about their life cycle and ecology.

Phytoplankton in the Tara Ocean

Photosynthesis evolved in the ocean more than 2 billion years ago and is now performed by a wide range of evolutionarily distinct organisms, including both prokaryotes and eukaryotes. Our appreciation of their abundance, distributions, and contributions to primary production in the ocean has been increasing since they were first discovered in the seventeenth century and has now been enhanced by data emerging from the Tara Oceans project, which performed a comprehensive worldwide sampling of plankton in the upper layers of the ocean between 2009 and 2013. Largely using recent data from Tara Oceans, here we review the geographic distributions of phytoplankton in the global ocean and their diversity, abundance, and standing stock biomass. We also discuss how omics-based information can be incorporated into studies of photosynthesis in the ocean and show the likely importance of mixotrophs and photosymbionts.

Diversity and biogeography of SAR11 bacteria from the Arctic Ocean

The Arctic Ocean is relatively isolated from other oceans and consists of strongly stratified water masses with distinct histories, nutrient, temperature, and salinity characteristics, therefore providing an optimal environment to investigate local adaptation. The globally distributed SAR11 bacterial group consists of multiple ecotypes that are associated with particular marine environments, yet relatively little is known about Arctic SAR11 diversity. Here, we examined SAR11 diversity using ITS analysis and metagenome-assembled genomes (MAGs). Arctic SAR11 assemblages were comprised of the S1a, S1b, S2, and S3 clades, and structured by water mass and depth. The fresher surface layer was dominated by an ecotype (S3-derived P3.2) previously associated with Arctic and brackish water. In contrast, deeper waters of Pacific origin were dominated by the P2.3 ecotype of the S2 clade, within which we identified a novel subdivision (P2.3s1) that was rare outside the Arctic Ocean. Arctic S2-derived SAR11 MAGs were restricted to high latitudes and included MAGs related to the recently defined S2b subclade, a finding consistent with bi-polar ecotypes and Arctic endemism. These results place the stratified Arctic Ocean into the SAR11 global biogeography and have identified SAR11 lineages for future investigation of adaptive evolution in the Arctic Ocean.

Into the bloom: Molecular response of pelagic tunicates to fluctuating food availability

The planktonic tunicates appendicularians and thaliaceans are highly efficient filter feeders on a wide range of prey size including bacteria and have shorter generation times than any other marine grazers. These traits allow some tunicate species to reach high population densities and ensure their success in a favourable environment. However, there are still few studies focusing on which genes and gene pathways are associated with responses of pelagic tunicates to environmental variability. Herein, we present the effect of food availability increase on tunicate community and gene expression at the Marquesas Islands (South-East Pacific Ocean). By using data from the Tara Oceans expedition, we show that changes in phytoplankton density and composition trigger the success of a dominant larvacean species (an undescribed appendicularian). Transcriptional signature to the autotroph bloom suggests key functions in specific physiological processes, i.e., energy metabolism, muscle contraction, membrane trafficking, and proteostasis. The relative abundance of reverse transcription-related Pfams was lower at bloom conditions, suggesting a link with adaptive genetic diversity in tunicates in natural ecosystems. Downstream of the bloom, pelagic tunicates were outcompeted by copepods. Our work represents the first metaomics study of the biological effects of phytoplankton bloom on a key zooplankton taxon.

HOPS: automated detection and authentication of pathogen DNA in archaeological remains

High-throughput DNA sequencing enables large-scale metagenomic analyses of complex biological systems. Such analyses are not restricted to present-day samples and can also be applied to molecular data from archaeological remains. Investigations of ancient microbes can provide valuable information on past bacterial commensals and pathogens, but their molecular detection remains a challenge. Here, we present HOPS (Heuristic Operations for Pathogen Screening), an automated bacterial screening pipeline for ancient DNA sequences that provides detailed information on species identification and authenticity. HOPS is a versatile tool for high-throughput screening of DNA from archaeological material to identify candidates for genome-level analyses.

Reconstruction of protein domain evolution using single-cell amplified genomes of uncultured choanoflagellates sheds light on the origin of animals

Understanding the origins of animal multicellularity is a fundamental biological question. Recent genome data have unravelled the role that co-option of pre-existing genes played in the origin of animals. However, there were also some important genetic novelties at the onset of Metazoa. To have a clear understanding of the specific genetic innovations and how they appeared, we need the broadest taxon sampling possible, especially among early-branching animals and their unicellular relatives. Here, we take advantage of single-cell genomics to expand our understanding of the genomic diversity of choanoflagellates, the sister-group to animals. With these genomes, we have performed an updated and taxon-rich reconstruction of protein evolution from the Last Eukaryotic Common Ancestor (LECA) to animals. Our novel data re-defines the origin of some genes previously thought to be metazoan-specific, like the POU transcription factor, which we show appeared earlier in evolution. Moreover, our data indicate that the acquisition of new genes at the stem of Metazoa was mainly driven by duplications and protein domain rearrangement processes at the stem of Metazoa. Furthermore, our analysis allowed us to reveal protein domains that are essential to the maintenance of animal multicellularity. Our analyses also demonstrate the utility of single-cell genomics from uncultured taxa to address evolutionary questions. This article is part of a discussion meeting issue 'Single cell ecology'.

Ecological and functional capabilities of an uncultured Kordia sp

Cultivable bacteria represent only a fraction of the diversity in microbial communities. However, the official procedures for classification and characterization of a novel prokaryotic species still rely on isolates. Nevertheless, due to single cell genomics, it is possible to retrieve genomes from environmental samples by sequencing them individually, and to assign specific genes to a specific taxon, regardless of their ability to grow in culture. In this study, a complete description was performed for uncultured Kordia sp. TARA_039_SRF, a proposed novel species within the genus Kordia, using culture-independent techniques. The type material was a high-quality draft genome (94.97% complete, 4.65% gene redundancy) co-assembled using ten nearly identical single amplified genomes (SAGs) from surface seawater in the North Indian Ocean during the Tara Oceans Expedition. The assembly process was optimized to obtain the best possible assembly metrics and a less fragmented genome. The closest relative of the species was Kordia periserrulae, which shared 97.56% similarity of the 16S rRNA gene, 75% orthologs and 89.13% average nucleotide identity. The functional potential of the proposed novel species included proteorhodopsin, the ability to incorporate nitrate, cytochrome oxidases with high affinity for oxygen, and CAZymes that were unique features within the genus. Its abundance at different depths and size fractions was also evaluated together with its functional annotation, revealing that its putative ecological niche could be particles of phytoplanktonic origin. It could putatively attach to these particles and consume them while sinking to the deeper and oxygen depleted layers of the North Indian Ocean.

Gene Expression Changes and Community Turnover Differentially Shape the Global Ocean Metatranscriptome

Ocean microbial communities strongly influence the biogeochemistry, food webs, and climate of our planet. Despite recent advances in understanding their taxonomic and genomic compositions, little is known about how their transcriptomes vary globally. Here, we present a dataset of 187 metatranscriptomes and 370 metagenomes from 126 globally distributed sampling stations and establish a resource of 47 million genes to study community-level transcriptomes across depth layers from pole-to-pole. We examine gene expression changes and community turnover as the underlying mechanisms shaping community transcriptomes along these axes of environmental variation and show how their individual contributions differ for multiple biogeochemically relevant processes. Furthermore, we find the relative contribution of gene expression changes to be significantly lower in polar than in non-polar waters and hypothesize that in polar regions, alterations in community activity in response to ocean warming will be driven more strongly by changes in organismal composition than by gene regulatory mechanisms. VIDEO ABSTRACT.

Fueling ab initio folding with marine metagenomics enables structure and function predictions of new protein families

INTRODUCTION

The ocean microbiome represents one of the largest microbiomes and produces nearly half of the primary energy on the planet through photosynthesis or chemosynthesis. Using recent advances in marine genomics, we explore new applications of oceanic metagenomes for protein structure and function prediction.

RESULTS

By processing 1.3 TB of high-quality reads from the Tara Oceans data, we obtain 97 million non-redundant genes. Of the 5721 Pfam families that lack experimental structures, 2801 have at least one member associated with the oceanic metagenomics dataset. We apply C-QUARK, a deep-learning contact-guided ab initio structure prediction pipeline, to model 27 families, where 20 are predicted to have a reliable fold with estimated template modeling score (TM-score) at least 0.5. Detailed analyses reveal that the abundance of microbial genera in the ocean is highly correlated to the frequency of occurrence in the modeled Pfam families, suggesting the significant role of the Tara Oceans genomes in the contact-map prediction and subsequent ab initio folding simulations. Of interesting note, PF15461, which has a majority of members coming from ocean-related bacteria, is identified as an important photosynthetic protein by structure-based function annotations. The pipeline is extended to a set of 417 Pfam families, built on the combination of Tara with other metagenomics datasets, which results in 235 families with an estimated TM-score over 0.5.

CONCLUSIONS

These results demonstrate a new avenue to improve the capacity of protein structure and function modeling through marine metagenomics, especially for difficult proteins with few homologous sequences.

Density-dependent oxylipin production in natural diatom communities: possible implications for plankton dynamics

Oxylipins are important signal transduction lipoxygenase-derived products of fatty acids that regulate a variety of physiological and pathological processes in plants and animals. In marine diatoms, these molecules can be highly bioactive, impacting zooplankton grazers, bacteria and other phytoplankton. However, the ultimate cause for oxylipin production in diatoms is still poorly understood, from an evolutionary perspective. Here we analysed production of particulate linear oxygenated fatty acids (LOFAs, previously named non-volatile oxylipins) from natural phytoplankton collected weekly for 1 year. We demonstrate for the first time that diatoms are the main LOFA producers in natural phytoplankton assemblages. Interestingly, LOFA-per-cell production decreased with increasing diatom density and was not due to major changes in diatom community composition. An inverse relation was confirmed at a global scale by analysing diatom lipoxygenase unigenes and metagenomes from Tara Oceans datasets. A network analysis suggested that different LOFAs could contribute to modulate co-variations of different diatom taxa. Overall, we offer new insights in diatom chemical ecology, possibly explaining the evolution of oxylipin synthesis in diatoms.

Single cell ecogenomics reveals mating types of individual cells and ssDNA viral infections in the smallest photosynthetic eukaryotes

Planktonic photosynthetic organisms of the class Mamiellophyceae include the smallest eukaryotes (less than 2 µm), are globally distributed and form the basis of coastal marine ecosystems. Eight complete fully annotated 13-22 Mb genomes from three genera, Ostreococcus, Bathycoccus and Micromonas, are available from previously isolated clonal cultured strains and provide an ideal resource to explore the scope and challenges of analysing single cell amplified genomes (SAGs) isolated from a natural environment. We assembled data from 12 SAGs sampled during the Tara Oceans expedition to gain biological insights about their in situ ecology, which might be lost by isolation and strain culture. Although the assembled nuclear genomes were incomplete, they were large enough to infer the mating types of four Ostreococcus SAGs. The systematic occurrence of sequences from the mitochondria and chloroplast, representing less than 3% of the total cell's DNA, intimates that SAGs provide suitable substrates for detection of non-target sequences, such as those of virions. Analysis of the non-Mamiellophyceae assemblies, following filtering out cross-contaminations during the sequencing process, revealed two novel 1.6 and 1.8 kb circular DNA viruses, and the presence of specific Bacterial and Oomycete sequences suggests that these organisms might co-occur with the Mamiellales. This article is part of a discussion meeting issue 'Single cell ecology'.

The Tara Pacific expedition-A pan-ecosystemic approach of the "-omics" complexity of coral reef holobionts across the Pacific Ocean

Coral reefs are the most diverse habitats in the marine realm. Their productivity, structural complexity, and biodiversity critically depend on ecosystem services provided by corals that are threatened because of climate change effects-in particular, ocean warming and acidification. The coral holobiont is composed of the coral animal host, endosymbiotic dinoflagellates, associated viruses, bacteria, and other microeukaryotes. In particular, the mandatory photosymbiosis with microalgae of the family Symbiodiniaceae and its consequences on the evolution, physiology, and stress resilience of the coral holobiont have yet to be fully elucidated. The functioning of the holobiont as a whole is largely unknown, although bacteria and viruses are presumed to play roles in metabolic interactions, immunity, and stress tolerance. In the context of climate change and anthropogenic threats on coral reef ecosystems, the Tara Pacific project aims to provide a baseline of the "-omics" complexity of the coral holobiont and its ecosystem across the Pacific Ocean and for various oceanographically distinct defined areas. Inspired by the previous Tara Oceans expeditions, the Tara Pacific expedition (2016-2018) has applied a pan-ecosystemic approach on coral reefs throughout the Pacific Ocean, drawing an east-west transect from Panama to Papua New Guinea and a south-north transect from Australia to Japan, sampling corals throughout 32 island systems with local replicates. Tara Pacific has developed and applied state-of-the-art technologies in very-high-throughput genetic sequencing and molecular analysis to reveal the entire microbial and chemical diversity as well as functional traits associated with coral holobionts, together with various measures on environmental forcing. This ambitious project aims at revealing a massive amount of novel biodiversity, shedding light on the complex links between genomes, transcriptomes, metabolomes, organisms, and ecosystem functions in coral reefs and providing a reference of the biological state of modern coral reefs in the Anthropocene.