Transcriptome reconstruction and functional analysis of eukaryotic marine plankton communities via high-throughput metagenomics and metatranscriptomics

DeepMicroClass sorts metagenomic contigs into prokaryotes, eukaryotes and viruses

Sequence classification facilitates a fundamental understanding of the structure of microbial communities. Binary metagenomic sequence classifiers are insufficient because environmental metagenomes are typically derived from multiple sequence sources. Here we introduce a deep-learning based sequence classifier, DeepMicroClass, that classifies metagenomic contigs into five sequence classes, i.e. viruses infecting prokaryotic or eukaryotic hosts, eukaryotic or prokaryotic chromosomes, and prokaryotic plasmids. DeepMicroClass achieved high performance for all sequence classes at various tested sequence lengths ranging from 500 bp to 100 kbps. By benchmarking on a synthetic dataset with variable sequence class composition, we showed that DeepMicroClass obtained better performance for eukaryotic, plasmid and viral contig classification than other state-of-the-art predictors. DeepMicroClass achieved comparable performance on viral sequence classification with geNomad and VirSorter2 when benchmarked on the CAMI II marine dataset. Using a coastal daily time-series metagenomic dataset as a case study, we showed that microbial eukaryotes and prokaryotic viruses are integral to microbial communities. By analyzing monthly metagenomes collected at HOT and BATS, we found relatively higher viral read proportions in the subsurface layer in late summer, consistent with the seasonal viral infection patterns prevalent in these areas. We expect DeepMicroClass will promote metagenomic studies of under-appreciated sequence types.

DMSOP-cleaving enzymes are diverse and widely distributed in marine microorganisms

Dimethylsulfoxonium propionate (DMSOP) is a recently identified and abundant marine organosulfur compound with roles in oxidative stress protection, global carbon and sulfur cycling and, as shown here, potentially in osmotolerance. Microbial DMSOP cleavage yields dimethyl sulfoxide, a ubiquitous marine metabolite, and acrylate, but the enzymes responsible, and their environmental importance, were unknown. Here we report DMSOP cleavage mechanisms in diverse heterotrophic bacteria, fungi and phototrophic algae not previously known to have this activity, and highlight the unappreciated importance of this process in marine sediment environments. These diverse organisms, including Roseobacter, SAR11 bacteria and Emiliania huxleyi, utilized their dimethylsulfoniopropionate lyase 'Ddd' or 'Alma' enzymes to cleave DMSOP via similar catalytic mechanisms to those for dimethylsulfoniopropionate. Given the annual teragram predictions for DMSOP production and its prevalence in marine sediments, our results highlight that DMSOP cleavage is likely a globally significant process influencing carbon and sulfur fluxes and ecological interactions.

Environmental transcriptomics under heat stress: Can environmental RNA reveal changes in gene expression of aquatic organisms?

To safeguard biodiversity in a changing climate, taxonomic information about species turnover and insights into the health of organisms are required. Environmental DNA approaches are increasingly used for species identification, but cannot provide functional insights. Transcriptomic methods reveal the physiological states of macroorganisms, but are currently species-specific and require tissue sampling or animal sacrifice, making community-wide assessments challenging. Here, we test whether broad functional information (expression level of the transcribed genes) can be harnessed from environmental RNA (eRNA), which includes extra-organismal RNA from macroorganisms along with whole microorganisms. We exposed Daphnia pulex as well as phytoplankton prey and microorganism colonizers to control (20°C) and heat stress (28°C) conditions for 7 days. We sequenced eRNA from tank water (after complete removal of Daphnia) as well as RNA from Daphnia tissue, enabling comparisons of extra-organismal and organismal RNA-based gene expression profiles. Both RNA types detected similar heat stress responses of Daphnia. Using eRNA, we identified 32 Daphnia genes to be differentially expressed following heat stress. Of these, 17 were also differentially expressed and exhibited similar levels of relative expression in organismal RNA. In addition to the extra-organismal Daphnia response, eRNA detected community-wide heat stress responses consisting of distinct functional profiles and 121 differentially expressed genes across eight taxa. Our study demonstrates that environmental transcriptomics based on extra-organismal eRNA can noninvasively reveal gene expression responses of macroorganisms following environmental changes, with broad potential implications for the biomonitoring of health across the trophic chain.

Holistic view of the seascape dynamics and environment impact on macro-scale genetic connectivity of marine plankton populations

