Localized high abundance of Marine Group II archaea in the subtropical Pearl River Estuary: implications for their niche adaptation

Distribution of protoporphyrin IX during Prorocentrum donghaiense blooms and its relationship with particle-attached and free-living bacterial communities

Particle-attached (PA) and free-living (FL) bacterial communities are essential for nutrient cycles and metabolite production and serve as a food source in aquatic systems. However, our understanding of how biotic factors influence community interactions, co-occurrence patterns, and niche occupancy remains limited. This study investigated the influence of protoporphyrin IX (PPIX) on bacteria with different lifestyles during Prorocentrum donghaiense bloom. The findings revealed that PPIX distribution responded variably to changes in physicochemical parameters induced by red tide bloom. Large-sized or particle-attached (PA) phytoplankton (cell size >3 μm) were identified as the primary contributors to environmental PPIX, while small-sized plankton or free-living (FL) microorganisms (<3 μm) contributed less. In red tide-affected areas, PPIX and its derivatives were significantly more abundant than in non-red tide areas, indicating an increased demand for porphyrins by plankton during red tides. Additionally, the red tide also significantly influenced the preference of bacterial lineages for PA or FL lifestyles, highlighting a close interaction between bacteria with different lifestyles and PPIX levels. This study quantitatively analyzed the distribution of PPIX across different cell sizes in red tide and non-red tide marine environments, providing insights into microbial interactions and dynamics in changing ecosystems and offering a reference for using PPIX to predict red tide ecological disasters.

Metagenomic insights into microbial adaptation to the salinity gradient of a typical short residence-time estuary

BACKGROUND

Microbial adaptation to salinity has been a classic inquiry in the field of microbiology. It has been demonstrated that microorganisms can endure salinity stress via either the "salt-in" strategy, involving inorganic ion uptake, or the "salt-out" strategy, relying on compatible solutes. While these insights are mostly based on laboratory-cultured isolates, exploring the adaptive mechanisms of microorganisms within natural salinity gradient is crucial for gaining a deeper understanding of microbial adaptation in the estuarine ecosystem.

RESULTS

Here, we conducted metagenomic analyses on filtered surface water samples collected from a typical subtropical short residence-time estuary and categorized them by salinity into low-, intermediate-, and high-salinity metagenomes. Our findings highlighted salinity-driven variations in microbial community composition and function, as revealed through taxonomic and Clusters of Orthologous Group (COG) functional annotations. Through metagenomic binning, 127 bacterial and archaeal metagenome-assembled genomes (MAGs) were reconstructed. These MAGs were categorized as stenohaline-specific to low-, intermediate-, or high-salinity-based on the average relative abundance in one salinity category significantly exceeding those in the other two categories by an order of magnitude. Those that did not meet this criterion were classified as euryhaline, indicating a broader range of salinity tolerance. Applying the Boruta algorithm, a machine learning-based feature selection method, we discerned important genomic features from the stenohaline bacterial MAGs. Of the total 12,162 COGs obtained, 40 were identified as important features, with the "inorganic ion transport and metabolism" COG category emerging as the most prominent. Furthermore, eight COGs were implicated in microbial osmoregulation, of which four were related to the "salt-in" strategy, three to the "salt-out" strategy, and one to the regulation of water channel activity. COG0168, annotated as the Trk-type K+ transporter related to the "salt-in" strategy, was ranked as the most important feature. The relative abundance of COG0168 was observed to increase with rising salinity across metagenomes, the stenohaline strains, and the dominant Actinobacteriota and Proteobacteria phyla.

CONCLUSIONS

We demonstrated that salinity exerts influences on both the taxonomic and functional profiles of the microbial communities inhabiting the estuarine ecosystem. Our findings shed light on diverse salinity adaptation strategies employed by the estuarine microbial communities, highlighting the crucial role of the "salt-in" strategy mediated by Trk-type K+ transporters for microorganisms thriving under osmotic stress in the short residence-time estuary. Video Abstract.

Gene inversion led to the emergence of brackish archaeal heterotrophs in the aftermath of the Cryogenian Snowball Earth

Land-ocean interactions greatly impact the evolution of coastal life on earth. However, the ancient geological forces and genetic mechanisms that shaped evolutionary adaptations and allowed microorganisms to inhabit coastal brackish waters remain largely unexplored. In this study, we infer the evolutionary trajectory of the ubiquitous heterotrophic archaea Poseidoniales (Marine Group II archaea) presently occurring across global aquatic habitats. Our results show that their brackish subgroups had a single origination, dated to over 600 million years ago, through the inversion of the magnesium transport gene corA that conferred osmotic-stress tolerance. The subsequent loss and gain of corA were followed by genome-wide adjustment, characterized by a general two-step mode of selection in microbial speciation. The coastal family of Poseidoniales showed a rapid increase in the evolutionary rate during and in the aftermath of the Cryogenian Snowball Earth (∼700 million years ago), possibly in response to the enhanced phosphorus supply and the rise of algae. Our study highlights the close interplay between genetic changes and ecosystem evolution that boosted microbial diversification in the Neoproterozoic continental margins, where the Cambrian explosion of animals soon followed.

Archaeal blooms and busts in an estuarine time series

Coastal bays, such as Delaware Bay, are highly productive, ecologically important transitions between rivers and the coastal ocean. They offer opportunities to investigate archaeal assemblages across seasons, with the exchange of water masses that occurs with tidal cycles, and in the context of variable organic matter quality. For a year-long estuarine, size-fractionated time series, we used amplicon sequencing, chemical measurements, and qPCR to follow archaeal groups through the seasons. We detected seasonally high abundances of Marine Group II archaea in summer months which correlate with indicators of phytoplankton production, although not phytoplankton biomass. Although previous studies have reported associations between Marine Group II archaea and particles, here they are almost entirely found in very small particles (0.22-0.7 μm), suggesting they are free-living cells. Populations of Nitrososphaeria did not vary with particle size or environmental conditions. Methanogens were significant fractions of archaeal sequences in large particles at low tide during winter months. Contrary to expectations, Nanoarchaeia were found predominantly in the free-living fraction despite the previous observation that they require an association with hosts. These results underscore the utility of time series studies in shallow, tidally mixed estuarine environments that capture variable conditions for understanding the ecology and biogeochemistry of planktic archaea.

Oxygen-driven divergence of marine group II archaea reflected by transitions of superoxide dismutases

IMPORTANCE

