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

Non-cyanobacterial diazotrophs: global diversity, distribution, ecophysiology, and activity in marine waters

Biological dinitrogen (N2) fixation supplies nitrogen to the oceans, supporting primary productivity, and is carried out by some bacteria and archaea referred to as diazotrophs. Cyanobacteria are conventionally considered to be the major contributors to marine N2 fixation, but non-cyanobacterial diazotrophs (NCDs) have been shown to be distributed throughout ocean ecosystems. However, the biogeochemical significance of marine NCDs has not been demonstrated. This review synthesizes multiple datasets, drawing from cultivation-independent molecular techniques and data from extensive oceanic expeditions, to provide a comprehensive view into the diversity, biogeography, ecophysiology, and activity of marine NCDs. A NCD nifH gene catalog was compiled containing sequences from both PCR-based and PCR-free methods, identifying taxa for future studies. NCD abundances from a novel database of NCD nifH-based abundances were colocalized with environmental data, unveiling distinct distributions and environmental drivers of individual taxa. Mechanisms that NCDs may use to fuel and regulate N2 fixation in response to oxygen and fixed nitrogen availability are discussed, based on a metabolic analysis of recently available Tara Oceans expedition data. The integration of multiple datasets provides a new perspective that enhances understanding of the biology, ecology, and biogeography of marine NCDs and provides tools and directions for future research.

Species' functional traits and interactions drive nitrate-mediated sulfur-oxidizing community structure and functioning

IMPORTANCE

Understanding the processes and mechanisms governing microbial community assembly and their linkages to ecosystem functioning has long been a core issue in microbial ecology. An in-depth insight still requires combining with analyses of species' functional traits and microbial interactions. Our study showed how species' functional traits and interactions determined microbial community structure and functions by a well-controlled laboratory experiment with nitrate-mediated sulfur oxidation systems using high-throughput sequencing and culture-dependent technologies. The results provided solid evidences that species' functional traits and interactions were the intrinsic factors determining community structure and function. More importantly, our study established quantitative links between community structure and function based on species' functional traits and interactions, which would have important implications for the design and synthesis of microbiomes with expected functions.

Hydrodynamic regimes modulate nitrogen fixation and the mode of diazotrophy in Lake Tanganyika

The factors that govern the geographical distribution of nitrogen fixation are fundamental to providing accurate nitrogen budgets in aquatic environments. Model-based insights have demonstrated that regional hydrodynamics strongly impact nitrogen fixation. However, the mechanisms establishing this physical-biological coupling have yet to be constrained in field surveys. Here, we examine the distribution of nitrogen fixation in Lake Tanganyika - a model system with well-defined hydrodynamic regimes. We report that nitrogen fixation is five times higher under stratified than under upwelling conditions. Under stratified conditions, the limited resupply of inorganic nitrogen to surface waters, combined with greater light penetration, promotes the activity of bloom-forming photoautotrophic diazotrophs. In contrast, upwelling conditions support predominantly heterotrophic diazotrophs, which are uniquely suited to chemotactic foraging in a more dynamic nutrient landscape. We suggest that these hydrodynamic regimes (stratification versus mixing) play an important role in governing both the rates and the mode of nitrogen fixation.

Globally distributed Myxococcota with photosynthesis gene clusters illuminate the origin and evolution of a potentially chimeric lifestyle

Photosynthesis is a fundamental biogeochemical process, thought to be restricted to a few bacterial and eukaryotic phyla. However, understanding the origin and evolution of phototrophic organisms can be impeded and biased by the difficulties of cultivation. Here, we analyzed metagenomic datasets and found potential photosynthetic abilities encoded in the genomes of uncultivated bacteria within the phylum Myxococcota. A putative photosynthesis gene cluster encoding a type-II reaction center appears in at least six Myxococcota families from three classes, suggesting vertical inheritance of these genes from an early common ancestor, with multiple independent losses in other lineages. Analysis of metatranscriptomic datasets indicate that the putative myxococcotal photosynthesis genes are actively expressed in various natural environments. Furthermore, heterologous expression of myxococcotal pigment biosynthesis genes in a purple bacterium supports that the genes can drive photosynthetic processes. Given that predatory abilities are thought to be widespread across Myxococcota, our results suggest the intriguing possibility of a chimeric lifestyle (combining predatory and photosynthetic abilities) in members of this phylum.

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.

Ecological mechanisms and current systems shape the modular structure of the global oceans' prokaryotic seascape

Major biogeographic features of the microbial seascape in the oceans have been established and their underlying ecological mechanisms in the (sub)tropical oceans and the Pacific Ocean identified. However, we still lack a unifying understanding of how prokaryotic communities and biogeographic patterns are affected by large-scale current systems in distinct ocean basins and how they are globally shaped in line with ecological mechanisms. Here we show that prokaryotic communities in the epipelagic Pacific and Atlantic Ocean, in the southern Indian Ocean, and the Mediterranean Sea are composed of modules of co-occurring taxa with similar environmental preferences. The relative partitioning of these modules varies along latitudinal and longitudinal gradients and are related to different hydrographic and biotic conditions. Homogeneous selection and dispersal limitation were identified as the major ecological mechanisms shaping these communities and their free-living (FL) and particle-associated (PA) fractions. Large-scale current systems govern the dispersal of prokaryotic modules leading to the highest diversity near subtropical fronts.

Exploring environmental intra-species diversity through non-redundant pangenome assemblies

At the genome level, microorganisms are highly adaptable both in terms of allele and gene composition. Such heritable traits emerge in response to different environmental niches and can have a profound influence on microbial community dynamics. As a consequence, any individual genome or population will contain merely a fraction of the total genetic diversity of any operationally defined "species", whose ecological potential can thus be only fully understood by studying all of their genomes and the genes therein. This concept, known as the pangenome, is valuable for studying microbial ecology and evolution, as it partitions genomes into core (present in all the genomes from a species, and responsible for housekeeping and species-level niche adaptation among others) and accessory regions (present only in some, and responsible for intra-species differentiation). Here we present SuperPang, an algorithm producing pangenome assemblies from a set of input genomes of varying quality, including metagenome-assembled genomes (MAGs). SuperPang runs in linear time and its results are complete, non-redundant, preserve gene ordering and contain both coding and non-coding regions. Our approach provides a modular view of the pangenome, identifying operons and genomic islands, and allowing to track their prevalence in different populations. We illustrate this by analysing intra-species diversity in Polynucleobacter, a bacterial genus ubiquitous in freshwater ecosystems, characterized by their streamlined genomes and their ecological versatility. We show how SuperPang facilitates the simultaneous analysis of allelic and gene content variation under different environmental pressures, allowing us to study the drivers of microbial diversification at unprecedented resolution.

A global atlas of marine antibiotic resistance genes and their expression

Oceans serve as global reservoirs of antibiotic-resistant bacteria and antibiotic resistance genes (ARGs). However, little is known about the traits and expression of ARGs in response to environmental factors. We analyzed 347 metagenomes and 182 metatranscriptomes to determine the distribution, hosts, and expression of ARGs in oceans. Our study found that the diversity and abundance of ARGs varied with latitude and depth. The core marine resistome mainly conferred glycopeptide and multidrug resistance. The hosts of this resistome were mainly limited to the core marine microbiome, with phylogenetic barriers to the horizontal transfer of ARGs, transfers being more frequent within species than between species. Sixty-five percent of the marine ARGs identified were expressed. More than 90% of high-risk ARGs were more likely to be expressed. Anthropogenic activity might affect the expression of ARGs by altering nitrate and phosphate concentrations and ocean temperature. Machine-learning models predict >97% of marine ARGs will change expression by 2100. High-risk ARGs will shift to low latitudes and regions with high anthropogenic activity, such as the Pacific and Atlantic Oceans. Certain ARGs serve a dual role in antibiotic resistance and potentially participate in element cycling, along with other unknown functions. Determining whether changes in ARG expression are beneficial to ecosystems and human health is challenging without comprehensive understanding of their functions. Our study identified a core resistome in the oceans and quantified the expression of ARGs for the development of future control strategies under global change.