BACKGROUND

Plankton seascape genomics studies have revealed different trends from large-scale weak differentiation to microscale structures. Previous studies have underlined the influence of the environment and seascape on species differentiation and adaptation. However, these studies have generally focused on a few single species, sparse molecular markers, or local scales. Here, we investigated the genomic differentiation of plankton at the macro-scale in a holistic approach using Tara Oceans metagenomic data together with a reference-free computational method.

RESULTS

We reconstructed the FST-based genomic differentiation of 113 marine planktonic taxa occurring in the North and South Atlantic Oceans, Southern Ocean, and Mediterranean Sea. These taxa belong to various taxonomic clades spanning Metazoa, Chromista, Chlorophyta, Bacteria, and viruses. Globally, population genetic connectivity was significantly higher within oceanic basins and lower in bacteria and unicellular eukaryotes than in zooplankton. Using mixed linear models, we tested six abiotic factors influencing connectivity, including Lagrangian travel time, as proxies of oceanic current effects. We found that oceanic currents were the main population genetic connectivity drivers, together with temperature and salinity. Finally, we classified the 113 taxa into parameter-driven groups and showed that plankton taxa belonging to the same taxonomic rank such as phylum, class or order presented genomic differentiation driven by different environmental factors.

CONCLUSION

Our results validate the isolation-by-current hypothesis for a non-negligible proportion of taxa and highlight the role of other physicochemical parameters in large-scale plankton genetic connectivity. The reference-free approach used in this study offers a new systematic framework to analyse the population genomics of non-model and undocumented marine organisms from a large-scale and holistic point of view.

Expanding the range of the respiratory infectome in Australian feedlot cattle with and without respiratory disease using metatranscriptomics

BACKGROUND

Bovine respiratory disease (BRD) is one of the most common diseases in intensively managed cattle, often resulting in high morbidity and mortality. Although several pathogens have been isolated and extensively studied, the complete infectome of the respiratory complex consists of a more extensive range unrecognised species. Here, we used total RNA sequencing (i.e., metatranscriptomics) of nasal and nasopharyngeal swabs collected from animals with and without BRD from two cattle feedlots in Australia.

RESULTS

A high abundance of bovine nidovirus, influenza D, bovine rhinitis A and bovine coronavirus was found in the samples. Additionally, we obtained the complete or near-complete genome of bovine rhinitis B, enterovirus E1, bovine viral diarrhea virus (sub-genotypes 1a and 1c) and bovine respiratory syncytial virus, and partial sequences of other viruses. A new species of paramyxovirus was also identified. Overall, the most abundant RNA virus, was the bovine nidovirus. Characterisation of bacterial species from the transcriptome revealed a high abundance and diversity of Mollicutes in BRD cases and unaffected control animals. Of the non-Mollicutes species, Histophilus somni was detected, whereas there was a low abundance of Mannheimia haemolytica.

CONCLUSION

This study highlights the use of untargeted sequencing approaches to study the unrecognised range of microorganisms present in healthy or diseased animals and the need to study previously uncultured viral species that may have an important role in cattle respiratory disease. Video Abstract.

Molecular discoveries in microbial DMSP synthesis

Dimethylsulfoniopropionate (DMSP) is one of the Earth's most abundant organosulfur compounds because many marine algae, bacteria, corals and some plants produce it to high mM intracellular concentrations. In these organisms, DMSP acts an anti-stress molecule with purported roles to protect against salinity, temperature, oxidative stress and hydrostatic pressure, amongst many other reported functions. However, DMSP is best known for being a major precursor of the climate-active gases and signalling molecules dimethylsulfide (DMS), methanethiol (MeSH) and, potentially, methane, through microbial DMSP catabolism. DMSP catabolism has been extensively studied and the microbes, pathways and enzymes involved have largely been elucidated through the application of molecular research over the last 17 years. In contrast, the molecular biology of DMSP synthesis is a much newer field, with the first DMSP synthesis enzymes only being identified in the last 5 years. In this review, we discuss how the elucidation of key DMSP synthesis enzymes has greatly expanded our knowledge of the diversity of DMSP-producing organisms, the pathways used, and what environmental factors regulate production, as well as to inform on the physiological roles of DMSP. Importantly, the identification of key DMSP synthesis enzymes in the major groups of DMSP producers has allowed scientists to study the distribution and predict the importance of different DMSP-producing organisms to global DMSP production in diverse marine and sediment environments. Finally, we highlight key challenges for future molecular research into DMSP synthesis that need addressing to better understand the cycling of this important marine organosulfur compound, and its magnitude in the environment.