Reactive oxygen species (ROS), including superoxide anion, is a series of substances that cause oxidative stress for all organisms. Marine group II (MGII) archaea are mainly live in the surface seawater and exposed to considerable ROS. Therefore, it is important to understand the antioxidant capacity of MGII. Our research found that Fe/Mn- superoxide dismutase (Fe/MnSOD) may be more suitable for MGII to resist oxidative damage, and the changes in oxygen concentrations and SOD metallic cofactors play an important role in the selection of SOD by the 17 clades of MGII, which in turn affects the species differentiation of MGII. Overall, this study provides insight into the co-evolutionary history of these uncultivated marine archaea with the earth system.

Exploring the Archaeal Virosphere by Metagenomics

During the past decade, environmental research has demonstrated that archaea are abundant and widespread in nature and play important ecological roles at a global scale. Currently, however, the majority of archaeal lineages cannot be cultivated under laboratory conditions and are known exclusively or nearly exclusively through metagenomics. A similar trend extends to the archaeal virosphere, where isolated representatives are available for a handful of model archaeal virus-host systems. Viral metagenomics provides an alternative way to circumvent the limitations of culture-based virus discovery and offers insight into the diversity, distribution, and environmental impact of uncultured archaeal viruses. Presently, metagenomics approaches have been successfully applied to explore the viromes associated with various lineages of extremophilic and mesophilic archaea, including Asgard archaea (Asgardarchaeota), ANME-1 archaea (Methanophagales), thaumarchaea (Nitrososphaeria), altiarchaea (Altiarchaeota), and marine group II archaea (Poseidoniales). Here, we provide an overview of methods widely used in archaeal virus metagenomics, covering metavirome preparation, genome annotation, phylogenetic and phylogenomic analyses, and archaeal host assignment. We hope that this summary will contribute to further exploration and characterization of the enigmatic archaeal virome lurking in diverse environments.

Seasonal dynamics of bacterial community and co-occurrence with eukaryotic phytoplankton in the Pearl River Estuary

In this study, we investigated the taxonomic composition of the bacteria and phytoplankton communities in the Pearl River Estuary (PRE) through Illumina sequencing of the V3-V4 region of the 16 S rRNA gene. Furthermore, their relationships as well as recorded environmental variables were explored by co-occurrence networks. Bacterial community composition was different in two size fractions, as well as along the salinity gradient across two seasons. Free-living (FL) communities were dominated by pico-sized Cyanobacteria (Synechococcus CC9902) while Exiguobacterium, Halomonas and Pseudomonas were predominantly associated with particle-associated (PA) lifestyle, and Cyanobium PCC-6307 exhibited seasonal shifts in lifestyles in different seasons. In wet season, bacterial community composition was characterized by abundance of Cyanobacteria, Actinobacteria, and Bacteroidetes, which were tightly linked with high riverine inflow. While in dry season, Proteobacteria increased in prevalence, especially for Psychrobacter, NOR5/OM60 clade and Pseudomonas, which were thrived in lower water temperature and higher salinity. Moreover, we discovered that differences between PA and FL composition were more significant in the wet season than in the dry season, which may be due to better nutritional conditions of particles (indicated by POC%) in the wet season and then attract more diverse PA populations. Based on the analysis of plastidial 16 S rRNA genes, abundant small-sized mixotrophic phytoplankton (Dinophyceae, Euglenida and Haptophyta) were identified in the PRE. The complexity of co-occurrence network increased from FL to PA fractions in both seasons, which suggested that suspended particles can provide ecological niches for particle-associated colonizers contributing to the maintenance of a more stable community structure. In addition, the majority of phytoplankton species exhibited positive co-occurrences with both other phytoplankton species and bacterial counterparts, indicating the mutual cooperation between phytoplankton assemblages and specific bacterial populations e likely benefited from phytoplankton-derived organic compounds. This study enhances our understanding of the seasonal and spatial dynamics of bacterial communities and their potential relationship with phytoplankton assembly in estuarine waters.

Diversity and salinity adaptations of ammonia oxidizing archaea in three estuaries of China

Ammonia-oxidizing archaea (AOA) are ubiquitously found in diverse habitats and play pivotal roles in the nitrogen and carbon cycle, especially in estuarine and coastal environments. Despite the fact that the diversity and distribution of AOA are thought to be tightly linked to habitats, little is known about the relationship that underpins their genomic traits, adaptive potentials, and ecological niches. Here, we have characterized and compared the AOA community in three estuaries of China using metagenomics. AOA were the dominant ammonia oxidizers in the three estuaries. Through phylogenetic analyses, five major AOA groups were identified, including the Nitrosomarinus-like, Nitrosopumilus-like, Aestuariumsis-like, Nitrosarchaeum-like, and Nitrosopelagicus-like groups. Statistical analyses showed that the aquatic and sedimentary AOA communities were mainly influenced by spatial factors (latitude and water depth) and environmental factors (salinity, pH, and dissolved oxygen) in estuaries, respectively. Compared to AOA dwelling in terrestrial and marine habitats, estuarine AOA encoded more genes involved in glucose and amino acid metabolism, transport systems, osmotic control, and cell motility. The low proteome isoelectric points (pI), high content of acidic amino acids, and the presence of potassium ion and mechanosensitive channels suggest a "salt-in" strategy for estuarine AOA to counteract high osmolarity in their surroundings. Our findings have indicated potential adaptation strategies and highlighted their importance in the estuarine nitrogen and carbon cycles. KEY POINTS: • Spatial and environmental factors influence water and sediment AOA respectively. • Estuarine AOA share low proteome isoelectric value and high acid amino acids content. • AOA adaptation to estuaries is likely resulted from their unique genomic features.

Distinct patterns of distribution, community assembly and cross-domain co-occurrence of planktonic archaea in four major estuaries of China

BACKGROUND

Archaea are key mediators of estuarine biogeochemical cycles, but comprehensive studies comparing archaeal communities among multiple estuaries with unified experimental protocols during the same sampling periods are scarce. Here, we investigated the distribution, community assembly, and cross-domain microbial co-occurrence of archaea in surface waters across four major estuaries (Yellow River, Yangtze River, Qiantang River, and Pearl River) of China cross climatic zones (~ 1,800 km) during the winter and summer cruises.

RESULTS

