Metabolic activities of marine ammonia-oxidizing archaea orchestrated by quorum sensing.

Ammonia-oxidizing archaea (AOA) play crucial roles in marine carbon and nitrogen cycles by fixing inorganic carbon and performing the initial step of nitrification. Evaluation of carbon and nitrogen metabolism popularly relies on functional genes such as amoA and accA. Increasing studies suggest that quorum sensing (QS) mainly studied in biofilms for bacteria may serve as a universal communication and regulatory mechanism among prokaryotes; however, this has yet to be demonstrated in marine planktonic archaea. To bridge this knowledge gap, we employed a combination of metabolic activity markers (amoA, accA, and grs) to elucidate the regulation of AOA-mediated nitrogen, carbon processes, and their interactions with the surrounding heterotrophic population. Through co-transcription investigations linking metabolic markers to potential key QS genes, we discovered that QS molecules could regulate AOA's carbon, nitrogen, and lipid metabolisms under different conditions. Interestingly, specific AOA ecotypes showed a preference for employing distinct QS systems and a distinct QS circuit involving a typical population. Overall, our data demonstrate that QS orchestrates nitrogen and carbon metabolism, including the exchange of organic metabolites between AOA and surrounding heterotrophic bacteria, which has been previously overlooked in marine AOA research.

Isolation of a methyl-reducing methanogen outside the Euryarchaeota.

Methanogenic archaea are main contributors to methane emissions, and have a crucial role in carbon cycling and global warming. Until recently, methanogens were confined to Euryarchaeota, but metagenomic studies revealed the presence of genes encoding the methyl coenzyme M reductase complex in other archaeal clades1-4, thereby opening up the premise that methanogenesis is taxonomically more widespread. Nevertheless, laboratory cultivation of these non-euryarchaeal methanogens was lacking to corroborate their potential methanogenic ability and physiology. Here we report the isolation of a thermophilic archaeon LWZ-6 from an oil field. This archaeon belongs to the class Methanosuratincolia (originally affiliated with 'Candidatus Verstraetearchaeota') in the phylum Thermoproteota. Methanosuratincola petrocarbonis LWZ-6 is a strict hydrogen-dependent methylotrophic methanogen. Although previous metagenomic studies speculated on the fermentative potential of Methanosuratincolia members, strain LWZ-6 does not ferment sugars, peptides or amino acids. Its energy metabolism is linked only to methanogenesis, with methanol and monomethylamine as electron acceptors and hydrogen as an electron donor. Comparative (meta)genome analysis confirmed that hydrogen-dependent methylotrophic methanogenesis is a widespread trait among Methanosuratincolia. Our findings confirm that the diversity of methanogens expands beyond the classical Euryarchaeota and imply the importance of hydrogen-dependent methylotrophic methanogenesis in global methane emissions and carbon cycle.

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.

Application of single-cell Raman-deuterium isotope probing to reveal the resistance of marine ammonia-oxidizing archaea SCM1 against common antibiotics.

Antimicrobial resistance (AMR) in oceans poses a significant threat to human health through the seafood supply chain. Ammonia-oxidizing archaea (AOA) are important marine microorganisms and play a key role in the biogeochemical nitrogen cycle around the world. However, the AMR of marine AOA to aquicultural antibiotics is poorly explored. Here, Raman-deuterium isotope probing (Raman-DIP), a single-cell tool, was developed to reveal the AMR of a typical marine species of AOA, Nitrosopumilus maritimus (designated SCM1), against six antibiotics, including erythromycin, tetracycline, novobiocin, neomycin, bacitracin, and vancomycin. The D2O concentration (30% v/v) and culture period (9 days) were optimized for the precise detection of metabolic activity in SCM1 cells through Raman-DIP. The relative metabolic activity of SCM1 upon exposure to antibiotics was semi-quantitatively calculated based on single-cell Raman spectra. SCM1 exhibited high resistance to erythromycin, tetracycline, novobiocin, neomycin, and vancomycin, with minimum inhibitory concentration (MIC) values between 100 and 400 mg/L, while SCM1 is very sensitive to bacitracin (MIC: 0.8 mg/L). Notably, SCM1 cells were completely inactive under the metabolic activity minimum inhibitory concentration conditions (MA-MIC: 1.6-800 mg/L) for the six antibiotics. Further genomic analysis revealed the antibiotic resistance genes (ARGs) of SCM1, including 14 types categorized into 33 subtypes. This work increases our knowledge of the AMR of marine AOA by linking the resistant phenome to the genome, contributing to the risk assessment of AMR in the underexplored ocean environment. As antibiotic resistance in marine microorganisms is significantly affected by the concentration of antibiotics in coastal environments, we encourage more studies concentrating on both the phenotypic and genotypic antibiotic resistance of marine archaea. This may facilitate a comprehensive evaluation of the capacity of marine microorganisms to spread AMR and the implementation of suitable control measures to protect environmental safety and human health.