Reverse engineering environmental metatranscriptomes clarifies best practices for eukaryotic assembly

BACKGROUND

Diverse communities of microbial eukaryotes in the global ocean provide a variety of essential ecosystem services, from primary production and carbon flow through trophic transfer to cooperation via symbioses. Increasingly, these communities are being understood through the lens of omics tools, which enable high-throughput processing of diverse communities. Metatranscriptomics offers an understanding of near real-time gene expression in microbial eukaryotic communities, providing a window into community metabolic activity.

RESULTS

Here we present a workflow for eukaryotic metatranscriptome assembly, and validate the ability of the pipeline to recapitulate real and manufactured eukaryotic community-level expression data. We also include an open-source tool for simulating environmental metatranscriptomes for testing and validation purposes. We reanalyze previously published metatranscriptomic datasets using our metatranscriptome analysis approach.

CONCLUSION

We determined that a multi-assembler approach improves eukaryotic metatranscriptome assembly based on recapitulated taxonomic and functional annotations from an in-silico mock community. The systematic validation of metatranscriptome assembly and annotation methods provided here is a necessary step to assess the fidelity of our community composition measurements and functional content assignments from eukaryotic metatranscriptomes.

Environmental DNA sequencing reveals the regional difference in diversity and community assembly mechanisms of eukaryotic plankton in coastal waters

The diversity and community assembly mechanisms of eukaryotic plankton in coastal waters is so far not clear. In this study, we selected the coastal waters of Guangdong-Hong Kong-Macao Greater Bay Area, which is a highly developed region in China, as the research area. By use of high-throughput sequencing technologies, the diversity and community assembly mechanisms of eukaryotic marine plankton were studied in which a total of 7,295 OTUs were obtained, and 2,307 species were annotated by doing environmental DNA survey of 17 sites consist of surface and bottom layer. Ultimately, the analysis reveals that the species abundance of bottom layer is, by and large, higher than that in the surface layer. In the bottom, Arthropoda is the first largest group, accounting for more than 20% while Arthropoda and Bacillariophyta are dominant groups in surface waters accounting for more than 40%. It is significant of the variance in alpha-diversity between sampling sites, and the difference of alpha-diversity between bottom sites is greater than that of surface sites. The result suggests that the environmental factors that have significant influence on alpha-diversity are total alkalinity and offshore distance for surface sites, and water depth and turbidity for bottom sites. Likewise, the plankton communities obey the typical distance-decay pattern. Analysis about community assembly mechanisms reveals that, overall, dispersal limitation is the major pattern of community formation, which accounts for more than 83% of the community formation processes, suggesting that stochastic processes are the crucial assembly mechanism of the eukaryotic plankton community in the study area.

Primer design for the amplification of the ammonium transporter genes from the uncultured haptophyte algal species symbiotic with the marine nitrogen-fixing cyanobacterium UCYN-A1

The multiple symbiotic partnerships between closely related species of the haptophyte algae Braarudosphaera bigelowii and the nitrogen-fixing cyanobacteria Candidatus Atelocyanobacterium thalassa (UCYN-A) contribute importantly to the nitrogen and carbon cycles in vast areas of the ocean. The diversity of the eukaryotic 18S rDNA phylogenetic gene marker has helped to identify some of these symbiotic haptophyte species, yet we still lack a genetic marker to assess its diversity at a finer scale. One of such genes is the ammonium transporter (amt) gene, which encodes the protein that might be involved in the uptake of ammonium from UCYN-A in these symbiotic haptophytes. Here, we designed three specific PCR primer sets targeting the amt gene of the haptophyte species (A1-Host) symbiotic with the open ocean UCYN-A1 sublineage, and tested them in samples collected from open ocean and near-shore environments. Regardless of the primer pair used at Station ALOHA, which is where UCYN-A1 is the pre-dominant UCYN-A sublineage, the most abundant amt amplicon sequence variant (ASV) was taxonomically classified as A1-Host. In addition, two out of the three PCR primer sets revealed the existence of closely-related divergent haptophyte amt ASVs (>95% nucleotide identity). These divergent amt ASVs had higher relative abundances than the haptophyte typically associated with UCYN-A1 in the Bering Sea, or co-occurred with the previously identified A1-Host in the Coral Sea, suggesting the presence of new diversity of closely-related A1-Hosts in polar and temperate waters. Therefore, our study reveals an overlooked diversity of haptophytes species with distinct biogeographic distributions partnering with UCYN-A, and provides new primers that will help to gain new knowledge of the UCYN-A/haptophyte symbiosis.