The relative abundance of archaea in the prokaryotic community and archaeal community composition varied with estuaries, seasons, and stations (reflecting local environmental changes such as salinity). Archaeal communities in four estuaries were overall predominated by ammonia-oxidizing archaea (AOA) (aka. Marine Group (MG) I; primarily Nitrosopumilus), while the genus Poseidonia of Poseidoniales (aka. MGII) was occasionally predominant in Pearl River estuary. The cross-estuary dispersal of archaea was largely limited and the assembly mechanism of archaea varied with estuaries in the winter cruise, while selection governed archaeal assembly in all estuaries in the summer cruise. Although the majority of archaea taxa in microbial networks were peripherals and/or connectors, extensive and distinct cross-domain associations of archaea with bacteria were found across the estuaries, with AOA as the most crucial archaeal group. Furthermore, the expanded associations of MGII taxa with heterotrophic bacteria were observed, speculatively indicating the endogenous demand for co-processing high amount and diversity of organic matters in the estuarine ecosystem highly impacted by terrestrial/anthropogenic input, which is worthy of further study.

CONCLUSIONS

Our results highlight the lack of common patterns in the dynamics of estuarine archaeal communities along the geographic gradient, expanding the understanding of roles of archaea in microbial networks of this highly dynamic ecosystem.

A moderately thermophilic origin of a novel family of marine group II euryarchaeota from deep ocean

Marine group II (MGII) is the most abundant planktonic heterotrophic archaea in the ocean. The evolutionary history of MGII archaea is elusive. In this study, 13 new MGII metagenome-assembled genomes were recovered from surface to the hadal zone in Challenger Deep of the Mariana Trench; four of them from the deep ocean represent a novel group. The optimal growth temperature (OGT) of the common ancestor of MGII has been estimated to be at about 60°C and OGTs of MGIIc, MGIIb, and MGIIa at 47°C-50ºC, 37°C-44ºC, and 30°C-37ºC, respectively, suggesting the adaptation of these species to different temperatures during evolution. The estimated OGT range of MGIIc was supported by experimental measurements of cloned β-galactosidase that showed optimal enzyme activity around 50°C. These results indicate that MGIIc may have originated from a common ancestor that lived in warm or even hot marine environment, such as hydrothermal vents.

Deterministic processes dominate archaeal community assembly from the Pearl River to the northern South China Sea

Archaea play a significant role in the biogeochemical cycling of nutrients in estuaries. However, comprehensive researches about their assembly processes remain notably insufficient. In this study, we systematically examined archaeal community dynamics distinguished between low-salinity and high-salinity groups in water and surface sediments over a 600-kilometer range from the upper Pearl River (PR) to the northern South China Sea (NSCS). Neutral community model analysis together with null model analysis showed that their C-score values were greater than 2, suggesting that deterministic processes could dominate the assembly of those planktonic or benthic archaeal communities at both the low-salinity and high-salinity sites. And deterministic processes contributed more in the low-salinity than high-salinity environments from the PR to the NSCS. Furthermore, through the co-occurrence network analysis, we found that the archaeal communities in the low-salinity groups possessed closer interactions and higher proportions of negative interactions than those in the high-salinity groups, which might be due to the larger environmental heterogeneities reflected by the nutrient concentrations of those low-salinity samples. Collectively, our work systematically investigated the composition and co-occurrence networks of archaeal communities in water as well as sediments from the PR to the NSCS, yielding new insights into the estuary's archaeal community assembly mechanisms.

Diversity, distribution, and functional potentials of magroviruses from marine and brackish waters

Marine group II (MGII) archaea (Ca. Poseidoniales) are among the most abundant microbes in global oceanic surface waters and play an important role in driving marine biogeochemical cycles. Magroviruses - the viruses of MGII archaea have been recently found to occur ubiquitously in surface ocean. However, their diversity, distribution, and potential ecological functions in coastal zones especially brackish waters are unknown. Here we obtained 234 non-redundant magroviral genomes from brackish surface waters by using homology searches for viral signature proteins highlighting the uncovered vast diversity of this novel viral group. Phylogenetic analysis based on these brackish magroviruses along with previously reported marine ones identified six taxonomic groups with close evolutionary connection to both haloviruses and the viruses of Marine Group I archaea. Magroviruses were present abundantly both in brackish and open ocean samples with some showing habitat specification and others having broad spectrums of distribution between different habitats. Genome annotation suggests they may be involved in regulating multiple metabolic pathways of MGII archaea. Our results uncover the previously overlooked diversity and ecological potentials of a major archaeal virial group in global ocean and brackish waters and shed light on the cryptic evolutionary history of archaeal viruses.

Non-negligible roles of archaea in coastal carbon biogeochemical cycling

Coastal zones are among the world's most productive ecosystems. They store vast amounts of organic carbon, as 'blue carbon' reservoirs, and impact global climate change. Archaeal communities are integral components of coastal microbiomes but their ecological roles are often overlooked. However, archaeal diversity, metabolism, evolution, and interactions, revealed by recent studies using rapidly developing cutting-edge technologies, place archaea as important players in coastal carbon biogeochemical cycling. We here summarize the latest advances in the understanding of archaeal carbon cycling processes in coastal ecosystems, specifically, archaeal involvement in CO₂ fixation, organic biopolymer transformation, and methane metabolism. We also showcase the potential to use of archaeal communities to increase carbon sequestration and reduce methane production, with implications for mitigating climate change.

Seasonal microbial dynamics in the ocean inferred from assembled and unassembled data: a view on the unknown biosphere

In environmental metagenomic experiments, a very high proportion of the microbial sequencing data (> 70%) remains largely unexploited because rare and closely related genomes are missed in short-read assemblies. The identity and the potential metabolisms of a large fraction of natural microbial communities thus remain inaccessible to researchers. The purpose of this study was to explore the genomic content of unassembled metagenomic data and test their level of novelty. We used data from a three-year microbial metagenomic time series of the NW Mediterranean Sea, and conducted reference-free and database-guided analysis. The results revealed a significant genomic difference between the assembled and unassembled reads. The unassembled reads had a lower mean identity against public databases, and fewer metabolic pathways could be reconstructed. In addition, the unassembled fraction presented a clear temporal pattern, unlike the assembled ones, and a specific community composition that was similar to the rare communities defined by metabarcoding using the 16S rRNA gene. The rare gene pool was characterised by keystone bacterial taxa, and the presence of viruses, suggesting that viral lysis could maintain some taxa in a state of rarity. Our study demonstrates that unassembled metagenomic data can provide important information on the structure and functioning of microbial communities.

Comparative Genomics of Thaumarchaeota From Deep-Sea Sponges Reveal Their Niche Adaptation