Clostridium scindens secretome suppresses virulence gene expression of Clostridioides difficile in a bile acid-independent manner

Clostridioides difficile infection (CDI) is a major health concern and one of the leading causes of hospital-acquired diarrhea in many countries. C. difficile infection is challenging to treat as C. difficile is resistant to multiple antibiotics. Alternative solutions are needed as conventional treatment with broad-spectrum antibiotics often leads to recurrent CDI. Recent studies have shown that specific microbiota-based therapeutics such as bile acids (BAs) are promising approaches to treat CDI. Clostridium scindens encodes the bile acid-induced (bai) operon that carries out 7-alpha-dehydroxylation of liver-derived primary BAs to secondary BAs. This biotransformation is thought to increase the antibacterial effects of BAs on C. difficile. Here, we used an automated multistage fermentor to study the antibacterial actions of C. scindens and BAs on C. difficile in the presence/absence of a gut microbial community derived from healthy human donor fecal microbiota. We observed that C. scindens inhibited C. difficile growth when the medium was supplemented with primary BAs. Transcriptomic analysis indicated upregulation of C. scindens bai operon and suppressed expression of C. difficile exotoxins that mediate CDI. We also observed BA-independent antibacterial activity of the secretome from C. scindens cultured overnight in a medium without supplementary primary BAs, which suppressed growth and exotoxin expression in C. difficile mono-culture. Further investigation of the molecular basis of our observation could lead to a more specific treatment for CDI than current approaches. IMPORTANCE There is an urgent need for new approaches to replace the available treatment options against Clostridioides difficile infection (CDI). Our novel work reports a bile acid-independent reduction of C. difficile growth and virulence gene expression by the secretome of Clostridium scindens. This potential treatment combined with other antimicrobial strategies could facilitate the development of alternative therapies in anticipation of CDI and in turn reduce the risk of antimicrobial resistance.

Consistency across multi-omics layers in a drug-perturbed gut microbial community

Multi-omics analyses are used in microbiome studies to understand molecular changes in microbial communities exposed to different conditions. However, it is not always clear how much each omics data type contributes to our understanding and whether they are concordant with each other. Here, we map the molecular response of a synthetic community of 32 human gut bacteria to three non-antibiotic drugs by using five omics layers (16S rRNA gene profiling, metagenomics, metatranscriptomics, metaproteomics and metabolomics). We find that all the omics methods with species resolution are highly consistent in estimating relative species abundances. Furthermore, different omics methods complement each other for capturing functional changes. For example, while nearly all the omics data types captured that the antipsychotic drug chlorpromazine selectively inhibits Bacteroidota representatives in the community, the metatranscriptome and metaproteome suggested that the drug induces stress responses related to protein quality control. Metabolomics revealed a decrease in oligosaccharide uptake, likely caused by Bacteroidota depletion. Our study highlights how multi-omics datasets can be utilized to reveal complex molecular responses to external perturbations in microbial communities.

A comprehensive genomic catalog from global cold seeps

Cold seeps harbor abundant and diverse microbes with tremendous potential for biological applications and that have a significant influence on biogeochemical cycles. Although recent metagenomic studies have expanded our understanding of the community and function of seep microorganisms, knowledge of the diversity and genetic repertoire of global seep microbes is lacking. Here, we collected a compilation of 165 metagenomic datasets from 16 cold seep sites across the globe to construct a comprehensive gene and genome catalog. The non-redundant gene catalog comprised 147 million genes, and 36% of them could not be assigned to a function with the currently available databases. A total of 3,164 species-level representative metagenome-assembled genomes (MAGs) were obtained, most of which (94%) belonged to novel species. Of them, 81 ANME species were identified that cover all subclades except ANME-2d, and 23 syntrophic SRB species spanned the Seep-SRB1a, Seep-SRB1g, and Seep-SRB2 clades. The non-redundant gene and MAG catalog is a valuable resource that will aid in deepening our understanding of the functions of cold seep microbiomes.

Human-induced salinity changes impact marine organisms and ecosystems

Climate change is fundamentally altering marine and coastal ecosystems on a global scale. While the effects of ocean warming and acidification on ecology and ecosystem functions and services are being comprehensively researched, less attention is directed toward understanding the impacts of human-driven ocean salinity changes. The global water cycle operates through water fluxes expressed as precipitation, evaporation, and freshwater runoff from land. Changes to these in turn modulate ocean salinity and shape the marine and coastal environment by affecting ocean currents, stratification, oxygen saturation, and sea level rise. Besides the direct impact on ocean physical processes, salinity changes impact ocean biological functions with the ecophysiological consequences are being poorly understood. This is surprising as salinity changes may impact diversity, ecosystem and habitat structure loss, and community shifts including trophic cascades. Climate model future projections (of end of the century salinity changes) indicate magnitudes that lead to modification of open ocean plankton community structure and habitat suitability of coral reef communities. Such salinity changes are also capable of affecting the diversity and metabolic capacity of coastal microorganisms and impairing the photosynthetic capacity of (coastal and open ocean) phytoplankton, macroalgae, and seagrass, with downstream ramifications on global biogeochemical cycling. The scarcity of comprehensive salinity data in dynamic coastal regions warrants additional attention. Such datasets are crucial to quantify salinity-based ecosystem function relationships and project such changes that ultimately link into carbon sequestration and freshwater as well as food availability to human populations around the globe. It is critical to integrate vigorous high-quality salinity data with interacting key environmental parameters (e.g., temperature, nutrients, oxygen) for a comprehensive understanding of anthropogenically induced marine changes and its impact on human health and the global economy.

Ocean-wide comparisons of mesopelagic planktonic community structures

For decades, marine plankton have been investigated for their capacity to modulate biogeochemical cycles and provide fishery resources. Between the sunlit (epipelagic) layer and the deep dark waters, lies a vast and heterogeneous part of the ocean: the mesopelagic zone. How plankton composition is shaped by environment has been well-explored in the epipelagic but much less in the mesopelagic ocean. Here, we conducted comparative analyses of trans-kingdom community assemblages thriving in the mesopelagic oxygen minimum zone (OMZ), mesopelagic oxic, and their epipelagic counterparts. We identified nine distinct types of intermediate water masses that correlate with variation in mesopelagic community composition. Furthermore, oxygen, NO3- and particle flux together appeared as the main drivers governing these communities. Novel taxonomic signatures emerged from OMZ while a global co-occurrence network analysis showed that about 70% of the abundance of mesopelagic plankton groups is organized into three community modules. One module gathers prokaryotes, pico-eukaryotes and Nucleo-Cytoplasmic Large DNA Viruses (NCLDV) from oxic regions, and the two other modules are enriched in OMZ prokaryotes and OMZ pico-eukaryotes, respectively. We hypothesize that OMZ conditions led to a diversification of ecological niches, and thus communities, due to selective pressure from limited resources. Our study further clarifies the interplay between environmental factors in the mesopelagic oxic and OMZ, and the compositional features of communities.

Glacial meltwater and seasonality influence community composition of diazotrophs in Arctic coastal and open waters