Light-driven processes: key players of the functional biodiversity in microalgae

Microalgae are prominent aquatic organisms, responsible for about half of the photosynthetic activity on Earth. Over the past two decades, breakthroughs in genomics and ecosystem biology, as well as the development of genetic resources in model species, have redrawn the boundaries of our knowledge on the relevance of these microbes in global ecosystems. However, considering their vast biodiversity and complex evolutionary history, our comprehension of algal biology remains limited. As algae rely on light, both as their main source of energy and for information about their environment, we focus here on photosynthesis, photoperception, and chloroplast biogenesis in the green alga Chlamydomonas reinhardtii and marine diatoms. We describe how the studies of light-driven processes are key to assessing functional biodiversity in evolutionary distant microalgae. We also emphasize that integration of laboratory and environmental studies, and dialogues between different scientific communities are both timely and essential to understand the life of phototrophs in complex ecosystems and to properly assess the consequences of environmental changes on aquatic environments globally.

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.

Strong population genomic structure of the toxic dinoflagellate Alexandrium minutum inferred from meta-transcriptome samples

Despite theoretical expectations, marine microeukaryote population are often highly structured and the mechanisms behind such patterns remain to be elucidated. These organisms display huge census population sizes, yet genotyping usually requires clonal strains originating from single cells, hindering proper population sampling. Estimating allelic frequency directly from population wide samples, without any isolation step, offers an interesting alternative. Here, we validate the use of meta-transcriptome environmental samples to determine the population genetic structure of the dinoflagellate Alexandrium minutum. Strain and meta-transcriptome based results both indicated a strong genetic structure for A. minutum in Western Europe, to the level expected between cryptic species. The presence of numerous private alleles, and even fixed polymorphism, would indicate ancient divergence and absence of gene flow between populations. Single nucleotide polymorphisms (SNPs) displaying strong allele frequency differences were distributed throughout the genome, which might indicate pervasive selection from standing genetic variation (soft selective sweeps). However, a few genomic regions displayed extremely low diversity that could result from the fixation of adaptive de novo mutations (hard selective sweeps) within the populations.

SEQ2MGS: an effective tool for generating realistic artificial metagenomes from the existing sequencing data

Assessment of bioinformatics tools for the metagenomics analysis from the whole genome sequencing data requires realistic benchmark sets. We developed an effective and simple generator of artificial metagenomes from real sequencing experiments. The tool (SEQ2MGS) analyzes the input FASTQ files, precomputes genomic content, and blends shotgun reads from different sequenced isolates, or spike isolate(s) in real metagenome, in desired proportions. SEQ2MGS eliminates the need for simulation of sequencing platform variations, reads distributions, presence of plasmids, viruses, and contamination. The tool is especially useful for a quick generation of multiple complex samples that include new or understudied organisms, even without assembled genomes. For illustration, we first demonstrated the ease of SEQ2MGS use for the simulation of altered Schaedler flora (ASF) in comparison with de novo metagenomics generators Grinder and CAMISIM. Next, we emulated the emergence of a pathogen in the human gut microbiome and observed that Kraken, Centrifuge, and MetaPhlAn, while correctly identified Klebsiella pneumoniae, produced inconsistent results for the rest of real metagenome. Finally, using the MG-RAST platform, we affirmed that SEQ2MGS properly transfers genomic information from an isolate into the simulated metagenome by the correct identification of antimicrobial resistance genes anticipated to appear compared to the original metagenome.

The Ocean Gene Atlas v2.0: online exploration of the biogeography and phylogeny of plankton genes