Ultrastructural insights into cellular organization, energy storage and ribosomal dynamics of an ammonia-oxidizing archaeon from oligotrophic oceans.

Introduction: Nitrososphaeria, formerly known as Thaumarchaeota, constitute a diverse and widespread group of ammonia-oxidizing archaea (AOA) inhabiting ubiquitously in marine and terrestrial environments, playing a pivotal role in global nitrogen cycling. Despite their importance in Earth's ecosystems, the cellular organization of AOA remains largely unexplored, leading to a significant unanswered question of how the machinery of these organisms underpins metabolic functions. Methods: In this study, we combined spherical-chromatic-aberration-corrected cryo-electron tomography (cryo-ET), scanning transmission electron microscopy (STEM), and energy dispersive X-ray spectroscopy (EDS) to unveil the cellular organization and elemental composition of Nitrosopumilus maritimus SCM1, a representative member of marine Nitrososphaeria. Results and Discussion: Our tomograms show the native ultrastructural morphology of SCM1 and one to several dense storage granules in the cytoplasm. STEM-EDS analysis identifies two types of storage granules: one type is possibly composed of polyphosphate and the other polyhydroxyalkanoate. With precise measurements using cryo-ET, we observed low quantity and density of ribosomes in SCM1 cells, which are in alignment with the documented slow growth of AOA in laboratory cultures. Collectively, these findings provide visual evidence supporting the resilience of AOA in the vast oligotrophic marine environment.

Lowering the Schottky Barrier Height by Quasi-van der Waals Contacts for High-Performance p-Type MoTe2 Field-Effect Transistors.

Two-dimensional (2D) transition-metal dichalcogenides (TMDs) offer advantages over traditional silicon in future electronics but are hampered by the prominent high contact resistance of metal-TMD interfaces, especially for p-type TMDs. Here, we present high-performance p-type MoTe2 field-effect transistors via a nondestructive van der Waals (vdW) transfer process, establishing low contact resistance between the 2D MoTe2 semiconductor and the PtTe2 semimetal. The integration of PtTe2 as contacts in MoTe2 field-effect transistors leads to significantly improved electrical characteristics compared to conventional metal contacts, evidenced by a mobility increase to 80 cm2 V-1 s-1, an on-state current rise to 5.0 µA/µm, and a reduction in Schottky barrier height (SBH) to 48 meV. Such a low SBH in quasi-van der Waals contacts can be assigned to the low electrical resistivity of PtTe2 and the high efficiency of carrier injection at the 2D semimetal/2D semiconductor interfaces. Imaging via transmission electron microscopy reveals that the 2D semimetal/two-dimensional semiconductor interfaces are atomically flat and exceptionally clean. This interface engineering strategy could enable low-resistance contacts based on vdW architectures in a facile manner, providing opportunities for 2D materials for next-generation optoelectronics and electronics.

Strong linkage between benthic oxygen uptake and bacterial tetraether lipids in deep-sea trench regions.

Oxygen in marine sediments regulates many key biogeochemical processes, playing a crucial role in shaping Earth's climate and benthic ecosystems. In this context, branched glycerol dialkyl glycerol tetraethers (brGDGTs), essential biomarkers in paleoenvironmental research, exhibit an as-yet-unresolved association with sediment oxygen conditions. Here, we investigated brGDGTs in sediments from three deep-sea regions (4045 to 10,100 m water depth) dominated by three respective trench systems and integrated the results with in situ oxygen microprofile data. Our results demonstrate robust correlations between diffusive oxygen uptake (DOU) obtained from microprofiles and brGDGT methylation and isomerization degrees, indicating their primary production within sediments and their strong linkage with microbial diagenetic activity. We establish a quantitative relationship between the Isomerization and Methylation index of Branched Tetraethers (IMBT) and DOU, suggesting its potential validity across deep-sea environments. Increased brGDGT methylation and isomerization likely enhance the fitness of source organisms in deep-sea habitats. Our study positions brGDGTs as a promising tool for quantifying benthic DOU in deep-sea settings, where DOU is a key metric for assessing sedimentary organic carbon degradation and microbial activity.

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.

RSPSSL: A novel high-fidelity Raman spectral preprocessing scheme to enhance biomedical applications and chemical resolution visualization.