The Arctic Ocean is particularly affected by climate change with unknown consequences for primary productivity. Diazotrophs-prokaryotes capable of converting atmospheric nitrogen to ammonia-have been detected in the often nitrogen-limited Arctic Ocean but distribution and community composition dynamics are largely unknown. We performed amplicon sequencing of the diazotroph marker gene nifH from glacial rivers, coastal, and open ocean regions and identified regionally distinct Arctic communities. Proteobacterial diazotrophs dominated all seasons, epi- to mesopelagic depths and rivers to open waters and, surprisingly, Cyanobacteria were only sporadically identified in coastal and freshwaters. The upstream environment of glacial rivers influenced diazotroph diversity, and in marine samples putative anaerobic sulphate-reducers showed seasonal succession with highest prevalence in summer to polar night. Betaproteobacteria (Burkholderiales, Nitrosomonadales, and Rhodocyclales) were typically found in rivers and freshwater-influenced waters, and Delta- (Desulfuromonadales, Desulfobacterales, and Desulfovibrionales) and Gammaproteobacteria in marine waters. The identified community composition dynamics, likely driven by runoff, inorganic nutrients, particulate organic carbon, and seasonality, imply diazotrophy a phenotype of ecological relevance with expected responsiveness to ongoing climate change. Our study largely expands baseline knowledge of Arctic diazotrophs-a prerequisite to understand underpinning of nitrogen fixation-and supports nitrogen fixation as a contributor of new nitrogen in the rapidly changing Arctic Ocean.

Gene expression dynamics of natural assemblages of heterotrophic flagellates during bacterivory

BACKGROUND

Marine heterotrophic flagellates (HF) are dominant bacterivores in the ocean, where they represent the trophic link between bacteria and higher trophic levels and participate in the recycling of inorganic nutrients for regenerated primary production. Studying their activity and function in the ecosystem is challenging since most of the HFs in the ocean are still uncultured. In the present work, we investigated gene expression of natural HF communities during bacterivory in four unamended seawater incubations.

RESULTS

The most abundant species growing in our incubations belonged to the taxonomic groups MAST-4, MAST-7, Chrysophyceae, and Telonemia. Gene expression dynamics were similar between incubations and could be divided into three states based on microbial counts, each state displaying distinct expression patterns. The analysis of samples where HF growth was highest revealed some highly expressed genes that could be related to bacterivory. Using available genomic and transcriptomic references, we identified 25 species growing in our incubations and used those to compare the expression levels of these specific genes. Video Abstract CONCLUSIONS: Our results indicate that several peptidases, together with some glycoside hydrolases and glycosyltransferases, are more expressed in phagotrophic than in phototrophic species, and thus could be used to infer the process of bacterivory in natural assemblages.

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.

A ubiquitous gammaproteobacterial clade dominates expression of sulfur oxidation genes across the mesopelagic ocean

The deep ocean (>200 m depth) is the largest habitat on Earth. Recent evidence suggests sulfur oxidation could be a major energy source for deep ocean microbes. However, the global relevance and the identity of the major players in sulfur oxidation in the oxygenated deep-water column remain elusive. Here we combined single-cell genomics, community metagenomics, metatranscriptomics and single-cell activity measurements on samples collected beneath the Ross Ice Shelf in Antarctica to characterize a ubiquitous mixotrophic bacterial group (UBA868) that dominates expression of RuBisCO genes and of key sulfur oxidation genes. Further analyses of the gene libraries from the 'Tara Oceans' and 'Malaspina' expeditions confirmed the ubiquitous distribution and global relevance of this enigmatic group in the expression of sulfur oxidation and dissolved inorganic carbon fixation genes across the global mesopelagic ocean. Our study also underscores the unrecognized importance of mixotrophic microbes in the biogeochemical cycles of the deep ocean.

Profiling the microbiome of oral and genital mucosal surfaces in Behçet's disease

Almost 90% of Behçet's patients present with oral and/or genital ulcers which influence the disease outcome. We hypothesised that the dysregulation of the oral and genital microbiome, coupled with dysregulation of the immune response, contributes to the aetiopathogenesis of Behçet's Disease (BD) and drives disease activation.

METHOD

152 BD patient samples; 70 matched oral and genital samples plus 12 unmatched samples (Female: Male, 58:12; mean age, 42 ± 13.9: 39.3 ± 10.3) to profile microbial community high-throughput sequencing of the microbiome using 16 s rRNA sequencing targeting the V1/V2 and V3/V4 hyper variable regions were used and results reviewed in relation to disease severity, Work and Social Adjustment Scale (WSAS) outcomes and medication.

RESULTS

Alpha and beta diversity were significantly decreased in genital compared to oral samples; p value<0.05. However, grouping the samples as to whether ulceration was present was not significant. Escherichia-Shigella was the only Amplicon Sequence Variants (ASVs) in the V1/V2 region that was shared between the oral mucosa with ulcer (O_U) and genital mucosa with ulcer (G_U) groups. This was in contrast to the V3/V4 region which indicated that Lachnospiraceae, Saccharimonadales, and Coriobacteriales were shared between the O_U and G_U groups. In addition, gender had no impact on the bacterial abundance in V1/V2 analysis of the oral and genital samples. V3/V4 analysis of genital samples demonstrated that Lactobacilli and Gardnerella were significantly increased in females (20 times) compared to the males in samples; p-adj <0.05. Interestingly in BD patients, Rothia which is commonly found in the mouth was present in both oral and genital samples. Streptococci were significantly increased while Veillonella significantly decreased in the presence of oral ulceration in the BD cohort. The clinical phenotype had no effect on V1/V2 and V3/V4 on the bacterial abundance of oral samples. However, medication e.g. colchicine had a significant effect on the oral microbial abundance (V1/V2; P = 0.020, V3/V4; P = 0.003). There was no relationship between colchicine and the presence/absence of genital ulcers. BD patients with active disease had higher WSAS scores, and their bacterial abundance differed significantly from the non-active BD patients (ADONIS, R2 = 0.05, p value =0.029).

CONCLUSION

The presence of the microbes Streptococcus, Veillonella, Gardnerella, Lactobacillus, Atopobium, Peptoniphilus, Corynebacterium and Staphylococcus may provide early evidence of BD patients are with active disease.

Distribution and survival strategies of endemic and cosmopolitan diazotrophs in the Arctic Ocean

Dinitrogen (N₂) fixation is the major source of reactive nitrogen in the ocean and has been considered to occur specifically in low-latitude oligotrophic oceans. Recent studies have shown that N₂ fixation also occurs in the polar regions and thus is a global process, although the physiological and ecological characteristics of polar diazotrophs are not yet known. Here, we successfully reconstructed diazotroph genomes, including that of cyanobacterium UCYN-A (Candidatus 'Atelocyanobacterium thalassa'), from metagenome data corresponding to 111 samples isolated from the Arctic Ocean. These diazotrophs were highly abundant in the Arctic Ocean (max., 1.28% of the total microbial community), suggesting that they have important roles in the Arctic ecosystem and biogeochemical cycles. Further, we show that diazotrophs within genera Arcobacter, Psychromonas, and Oceanobacter are prevalent in the <0.2 µm fraction in the Arctic Ocean, indicating that current methods cannot capture their N₂ fixation. Diazotrophs in the Arctic Ocean were either Arctic-endemic or cosmopolitan species from their global distribution patterns. Arctic-endemic diazotrophs, including Arctic UCYN-A, were similar to low-latitude-endemic and cosmopolitan diazotrophs in genome-wide function, however, they had unique gene sets (e.g., diverse aromatics degradation genes), suggesting adaptations to Arctic-specific conditions. Cosmopolitan diazotrophs were generally non-cyanobacteria and commonly had the gene that encodes the cold-inducible RNA chaperone, which presumably makes their survival possible even in deep, cold waters of global ocean and polar surface waters. This study shows global distribution pattern of diazotrophs with their genomes and provides clues to answering the question of how diazotrophs can inhabit polar waters.

Aerobic bacteria produce nitric oxide via denitrification and promote algal population collapse