Testing hypothesis about the biogeography of genes using large data resources such as Tara Oceans marine metagenomes and metatranscriptomes requires significant hardware resources and programming skills. The new release of the ‘Ocean Gene Atlas’ (OGA2) is a freely available intuitive online service to mine large and complex marine environmental genomic databases. OGA2 datasets available have been extended and now include, from the Tara Oceans portfolio: (i) eukaryotic Metagenome-Assembled-Genomes (MAGs) and Single-cell Assembled Genomes (SAGs) (10.2E+6 coding genes), (ii) version 2 of Ocean Microbial Reference Gene Catalogue (46.8E+6 non-redundant genes), (iii) 924 MetaGenomic Transcriptomes (7E+6 unigenes), (iv) 530 MAGs from an Arctic MAG catalogue (1E+6 genes) and (v) 1888 Bacterial and Archaeal Genomes (4.5E+6 genes), and an additional dataset from the Malaspina 2010 global circumnavigation: (vi) 317 Malaspina Deep Metagenome Assembled Genomes (0.9E+6 genes). Novel analyses enabled by OGA2 include phylogenetic tree inference to visualize user queries within their context of sequence homologues from both the marine environmental dataset and the RefSeq database. An Application Programming Interface (API) now allows users to query OGA2 using command-line tools, hence providing local workflow integration. Finally, gene abundance can be interactively filtered directly on map displays using any of the available environmental variables. Ocean Gene Atlas v2.0 is freely-available at: https://tara-oceans.mio.osupytheas.fr/ocean-gene-atlas/.

Development of transcriptomics-based growth rate indices in two model eukaryotes and relevance to metatranscriptomic datasets

Growth rate estimation is important to understand the flow of energy and nutrient elements in an ecosystem, but it has remained challenging, especially on microscopic organisms. In this study, we propose four growth rate indices that use mRNA abundance ratios between nuclear and mitochondrial genes: (1) total nuclear and mitochondrial mRNA ratio (Nuc:Mito‐TmRNA); (2) nuclear and mitochondrial ribosomal protein mRNA ratio (Nuc:Mito‐RPmRNA); (3) gene ontology (GO) terms and total mitochondrial mRNA ratios; and (4) nuclear and mitochondrial specific gene mRNA ratio. We examine these proposed ratios using RNA‐Seq datasets of Daphnia magna, and Saccharomyces cerevisiae retrieved from the NCBI Short Read Archive. The results showed that both Nuc:Mito‐TmRNA and Nuc:Mito‐RPmRNA ratio indices showed significant correlations with the growth rate for both species. A large number of GO terms mRNA ratios showed significant correlations with the growth rate of S. cerevisiae. Lastly, we identified mRNA ratios of several specific nuclear and mitochondrial gene pairs that showed significant correlations. We foresee future implications for the proposed mRNA ratios used in metatranscriptome analyses to estimate the growth rate of communities and species.

Functional repertoire convergence of distantly related eukaryotic plankton lineages abundant in the sunlit ocean

Marine planktonic eukaryotes play critical roles in global biogeochemical cycles and climate. However, their poor representation in culture collections limits our understanding of the evolutionary history and genomic underpinnings of planktonic ecosystems. Here, we used 280 billion Tara Oceans metagenomic reads from polar, temperate, and tropical sunlit oceans to reconstruct and manually curate more than 700 abundant and widespread eukaryotic environmental genomes ranging from 10 Mbp to 1.3 Gbp. This genomic resource covers a wide range of poorly characterized eukaryotic lineages that complement long-standing contributions from culture collections while better representing plankton in the upper layer of the oceans. We performed the first, to our knowledge, comprehensive genome-wide functional classification of abundant unicellular eukaryotic plankton, revealing four major groups connecting distantly related lineages. Neither trophic modes of plankton nor its vertical evolutionary history could completely explain the functional repertoire convergence of major eukaryotic lineages that coexisted within oceanic currents for millions of years.

Unifying the known and unknown microbial coding sequence space

Genes of unknown function are among the biggest challenges in molecular biology, especially in microbial systems, where 40–60% of the predicted genes are unknown. Despite previous attempts, systematic approaches to include the unknown fraction into analytical workflows are still lacking. Here, we present a conceptual framework, its translation into the computational workflow AGNOSTOS and a demonstration on how we can bridge the known-unknown gap in genomes and metagenomes. By analyzing 415,971,742 genes predicted from 1749 metagenomes and 28,941 bacterial and archaeal genomes, we quantify the extent of the unknown fraction, its diversity, and its relevance across multiple organisms and environments. The unknown sequence space is exceptionally diverse, phylogenetically more conserved than the known fraction and predominantly taxonomically restricted at the species level. From the 71 M genes identified to be of unknown function, we compiled a collection of 283,874 lineage-specific genes of unknown function for Cand. Patescibacteria (also known as Candidate Phyla Radiation, CPR), which provides a significant resource to expand our understanding of their unusual biology. Finally, by identifying a target gene of unknown function for antibiotic resistance, we demonstrate how we can enable the generation of hypotheses that can be used to augment experimental data.