Thaumarchaeota account for a large portion of microbial symbionts in deep-sea sponges and are even dominant in some cases. In this study, we investigated three new sponge-associated Thaumarchaeota from the deep West Pacific Ocean. Thaumarchaeota were found to be the most dominant phylum in this sponge by both prokaryotic 16S rRNA amplicons and metagenomic sequencing. Fifty-seven published Thaumarchaeota genomes from sponges and other habitats were included for genomic comparison. Similar to shallow sponge-associated Thaumarchaeota, those Thaumarchaeota in deep-sea sponges have extended genome sizes and lower coding density compared with their free-living lineages. Thaumarchaeota in deep-sea sponges were specifically enriched in genes related to stress adapting, symbiotic adhesion and stability, host–microbe interaction and protein transportation. The genes involved in defense mechanisms, such as the restriction-modification system, clustered regularly interspaced short palindromic repeats (CRISPR)/Cas system, and toxin-antitoxin system were commonly enriched in both shallow and deep sponge-associated Thaumarchaeota. Our study demonstrates the significant effects of both depth and symbiosis on forming genomic characteristics of Thaumarchaeota, and provides novel insights into their niche adaptation in deep-sea sponges.

Variations in Marine Bacterial and Archaeal Communities during an Ulva prolifera Green Tide in Coastal Qingdao Areas

Green tides caused by Ulva prolifera occur annually in the Yellow Sea, potentially influencing the marine microorganisms. Here, we focused on the variations in marine bacterial and archaeal communities during an U. prolifera green tide in coastal Qingdao areas with Illumina high-throughput sequencing analysis. Our results revealed that the diversity and structure of bacterial and archaeal communities, as well as the organization and structure of microbial co-occurrence networks, varied during the green tide. The decline phase may be favorable to the bacterial and archaeal diversity and richness. The bacterial community, as well as the archaeal community, showed clear variations between the outbreak and decline phases. A simpler and less connected microbial co-occurrence network was observed during the outbreak phase compared with the decline phase. Flavobacteriales and Rhodobacterales separately dominated the bacterial community during the outbreak and decline phase, and Marine Group II (MGII) dominated the archaeal community during the green tide. Combined with microbial co-occurrence network analysis, Flavobacteriales, Rhodobacterales and MGII may be important organisms during the green tide. Temperature, chlorophyll a content and salinity may have an important impact on the variations in bacterial and archaeal communities during the green tide.

A holistic genome dataset of bacteria, archaea and viruses of the Pearl River estuary

Estuaries are one of the most important coastal ecosystems. While microbiomes and viromes have been separately investigated in some estuaries, few studies holistically deciphered the genomes and connections of viruses and their microbial hosts along an estuarine salinity gradient. Here we applied deep metagenomic sequencing on microbial and viral communities in surface waters of the Pearl River estuary, one of China’s largest estuaries with strong anthropogenic impacts. Overall, 1,205 non-redundant prokaryotic genomes with ≥50% completeness and ≤10% contamination, and 78,502 non-redundant viral-like genomes were generated from samples of three size fractions and five salinity levels. Phylogenomic analysis and taxonomy classification show that majority of these estuarine prokaryotic and viral genomes are novel at species level according to public databases. Potential connections between the microbial and viral populations were further investigated by host-virus matching. These combined microbial and viral genomes provide an important complement of global marine genome datasets and should greatly facilitate our understanding of microbe-virus interactions, evolution and their implications in estuarine ecosystems.

Measurement(s)	bacteria • Archaea • viruses
Technology Type(s)	Shotgun Sequencing
Sample Characteristic - Organism	estuary metagenome
Sample Characteristic - Environment	subtropical estuarine
Sample Characteristic - Location	Pearl river estuary

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.17139998

Interactions Between Marine Group II Archaea and Phytoplankton Revealed by Population Correlations in the Northern Coast of South China Sea

Marine Group II (MGII) archaea (Poseidoniales) are the most abundant surface marine planktonic archaea and are widely distributed in both coastal and pelagic waters. The factors affecting their distribution and activity are poorly understood. MGII archaea have the metabolic potential to utilize algae-derived organic matter and are frequently observed in high abundance during or following phytoplankton blooms, suggesting that they are key players of the marine food web. In this study, we studied interactions between MGII archaea and the diverse taxa of phytoplankton in the northern coast of South China Sea. Non-metric multidimensional scaling and cluster analyses demonstrated distinct MGII community patterns in the Pearl River plume (PRP) and the open regions of the northern South China Sea (ONSCS), with MGIIb dominating the former and MGIIa and MGIIb showing remarkable variations in the latter for the same sampling season. Nevertheless, positive correlations (Pearson correlation: R > 0.8 and P < 0.01) in absolute abundances of ribosomal RNA (rRNA)-derived complementary DNA and rRNA genes from network analyses were found between MGII archaea and phytoplankton (cyanobacteria, haptophytes, and stramenopiles in both PRP and ONSCS) among different particle size fractions, indicating their intrinsic relationships under changing environmental conditions. The results of this study may shed light on the multiple interactions between co-existing species in the micro-niches of different oceanic regions.

Genome-Resolved Metagenomic Insights into Massive Seasonal Ammonia-Oxidizing Archaea Blooms in San Francisco Bay

Ammonia-oxidizing archaea (AOA) are key for the transformation of ammonia to oxidized forms of nitrogen in aquatic environments around the globe, including nutrient-rich coastal and estuarine waters such as San Francisco Bay (SFB). Using metagenomics and 16S rRNA gene amplicon libraries, we found that AOA are more abundant than ammonia-oxidizing bacteria (AOB) and nitrite-oxidizing bacteria (NOB), except in the freshwater stations in SFB. In South SFB, we observed recurrent AOA blooms of “Candidatus Nitrosomarinus catalina” SPOT01-like organisms, which account for over 20% of 16S rRNA gene amplicons in both surface and bottom waters and co-occur with weeks of high nitrite concentrations (>10 μM) in the oxic water column. We observed pronounced nitrite peaks occurring in the autumn for 7 of the last 9 years (2012 to 2020), suggesting that seasonal AOA blooms are common in South SFB. We recovered two high-quality AOA metagenome-assembled genomes (MAGs), including a Nitrosomarinus-like genome from the South SFB bloom and another Nitrosopumilus genome originating from Suisun Bay in North SFB. Both MAGs cluster with genomes from other estuarine/coastal sites. Analysis of Nitrosomarinus-like genomes show that they are streamlined, with low GC content and high coding density, and harbor urease genes. Our findings support the unique niche of Nitrosomarinus-like organisms which dominate coastal/estuarine waters and provide insights into recurring AOA blooms in SFB.