Microbial interactions govern marine biogeochemistry. These interactions are generally considered to rely on exchange of organic molecules. Here we report on a novel inorganic route of microbial communication, showing that algal-bacterial interactions between Phaeobacter inhibens bacteria and Gephyrocapsa huxleyi algae are mediated through inorganic nitrogen exchange. Under oxygen-rich conditions, aerobic bacteria reduce algal-secreted nitrite to nitric oxide (NO) through denitrification, a well-studied anaerobic respiratory mechanism. The bacterial NO is involved in triggering a cascade in algae akin to programmed cell death. During death, algae further generate NO, thereby propagating the signal in the algal population. Eventually, the algal population collapses, similar to the sudden demise of oceanic algal blooms. Our study suggests that the exchange of inorganic nitrogen species in oxygenated environments is a potentially significant route of microbial communication within and across kingdoms.

Rhodopsin-mediated nutrient uptake by cultivated photoheterotrophic Verrucomicrobiota

Rhodopsin photosystems convert light energy into electrochemical gradients used by the cell to produce ATP, or for other energy-demanding processes. While these photosystems are widespread in the ocean and have been identified in diverse microbial taxonomic groups, their physiological role in vivo has only been studied in few marine bacterial strains. Recent metagenomic studies revealed the presence of rhodopsin genes in the understudied Verrucomicrobiota phylum, yet their distribution within different Verrucomicrobiota lineages, their diversity, and function remain unknown. In this study, we show that more than 7% of Verrucomicrobiota genomes (n = 2916) harbor rhodopsins of different types. Furthermore, we describe the first two cultivated rhodopsin-containing strains, one harboring a proteorhodopsin gene and the other a xanthorhodopsin gene, allowing us to characterize their physiology under laboratory-controlled conditions. The strains were isolated in a previous study from the Eastern Mediterranean Sea and read mapping of 16S rRNA gene amplicons showed the highest abundances of these strains at the deep chlorophyll maximum (source of their inoculum) in winter and spring, with a substantial decrease in summer. Genomic analysis of the isolates suggests that motility and degradation of organic material, both energy demanding functions, may be supported by rhodopsin phototrophy in Verrucomicrobiota. Under culture conditions, we show that rhodopsin phototrophy occurs under carbon starvation, with light-mediated energy generation supporting sugar transport into the cells. Overall, this study suggests that photoheterotrophic Verrucomicrobiota may occupy an ecological niche where energy harvested from light enables bacterial motility toward organic matter and supports nutrient uptake.

Functional annotation and importance of marine bacterial transporters of plankton exometabolites

Metabolite exchange within marine microbial communities transfers carbon and other major elements through global cycles and forms the basis of microbial interactions. Yet lack of gene annotations and concern about the quality of existing ones remain major impediments to revealing currencies of carbon flux. We employed an arrayed mutant library of the marine bacterium Ruegeria pomeroyi DSS-3 to experimentally annotate substrates of organic compound transporter systems, using mutant growth and compound drawdown analyses to link transporters to their cognate substrates. Mutant experiments verified substrates for thirteen R. pomeroyi transporters. Four were previously hypothesized based on gene expression data (taurine, glucose/xylose, isethionate, and cadaverine/putrescine/spermidine); five were previously hypothesized based on homology to experimentally annotated transporters in other bacteria (citrate, glycerol, N-acetylglucosamine, fumarate/malate/succinate, and dimethylsulfoniopropionate); and four had no previous annotations (thymidine, carnitine, cysteate, and 3-hydroxybutyrate). These bring the total number of experimentally-verified organic carbon influx transporters to 18 of 126 in the R. pomeroyi genome. In a longitudinal study of a coastal phytoplankton bloom, expression patterns of the experimentally annotated transporters linked them to different stages of the bloom, and also led to the hypothesis that citrate and 3-hydroxybutyrate were among the most highly available bacterial substrates. Improved functional annotation of the gatekeepers of organic carbon uptake is critical for deciphering carbon flux and fate in microbial ecosystems.

Transcriptional Mechanisms of Thermal Acclimation in Prochlorococcus

Low temperature limits the growth and the distribution of the key oceanic primary producer Prochlorococcus, which does not proliferate above a latitude of ca. 40°. Yet, the molecular basis of thermal acclimation in this cyanobacterium remains unexplored. We analyzed the transcriptional response of the Prochlorococcus marinus strain MIT9301 in long-term acclimations and in natural Prochlorococcus populations along a temperature range enabling its growth (17 to 30°C). MIT9301 upregulated mechanisms of the global stress response at the temperature minimum (17°C) but maintained the expression levels of genes involved in essential metabolic pathways (e.g., ATP synthesis and carbon fixation) along the whole thermal niche. Notably, the declining growth of MIT9301 from the optimum to the minimum temperature was coincident with a transcriptional suppression of the photosynthetic apparatus and a dampening of its circadian expression patterns, indicating a loss in their regulatory capacity under cold conditions. Under warm conditions, the cellular transcript inventory of MIT9301 was strongly streamlined, which may also induce regulatory imbalances due to stochasticity in gene expression. The daytime transcriptional suppression of photosynthetic genes at low temperature was also observed in metatranscriptomic reads mapping to MIT9301 across the global ocean, implying that this molecular mechanism may be associated with the restricted distribution of Prochlorococcus to temperate zones. IMPORTANCE Prochlorococcus is a major marine primary producer with a global impact on atmospheric CO₂ fixation. This cyanobacterium is widely distributed across the temperate ocean, but virtually absent at latitudes above 40° for yet unknown reasons. Temperature has been suggested as a major limiting factor, but the exact mechanisms behind Prochlorococcus thermal growth restriction remain unexplored. This study brings us closer to understanding how Prochlorococcus functions under challenging temperature conditions, by focusing on its transcriptional response after long-term acclimation from its optimum to its thermal thresholds. Our results show that the drop in Prochlorococcus growth rate under cold conditions was paralleled by a transcriptional suppression of the photosynthetic machinery during daytime and a loss in the organism's regulatory capacity to maintain circadian expression patterns. Notably, warm temperature induced a marked shrinkage of the organism's cellular transcript inventory, which may also induce regulatory imbalances in the future functioning of this cyanobacterium.

Anaerobic thiosulfate oxidation by the Roseobacter group is prevalent in marine biofilms

Thiosulfate oxidation by microbes has a major impact on global sulfur cycling. Here, we provide evidence that bacteria within various Roseobacter lineages are important for thiosulfate oxidation in marine biofilms. We isolate and sequence the genomes of 54 biofilm-associated Roseobacter strains, finding conserved sox gene clusters for thiosulfate oxidation and plasmids, pointing to a niche-specific lifestyle. Analysis of global ocean metagenomic data suggests that Roseobacter strains are abundant in biofilms and mats on various substrates, including stones, artificial surfaces, plant roots, and hydrothermal vent chimneys. Metatranscriptomic analysis indicates that the majority of active sox genes in biofilms belong to Roseobacter strains. Furthermore, we show that Roseobacter strains can grow and oxidize thiosulfate to sulfate under both aerobic and anaerobic conditions. Transcriptomic and membrane proteomic analyses of biofilms formed by a representative strain indicate that thiosulfate induces sox gene expression and alterations in cell membrane protein composition, and promotes biofilm formation and anaerobic respiration. We propose that bacteria of the Roseobacter group are major thiosulfate-oxidizers in marine biofilms, where anaerobic thiosulfate metabolism is preferred.

The Landscape of Global Ocean Microbiome: From Bacterioplankton to Biofilms

The development of metagenomics has opened up a new era in the study of marine microbiota, which play important roles in biogeochemical cycles. In recent years, the global ocean sampling expeditions have spurred this research field toward a deeper understanding of the microbial diversities and functions spanning various lifestyles, planktonic (free-living) or sessile (biofilm-associated). In this review, we deliver a comprehensive summary of marine microbiome datasets generated in global ocean expeditions conducted over the last 20 years, including the Sorcerer II GOS Expedition, the Tara Oceans project, the bioGEOTRACES project, the Micro B3 project, the Bio-GO-SHIP project, and the Marine Biofilms. These datasets have revealed unprecedented insights into the microscopic life in our oceans and led to the publication of world-leading research. We also note the progress of metatranscriptomics and metaproteomics, which are confined to local marine microbiota. Furthermore, approaches to transforming the global ocean microbiome datasets are highlighted, and the state-of-the-art techniques that can be combined with data analyses, which can present fresh perspectives on marine molecular ecology and microbiology, are proposed.

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.