It is estimated that scientists do not know what half of microbial genes actually do. When these genes are discovered in microorganisms grown in the lab or found in environmental samples, it is not possible to identify what their roles are. Many of these genes are excluded from further analyses for these reasons, meaning that the study of microbial genes tends to be limited to genes that have already been described.

These limitations hinder research into microbiology, because information from newly discovered genes cannot be integrated to better understand how these organisms work. Experiments to understand what role these genes have in the microorganisms are labor-intensive, so new analytical strategies are needed.

To do this, Vanni et al. developed a new framework to categorize genes with unknown roles, and a computational workflow to integrate them into traditional analyses. When this approach was applied to over 400 million microbial genes (both with known and unknown roles), it showed that the share of genes with unknown functions is only about 30 per cent, smaller than previously thought. The analysis also showed that these genes are very diverse, revealing a huge space for future research and potential applications. Combining their approach with experimental data, Vanni et al. were able to identify a gene with a previously unknown purpose that could be involved in antibiotic resistance.

This system could be useful for other scientists studying microorganisms to get a more complete view of microbial systems. In future, it may also be used to analyze the genetics of other organisms, such as plants and animals.

Evolution of Phytoplankton as Estimated from Genetic Diversity

Phytoplankton are photosynthetic, single-celled organisms producing almost half of all oxygen on Earth and play a central role as prey for higher organisms, making them irreplaceable in the marine food web. As Global Change proceeds, imposing rapidly intensifying selection pressures, phytoplankton are forced to undergo evolution, local extinction, or redistribution, with potentially cascading effects throughout the marine ecosystem. Recent results from the field of population genetics display high levels of standing genetic diversity in natural phytoplankton populations, providing ample ‘evolutionary options’ and implying high adaptive potential to changing conditions. This potential for adaptive evolution is realized in several studies of experimental evolution, even though most of these studies investigate the evolution of only single strains. This, however, shows that phytoplankton not only evolve from standing genetic diversity, but also rely on de novo mutations. Recent global sampling campaigns show that the immense intraspecific diversity of phytoplankton in the marine ecosystem has been significantly underestimated, meaning we are only studying a minor portion of the relevant variability in the context of Global Change and evolution. An increased understanding of genomic diversity is primarily hampered by the low number of ecologically representative reference genomes of eukaryotic phytoplankton and the functional annotation of these. However, emerging technologies relying on metagenome and transcriptome data may offer a more realistic understanding of phytoplankton diversity.

Spatiotemporal distribution of bacterial dimethylsulfoniopropionate producing and catabolic genes in the Changjiang Estuary

The osmolyte dimethylsulfoniopropionate (DMSP) is produced in petagram amounts by marine microorganisms. Estuaries provide natural gradients in salinity and nutrients, factors known to regulate DMSP production; yet there have been no molecular studies of DMSP production and cycling across these gradients. Here, we study the abundance, distribution and transcription of key DMSP synthesis (e.g. dsyB and mmtN) and catabolic (e.g. dddP and dmdA) genes along the salinity gradient of the Changjiang Estuary. DMSP levels did not correlate with Chl a across the salinity gradient. In contrast, DMSP concentration, abundance of bacterial DMSP producers and their dsyB and mmtN transcripts were lowest in the freshwater samples and increased abruptly with salinity in the transitional and seawater samples. Metagenomics analysis suggests bacterial DMSP-producers were more abundant than their algal equivalents and were more prominent in summer than winter samples. Bacterial DMSP catabolic genes and their transcripts followed the same trend of being greatly enhanced in transitional and seawater samples with higher DMSP levels than freshwater samples. DMSP cleavage was likely the dominant catabolic pathway, with DMSP lyase genes being ~4.3-fold more abundant than the demethylase gene dmdA. This is an exemplar study for future research on microbial DMSP cycling in estuary environments.

Metatranscriptomics for the Human Microbiome and Microbial Community Functional Profiling