IMPORTANCE Ammonia-oxidizing archaea (AOA) carry out key transformations of ammonia in estuarine systems such as San Francisco Bay (SFB)—the largest estuary on the west coast of North America—and play a significant role in both local and global nitrogen cycling. Using metagenomics and 16S rRNA gene amplicon libraries, we document a massive, recurrent AOA bloom in South SFB that co-occurs with months of high nitrite concentrations in the oxic water column. Our study is the first to generate metagenome-assembled genomes (MAGs) from SFB, and through this process we recovered two high-quality AOA MAGs, one of which originated from bloom samples. These AOA MAGs yield new insight into the Nitrosopumilus and Nitrosomarinus-like lineages and their potential niches in coastal and estuarine systems. Nitrosomarinus-like AOA are abundant in coastal regions around the globe, and we highlight the common occurrence of urease genes, low GC content, and range of salinity tolerances within this lineage.

Diversity Distribution, Driving Factors and Assembly Mechanisms of Free-Living and Particle-Associated Bacterial Communities at a Subtropical Marginal Sea

Free-living (FL) and particle-associated (PA) bacterioplankton communities play critical roles in biogeochemical cycles in the ocean. However, their community composition, assembly process and functions in the continental shelf and slope regions are poorly understood. Based on 16S rRNA gene amplicon sequencing, we investigated bacterial communities’ driving factors, assembly processes and functional potentials at a subtropical marginal sea. The bacterioplankton community showed specific distribution patterns with respect to lifestyle (free living vs. particle associated), habitat (slope vs. shelf) and depth (surface vs. DCM and Bottom). Salinity and water temperature were the key factors modulating turnover in the FL community, whereas nitrite, silicate and phosphate were the key factors for the PA community. Model analyses revealed that stochastic processes outweighed deterministic processes and had stronger influences on PA than FL. Homogeneous selection (Hos) was more responsible for the assembly and turnover of FL, while drift and dispersal limitation contributed more to the assembly of PA. Importantly, the primary contributor to Hos in PA was Gammaproteobacteria:Others, whereas that in FL was Cyanobacteria:Bin6. Finally, the PICRUSt2 analysis indicated that the potential metabolisms of carbohydrates, cofactors, amino acids, terpenoids, polyketides, lipids and antibiotic resistance were markedly enriched in PA than FL.

Characterization of the Exometabolome of Nitrosopumilus maritimus SCM1 by Liquid Chromatography-Ion Mobility Mass Spectrometry

Marine Thaumarchaeota (formerly known as the marine group I archaea) have received much research interest in recent years since these chemolithoautotrophic organisms are abundant in the subsurface ocean and oxidize ammonium to nitrite, which makes them a major contributor to the marine carbon and nitrogen cycles. However, few studies have investigated the chemical composition of their exometabolome and their contributions to the pool of dissolved organic matter (DOM) in seawater. This study exploits the recent advances in ion mobility mass spectrometry (IM-MS) and integrates this instrumental capability with bioinformatics to reassess the exometabolome of a model ammonia-oxidizing archaeon, Nitrosopumilus maritimus strain SCM1. Our method has several advantages over the conventional approach using an Orbitrap or ion cyclotron resonance mass analyzer and allows assignments or annotations of spectral features to known metabolites confidently and indiscriminately, as well as distinction of biological molecules from background organics. Consistent with the results of a previous report, the SPE-extracted exometabolome of N. maritimus is dominated by biologically active nitrogen-containing metabolites, in addition to peptides secreted extracellularly. Cobalamin and associated intermediates, including α-ribazole and α-ribazole 5'-phosphate, are major components of the SPE-extracted exometabolome of N. maritimus. This supports the proposition that Thaumarchaeota have the capacity of de novo biosynthesizing cobalamin. Other biologically significant metabolites, such as agmatidine and medicagenate, predicted by genome screening are also detected, which indicates that Thaumarchaeota have remarkable metabolic potentials, underlining their importance in driving elemental cycles critical to biological processes in the ocean.

Network analysis methods for studying microbial communities: A mini review

Microorganisms including bacteria, fungi, viruses, protists and archaea live as communities in complex and contiguous environments. They engage in numerous inter- and intra- kingdom interactions which can be inferred from microbiome profiling data. In particular, network-based approaches have proven helpful in deciphering complex microbial interaction patterns. Here we give an overview of state-of-the-art methods to infer intra-kingdom interactions ranging from simple correlation- to complex conditional dependence-based methods. We highlight common biases encountered in microbial profiles and discuss mitigation strategies employed by different tools and their trade-off with increased computational complexity. Finally, we discuss current limitations that motivate further method development to infer inter-kingdom interactions and to robustly and comprehensively characterize microbial environments in the future.

Seasonality of archaeal proteorhodopsin and associated Marine Group IIb ecotypes (Ca. Poseidoniales) in the North Western Mediterranean Sea

The Archaea Marine Group II (MGII) is widespread in the world's ocean where it plays an important role in the carbon cycle. Despite recent discoveries on the group's metabolisms, the ecology of this newly proposed order (Candidatus Poseidoniales) remains poorly understood. Here we used a combination of time-series metagenome-assembled genomes (MAGs) and high-frequency 16S rRNA data from the NW Mediterranean Sea to test if the taxonomic diversity within the MGIIb family (Candidatus Thalassarchaeaceae) reflects the presence of different ecotypes. The MAGs' seasonality revealed a MGIIb family composed of different subclades that have distinct lifestyles and physiologies. The vitamin metabolisms were notably different between ecotypes with, in some, a possible link to sunlight's energy. Diverse archaeal proteorhodopsin variants, with unusual signature in key amino acid residues, had distinct seasonal patterns corresponding to changing day length. In addition, we show that in summer, archaea, as opposed to bacteria, disappeared completely from surface waters. Our results shed light on the diversity and the distribution of the euryarchaeotal proteorhodopsin, and highlight that MGIIb is a diverse ecological group. The work shows that time-series based studies of the taxonomy, seasonality, and metabolisms of marine prokaryotes is critical to uncover their diverse role in the ocean.

In-depth Spatiotemporal Characterization of Planktonic Archaeal and Bacterial Communities in North and South San Francisco Bay