Phototrophy by antenna-containing rhodopsin pumps in aquatic environments

Energy transfer from light-harvesting ketocarotenoids to the light-driven proton pump xanthorhodopsins has been previously demonstrated in two unique cases: an extreme halophilic bacterium¹ and a terrestrial cyanobacterium². Attempts to find carotenoids that bind and transfer energy to abundant rhodopsin proton pumps³ from marine photoheterotrophs have thus far failed^4-6. Here we detected light energy transfer from the widespread hydroxylated carotenoids zeaxanthin and lutein to the retinal moiety of xanthorhodopsins and proteorhodopsins using functional metagenomics combined with chromophore extraction from the environment. The light-harvesting carotenoids transfer up to 42% of the harvested energy in the violet- or blue-light range to the green-light absorbing retinal chromophore. Our data suggest that these antennas may have a substantial effect on rhodopsin phototrophy in the world's lakes, seas and oceans. However, the functional implications of our findings are yet to be discovered.

Structure-informed microbial population genetics elucidate selective pressures that shape protein evolution

Comprehensive sampling of natural genetic diversity with metagenomics enables highly resolved insights into the interplay between ecology and evolution. However, resolving adaptive, neutral, or purifying processes of evolution from intrapopulation genomic variation remains a challenge, partly due to the sole reliance on gene sequences to interpret variants. Here, we describe an approach to analyze genetic variation in the context of predicted protein structures and apply it to a marine microbial population within the SAR11 subclade 1a.3.V, which dominates low-latitude surface oceans. Our analyses reveal a tight association between genetic variation and protein structure. In a central gene in nitrogen metabolism, we observe decreased occurrence of nonsynonymous variants from ligand-binding sites as a function of nitrate concentrations, revealing genetic targets of distinct evolutionary pressures maintained by nutrient availability. Our work yields insights into the governing principles of evolution and enables structure-aware investigations of microbial population genetics.

Molecular hydrogen in seawater supports growth of diverse marine bacteria

Molecular hydrogen (H₂) is an abundant and readily accessible energy source in marine systems, but it remains unknown whether marine microbial communities consume this gas. Here we use a suite of approaches to show that marine bacteria consume H₂ to support growth. Genes for H₂-uptake hydrogenases are prevalent in global ocean metagenomes, highly expressed in metatranscriptomes and found across eight bacterial phyla. Capacity for H₂ oxidation increases with depth and decreases with oxygen concentration, suggesting that H₂ is important in environments with low primary production. Biogeochemical measurements of tropical, temperate and subantarctic waters, and axenic cultures show that marine microbes consume H₂ supplied at environmentally relevant concentrations, yielding enough cell-specific power to support growth in bacteria with low energy requirements. Conversely, our results indicate that oxidation of carbon monoxide (CO) primarily supports survival. Altogether, H₂ is a notable energy source for marine bacteria and may influence oceanic ecology and biogeochemistry.

Postglacial adaptations enabled colonization and quasi-clonal dispersal of ammonia-oxidizing archaea in modern European large lakes

Ammonia-oxidizing archaea (AOA) play a key role in the aquatic nitrogen cycle. Their genetic diversity is viewed as the outcome of evolutionary processes that shaped ancestral transition from terrestrial to marine habitats. However, current genome-wide insights into AOA evolution rarely consider brackish and freshwater representatives or provide their divergence timeline in lacustrine systems. An unbiased global assessment of lacustrine AOA diversity is critical for understanding their origins, dispersal mechanisms, and ecosystem roles. Here, we leveraged continental-scale metagenomics to document that AOA species diversity in freshwater systems is remarkably low compared to marine environments. We show that the uncultured freshwater AOA, "Candidatus Nitrosopumilus limneticus," is ubiquitous and genotypically static in various large European lakes where it evolved 13 million years ago. We find that extensive proteome remodeling was a key innovation for freshwater colonization of AOA. These findings reveal the genetic diversity and adaptive mechanisms of a keystone species that has survived clonally in lakes for millennia.

Metatranscriptomes reveal the diverse responses of Thaumarchaeota ecotypes to environmental variations in the northern slope of the South China Sea

Thaumarchaeota are among the most abundant prokaryotes in the ocean, playing important roles in carbon and nitrogen cycling. Marine Thaumarchaeota ecotypes exhibit depth-related diversification and seasonal changes. However, transcriptomic activities concerning niche partitioning among thaumarchaeal ecotypes remain unclear. Here, we examined the variations in the distribution and transcriptomic activity of marine Thaumarchaeota ecotypes. Three primary ecotypes were identified: a Nitrosopumilus-like clade; a Nitrosopelagicus-like water column A (WCA) clade, thriving in epipelagic water; and a water column B (WCB) clade, dominant in deep water. Depth-related partitioning of the three ecotypes and the seasonal variability of the WCA and WCB ecotypes were observed. Nutrient concentrations, chlorophyll α and salinity were the primary environmental factors. The relative abundance of the WCA ecotype and its transcript abundance of amoA gene were positively correlated with chlorophyll α and salinity, while the WCB ecotype was positively correlated with nitrate and phosphate. Based on high-quality metagenome-assembled genomes, transcriptomic analysis revealed that the three ecotypes exhibited various co-occurring expression patterns of the elemental cycling genes in the nitrogen, carbon, phosphorus, and sulfur cycles. Our results provide transcriptomic evidence of the niche differentiation of marine Thaumarchaeota ecotypes, highlighting the diverse roles of ecotypes and WCA subclades in biogeochemical cycles.

Metagenomic Views of Microbial Communities in Sand Sediments Associated with Coral Reefs

Reef sediments, the home for microbes with high abundances, provide an important source of carbonates and nutrients for the growth and maintenance of coral reefs. However, there is a lack of systematic research on the composition of microbial community in sediments of different geographic sites and their potential effect on nutrient recycling and health of the coral reef ecosystem. In combination of biogeochemical measurements with gene- and genome-centric metagenomics, we assessed microbial community compositions and functional diversity, as well as profiles of antibiotic resistance genes in surface sediments of 16 coral reef sites at different depths from the Xisha islands in the South China Sea. Reef sediment microbiomes are diverse and novel at lower taxonomic ranks, dominated by Proteobacteria and Planctomycetota. Most reef sediment bacteria potentially participate in biogeochemical cycling via oxidizing various organic and inorganic compounds as energy sources. High abundances of Proteobacteria (mostly Rhizobiales and Woeseiales) are metabolically flexible and contain rhodopsin genes. Various classes of antibiotic resistance genes, hosted by diverse bacterial lineages, were identified to confer resistance to multidrug, aminoglycoside, and other antibiotics. Overall, our findings expanded the understanding of reef sediment microbial ecology and provided insights for their link to the coral reef ecosystem health.

The hitchhikers' guide to RNA sequencing and functional analysis