Raman spectroscopy has tremendous potential for material analysis with its molecular fingerprinting capability in many branches of science and technology. It is also an emerging omics technique for metabolic profiling to shape precision medicine. However, precisely attributing vibration peaks coupled with specific environmental, instrumental, and specimen noise is problematic. Intelligent Raman spectral preprocessing to remove statistical bias noise and sample-related errors should provide a powerful tool for valuable information extraction. Here, we propose a novel Raman spectral preprocessing scheme based on self-supervised learning (RSPSSL) with high capacity and spectral fidelity. It can preprocess arbitrary Raman spectra without further training at a speed of ~1 900 spectra per second without human interference. The experimental data preprocessing trial demonstrated its excellent capacity and signal fidelity with an 88% reduction in root mean square error and a 60% reduction in infinite norm ([Formula: see text]) compared to established techniques. With this advantage, it remarkably enhanced various biomedical applications with a 400% accuracy elevation (ΔAUC) in cancer diagnosis, an average 38% (few-shot) and 242% accuracy improvement in paraquat concentration prediction, and unsealed the chemical resolution of biomedical hyperspectral images, especially in the spectral fingerprint region. It precisely preprocessed various Raman spectra from different spectroscopy devices, laboratories, and diverse applications. This scheme will enable biomedical mechanism screening with the label-free volumetric molecular imaging tool on organism and disease metabolomics profiling with a scenario of high throughput, cross-device, various analyte complexity, and diverse applications.

Novel order-level lineage of ammonia-oxidizing archaea widespread in marine and terrestrial environments.

Ammonia-oxidizing archaea (AOA) are among the most ubiquitous and abundant archaea on Earth, widely distributed in marine, terrestrial, and geothermal ecosystems. However, the genomic diversity, biogeography, and evolutionary process of AOA populations in subsurface environments are vastly understudied compared to those in marine and soil systems. Here, we report a novel AOA order Candidatus (Ca.) Nitrosomirales which forms a sister lineage to the thermophilic Ca. Nitrosocaldales. Metagenomic and 16S rRNA gene-read mapping demonstrates the abundant presence of Nitrosomirales AOA in various groundwater environments and their widespread distribution across a range of geothermal, terrestrial, and marine habitats. Terrestrial Nitrosomirales AOA show the genetic capacity of using formate as a source of reductant and using nitrate as an alternative electron acceptor. Nitrosomirales AOA appear to have acquired key metabolic genes and operons from other mesophilic populations via horizontal gene transfer, including genes encoding urease, nitrite reductase, and V-type ATPase. The additional metabolic versatility conferred by acquired functions may have facilitated their radiation into a variety of subsurface, marine, and soil environments. We also provide evidence that each of the four AOA orders spans both marine and terrestrial habitats, which suggests a more complex evolutionary history for major AOA lineages than previously proposed. Together, these findings establish a robust phylogenomic framework of AOA and provide new insights into the ecology and adaptation of this globally abundant functional guild.

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.

Lipidomic chemotaxonomy aligned with phylogeny of Halobacteria.

Archaea play an important role in global biogeochemical cycles and are considered ancestral to eukaryotes. The unique lipid composition of archaea, characterized by isoprenoid alkyl chains and ether linkage to glycerol-1-phosphate, offers valuable insights into archaeal phylogeny and evolution. However, comprehensive studies focusing on archaeal lipidomes, especially at the intact polar lipid level, are currently limited. Here, we built an in-house library of archaeal lipids by using high-performance liquid chromatography coupled with mass-spectrometry, which was integrated with bioinformatics and molecular network analyses. Seven halobacterial strains, representing three distinct orders, were cultured under identical conditions to investigate their lipidomes. A total of 162 features were identified, corresponding to 107 lipids that could be assigned to different strains. Clustering analyses of both core lipids and total lipids matched the phylogeny of Halobacteria at the order level. Notably, lipids such as triglycosyl diether-phosphatidyl acid and bis-sulfate glycosyl lipids were specific to particular groups and could serve as diagnostic intact lipid biomarkers for Halobacteria. Furthermore, the analysis of network-coordinated features facilitated the linkage of unknown lipid compounds to phylogeny, which promotes a lipidome to phylogeny matchup among three Haloferax strains, thereby expanding the knowledge of the halobacterial lipidome. Our study provides a comprehensive view of the lipidomes of the seven strains of Halobacteria and highlights the potential of lipidomics for studying archaeal phylogeny.

Lipidomic diversity and proxy implications of archaea from cold seep sediments of the South China Sea.