Despite being the largest estuary on the west coast of North America, no in-depth survey of microbial communities in San Francisco Bay (SFB) waters currently exists. In this study, we analyze bacterioplankton and archaeoplankton communities at several taxonomic levels and spatial extents (i.e., North versus South Bay) to reveal patterns in alpha and beta diversity. We assess communities using high-throughput sequencing of the 16S rRNA gene in 177 water column samples collected along a 150-km transect over a 2-year monthly time-series. In North Bay, the microbial community is strongly structured by spatial salinity changes while in South Bay seasonal variations dominate community dynamics. Along the steep salinity gradient in North Bay, we find that operational taxonomic units (OTUs; 97% identity) have higher site specificity than at coarser taxonomic levels and turnover ("species" replacement) is high, revealing a distinct brackish community (in oligo-, meso-, and polyhaline samples) from fresh and marine end-members. At coarser taxonomic levels (e.g., phylum, class), taxa are broadly distributed across salinity zones (i.e., present/abundant in a large number of samples) and brackish communities appear to be a mix of fresh and marine communities. We also observe variations in brackish communities between samples with similar salinities, likely related to differences in water residence times between North and South Bay. Throughout SFB, suspended particulate matter is positively correlated with richness and influences changes in beta diversity. Within several abundant groups, including the SAR11 clade (comprising up to 30% of reads in a sample), OTUs appear to be specialized to a specific salinity range. Some other organisms also showed pronounced seasonal abundance, including Synechococcus, Ca. Actinomarina, and Nitrosopumilus-like OTUs. Overall, this study represents the first in-depth spatiotemporal survey of SFB microbial communities and provides insight into how planktonic microorganisms have specialized to different niches along the salinity gradient.

Transcriptional activity differentiates families of Marine Group II Euryarchaeota in the coastal ocean

Marine Group II Euryarchaeota (Candidatus Poseidoniales), abundant but yet-uncultivated members of marine microbial communities, are thought to be (photo)heterotrophs that metabolize dissolved organic matter (DOM), such as lipids and peptides. However, little is known about their transcriptional activity. We mapped reads from a metatranscriptomic time series collected at Sapelo Island (GA, USA) to metagenome-assembled genomes to determine the diversity of transcriptionally active Ca. Poseidoniales. Summer metatranscriptomes had the highest abundance of Ca. Poseidoniales transcripts, mostly from the O1 and O3 genera within Ca. Thalassarchaeaceae (MGIIb). In contrast, transcripts from fall and winter samples were predominantly from Ca. Poseidoniaceae (MGIIa). Genes encoding proteorhodopsin, membrane-bound pyrophosphatase, peptidase/proteases, and part of the ß-oxidation pathway were highly transcribed across abundant genera. Highly transcribed genes specific to Ca. Thalassarchaeaceae included xanthine/uracil permease and receptors for amino acid transporters. Enrichment of Ca. Thalassarchaeaceae transcript reads related to protein/peptide, nucleic acid, and amino acid transport and metabolism, as well as transcript depletion during dark incubations, provided further evidence of heterotrophic metabolism. Quantitative PCR analysis of South Atlantic Bight samples indicated consistently abundant Ca. Poseidoniales in nearshore and inshore waters. Together, our data suggest that Ca. Thalassarchaeaceae are important photoheterotrophs potentially linking DOM and nitrogen cycling in coastal waters.

Marine Group-II archaea dominate particle-attached as well as free-living archaeal assemblages in the surface waters of Kongsfjorden, Svalbard, Arctic Ocean

Marine archaea are a significant component of the global oceanic ecosystems, including the polar oceans. However, only a few attempts have been made to study archaea in the high Arctic fjords. Given the importance of Archaea in carbon and nitrogen cycling, it is imperative to explore their diversity and community composition in the high Arctic fjords, such as Kongsfjorden (Svalbard). In the present study, we evaluated archaeal diversity and community composition in the size-fractionated microbial population, viz-a-viz free-living (FL; 0.2-3 μm) and particle-attached (PA; > 3 μm) using archaeal V3-V4 16S rRNA gene amplicon sequencing. Our results indicate that the overall archaeal community in the surface water of Kongsfjorden was dominated by the members of the marine group-II (MGII) archaea, followed by the MGI group members, including Nitrosopumilaceae and Nitrososphaeraceae. Although a clear niche partitioning between PA and FL archaeal communities was not observed, 2 OTUs among 682 OTUs, and 3 ASVs out of 1932 ASVs were differentially abundant among the fractions. OTU001/ASV0002, classified as MGIIa, was differentially abundant in the PA fraction. OTU006/ASV0006/ASV0010 affiliated with MGIIb were differentially abundant in the FL fraction. Particulate organic nitrogen and C:N ratio were the most significant variables (P < 0.05) explaining the observed variation in the FL and PA archaeal communities, respectively. These results indicate an exchange between archaeal communities or a generalist lifestyle switching between FL and PA fractions. Besides, the particles' elemental composition (carbon and nitrogen) seems to play an essential role in shaping the PA archaeal communities in the surface waters of Kongsfjorden.

Spatial and seasonal variations of sediment bacterial communities in a river-bay system in South China

River-bay systems are transitional areas that hold important roles in biogeochemical processes between continents and oceans. However, composition and structure of microbial communities shaped by such environments have not been clear yet. In this study, we used high-throughput sequencing of 16S rRNA genes to analyze the diversity and composition of sediment bacterial communities from the Shenzhen river-bay system during dry and wet seasons. The results showed that sediment bacterial community structure was varied according to habitats (river vs. estuary) and seasons (wet season vs. dry season). The alpha diversity of sediment bacterial community was significantly higher in the dry season than in the wet season, while no significant difference of alpha diversity was found between river and estuary. Neutral community model revealed a significant influence of stochastic processes on sediment bacterial community assembly, especially in the wet season. However, the beta nearest-taxon index indicated that deterministic processes were more responsible for the assembly of sediment bacterial community. Additionally, redundancy analysis suggested strong links between sediment bacterial communities and environmental factors in Shenzhen river-bay system, with the environmental factors explaining 63.5% of the bacterial community variation. Specifically, NH₄⁺, pH, and salinity were the three most important contributing factors that shaped the sediment bacterial communities. Overall, this study provides a valuable reference to get insights into the spatiotemporal pattern of sediment bacterial communities in a typical river-bay system. KEY POINTS: • Stochastic processes contribute sediment bacterial community assembly. • Deterministic processes dominate sediment bacterial community assembly. • Environmental factors shape sediment bacterial communities.

Distinct Distribution of Archaea From Soil to Freshwater to Estuary: Implications of Archaeal Composition and Function in Different Environments