Shotgun metatranscriptomics (MTX) is an increasingly practical way to survey microbial community gene function and regulation at scale. This review begins by summarizing the motivations for community transcriptomics and the history of the field. We then explore the principles, best practices, and challenges of contemporary MTX workflows: beginning with laboratory methods for isolation and sequencing of community RNA, followed by informatics methods for quantifying RNA features, and finally statistical methods for detecting differential expression in a community context. In thesecond half of the review, we survey important biological findings from the MTX literature, drawing examples from the human microbiome, other (nonhuman) host-associated microbiomes, and the environment. Across these examples, MTX methods prove invaluable for probing microbe-microbe and host-microbe interactions, the dynamics of energy harvest and chemical cycling, and responses to environmental stresses. We conclude with a review of open challenges in the MTX field, including making assays and analyses more robust, accessible, and adaptable to new technologies; deciphering roles for millions of uncharacterized microbial transcripts; and solving applied problems such as biomarker discovery and development of microbial therapeutics.

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.

Statistical approaches for differential expression analysis in metatranscriptomics

Motivation

Metatranscriptomics (MTX) has become an increasingly practical way to profile the functional activity of microbial communities in situ. However, MTX remains underutilized due to experimental and computational limitations. The latter are complicated by non-independent changes in both RNA transcript levels and their underlying genomic DNA copies (as microbes simultaneously change their overall abundance in the population and regulate individual transcripts), genetic plasticity (as whole loci are frequently gained and lost in microbial lineages) and measurement compositionality and zero-inflation. Here, we present a systematic evaluation of and recommendations for differential expression (DE) analysis in MTX.

Results

We designed and assessed six statistical models for DE discovery in MTX that incorporate different combinations of DNA and RNA normalization and assumptions about the underlying changes of gene copies or species abundance within communities. We evaluated these models on multiple simulated and real multi-omic datasets. Models adjusting transcripts relative to their encoding gene copies as a covariate were significantly more accurate in identifying DE from MTX in both simulated and real datasets. Moreover, we show that when paired DNA measurements (metagenomic data) are not available, models normalizing MTX measurements within-species while also adjusting for total-species RNA balance sensitivity, specificity and interpretability of DE detection, as does filtering likely technical zeros. The efficiency and accuracy of these models pave the way for more effective MTX-based DE discovery in microbial communities.

Availability and implementation

The analysis code and synthetic datasets used in this evaluation are available online at http://huttenhower.sph.harvard.edu/mtx2021.

Supplementary information

Supplementary data are available at Bioinformatics online.

Global distribution patterns of marine nitrogen-fixers by imaging and molecular methods

Nitrogen fixation has a critical role in marine primary production, yet our understanding of marine nitrogen-fixers (diazotrophs) is hindered by limited observations. Here, we report a quantitative image analysis pipeline combined with mapping of molecular markers for mining >2,000,000 images and >1300 metagenomes from surface, deep chlorophyll maximum and mesopelagic seawater samples across 6 size fractions (<0.2–2000 μm). We use this approach to characterise the diversity, abundance, biovolume and distribution of symbiotic, colony-forming and particle-associated diazotrophs at a global scale. We show that imaging and PCR-free molecular data are congruent. Sequence reads indicate diazotrophs are detected from the ultrasmall bacterioplankton (<0.2 μm) to mesoplankton (180–2000 μm) communities, while images predict numerous symbiotic and colony-forming diazotrophs (>20 µm). Using imaging and molecular data, we estimate that polyploidy can substantially affect gene abundances of symbiotic versus colony-forming diazotrophs. Our results support the canonical view that larger diazotrophs (>10 μm) dominate the tropical belts, while unicellular cyanobacterial and non-cyanobacterial diazotrophs are globally distributed in surface and mesopelagic layers. We describe co-occurring diazotrophic lineages of different lifestyles and identify high-density regions of diazotrophs in the global ocean. Overall, we provide an update of marine diazotroph biogeographical diversity and present a new bioimaging-bioinformatic workflow.

Nitrogen fixation by diazotrophs is critical for marine primary production. Using Tara Oceans datasets, this study combines a quantitative image analysis pipeline with metagenomic mining to provide an improved global overview of diazotroph abundance, diversity and distribution.

SPX-related genes regulate phosphorus homeostasis in the marine phytoplankton, Phaeodactylum tricornutum