DNA and RNA sequencing technologies have revolutionized biology and biomedical sciences, sequencing full genomes and transcriptomes at very high speeds and reasonably low costs. RNA sequencing (RNA-Seq) enables transcript identification and quantification, but once sequencing has concluded researchers can be easily overwhelmed with questions such as how to go from raw data to differential expression (DE), pathway analysis and interpretation. Several pipelines and procedures have been developed to this effect. Even though there is no unique way to perform RNA-Seq analysis, it usually follows these steps: 1) raw reads quality check, 2) alignment of reads to a reference genome, 3) aligned reads' summarization according to an annotation file, 4) DE analysis and 5) gene set analysis and/or functional enrichment analysis. Each step requires researchers to make decisions, and the wide variety of options and resulting large volumes of data often lead to interpretation challenges. There also seems to be insufficient guidance on how best to obtain relevant information and derive actionable knowledge from transcription experiments. In this paper, we explain RNA-Seq steps in detail and outline differences and similarities of different popular options, as well as advantages and disadvantages. We also discuss non-coding RNA analysis, multi-omics, meta-transcriptomics and the use of artificial intelligence methods complementing the arsenal of tools available to researchers. Lastly, we perform a complete analysis from raw reads to DE and functional enrichment analysis, visually illustrating how results are not absolute truths and how algorithmic decisions can greatly impact results and interpretation.

Water masses shape pico-nano eukaryotic communities of the Weddell Sea

Polar oceans belong to the most productive and rapidly changing environments, yet our understanding of this fragile ecosystem remains limited. Here we present an analysis of a unique set of DNA metabarcoding samples from the western Weddell Sea sampled throughout the whole water column and across five water masses with different characteristics and different origin. We focus on factors affecting the distribution of planktonic pico-nano eukaryotes and observe an ecological succession of eukaryotic communities as the water masses move away from the surface and as oxygen becomes depleted with time. At the beginning of this succession, in the photic zone, algae, bacteriovores, and predators of small eukaryotes dominate the community, while another community develops as the water sinks deeper, mostly composed of parasitoids (syndinians), mesoplankton predators (radiolarians), and diplonemids. The strongly correlated distribution of syndinians and diplonemids along the depth and oxygen gradients suggests their close ecological link and moves us closer to understanding the biological role of the latter group in the ocean ecosystem.

Diazotrophs are overlooked contributors to carbon and nitrogen export to the deep ocean

Diazotrophs are widespread microorganisms that alleviate nitrogen limitation in 60% of our oceans, thereby regulating marine productivity. Yet, the group-specific contribution of diazotrophs to organic matter export has not been quantified, which so far has impeded an accurate assessment of their impact on the biological carbon pump. Here, we examine the fate of five groups of globally-distributed diazotrophs by using an original combination of mesopelagic particle sampling devices across the subtropical South Pacific Ocean. We demonstrate that cyanobacterial and non-cyanobacterial diazotrophs are exported down to 1000 m depth. Surprisingly, group-specific export turnover rates point to a more efficient export of small unicellular cyanobacterial diazotrophs (UCYN) relative to the larger and filamentous Trichodesmium. Phycoerythrin-containing UCYN-B and UCYN-C-like cells were recurrently found embedded in large (>50 µm) organic aggregates or organized into clusters of tens to hundreds of cells linked by an extracellular matrix, presumably facilitating their export. Beyond the South Pacific, our data are supported by analysis of the Tara Oceans metagenomes collected in other ocean basins, extending the scope of our results globally. We show that, when diazotrophs are found in the euphotic zone, they are also systematically present in mesopelagic waters, suggesting their transport to the deep ocean. We thus conclude that diazotrophs are a significant part of the carbon sequestered in the deep ocean and, therefore, they need to be accounted in regional and global estimates of export.

Assessing and removing the effect of unwanted technical variations in microbiome data

Varying technologies and experimental approaches used in microbiome studies often lead to irreproducible results due to unwanted technical variations. Such variations, often unaccounted for and of unknown source, may interfere with true biological signals, resulting in misleading biological conclusions. In this work, we aim to characterize the major sources of technical variations in microbiome data and demonstrate how in-silico approaches can minimize their impact. We analyzed 184 pig faecal metagenomes encompassing 21 specific combinations of deliberately introduced factors of technical and biological variations. Using the novel Removing Unwanted Variations-III-Negative Binomial (RUV-III-NB), we identified several known experimental factors, specifically storage conditions and freeze-thaw cycles, as likely major sources of unwanted variation in metagenomes. We also observed that these unwanted technical variations do not affect taxa uniformly, with freezing samples affecting taxa of class Bacteroidia the most, for example. Additionally, we benchmarked the performances of different correction methods, including ComBat, ComBat-seq, RUVg, RUVs, and RUV-III-NB. While RUV-III-NB performed consistently robust across our sensitivity and specificity metrics, most other methods did not remove unwanted variations optimally. Our analyses suggest that a careful consideration of possible technical confounders is critical during experimental design of microbiome studies, and that the inclusion of technical replicates is necessary to efficiently remove unwanted variations computationally.

Transcriptomic-Guided Phosphonate Utilization Analysis Unveils Evidence of Clathrin-Mediated Endocytosis and Phospholipid Synthesis in the Model Diatom, Phaeodactylum tricornutum

Phosphonates are important components of marine organic phosphorus, but their bioavailability and catabolism by eukaryotic phytoplankton remain enigmatic. Here, diatom Phaeodactylum tricornutum was used to investigate the bioavailability of phosphonates and describe the underlying molecular mechanism. The results showed that 2-aminoethylphosphonic acid (2-AEP) can be utilized as an alternative phosphorus source. Comparative transcriptomics revealed that the utilization of 2-AEP comprised 2 steps, including molecular uptake through clathrin-mediated endocytosis and incorporation into the membrane phospholipids in the form of diacylglyceryl-2-AEP (DAG-2-AEP). In the global ocean, we found the prevalence and dynamic expression pattern of key genes that are responsible for vesicle formation (CLTC, AP-2) and DAG-AEP synthesis (PCYT2, EPT1) in diatom assemblages. This study elucidates a distinctive mechanism of phosphonate utilization by diatoms, and discusses the ecological implications. IMPORTANCE Phosphonates contribute ~25% of total dissolved organic phosphorus in the ocean, and are found to be important for marine phosphorus biogeochemical cycle. As a type of biogenic phosphonate produced by microorganisms, 2-aminoethylphosphonic acid (2-AEP) widely exists in the ocean. It is well known that 2-AEP can be cleaved and utilized by prokaryotes, but its ability to support the growth of eukaryotic phytoplankton remains unclear. Our research identified the bioavailability of 2-AEP for the diatom Phaeodactylum tricornutum, and proposed a distinctive metabolic pathway of 2-AEP utilization. Different from the enzymatic hydrolysis of phosphonates, the results suggested that P. tricornutum utilizes 2-AEP by incorporating it into phospholipid instead of cleaving the C-P bond. Moreover, the ubiquitous distribution of associated representative gene transcripts in the environmental assemblages and the higher gene transcript abundance in the cold regions were observed, which suggests the possible environmental adaption of 2-AEP utilization by diatoms.

Unleashing the power within short-read RNA-seq for plant research: Beyond differential expression analysis and toward regulomics

RNA-seq has become a state-of-the-art technique for transcriptomic studies. Advances in both RNA-seq techniques and the corresponding analysis tools and pipelines have unprecedently shaped our understanding in almost every aspects of plant sciences. Notably, the integration of huge amount of RNA-seq with other omic data sets in the model plants and major crop species have facilitated plant regulomics, while the RNA-seq analysis has still been primarily used for differential expression analysis in many less-studied plant species. To unleash the analytical power of RNA-seq in plant species, especially less-studied species and biomass crops, we summarize recent achievements of RNA-seq analysis in the major plant species and representative tools in the four types of application: (1) transcriptome assembly, (2) construction of expression atlas, (3) network analysis, and (4) structural alteration. We emphasize the importance of expression atlas, coexpression networks and predictions of gene regulatory relationships in moving plant transcriptomes toward regulomics, an omic view of genome-wide transcription regulation. We highlight what can be achieved in plant research with RNA-seq by introducing a list of representative RNA-seq analysis tools and resources that are developed for certain minor species or suitable for the analysis without species limitation. In summary, we provide an updated digest on RNA-seq tools, resources and the diverse applications for plant research, and our perspective on the power and challenges of short-read RNA-seq analysis from a regulomic point view. A full utilization of these fruitful RNA-seq resources will promote plant omic research to a higher level, especially in those less studied species.