In addition to inhabiting extreme territories, Archaea are widely distributed in common environments spanning from terrestrial to aquatic environments. This study investigated and compared archaeal community structures from three different habitats (representing distinct environments): agriculture soils (from farming system trials FST, PA, United States), freshwater biofilms (from White Clay Creek, PA, United States), and estuary water (Chesapeake Bay, United States). High-throughput sequencing of 16S rRNA genes indicated that Thaumarchaeota, Euryarchaeota, Nanoarchaeota, Crenarchaeota, and Diapherotrites were the commonly found dominant phyla across these three environments. Similar to Bacteria, distinct community structure and distribution patterns for Archaea were observed in soils vs. freshwater vs. estuary. However, the abundance, richness, evenness, and diversity of archaeal communities were significantly greater in soils than it was in freshwater and estuarine environments. Indicator species (or amplicon sequence variants, ASVs) were identified from different nitrogen and carbon cycling archaeal groups in soils (Nitrososphaerales, Nitrosotaleales, Nitrosopumilales, Methanomassiliicoccales, Lainarchaeales), freshwater biofilms (Methanobacteria, Nitrososphaerales) and Chesapeake Bay (Marine Group II, Nitrosopumilales), suggesting the habitat-specificity of their biogeochemical contributions to different environments. Distinct functional aspects of Archaea were also confirmed by functional predictions (PICRUSt2 analysis). Further, co-occurrence network analysis indicated that only soil Archaea formed stable modules. Keystone species (ASVs) were identified mainly from Methanomassiliicoccales, Nitrososphaerales, Nitrosopumilales. Overall, these results indicate a strong habitat-dependent distribution of Archaea and their functional partitions within the local environments.

Simulation of Enhanced Growth of Marine Group II Euryarchaeota From the Deep Chlorophyll Maximum of the Western Pacific Ocean: Implication for Upwelling Impact on Microbial Functions in the Photic Zone

Mesoscale eddies can have a strong impact on regional biogeochemistry and primary productivity. To investigate the effect of the upwelling of seawater by western Pacific eddies on the composition of the active planktonic marine archaeal community composition of the deep chlorophyll maximum (DCM) layer, mesoscale cold-core eddies were simulated in situ by mixing western Pacific DCM layer water with mesopelagic layer (400 m) water. Illumina sequencing of the 16S rRNA gene and 16S rRNA transcripts indicated that the specific heterotrophic Marine Group IIb (MGIIb) taxonomic group of the DCM layer was rapidly stimulated after receiving fresh substrate from 400 m water, which was dominated by uncultured autotrophic Marine Group I (MGI) archaea. Furthermore, niche differentiation of autotrophic ammonia-oxidizing archaea (MGI) was demonstrated by deep sequencing of 16S rRNA, amoA, and accA genes, respectively. Similar distribution patterns of active Marine Group III (MGIII) were observed in the DCM layer with or without vertical mixing, indicating that they are inclined to utilize the substrates already present in the DCM layer. These findings underscore the importance of mesoscale cyclonic eddies in stimulating microbial processes involved in the regional carbon cycle.

Community, Distribution, and Ecological Roles of Estuarine Archaea

Archaea are diverse and ubiquitous prokaryotes present in both extreme and moderate environments. Estuaries, serving as links between the land and ocean, harbor numerous microbes that are relatively highly active because of massive terrigenous input of nutrients. Archaea account for a considerable portion of the estuarine microbial community. They are diverse and play key roles in the estuarine biogeochemical cycles. Ammonia-oxidizing archaea (AOA) are an abundant aquatic archaeal group in estuaries, greatly contributing estuarine ammonia oxidation. Bathyarchaeota are abundant in sediments, and they may involve in sedimentary organic matter degradation, acetogenesis, and, potentially, methane metabolism, based on genomics. Other archaeal groups are also commonly detected in estuaries worldwide. They include Euryarchaeota, and members of the DPANN and Asgard archaea. Based on biodiversity surveys of the 16S rRNA gene and some functional genes, the distribution and abundance of estuarine archaea are driven by physicochemical factors, such as salinity and oxygen concentration. Currently, increasing amount of genomic information for estuarine archaea is becoming available because of the advances in sequencing technologies, especially for AOA and Bathyarchaeota, leading to a better understanding of their functions and environmental adaptations. Here, we summarized the current knowledge on the community composition and major archaeal groups in estuaries, focusing on AOA and Bathyarchaeota. We also highlighted the unique genomic features and potential adaptation strategies of estuarine archaea, pointing out major unknowns in the field and scope for future research.

Marine Group II Euryarchaeota Contribute to the Archaeal Lipid Pool in Northwestern Pacific Ocean Surface Waters

Planktonic archaea include predominantly Marine Group I Thaumarchaeota (MG I) and Marine Group II Euryarchaeota (MG II), which play important roles in the oceanic carbon cycle. MG I produce specific lipids called isoprenoid glycerol dibiphytanyl glycerol tetraethers (GDGTs), which are being used in the sea surface temperature proxy named TEX₈₆. Although MG II may be the most abundant planktonic archaeal group in surface water, their lipid composition remains poorly characterized because of the lack of cultured representatives. Circumstantial evidence from previous studies of marine suspended particulate matter suggests that MG II may produce both GDGTs and archaeol-based lipids. In this study, integration of the 16S rRNA gene quantification and sequencing and lipid analysis demonstrated that MG II contributed significantly to the pool of archaeal tetraether lipids in samples collected from MG II-dominated surface waters of the Northwestern Pacific Ocean (NWPO). The archaeal lipid composition in MG II-dominated NWPO waters differed significantly from that of known MG I cultures, containing relatively more 2G-OH-, 2G- and 1G- GDGTs, especially in their acyclic form. Lipid composition in NWPO waters was also markedly different from MG I-dominated surface water samples collected in the East China Sea. GDGTs from MG II-dominated samples seemed to respond to temperature similarly to GDGTs from the MG I-dominated samples, which calls for further study using pure cultures to determine the exact impact of MG II on GDGT-based proxies.

Metabolic response of prokaryotic microbes to sporadic hypoxia in a eutrophic subtropical estuary