Phosphorus (P) is an essential nutrient for marine phytoplankton. Maintaining intracellular P homeostasis against environmental P variability is critical for phytoplankton, but how they achieve this is poorly understood. Here we identify a SPX gene and investigate its role in Phaeodactylum tricornutum. SPX knockout led to significant increases in the expression of phosphate transporters, alkaline phosphatases (the P acquisition machinery) and phospholipid hydrolases (a mechanism to reduce P demand). These demonstrate that SPX is a negative regulator of both P uptake and P-stress responses. Furthermore, we show that SPX regulation of P uptake and metabolism involves a phosphate starvation response regulator (PHR) as an intermediate. Additionally, we find the SPX related genes exist and operate across the phytoplankton phylogenetic spectrum and in the global oceans, indicating its universal importance in marine phytoplankton. This study lays a foundation for better understanding phytoplankton adaptation to P variability in the future changing oceans.

De novo transcriptome assembly of the green alga Ankistrodesmus falcatus

Ankistrodesmus falcatus is a globally distributed freshwater chlorophyte that is a candidate for biofuel production, is used to study the effects of toxins on aquatic communities, and is used as food in zooplankton research. Each of these research fields is transitioning to genomic tools. We created a reference transcriptome for of A. falcatus using NextGen sequencing and de novo assembly methods including Trinity, Velvet-Oases, and EvidentialGene. The assembled transcriptome has a total of 17,997 contigs, an N50 value of 2,462, and a GC content of 64.8%. BUSCO analysis recovered 83.3% of total chlorophyte BUSCOs and 82.5% of the eukaryotic BUSCOs. A portion (7.9%) of these supposedly single-copy genes were found to have transcriptionally active, distinct duplicates. We annotated the assembly using the dammit annotation pipeline, resulting in putative functional annotation for 68.89% of the assembly. Using available rbcL sequences from 16 strains (10 species) of Ankistrodesmus, we constructed a neighbor-joining phylogeny to illustrate genetic distances of our A. falcatus strain to other members of the genus. This assembly will be valuable for researchers seeking to identify Ankistrodesmus sequences in metatranscriptomic and metagenomic field studies and in experiments where separating expression responses of zooplankton and their algal food sources through bioinformatics is important.

Environmental RNA: A Revolution in Ecological Resolution?

Current advancements in environmental RNA (eRNA) exploit its relatively fast turnover rate relative to environmental DNA (eDNA) to assess 'metabolically active' or temporally/spatially recent community diversity. However, this focus significantly underutilizes the trove of potential ecological information encrypted in eRNA. Here, we argue for pushing beyond current species-level eDNA detection capabilities by using eRNA to detect any organisms with unique eRNA profiles, potentially including different life-history stages, sexes, or even specific phenotypes within a species. We also discuss the future of eRNA as a means of assessing the physiological status of organisms and the ecological health of populations and communities, reflecting ecosystem-level conditions. We posit that eRNA has the potential to significantly improve the resolution of organism detection, biological monitoring, and biomonitoring applications in ecology.

DMSP synthesis genes distinguish two types of DMSP producer phenotypes

Dimethylsulfoniopropionate (DMSP) is an important organic carbon and sulfur source in the surface ocean that fuels microbial activity and significantly impacts Earth's climate. After three decades of research, the cellular role(s) of DMSP and environmental drivers of production remain enigmatic. Recent work suggests that cellular DMSP concentrations, and changes in these concentrations in response to environmental stressors, define two major groups of DMSP producers: high DMSP producers that contain ≥ 50 mM intracellular DMSP and low DMSP producers that contain < 50 mM. Here we show that two recently described DMSP synthesis genes (DSYB and TpMT2) may differentiate these two DMSP phenotypes. A survey of prokaryotic and eukaryotic isolates found a significant correlation between the presence of DSYB and TpMT2 genes and previous measurements of high and low DMSP concentrations, respectively. Phylogenetic analysis demonstrated that DSYB and TpMT2 form two distinct clades. DSYB and TpMT2 were also found to be globally abundant in in situ surface communities, and their taxonomic annotations were similar to those observed for isolates. The strong correlation of the DSYB and TpMT2 synthesis genes with high and low producer phenotypes establishes a foundation for direct quantification of DMSP producers, enabling significantly improved predictions of DMSP in situ.