Cultivation-independent genomes greatly expand taxonomic-profiling capabilities of mOTUs across various environments

BACKGROUND

Taxonomic profiling is a fundamental task in microbiome research that aims to detect and quantify the relative abundance of microorganisms in biological samples. Available methods using shotgun metagenomic data generally depend on the deposition of sequenced and taxonomically annotated genomes, usually from cultures of isolated strains, in reference databases (reference genomes). However, the majority of microorganisms have not been cultured yet. Thus, a substantial fraction of microbial community members remains unaccounted for during taxonomic profiling, particularly in samples from underexplored environments. To address this issue, we developed the mOTU profiler, a tool that enables reference genome-independent species-level profiling of metagenomes. As such, it supports the identification and quantification of both "known" and "unknown" species based on a set of select marker genes.

RESULTS

We present mOTUs3, a command line tool that enables the profiling of metagenomes for >33,000 species-level operational taxonomic units. To achieve this, we leveraged the reconstruction of >600,000 draft genomes, most of which are metagenome-assembled genomes (MAGs), from diverse microbiomes, including soil, freshwater systems, and the gastrointestinal tract of ruminants and other animals, which we found to be underrepresented by reference genomes. Overall, two thirds of all species-level taxa lacked a reference genome. The cumulative relative abundance of these newly included taxa was low in well-studied microbiomes, such as the human body sites (6-11%). By contrast, they accounted for substantial proportions (ocean, freshwater, soil: 43-63%) or even the majority (pig, fish, cattle: 60-80%) of the relative abundance across diverse non-human-associated microbiomes. Using community-developed benchmarks and datasets, we found mOTUs3 to be more accurate than other methods and to be more congruent with 16S rRNA gene-based methods for taxonomic profiling. Furthermore, we demonstrate that mOTUs3 increases the resolution of well-known microbial groups into species-level taxa and helps identify new differentially abundant taxa in comparative metagenomic studies.

CONCLUSIONS

We developed mOTUs3 to enable accurate species-level profiling of metagenomes. Compared to other methods, it provides a more comprehensive view of prokaryotic community diversity, in particular for currently underexplored microbiomes. To facilitate comparative analyses by the research community, it is released with >11,000 precomputed profiles for publicly available metagenomes and is freely available at: https://github.com/motu-tool/mOTUs . Video Abstract.

Comparison of Metagenomics and Metatranscriptomics Tools: A Guide to Making the Right Choice

The study of microorganisms is a field of great interest due to their environmental (e.g., soil contamination) and biomedical (e.g., parasitic diseases, autism) importance. The advent of revolutionary next-generation sequencing techniques, and their application to the hypervariable regions of the 16S, 18S or 23S ribosomal subunits, have allowed the research of a large variety of organisms more in-depth, including bacteria, archaea, eukaryotes and fungi. Additionally, together with the development of analysis software, the creation of specific databases (e.g., SILVA or RDP) has boosted the enormous growth of these studies. As the cost of sequencing per sample has continuously decreased, new protocols have also emerged, such as shotgun sequencing, which allows the profiling of all taxonomic domains in a sample. The sequencing of hypervariable regions and shotgun sequencing are technologies that enable the taxonomic classification of microorganisms from the DNA present in microbial communities. However, they are not capable of measuring what is actively expressed. Conversely, we advocate that metatranscriptomics is a "new" technology that makes the identification of the mRNAs of a microbial community possible, quantifying gene expression levels and active biological pathways. Furthermore, it can be also used to characterise symbiotic interactions between the host and its microbiome. In this manuscript, we examine the three technologies above, and discuss the implementation of different software and databases, which greatly impact the obtaining of reliable results. Finally, we have developed two easy-to-use pipelines leveraging Nextflow technology. These aim to provide everything required for an average user to perform a metagenomic analysis of marker genes with QIMME2 and a metatranscriptomic study using Kraken2/Bracken.

High-quality genome of the basidiomycete yeast Dioszegia hungarica PDD-24b-2 isolated from cloud water

The genome of the basidiomycete yeast Dioszegia hungarica strain PDD-24b-2 isolated from cloud water at the summit of puy de Dôme (France) was sequenced using a hybrid PacBio and Illumina sequencing strategy. The obtained assembled genome of 20.98 Mb and a GC content of 57% is structured in 16 large-scale contigs ranging from 90 kb to 5.56 Mb, and another 27.2 kb contig representing the complete circular mitochondrial genome. In total, 8,234 proteins were predicted from the genome sequence. The mitochondrial genome shows 16.2% cgu codon usage for arginine but has no canonical cognate tRNA to translate this codon. Detected transposable element (TE)-related sequences account for about 0.63% of the assembled genome. A dataset of 2,068 hand-picked public environmental metagenomes, representing over 20 Tbp of raw reads, was probed for D. hungarica related ITS sequences, and revealed worldwide distribution of this species, particularly in aerial habitats. Growth experiments suggested a psychrophilic phenotype and the ability to disperse by producing ballistospores. The high-quality assembled genome obtained for this D. hungarica strain will help investigate the behavior and ecological functions of this species in the environment.

Discovering marine biodiversity in the 21st century

We review the current knowledge of the biodiversity of the ocean as well as the levels of decline and threat for species and habitats. The lack of understanding of the distribution of life in the ocean is identified as a significant barrier to restoring its biodiversity and health. We explore why the science of taxonomy has failed to deliver knowledge of what species are present in the ocean, how they are distributed and how they are responding to global and regional to local anthropogenic pressures. This failure prevents nations from meeting their international commitments to conserve marine biodiversity with the results that investment in taxonomy has declined in many countries. We explore a range of new technologies and approaches for discovery of marine species and their detection and monitoring. These include: imaging methods, molecular approaches, active and passive acoustics, the use of interconnected databases and citizen science. Whilst no one method is suitable for discovering or detecting all groups of organisms many are complementary and have been combined to give a more complete picture of biodiversity in marine ecosystems. We conclude that integrated approaches represent the best way forwards for accelerating species discovery, description and biodiversity assessment. Examples of integrated taxonomic approaches are identified from terrestrial ecosystems. Such integrated taxonomic approaches require the adoption of cybertaxonomy approaches and will be boosted by new autonomous sampling platforms and development of machine-speed exchange of digital information between databases.

Taxonomic and functional dynamics during chytrid epidemics in an aquatic ecosystem

Fungal parasitism is common in plankton communities and plays a crucial role in the ecosystem by balancing nutrient cycling in the food web. Previous studies of aquatic ecosystems revealed that zoosporic chytrid epidemics represent an important driving factor in phytoplankton seasonal successions. In this study, host-parasite dynamics in Lake Pavin (France) were investigated during the spring diatom bloom while following chytrid epidemics using next generation sequencing (NGS). Metabarcoding analyses were applied to study changes in the eukaryotic microbial community throughout diatom bloom-chytrid epidemics. Relative read abundances of metabarcoding data revealed potential "beneficiaries" and "victims" during the studied period. Subsequently, metatranscriptomic analyses on samples before and during the chytrid epidemic unveiled the active part of the community and functional/metabolic dynamics in association with the progress of chytrid infection. Diatom functions involving lipases, transporters, histones, vacuolar systems, the proteasome, proteases and DNA/RNA polymerases were more abundant during the diatom bloom. Chytrid functions related to a parasitic lifestyle including invasion, colonization and stress tolerance were up-regulated during the chytrid epidemic. In addition, functions related to the degradation/metabolism of proteins, lipids and chitin were in higher proportion in the community during the epidemic event. Results of NGS and bioinformatics analyses offered a panorama of dynamic biodiversity and biological functioning of the community.