Coastal eutrophication and consequent oxygen depletion (hypoxia) occurs worldwide due to increased human activity. The paucity of genomic information of microbes in hypoxia prone coastal waters have hindered our understanding of microorganism related causation and adaption to the environment. Here, using metagenomic approach, we investigated microbial metabolic capability in heavily polluted Pearl River estuary. Our results highlighted the possible roles of microbial metabolic activity in the formation of bottom water hypoxia by revealing enriched organic degradation related microbial genes in the bottom layer beneath surface phytoplankton bloom. Microbial nitrate reduction in hypoxia layer was low, possibly due to the low pH and fluctuating oxygen level. On contrary, high abundance of sulfate-reducing, and antibiotic and metal resistance related genes were detected in bottom and surface layers, respectively, indicating microbial adaptation to oxygen depletion and pollution. Our study provides gene level information on the interactive relations between microbial functions and environmental stress.

Distinct temporal dynamics of planktonic archaeal and bacterial assemblages in the bays of the Yellow Sea

The Yellow Sea features unique characteristics due to strong tides and nutrient-enriched freshwater outflows from China and Korea. The coupling of archaeal and bacterial assemblages associated with environmental factors at two bay areas in the Yellow Sea was investigated. Temporal variations of the archaeal and bacterial assemblages were shown to be greater than the spatial variations based on an analysis of the 16S rRNA gene sequences. Distinct temporal dynamics of both planktonic archaeal and bacterial assemblages was associated with temperature, NO₂^-, and chlorophyll a ([chl-a]) concentrations in the bays of the Yellow Sea. The [chl-a] was the prime predictor of bacterial abundance, and some taxa were clearly correlated with [chl-a]. Bacteroidetes and Alpha-proteobacteria dominated at high [chl-a] stations while Gamma-proteobacteria (esp. SAR86 clade) and Actinobacteria (Candidatus Actinomarina clade) were abundant at low [chl-a] stations. The archaeal abundance was comparable with the bacterial abundance in most of the October samples. Co-dominance of Marine Group II (MGII) and Candidatus Nitrosopumilus suggests that the assimilation of organic nitrogen by MGII could be coupled with nitrification by ammonia-oxidizing archaea. The distinct temporal dynamics of the archaeal and bacterial assemblages might be attributable to the strong tides and the inflow of nutrient-rich freshwater.

Genomic ecology of Marine Group II, the most common marine planktonic Archaea across the surface ocean

Planktonic Archaea have been detected in all the world's oceans and are found from surface waters to the deep sea. The two most common Archaea phyla are Thaumarchaeota and Euryarchaeota. Euryarchaeota are generally more common in surface waters, but very little is known about their ecology and their potential metabolisms. In this study, we explore the genomic ecology of the Marine Group II (MGII), the main marine planktonic Euryarchaeota, and test if it is composed of different ecologically relevant units. We re‐analyzed Tara Oceans metagenomes from the photic layer and the deep ocean by annotating sequences against a custom MGII database and by mapping gene co‐occurrences. Our data provide a global view of the distribution of Euryarchaeota, and more specifically of MGII subgroups, and reveal their association to a number of gene‐coding sequences. In particular, we show that MGII proteorhodopsins were detected in both the surface and the deep chlorophyll maximum layer and that different clusters of these light harvesting proteins were present. Our approach helped describing the set of genes found together with specific MGII subgroups. We could thus define genomic environments that could theoretically describe ecologically meaningful units and the ecological niche that they occupy.

Archaeal community variation in the Qinhuangdao coastal aquaculture zone revealed by high-throughput sequencing

The differences in archaeal diversity and community composition in the sediments and waters of the Qinhuangdao coastal aquaculture zone were investigated. Furthermore, the associations between dominant archaeal taxa with geographic and environmental variables were evaluated. High-throughput sequencing of archaeal 16S rRNA genes yielded a total of 176,211 quality-filtered reads and 1,178 operational taxonomic units (OTUs) overall. The most abundant phylum and class among all communities were Thaumarchaeota and Nitrososphaeria, respectively. Beta diversity analysis indicated that community composition was divided into two groups according to the habitat type (i.e., sediments or waters). Only 9.8% OTUs were shared by communities from the two habitats, while 73.9% and 16.3% of the OTUs were unique to sediment or water communities, respectively. Furthermore, the relative abundances of the dominant OTUs differed with habitat type. Investigations of relationships between dominant OTUs and environmental variables indicated that some dominant OTUs were more sensitive to variation in environmental factors, which could be due to individual taxonomic differences in lifestyles and biological processes. Overall, the investigation of archaeal community variation within the Qinhuangdao coastal aquaculture zone provides an important baseline understanding of the microbial ecology in this important ecosystem.

Metabolic diversity within the globally abundant Marine Group II Euryarchaea offers insight into ecological patterns

Despite their discovery over 25 years ago, the Marine Group II Euryarchaea (MGII) remain a difficult group of organisms to study, lacking cultured isolates and genome references. The MGII have been identified in marine samples from around the world, and evidence supports a photoheterotrophic lifestyle combining phototrophy via proteorhodopsins with the remineralization of high molecular weight organic matter. Divided between two clades, the MGII have distinct ecological patterns that are not understood based on the limited number of available genomes. Here, I present a comparative genomic analysis of 250 MGII genomes, providing a comprehensive investigation of these mesophilic archaea. This analysis identifies 17 distinct subclades including nine subclades that previously lacked reference genomes. The metabolic potential and distribution of the MGII genera reveals distinct roles in the environment, identifying algal-saccharide-degrading coastal subclades, protein-degrading oligotrophic surface ocean subclades, and mesopelagic subclades lacking proteorhodopsins, common in all other subclades.

Evolving paradigms in biological carbon cycling in the ocean

Carbon is a keystone element in global biogeochemical cycles. It plays a fundamental role in biotic and abiotic processes in the ocean, which intertwine to mediate the chemistry and redox status of carbon in the ocean and the atmosphere. The interactions between abiotic and biogenic carbon (e.g. CO2, CaCO3, organic matter) in the ocean are complex, and there is a half-century-old enigma about the existence of a huge reservoir of recalcitrant dissolved organic carbon (RDOC) that equates to the magnitude of the pool of atmospheric CO2. The concepts of the biological carbon pump (BCP) and the microbial loop (ML) shaped our understanding of the marine carbon cycle. The more recent concept of the microbial carbon pump (MCP), which is closely connected to those of the BCP and the ML, explicitly considers the significance of the ocean's RDOC reservoir and provides a mechanistic framework for the exploration of its formation and persistence. Understanding of the MCP has benefited from advanced ‘omics’ and novel research in biological oceanography and microbial biogeochemistry. The need to predict the ocean's response to climate change makes an integrative understanding of the BCP, ML and MCP a high priority. In this review, we summarize and discuss progress since the proposal of the MCP in 2010 and formulate research questions for the future.