Cold seeps on the continental margins are characterized by intense microbial activities that consume a large portion of methane by anaerobic methanotrophic archaea (ANME) through anaerobic oxidation of methane (AOM). Although ANMEs are known to contain unique ether lipids that may have an important function in marine carbon cycling, their full lipidomic profiles and functional distribution in particular cold-seep settings are still poorly characterized. Here, we combined the 16S rRNA gene sequencing and lipidomic approaches to analyze archaeal communities and their lipids in cold seep sediments with distinct methane supplies from the South China Sea. The archaeal community was dominated by ANME-1 in the moderate seepage area with strong methane emission. Low seepage area presented higher archaeal diversity covering Lokiarchaeia, Bathyarchaeia, and Thermoplasmata. A total of 55 core lipids (CLs) and intact polar lipids (IPLs) of archaea were identified, which included glycerol dialkyl glycerol tetraethers (GDGTs), hydroxy-GDGTs (OH-GDGTs), archaeol (AR), hydroxyarchaeol (OH-AR), and dihydroxyarchaeol (2OH-AR). Diverse polar headgroups constituted the archaeal IPLs. High concentrations of dissolved inorganic carbon (DIC) with depleted δ13CDIC and high methane index (MI) values based on both CLs (MICL) and IPLs (MIIPL) indicate that ANMEs were active in the moderate seepage area. The ANME-2 and ANME-3 clades were characterized by enhanced glycosidic and phosphoric diether lipids production, indicating their potential role in coupling carbon and phosphurus cycling in cold seep ecosystems. ANME-1, though representing a smaller proportion of total archaea than ANME-2 and ANME-3 in the low seepage area, showed a positive correlation with MIIPL, indicating a different mechanism contributing to the IPL-GDGT pool. This also suggests that MIIPL could be a sensitive index to trace AOM activities performed by ANME-1. Overall, our study expands the understanding of the archaeal lipid composition in the cold seep and improves the application of MI using intact polar lipids that potentially link to extent ANME activities.

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.

Changing microbiome community structure and functional potential during permafrost thawing on the Tibetan Plateau.

Large amounts of carbon sequestered in permafrost on the Tibetan Plateau (TP) are becoming vulnerable to microbial decomposition in a warming world. However, knowledge about how the responsible microbial community responds to warming-induced permafrost thaw on the TP is still limited. This study aimed to conduct a comprehensive comparison of the microbial communities and their functional potential in the active layer of thawing permafrost on the TP. We found that the microbial communities were diverse and varied across soil profiles. The microbial diversity declined and the relative abundance of Chloroflexi, Bacteroidetes, Euryarchaeota, and Bathyarchaeota significantly increased with permafrost thawing. Moreover, warming reduced the similarity and stability of active layer microbial communities. The high-throughput qPCR results showed that the abundance of functional genes involved in liable carbon degradation and methanogenesis increased with permafrost thawing. Notably, the significantly increased mcrA gene abundance and the higher methanogens to methanotrophs ratio implied enhanced methanogenic activities during permafrost thawing. Overall, the composition and functional potentials of the active layer microbial community in the Tibetan permafrost region are susceptible to warming. These changes in the responsible microbial community may accelerate carbon degradation, particularly in the methane releases from alpine permafrost ecosystems on the TP.

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.

Virus impacted community adaptation in oligotrophic groundwater environment revealed by Hi-C coupled metagenomic and viromic study.

Viruses play a crucial role in microbial mortality, diversity and biogeochemical cycles. Groundwater is the largest global freshwater and one of the most oligotrophic aquatic systems on Earth, but how microbial and viral communities are shaped in this special habitat is largely unexplored. In this study, we collected groundwater samples from 23 to 60 m aquifers at Yinchuan Plain, China. In total, 1920 non-reductant viral contigs were retrieved from metagenomes and viromes constructed by Illumina and Nanopore hybrid sequencing. Only 3% of them could be clustered with known viruses, most of which were Caudoviricetes. Coupling 1.2 Tb Hi-C sequencing with CRISPR matching and homology search, we connected 469 viruses with their hosts while some viral clusters presented a broad-host-range trait. Meanwhile, a large proportion of biosynthesis related auxiliary metabolism genes were identified. Those characteristics might benefit viruses for a better survival in this special oligotrophic environment. Additionally, the groundwater virome showed genomic features distinct from those of the open ocean and wastewater treatment facilities in GC distribution and unannotated gene compositions. This paper expands the current knowledge of the global viromic records and serves as a foundation for a more thorough understanding of viruses in groundwater.

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.