Chemotaxis may assist marine heterotrophic bacterial diazotrophs to find microzones suitable for N2 fixation in the pelagic ocean

Heterotrophic bacterial diazotrophs (HBDs) are ubiquitous in the pelagic ocean, where they have been predicted to carry out the anaerobic process of nitrogen fixation within low-oxygen microenvironments associated with marine pelagic particles. However, the mechanisms enabling particle colonization by HBDs are unknown. We hypothesized that HBDs use chemotaxis to locate and colonize suitable microenvironments, and showed that a cultivated marine HBD is chemotactic toward amino acids and phytoplankton-derived DOM. Using an in situ chemotaxis assay, we also discovered that diverse HBDs at a coastal site are motile and chemotactic toward DOM from various phytoplankton taxa and, indeed, that the proportion of diazotrophs was up to seven times higher among the motile fraction of the bacterial community compared to the bulk seawater community. Finally, three of four HBD isolates and 16 of 17 HBD metagenome assembled genomes, recovered from major ocean basins and locations along the Australian coast, each encoded >85% of proteins affiliated with the bacterial chemotaxis pathway. These results document the widespread capacity for chemotaxis in diverse and globally relevant marine HBDs. We suggest that HBDs could use chemotaxis to seek out and colonize low-oxygen microenvironments suitable for nitrogen fixation, such as those formed on marine particles. Chemotaxis in HBDs could therefore affect marine nitrogen and carbon biogeochemistry by facilitating nitrogen fixation within otherwise oxic waters, while also altering particle degradation and the efficiency of the biological pump.

Expansion of the global RNA virome reveals diverse clades of bacteriophages

High-throughput RNA sequencing offers broad opportunities to explore the Earth RNA virome. Mining 5,150 diverse metatranscriptomes uncovered >2.5 million RNA virus contigs. Analysis of >330,000 RNA-dependent RNA polymerases (RdRPs) shows that this expansion corresponds to a 5-fold increase of the known RNA virus diversity. Gene content analysis revealed multiple protein domains previously not found in RNA viruses and implicated in virus-host interactions. Extended RdRP phylogeny supports the monophyly of the five established phyla and reveals two putative additional bacteriophage phyla and numerous putative additional classes and orders. The dramatically expanded phylum Lenarviricota, consisting of bacterial and related eukaryotic viruses, now accounts for a third of the RNA virome. Identification of CRISPR spacer matches and bacteriolytic proteins suggests that subsets of picobirnaviruses and partitiviruses, previously associated with eukaryotes, infect prokaryotic hosts.

Integrating phylogenetic and functional data in microbiome studies

MOTIVATION

Microbiome functional data are frequently analyzed to identify associations between microbial functions (e.g. genes) and sample groups of interest. However, it is challenging to distinguish between different possible explanations for variation in community-wide functional profiles by considering functions alone. To help address this problem, we have developed POMS, a package that implements multiple phylogeny-aware frameworks to more robustly identify enriched functions.

RESULTS

The key contribution is an extended balance-tree workflow that incorporates functional and taxonomic information to identify functions that are consistently enriched in sample groups across independent taxonomic lineages. Our package also includes a workflow for running phylogenetic regression. Based on simulated data we demonstrate that these approaches more accurately identify gene families that confer a selective advantage compared with commonly used tools. We also show that POMS in particular can identify enriched functions in real-world metagenomics datasets that are potential targets of strong selection on multiple members of the microbiome.

AVAILABILITY AND IMPLEMENTATION

These workflows are freely available in the POMS R package at https://github.com/gavinmdouglas/POMS.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Microbial functional diversity across biogeochemical provinces in the central Pacific Ocean

Enzymes catalyze key reactions within Earth's life-sustaining biogeochemical cycles. Here, we use metaproteomics to examine the enzymatic capabilities of the microbial community (0.2 to 3 µm) along a 5,000-km-long, 1-km-deep transect in the central Pacific Ocean. Eighty-five percent of total protein abundance was of bacterial origin, with Archaea contributing 1.6%. Over 2,000 functional KEGG Ontology (KO) groups were identified, yet only 25 KO groups contributed over half of the protein abundance, simultaneously indicating abundant key functions and a long tail of diverse functions. Vertical attenuation of individual proteins displayed stratification of nutrient transport, carbon utilization, and environmental stress. The microbial community also varied along horizontal scales, shaped by environmental features specific to the oligotrophic North Pacific Subtropical Gyre, the oxygen-depleted Eastern Tropical North Pacific, and nutrient-rich equatorial upwelling. Some of the most abundant proteins were associated with nitrification and C1 metabolisms, with observed interactions between these pathways. The oxidoreductases nitrite oxidoreductase (NxrAB), nitrite reductase (NirK), ammonia monooxygenase (AmoABC), manganese oxidase (MnxG), formate dehydrogenase (FdoGH and FDH), and carbon monoxide dehydrogenase (CoxLM) displayed distributions indicative of biogeochemical status such as oxidative or nutritional stress, with the potential to be more sensitive than chemical sensors. Enzymes that mediate transformations of atmospheric gases like CO, CO₂, NO, methanethiol, and methylamines were most abundant in the upwelling region. We identified hot spots of biochemical transformation in the central Pacific Ocean, highlighted previously understudied metabolic pathways in the environment, and provided rich empirical data for biogeochemical models critical for forecasting ecosystem response to climate change.

Phylogenetically and catabolically diverse diazotrophs reside in deep-sea cold seep sediments

Microbially mediated nitrogen cycling in carbon-dominated cold seep environments remains poorly understood. So far anaerobic methanotrophic archaea (ANME-2) and their sulfate-reducing bacterial partners (SEEP-SRB1 clade) have been identified as diazotrophs in deep sea cold seep sediments. However, it is unclear whether other microbial groups can perform nitrogen fixation in such ecosystems. To fill this gap, we analyzed 61 metagenomes, 1428 metagenome-assembled genomes, and six metatranscriptomes derived from 11 globally distributed cold seeps. These sediments contain phylogenetically diverse nitrogenase genes corresponding to an expanded diversity of diazotrophic lineages. Diverse catabolic pathways were predicted to provide ATP for nitrogen fixation, suggesting diazotrophy in cold seeps is not necessarily associated with sulfate-dependent anaerobic oxidation of methane. Nitrogen fixation genes among various diazotrophic groups in cold seeps were inferred to be genetically mobile and subject to purifying selection. Our findings extend the capacity for diazotrophy to five candidate phyla (Altarchaeia, Omnitrophota, FCPU426, Caldatribacteriota and UBA6262), and suggest that cold seep diazotrophs might contribute substantially to the global nitrogen balance.

A consensus protocol for the recovery of mercury methylation genes from metagenomes

Mercury (Hg) methylation genes (hgcAB) mediate the formation of the toxic methylmercury and have been identified from diverse environments, including freshwater and marine ecosystems, Arctic permafrost, forest and paddy soils, coal-ash amended sediments, chlor-alkali plants discharges and geothermal springs. Here we present the first attempt at a standardized protocol for the detection, identification and quantification of hgc genes from metagenomes. Our Hg-MATE (Hg-cycling Microorganisms in Aquatic and Terrestrial Ecosystems) database, a catalogue of hgc genes, provides the most accurate information to date on the taxonomic identity and functional/metabolic attributes of microorganisms responsible for Hg methylation in the environment. Furthermore, we introduce "marky-coco", a ready-to-use bioinformatic pipeline based on de novo single-metagenome assembly, for easy and accurate characterization of hgc genes from environmental samples. We compared the recovery of hgc genes from environmental metagenomes using the marky-coco pipeline with an approach based on co-assembly of multiple metagenomes. Our data show similar efficiency in both approaches for most environments except those with high diversity (i.e., paddy soils) for which a co-assembly approach was preferred. Finally, we discuss the definition of true hgc genes and methods to normalize hgc gene counts from metagenomes.