Hindsight in the relative abundance, metabolic potential and genome dynamics of uncultivated marine archaea from comparative metagenomic analyses of bathypelagic plankton of different oceanic regions

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.

Influence of host genotype in establishing root associated microbiome of indica rice cultivars for plant growth promotion

Rice plants display a unique root ecosystem comprising oxic-anoxic zones, harboring a plethora of metabolic interactions mediated by its root microbiome. Since agricultural land is limited, an increase in rice production will rely on novel methods of yield enhancement. The nascent concept of tailoring plant phenotype through the intervention of synthetic microbial communities (SynComs) is inspired by the genetics and ecology of core rhizobiome. In this direction, we have studied structural and functional variations in the root microbiome of 10 indica rice varieties. The studies on α and β-diversity indices of rhizospheric root microbiome with the host genotypes revealed variations in the structuring of root microbiome as well as a strong association with the host genotypes. Biomarker discovery, using machine learning, highlighted members of class Anaerolineae, α-Proteobacteria, and bacterial genera like Desulfobacteria, Ca. Entotheonella, Algoriphagus, etc. as the most important features of indica rice microbiota having a role in improving the plant's fitness. Metabolically, rice rhizobiomes showed an abundance of genes related to sulfur oxidation and reduction, biofilm production, nitrogen fixation, denitrification, and phosphorus metabolism. This comparative study of rhizobiomes has outlined the taxonomic composition and functional diversification of rice rhizobiome, laying the foundation for the development of next-generation microbiome-based technologies for yield enhancement in rice and other crops.

Planktonic Archaeal Ether Lipid Origins in Surface Waters of the North Pacific Subtropical Gyre

Thaumarchaeota and Thermoplasmatota are the most abundant planktonic archaea in the sea. Thaumarchaeota contain tetraether lipids as their major membrane lipids, but the lipid composition of uncultured planktonic Thermoplasmatota representatives remains unknown. To address this knowledge gap, we quantified archaeal cells and ether lipids in open ocean depth profiles (0–200 m) of the North Pacific Subtropical Gyre. Planktonic archaeal community structure and ether lipid composition in the water column partitioned into two separate clusters: one above the deep chlorophyll maximum, the other within and below it. In surface waters, Thermoplasmatota densities ranged from 2.11 × 10⁶ to 6.02 × 10⁶ cells/L, while Thaumarchaeota were undetectable. As previously reported for Thaumarchaeota, potential homologs of archaeal tetraether ring synthases were present in planktonic Thermoplasmatota metagenomes. Despite the absence of Thaumarchaeota in surface waters, measurable amounts of intact polar ether lipids were found there. Based on cell abundance estimates, these surface water archaeal ether lipids contributed only 1.21 × 10^–9 ng lipid/Thermoplasmatota cell, about three orders of magnitude less than that reported for Thaumarchaeota cells. While these data indicate that even if some tetraether and diether lipids may be derived from Thermoplasmatota, they would only comprise a small fraction of Thermoplasmatota total biomass. Therefore, while both MGI Thaumarchaeota and MGII/III Thermoplasmatota are potential biological sources of archaeal GDGTs, the Thaumarchaeota appear to be the major contributors of archaeal tetraether lipids in planktonic marine habitats. These results extend and confirm previous reports of planktonic archaeal lipid sources, and further emphasize the need for Thermoplasmatota cultivation, to better characterize the membrane lipid constituents of marine planktonic Thermoplasmatota, and more precisely define the sources and patterns of archaeal tetraether lipid distributions in marine plankton.

Gulf of Mexico blue hole harbors high levels of novel microbial lineages

Exploration of oxygen-depleted marine environments has consistently revealed novel microbial taxa and metabolic capabilities that expand our understanding of microbial evolution and ecology. Marine blue holes are shallow karst formations characterized by low oxygen and high organic matter content. They are logistically challenging to sample, and thus our understanding of their biogeochemistry and microbial ecology is limited. We present a metagenomic and geochemical characterization of Amberjack Hole on the Florida continental shelf (Gulf of Mexico). Dissolved oxygen became depleted at the hole's rim (32 m water depth), remained low but detectable in an intermediate hypoxic zone (40-75 m), and then increased to a secondary peak before falling below detection in the bottom layer (80-110 m), concomitant with increases in nutrients, dissolved iron, and a series of sequentially more reduced sulfur species. Microbial communities in the bottom layer contained heretofore undocumented levels of the recently discovered phylum Woesearchaeota (up to 58% of the community), along with lineages in the bacterial Candidate Phyla Radiation (CPR). Thirty-one high-quality metagenome-assembled genomes (MAGs) showed extensive biochemical capabilities for sulfur and nitrogen cycling, as well as for resisting and respiring arsenic. One uncharacterized gene associated with a CPR lineage differentiated hypoxic from anoxic zone communities. Overall, microbial communities and geochemical profiles were stable across two sampling dates in the spring and fall of 2019. The blue hole habitat is a natural marine laboratory that provides opportunities for sampling taxa with under-characterized but potentially important roles in redox-stratified microbial processes.

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.

Ecogenomics and metabolic potential of the South Atlantic Ocean microbiome

The unique combination of depth, salinity, and water masses make the South Atlantic Ocean an ecosystem of special relevance within the global ocean. Yet, the microbiome of this ecosystem has received less attention than other regions of the global Ocean. This has hampered our understanding of the diversity and metabolic potential of the microorganisms that dwell in this habitat. To fill this knowledge gap, we analyzed a collection of 31 metagenomes from the Atlantic Ocean that spanned the epipelagic, mesopelagic and bathypelagic zones (surface to 4000 m). Read-centric and gene-centric analysis revealed the unique taxonomic and functional composition of metagenomes from each depth zone, which was driven by differences in physical and chemical parameters. In parallel, a total of 40 metagenome-assembled genomes were obtained, which recovered one third of the total community. Phylogenomic reconstruction revealed that many of these genomes are derived from poorly characterized taxa of Bacteria and Archaea. Genomes derived from heterotrophic bacteria of the aphotic zone displayed a large apparatus of genes suited for the utilization of recalcitrant organic compounds such as cellulose, chitin and alkanes. In addition, we found genomic evidence suggesting that mixotrophic bacteria from the bathypelagic zone could perform carbon fixation through the Calvin-Benson-Bassham cycle, fueled by sulfur oxidation. Finally, we found that the viral communities shifted throughout the water column regarding their targeted hosts and virus-to-microbe ratio, in response to shifts in the composition and functioning their microbial counterparts. Our findings shed light on the microbial and viral drivers of important biogeochemical processes that take place in the South Atlantic Ocean.

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.

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.

Heterotrophic Thaumarchaea with Small Genomes Are Widespread in the Dark Ocean

The Thaumarchaeota is a diverse archaeal phylum comprising numerous lineages that play key roles in global biogeochemical cycling, particularly in the ocean. To date, all genomically characterized marine thaumarchaea are reported to be chemolithoautotrophic ammonia oxidizers. In this study, we report a group of putatively heterotrophic marine thaumarchaea (HMT) with small genome sizes that is globally abundant in the mesopelagic, apparently lacking the ability to oxidize ammonia. We assembled five HMT genomes from metagenomic data and show that they form a deeply branching sister lineage to the ammonia-oxidizing archaea (AOA). We identify this group in metagenomes from mesopelagic waters in all major ocean basins, with abundances reaching up to 6% of that of AOA. Surprisingly, we predict the HMT have small genomes of ∼1 Mbp, and our ancestral state reconstruction indicates this lineage has undergone substantial genome reduction compared to other related archaea. The genomic repertoire of HMT indicates a versatile metabolism for aerobic chemoorganoheterotrophy that includes a divergent form III-a RuBisCO, a 2M respiratory complex I that has been hypothesized to increase energetic efficiency, and a three-subunit heme-copper oxidase complex IV that is absent from AOA. We also identify 21 pyrroloquinoline quinone (PQQ)-dependent dehydrogenases that are predicted to supply reducing equivalents to the electron transport chain and are among the most highly expressed HMT genes, suggesting these enzymes play an important role in the physiology of this group. Our results suggest that heterotrophic members of the Thaumarchaeota are widespread in the ocean and potentially play key roles in global chemical transformations.IMPORTANCE It has been known for many years that marine Thaumarchaeota are abundant constituents of dark ocean microbial communities, where their ability to couple ammonia oxidation and carbon fixation plays a critical role in nutrient dynamics. In this study, we describe an abundant group of putatively heterotrophic marine Thaumarchaeota (HMT) in the ocean with physiology distinct from those of their ammonia-oxidizing relatives. HMT lack the ability to oxidize ammonia and fix carbon via the 3-hydroxypropionate/4-hydroxybutyrate pathway but instead encode a form III-a RuBisCO and diverse PQQ-dependent dehydrogenases that are likely used to conserve energy in the dark ocean. Our work expands the scope of known diversity of Thaumarchaeota in the ocean and provides important insight into a widespread marine lineage.

Metagenome-assembled genomes reveal unique metabolic adaptations of a basal marine Thaumarchaeota lineage

Thaumarchaeota constitute an abundant and ubiquitous phylum of Archaea that play critical roles in the global nitrogen and carbon cycles. Most well-characterized members of the phylum are chemolithoautotrophic ammonia-oxidizing archaea (AOA), which comprise up to 5 and 20% of the total single-celled life in soil and marine systems, respectively. Using two high-quality metagenome-assembled genomes (MAGs), here we describe a divergent marine thaumarchaeal clade that is devoid of the ammonia-oxidation machinery and the AOA-specific carbon-fixation pathway. Phylogenomic analyses placed these genomes within the uncultivated and largely understudied marine pSL12-like thaumarchaeal clade. The predominant mode of nutrient acquisition appears to be aerobic heterotrophy, evidenced by the presence of respiratory complexes and various organic carbon degradation pathways. Both genomes encoded several pyrroloquinoline quinone (PQQ)-dependent alcohol dehydrogenases, as well as a form III RuBisCO. Metabolic reconstructions suggest anaplerotic CO₂ assimilation mediated by RuBisCO, which may be linked to the central carbon metabolism. We conclude that these genomes represent a hitherto unrecognized evolutionary link between predominantly anaerobic basal thaumarchaeal lineages and mesophilic marine AOA, with important implications for diversification within the phylum Thaumarchaeota.

Prokaryotic Population Dynamics and Viral Predation in a Marine Succession Experiment Using Metagenomics

We performed an incubation experiment of seawater confined in plastic bottles with samples collected at three depths (15, 60, and 90 m) after retrieval from a single offshore location in the Mediterranean Sea, from a late summer stratified water column. Two samples representative of each depth were collected and stored in opaque bottles after two periods of 7 h. We took advantage of the "bottle effect" to investigate changes in the natural microbial communities (abundant and rare). We recovered 94 metagenome-assembled genomes (MAGs) and 1089 metagenomic viral contigs and examined their abundance using metagenomic recruitment. We detected a significant fast growth of copiotrophic bacteria such as Alteromonas or Erythrobacter throughout the entire water column with different dynamics that we assign to "clonal," "polyclonal," or "multispecies" depending on the recruitment pattern. Results also showed a marked ecotype succession in the phototropic picocyanobacteria that were able to grow at all the depths in the absence of light, highlighting the importance of their mixotrophic potential. In addition, "wall-chain-reaction" hypothesis based on the study of phage-host dynamics showed the higher impact of viral predation on archaea in deeper waters, evidencing their prominent role during incubations. Our results provide a step forward in understanding the mechanisms underlying dynamic patterns and ecology of the marine microbiome and the importance of processing the samples immediately after collection to avoid changes in the community structure.

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 biogeography and interactions with microbial community across complex subtropical coastal waters

Marine Archaea are crucial in biogeochemical cycles, but their horizontal spatial variability, assembly processes, and microbial associations across complex coastal waters still lack characterizations at high coverage. Using a dense sampling strategy, we investigated horizontal variability in total archaeal, Thaumarchaeota Marine Group (MG) I, and Euryarchaeota MGII communities and associations of MGI/MGII with other microbes in surface waters with contrasting environmental characteristics across ~200 km by 16S rRNA gene amplicon sequencing. Total archaeal communities were extremely dominated by MGI and/or MGII (98.9% in average relative abundance). Niche partitioning between MGI and MGII or within each group was found across multiple environmental gradients. "Selection" was more important than "dispersal limitation" in governing biogeographic patterns of total archaeal, MGI, and MGII communities, and basic abiotic parameters (such as salinity) and inorganic/organic resources as a whole could be the main driver of "selection". While "homogenizing dispersal" also considerably governed their biogeography. MGI-Nitrospira assemblages were speculatively responsible for complete nitrification. MGI taxa commonly had negative correlations with members of Synechococcus but positive correlations with members of eukaryotic phytoplankton, suggesting that competition or synergy between MGI and phytoplankton depends on specific MGI-phytoplankton assemblages. MGII taxa showed common associations with presumed (photo)heterotrophs including members of SAR11, SAR86, SAR406, and Candidatus Actinomarina. This study sheds light on ecological processes and drivers shaping archaeal biogeography and many strong MGI/MGII-bacterial associations across complex subtropical coastal waters. Future efforts should be made on seasonality of archaeal biogeography and biological, environmental, or ecological mechanisms underlying these statistical microbial associations.

Extent of the annual Gulf of Mexico hypoxic zone influences microbial community structure

Rich geochemical datasets generated over the past 30 years have provided fine-scale resolution on the northern Gulf of Mexico (nGOM) coastal hypoxic (≤ 2 mg of O₂ L^-1) zone. In contrast, little is known about microbial community structure and activity in the hypoxic zone despite the implication that microbial respiration is responsible for forming low dissolved oxygen (DO) conditions. Here, we hypothesized that the extent of the hypoxic zone is a driver in determining microbial community structure, and in particular, the abundance of ammonia-oxidizing archaea (AOA). Samples collected across the shelf for two consecutive hypoxic seasons in July 2013 and 2014 were analyzed using 16S rRNA gene sequencing, oligotyping, microbial co-occurrence analysis, and quantification of thaumarchaeal 16S rRNA and archaeal ammonia-monooxygenase (amoA) genes. In 2014 Thaumarchaeota were enriched and inversely correlated with DO while Cyanobacteria, Acidimicrobiia, and Proteobacteria where more abundant in oxic samples compared to hypoxic. Oligotyping analysis of Nitrosopumilus 16S rRNA gene sequences revealed that one oligotype was significantly inversely correlated with DO in both years. Oligotyping analysis revealed single nucleotide variation among all Nitrosopumilaceae, including Nitrosopumilus 16S rRNA gene sequences, with one oligotype possibly being better adapted to hypoxia. We further provide evidence that in the hypoxic zone of both year 2013 and 2014, low DO concentrations and high Thaumarchaeota abundances influenced microbial co-occurrence patterns. Taken together, the data demonstrated that the extent of hypoxic conditions could potentially drive patterns in microbial community structure, with two years of data revealing the annual nGOM hypoxic zone to be emerging as a low DO adapted AOA hotspot.

Archaeal Sources of Intact Membrane Lipid Biomarkers in the Oxygen Deficient Zone of the Eastern Tropical South Pacific

Archaea are ubiquitous in the modern ocean where they are involved in the carbon and nitrogen biogeochemical cycles. However, the majority of Archaea remain uncultured. Archaeal specific membrane intact polar lipids (IPLs) are biomarkers of the presence and abundance of living cells. They comprise archaeol and glycerol dibiphytanyl glycerol tetraethers (GDGTs) attached to various polar headgroups. However, little is known of the IPLs of uncultured marine Archaea, complicating their use as biomarkers. Here, we analyzed suspended particulate matter (SPM) obtained in high depth resolution from a coastal and open ocean site in the eastern tropical South Pacific (ETSP) oxygen deficient zone (ODZ) with the aim of determining possible biological sources of archaeal IPL by comparing their composition by Ultra High Pressure Liquid Chromatography coupled to high resolution mass spectrometry with the archaeal diversity by 16S rRNA gene amplicon sequencing and their abundance by quantitative PCR. Thaumarchaeotal Marine Group I (MGI) closely related to Ca. Nitrosopelagicus and Nitrosopumilus dominated the oxic surface and upper ODZ water together with Marine Euryarchaeota Group II (MGII). High relative abundance of hexose phosphohexose- (HPH) crenarchaeol, the specific biomarker for living Thaumarchaeota, and HPH-GDGT-0, dihexose- (DH) GDGT-3 and -4 were detected in these water masses. Within the ODZ, DPANN (Diapherotrites, Parvarchaeota, Aenigmarchaeota, Nanoarchaeota, and Nanohaloarchaea) of the Woesearchaeota DHVE-6 group and Marine Euryarchaeota Group III (MGIII) were present together with a higher proportion of archaeol-based IPLs, which were likely made by MGIII, since DPANN archaea are supposedly unable to synthesize their own IPLs and possibly have a symbiotic or parasitic partnership with MGIII. Finally, in deep suboxic/oxic waters a different MGI population occurred with HPH-GDGT-1, -2 and DH-GDGT-0 and -crenarchaeol, indicating that here MGI synthesize membranes with IPLs in a different relative abundance which could be attributed to the different detected population or to an environmental adaptation. Our study sheds light on the complex archaeal community of one of the most prominent ODZs and on the IPL biomarkers they potentially synthesize.

Planktonic Marine Archaea

Archaea are ubiquitous and abundant members of the marine plankton. Once thought of as rare organisms found in exotic extremes of temperature, pressure, or salinity, archaea are now known in nearly every marine environment. Though frequently referred to collectively, the planktonic archaea actually comprise four major phylogenetic groups, each with its own distinct physiology and ecology. Only one group-the marine Thaumarchaeota-has cultivated representatives, making marine archaea an attractive focus point for the latest developments in cultivation-independent molecular methods. Here, we review the ecology, physiology, and biogeochemical impact of the four archaeal groups using recent insights from cultures and large-scale environmental sequencing studies. We highlight key gaps in our knowledge about the ecological roles of marine archaea in carbon flow and food web interactions. We emphasize the incredible uncultivated diversity within each of the four groups, suggesting there is much more to be done.

Fine metagenomic profile of the Mediterranean stratified and mixed water columns revealed by assembly and recruitment

BACKGROUND

The photic zone of aquatic habitats is subjected to strong physicochemical gradients. To analyze the fine-scale variations in the marine microbiome, we collected seven samples from a single offshore location in the Mediterranean at 15 m depth intervals during a period of strong stratification, as well as two more samples during the winter when the photic water column was mixed. We were able to recover 94 new metagenome-assembled genomes (MAGs) from these metagenomes and examine the distribution of key marine microbes within the photic zone using metagenomic recruitment.

RESULTS

Our results showed significant differences in the microbial composition of different layers within the stratified photic water column. The majority of microorganisms were confined to discreet horizontal layers of no more than 30 m (stenobathic). Only a few such as members of the SAR11 clade appeared at all depths (eurybathic). During the winter mixing period, only some groups of bloomers such as Pseudomonas were favored. Although most microbes appeared in both seasons, some groups like the SAR116 clade and some Bacteroidetes and Verrucomicrobia seemed to disappear during the mixing period. Furthermore, we found that some microbes previously considered seasonal (e.g., Archaea or Actinobacteria) were living in deeper layers within the photic zone during the stratification period. A strong depth-related specialization was detected, not only at the taxonomic level but also at the functional level, even within the different clades, for the manipulation and uptake of specific polysaccharides. Rhodopsin sequences (green or blue) also showed narrow depth distributions that correlated with the taxonomy of the microbe in which they were found but not with depth.

CONCLUSIONS

Although limited to a single location in the Mediterranean, this study has profound implications for our understanding of how marine microbial communities vary with depth within the photic zone when stratified. Our results highlight the importance of collecting samples at different depths in the water column when comparing seasonal variations and have important ramifications for global marine studies that most often take samples from only one single depth. Furthermore, our perspective and approaches (metagenomic assembly and recruitment) are broadly applicable to other metagenomic studies.

High Prevalence of Quorum-Sensing and Quorum-Quenching Activity among Cultivable Bacteria and Metagenomic Sequences in the Mediterranean Sea

There is increasing evidence being accumulated regarding the importance of N-acyl homoserine lactones (AHL)-mediated quorum-sensing (QS) and quorum-quenching (QQ) processes in the marine environment, but in most cases, data has been obtained from specific microhabitats, and subsequently little is known regarding these activities in free-living marine bacteria. The QS and QQ activities among 605 bacterial isolates obtained at 90 and 2000 m depths in the Mediterranean Sea were analyzed. Additionally, putative QS and QQ sequences were searched in metagenomic data obtained at different depths (15-2000 m) at the same sampling site. The number of AHL producers was higher in the 90 m sample (37.66%) than in the 2000 m sample (4.01%). However, the presence of QQ enzymatic activity was 1.63-fold higher in the 2000 m sample. The analysis of putative QQ enzymes in the metagenomes supports the relevance of QQ processes in the deepest samples, found in cultivable bacteria. Despite the unavoidable biases in the cultivation methods and biosensor assays and the possible promiscuous activity of the QQ enzymes retrieved in the metagenomic analysis, the results indicate that AHL-related QS and QQ processes could be common activity in the marine environment.

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

Marine Group II archaea are widely distributed in global oceans and dominate the total archaeal community within the upper euphotic zone of temperate waters. However, factors controlling the distribution of MGII are poorly delineated and the physiology and ecological functions of these still-uncultured organisms remain elusive. In this study, we investigated the planktonic MGII associated with particles and in free-living forms in the Pearl River Estuary (PRE) over a 10-month period. We detected high abundance of particle-associated MGII in PRE (up to ∼10⁸ 16S rRNA gene copies/l), which was around 10-fold higher than the free-living MGII in the same region, and an order of magnitude higher than previously reported in other marine environments. 10‰ salinity appeared to be a threshold value for these MGII because MGII abundance decreased sharply below it. Above 10‰ salinity, the abundance of MGII on the particles was positively correlated with phototrophs and MGII in the surface water was negatively correlated with irradiance. However, the abundances of those free-living MGII showed positive correlations with salinity and temperature, suggesting the different physiological characteristics between particle-attached and free-living MGIIs. A nearly completely assembled metagenome, MGIIa_P, was recovered using metagenome binning methods. Compared with the other two MGII genomes from surface ocean, MGIIa_P contained higher proportions of glycoside hydrolases, indicating the ability of MGIIa_P to hydrolyse glycosidic bonds in complex sugars in PRE. MGIIa_P is the first assembled MGII metagenome containing a catalase gene, which might be involved in scavenging reactive oxygen species generated by the abundant phototrophs in the eutrophic PRE. Our study presented the widespread and high abundance of MGII in the water columns of PRE, and characterized the determinant abiotic factors affecting their distribution. Their association with heterotrophs, preference for particles and resourceful metabolic traits indicate MGII might play a significant role in metabolising organic matters in the PRE and other temperate estuarine systems.

Genome diversity of marine phages recovered from Mediterranean metagenomes: Size matters

Marine viruses play a critical role not only in the global geochemical cycles but also in the biology and evolution of their hosts. Despite their importance, viral diversity remains underexplored mostly due to sampling and cultivation challenges. Direct sequencing approaches such as viromics has provided new insights into the marine viral world. As a complementary approach, we analysed 24 microbial metagenomes (>0.2 μm size range) obtained from six sites in the Mediterranean Sea that vary by depth, season and filter used to retrieve the fraction. Filter-size comparison showed a significant number of viral sequences that were retained on the larger-pore filters and were different from those found in the viral fraction from the same sample, indicating that some important viral information is missing using only assembly from viromes. Besides, we were able to describe 1,323 viral genomic fragments that were more than 10Kb in length, of which 36 represented complete viral genomes including some of them retrieved from a cross-assembly from different metagenomes. Host prediction based on sequence methods revealed new phage groups belonging to marine prokaryotes like SAR11, Cyanobacteria or SAR116. We also identified the first complete virophage from deep seawater and a new endemic clade of the recently discovered Marine group II Euryarchaeota virus. Furthermore, analysis of viral distribution using metagenomes and viromes indicated that most of the new phages were found exclusively in the Mediterranean Sea and some of them, mostly the ones recovered from deep metagenomes, do not recruit in any database probably indicating higher variability and endemicity in Mediterranean bathypelagic waters. Together these data provide the first detailed picture of genomic diversity, spatial and depth variations of viral communities within the Mediterranean Sea using metagenome assembly.

New insights into marine group III Euryarchaeota, from dark to light

Marine Euryarchaeota remain among the least understood major components of marine microbial communities. Marine group II Euryarchaeota (MG-II) are more abundant in surface waters (4-20% of the total prokaryotic community), whereas marine group III Euryarchaeota (MG-III) are generally considered low-abundance members of deep mesopelagic and bathypelagic communities. Using genome assembly from direct metagenome reads and metagenomic fosmid clones, we have identified six novel MG-III genome sequence bins from the photic zone (Epi1-6) and two novel bins from deep-sea samples (Bathy1-2). Genome completeness in those genome bins varies from 44% to 85%. Photic-zone MG-III bins corresponded to novel groups with no similarity, and significantly lower GC content, when compared with previously described deep-MG-III genome bins. As found in many other epipelagic microorganisms, photic-zone MG-III bins contained numerous photolyase and rhodopsin genes, as well as genes for peptide and lipid uptake and degradation, suggesting a photoheterotrophic lifestyle. Phylogenetic analysis of these photolyases and rhodopsins as well as their genomic context suggests that these genes are of bacterial origin, supporting the hypothesis of an MG-III ancestor that lived in the dark ocean. Epipelagic MG-III occur sporadically and in relatively small proportions in marine plankton, representing only up to 0.6% of the total microbial community reads in metagenomes. None of the reconstructed epipelagic MG-III genomes were present in metagenomes from aphotic zone depths or from high latitude regions. Most low-GC bins were highly enriched at the deep chlorophyll maximum zones, with the exception of Epi1, which appeared evenly distributed throughout the photic zone worldwide.

Comparative metagenomics unveils functions and genome features of microbialite-associated communities along a depth gradient

Modern microbialites are often used as analogs of Precambrian stromatolites; therefore, studying the metabolic interplay within their associated microbial communities can help formulating hypotheses on their formation and long-term preservation within the fossil record. We performed a comparative metagenomic analysis of microbialite samples collected at two sites and along a depth gradient in Lake Alchichica (Mexico). The community structure inferred from single-copy gene family identification and long-contig (>10 kb) assignation, consistently with previous rRNA gene surveys, showed a wide prokaryotic diversity dominated by Alphaproteobacteria, Gammaproteobacteria, Cyanobacteria, and Bacteroidetes, while eukaryotes were largely dominated by green algae or diatoms. Functional analyses based on RefSeq, COG and SEED assignations revealed the importance of housekeeping functions, with an overrepresentation of genes involved in carbohydrate metabolism, as compared with other metabolic capacities. The search for genes diagnostic of specific metabolic functions revealed the important involvement of Alphaproteobacteria in anoxygenic photosynthesis and sulfide oxidation, and Cyanobacteria in oxygenic photosynthesis and nitrogen fixation. Surprisingly, sulfate reduction appeared negligible. Comparative analyses suggested functional similarities among various microbial mat and microbialite metagenomes as compared with soil or oceans, but showed differences in microbial processes among microbialite types linked to local environmental conditions.

Genomes of Abundant and Widespread Viruses from the Deep Ocean

The deep sea is a massive, largely oligotrophic ecosystem, stretched over nearly 65% of the planet’s surface. Deep-sea planktonic communities are almost completely dependent upon organic carbon sinking from the productive surface, forming a vital component of global biogeochemical cycles. However, despite their importance, viruses from the deep ocean remain largely unknown. Here, we describe the first complete genomes of deep-sea viruses assembled from metagenomic fosmid libraries. “Candidatus Pelagibacter” (SAR11) phage HTVC010P and Puniceispirillum phage HMO-2011 are considered the most abundant cultured marine viruses known to date. Remarkably, some of the viruses described here recruited as many reads from deep waters as these viruses do in the photic zone, and, considering the gigantic scale of the bathypelagic habitat, these genomes provide information about what could be some of the most abundant viruses in the world at large. Their role in the viral shunt in the global ocean could be very significant. Despite the challenges encountered in inferring the identity of their hosts, we identified one virus predicted to infect members of the globally distributed SAR11 cluster. We also identified a number of putative proviruses from diverse taxa, including deltaproteobacteria, bacteroidetes, SAR11, and gammaproteobacteria. Moreover, our findings also indicate that lysogeny is the preferred mode of existence for deep-sea viruses inhabiting an energy-limited environment, in sharp contrast to the predominantly lytic lifestyle of their photic-zone counterparts. Some of the viruses show a widespread distribution, supporting the tenet “everything is everywhere” for the deep-ocean virome.

The deep sea is among the largest known habitats and a critical cog in biogeochemical cycling but remains underexplored in its microbiology. Even more than is the case for its prokaryotic community, our knowledge of its viral component has remained limited by the paucity of information provided by studies dependent upon short sequence fragments. In this work, we attempt to fill this existing gap by using a combination of classical fosmid libraries with next-generation sequencing and assembly to recover long viral genomic fragments. We have sequenced ca. 6,000 fosmids from two metagenomics libraries made from prokaryotic biomass from the deep Mediterranean Sea and recovered twenty-eight complete viral genomes, all of them novel and quite distinct from all previously described viral genomes. They are preferentially found in deeper waters and are widely distributed all over the oceans. To our knowledge, this is the first report on complete and cosmopolitan viral genomes from the bathypelagic habitat.

Identification of Free-Living and Particle-Associated Microbial Communities Present in Hadal Regions of the Mariana Trench

Relatively few studies have described the microbial populations present in ultra-deep hadal environments, largely as a result of difficulties associated with sampling. Here we report Illumina-tag V6 16S rRNA sequence-based analyses of the free-living and particle-associated microbial communities recovered from locations within two of the deepest hadal sites on Earth, the Challenger Deep (10,918 meters below surface-mbs) and the Sirena Deep (10,667 mbs) within the Mariana Trench, as well as one control site (Ulithi Atoll, 761 mbs). Seawater samples were collected using an autonomous lander positioned ~1 m above the seafloor. The bacterial populations within the Mariana Trench bottom water samples were dissimilar to other deep-sea microbial communities, though with overlap with those of diffuse flow hydrothermal vents and deep-subsurface locations. Distinct particle-associated and free-living bacterial communities were found to exist. The hadal bacterial populations were also markedly different from one another, indicating the likelihood of different chemical conditions at the two sites. In contrast to the bacteria, the hadal archaeal communities were more similar to other less deep datasets and to each other due to an abundance of cosmopolitan deep-sea taxa. The hadal communities were enriched in 34 bacterial and 4 archaeal operational taxonomic units (OTUs) including members of the Gammaproteobacteria, Epsilonproteobacteria, Marinimicrobia, Cyanobacteria, Deltaproteobacteria, Gemmatimonadetes, Atribacteria, Spirochaetes, and Euryarchaeota. Sequences matching cultivated piezophiles were notably enriched in the Challenger Deep, especially within the particle-associated fraction, and were found in higher abundances than in other hadal studies, where they were either far less prevalent or missing. Our results indicate the importance of heterotrophy, sulfur-cycling, and methane and hydrogen utilization within the bottom waters of the deeper regions of the Mariana Trench, and highlight novel community features of these extreme habitats.

Genomic and transcriptomic evidence for scavenging of diverse organic compounds by widespread deep-sea archaea

Microbial activity is one of the most important processes to mediate the flux of organic carbon from the ocean surface to the seafloor. However, little is known about the microorganisms that underpin this key step of the global carbon cycle in the deep oceans. Here we present genomic and transcriptomic evidence that five ubiquitous archaeal groups actively use proteins, carbohydrates, fatty acids and lipids as sources of carbon and energy at depths ranging from 800 to 4,950 m in hydrothermal vent plumes and pelagic background seawater across three different ocean basins. Genome-enabled metabolic reconstructions and gene expression patterns show that these marine archaea are motile heterotrophs with extensive mechanisms for scavenging organic matter. Our results shed light on the ecological and physiological properties of ubiquitous marine archaea and highlight their versatile metabolic strategies in deep oceans that might play a critical role in global carbon cycling.

Marine Group II Archaea, potentially important players in the global ocean carbon cycle

Marine Group (MG) I (currently known as Thaumarchaeota) and MG II Archaea were first reported over two decades ago. While significant progress has been made on MG I microbiology and ecology, the progress on MG II has been noticeably slower. The common understanding is that while MG I mainly function as chemolithoautotrophs and occur predominantly in the deep ocean, MG II reside mostly in the photic zone and live heterotrophically. Studies to date have shown that MG II are abundant in the marine aquatic environment and display great seasonal and spatial variation and phylogenetic diversity. They also show unique patterns of organic carbon degradation and their energy requirements may be augmented by light in the photic zone. However, no pure culture of MG II has been obtained and thus their precise ecological role remains elusive.

Environmental Variables Shaping the Ecological Niche of Thaumarchaeota in Soil: Direct and Indirect Causal Effects

To find environmental variables (EVs) shaping the ecological niche of the archaeal phylum Thaumarchaeota in terrestrial environments, we determined the abundance of Thaumarchaeota in various soil samples using real-time PCR targeting thaumarchaeotal 16S rRNA gene sequences. We employed our previously developed primer, THAUM-494, which had greater coverage for Thaumarchaeota and lower tolerance to nonthaumarchaeotal taxa than previous Thaumarchaeota-directed primers. The relative abundance estimates (RVs) of Thaumarchaeota (RTHAUM), Archaea (RARCH), and Bacteria (RBACT) were subjected to a series of statistical analyses. Redundancy analysis (RDA) showed a significant (p < 0.05) canonical relationship between RVs and EVs. Negative causal relationships between RTHAUM and nutrient level-related EVs were observed in an RDA biplot. These negative relationships were further confirmed by correlation and regression analyses. Total nitrogen content (TN) appeared to be the EV that affected RTHAUM most strongly, and total carbon content (TC), which reflected the content of organic matter (OM), appeared to be the EV that affected it least. However, in the path analysis, a path model indicated that TN might be a mediator EV that could be controlled directly by the OM. Additionally, another path model implied that water content (WC) might also indirectly affect RTHAUM by controlling ammonium nitrogen (NH4+-N) level through ammonification. Thus, although most directly affected by NH4+-N, RTHAUM could be ultimately determined by OM content, suggesting that Thaumarchaeota could prefer low-OM or low-WC conditions, because either of these EVs could subsequently result in low levels of NH4+-N in soil.

Bacterial gene import and mesophilic adaptation in archaea

It is widely believed that the archaeal ancestor was hyperthermophilic, but during archaeal evolution, several lineages - including haloarchaea and their sister methanogens, the Thaumarchaeota, and the uncultured Marine Group II and Marine Group III Euryarchaeota (MGII/III) - independently adapted to lower temperatures. Recent phylogenomic studies suggest that the ancestors of these lineages were recipients of massive horizontal gene transfer from bacteria. Many of the acquired genes, which are often involved in metabolism and cell envelope biogenesis, were convergently acquired by distant mesophilic archaea. In this Opinion article, we explore the intriguing hypothesis that the import of these bacterial genes was crucial for the adaptation of archaea to mesophilic lifestyles.

Shifts in the meso- and bathypelagic archaea communities composition during recovery and short-term handling of decompressed deep-sea samples

Dark ocean microbial communities are actively involved in chemoautotrophic and anaplerotic fixation of bicarbonate. Thus, aphotic pelagic realm of the ocean might represent a significant sink of CO2 and source of primary production. However, the estimated metabolic activities in the dark ocean are fraught with uncertainties. Typically, deep-sea samples are recovered to the sea surface for downstream processing on deck. Shifts in ambient settings, associated with such treatments, can likely change the metabolic activity and community structure of deep-sea adapted autochthonous microbial populations. To estimate influence of recovery and short-term handling of deep-sea samples, we monitored the succession of bathypelagic microbial community during its 3 days long on deck incubation. We demonstrated that at the end of exposition, the deep-sea archaeal population decreased threefold, whereas the bacterial fraction doubled in size. As revealed by phylogenetic analyses of amoA gene transcripts, dominance of the active ammonium-oxidizing bathypelagic Thaumarchaeota groups shifted over time very fast. These findings demonstrated the simultaneous existence of various 'deep-sea ecotypes', differentially reacting to the sampling and downstream handling. Our study supports the hypothesis that metabolically active members of meso- and bathypelagic Thaumarchaeota possess the habitat-specific distribution, metabolic complexity and genetic divergence at subpopulation level.

Ecophysiology of uncultivated marine euryarchaea is linked to particulate organic matter

Particles in aquatic environments host distinct communities of microbes, yet the evolution of particle-specialized taxa and the extent to which specialized microbial metabolism is associated with particles is largely unexplored. Here, we investigate the hypothesis that a widely distributed and uncultivated microbial group—the marine group II euryarchaea (MGII)—interacts with living and detrital particulate organic matter (POM) in the euphotic zone of the central California Current System. Using fluorescent in situ hybridization, we verified the association of euryarchaea with POM. We further quantified the abundance and distribution of MGII 16 S ribosomal RNA genes in size-fractionated seawater samples and compared MGII functional capacity in metagenomes from the same fractions. The abundance of MGII in free-living and >3 μm fractions decreased with increasing distance from the coast, whereas MGII abundance in the 0.8–3 μm fraction remained constant. At several offshore sites, MGII abundance was highest in particle fractions, indicating that particle-attached MGII can outnumber free-living MGII under oligotrophic conditions. Compared with free-living MGII, the genome content of MGII in particle-associated fractions exhibits an increased capacity for surface adhesion, transcriptional regulation and catabolism of high molecular weight substrates. Moreover, MGII populations in POM fractions are phylogenetically distinct from and more diverse than free-living MGII. Eukaryotic phytoplankton additions stimulated MGII growth in bottle incubations, providing the first MGII net growth rate measurements. These ranged from 0.47 to 0.54 d⁻¹. However, MGII were not recovered in whole-genome amplifications of flow-sorted picoeukaryotic phytoplankton and heterotrophic nanoflagellates, suggesting that MGII in particle fractions are not physically attached to living POM. Collectively, our results support a linkage between MGII ecophysiology and POM, implying that marine archaea have a role in elemental cycling through interactions with particles.

A new class of marine Euryarchaeota group II from the Mediterranean deep chlorophyll maximum

We have analyzed metagenomic fosmid clones from the deep chlorophyll maximum (DCM), which, by genomic parameters, correspond to the 16S ribosomal RNA (rRNA)-defined marine Euryarchaeota group IIB (MGIIB). The fosmid collections associated with this group add up to 4 Mb and correspond to at least two species within this group. From the proposed essential genes contained in the collections, we infer that large sections of the conserved regions of the genomes of these microbes have been recovered. The genomes indicate a photoheterotrophic lifestyle, similar to that of the available genome of MGIIA (assembled from an estuarine metagenome in Puget Sound, Washington Pacific coast), with a proton-pumping rhodopsin of the same kind. Several genomic features support an aerobic metabolism with diversified substrate degradation capabilities that include xenobiotics and agar. On the other hand, these MGIIB representatives are non-motile and possess similar genome size to the MGIIA-assembled genome, but with a lower GC content. The large phylogenomic gap with other known archaea indicates that this is a new class of marine Euryarchaeota for which we suggest the name Thalassoarchaea. The analysis of recruitment from available metagenomes indicates that the representatives of group IIB described here are largely found at the DCM (ca. 50 m deep), in which they are abundant (up to 0.5% of the reads), and at the surface mostly during the winter mixing, which explains formerly described 16S rRNA distribution patterns. Their uneven representation in environmental samples that are close in space and time might indicate sporadic blooms.

Pangenome evidence for extensive interdomain horizontal transfer affecting lineage core and shell genes in uncultured planktonic thaumarchaeota and euryarchaeota

Horizontal gene transfer (HGT) is an important force in evolution, which may lead, among other things, to the adaptation to new environments by the import of new metabolic functions. Recent studies based on phylogenetic analyses of a few genome fragments containing archaeal 16S rRNA genes and fosmid-end sequences from deep-sea metagenomic libraries have suggested that marine planktonic archaea could be affected by high HGT frequency. Likewise, a composite genome of an uncultured marine euryarchaeote showed high levels of gene sequence similarity to bacterial genes. In this work, we ask whether HGT is frequent and widespread in genomes of these marine archaea, and whether HGT is an ancient and/or recurrent phenomenon. To answer these questions, we sequenced 997 fosmid archaeal clones from metagenomic libraries of deep-Mediterranean waters (1,000 and 3,000 m depth) and built comprehensive pangenomes for planktonic Thaumarchaeota (Group I archaea) and Euryarchaeota belonging to the uncultured Groups II and III Euryarchaeota (GII/III-Euryarchaeota). Comparison with available reference genomes of Thaumarchaeota and a composite marine surface euryarchaeote genome allowed us to define sets of core, lineage-specific core, and shell gene ortholog clusters for the two archaeal lineages. Molecular phylogenetic analyses of all gene clusters showed that 23.9% of marine Thaumarchaeota genes and 29.7% of GII/III-Euryarchaeota genes had been horizontally acquired from bacteria. HGT is not only extensive and directional but also ongoing, with high HGT levels in lineage-specific core (ancient transfers) and shell (recent transfers) genes. Many of the acquired genes are related to metabolism and membrane biogenesis, suggesting an adaptive value for life in cold, oligotrophic oceans. We hypothesize that the acquisition of an important amount of foreign genes by the ancestors of these archaeal groups significantly contributed to their divergence and ecological success.

Novel PCR primers for the archaeal phylum Thaumarchaeota designed based on the comparative analysis of 16S rRNA gene sequences

Based on comparative phylogenetic analysis of 16S rRNA gene sequences deposited in an RDP database, we constructed a local database of thaumarchaeotal 16S rRNA gene sequences and developed a novel PCR primer specific for the archaeal phylum Thaumarchaeota. Among 9,727 quality-filtered (chimeral-checked, size >1.2 kb) archaeal sequences downloaded from the RDP database, 1,549 thaumarchaeotal sequences were identified and included in our local database. In our study, Thaumarchaeota included archaeal groups MG-I, SAGMCG-I, SCG, FSCG, RC, and HWCG-III, forming a monophyletic group in the phylogenetic tree. Cluster analysis revealed 114 phylotypes for Thaumarchaeota. The majority of the phylotypes (66.7%) belonged to the MG-I and SCG, which together contained most (93.9%) of the thaumarchaeotal sequences in our local database. A phylum-directed primer was designed from a consensus sequence of the phylotype sequences, and the primer's specificity was evaluated for coverage and tolerance both in silico and empirically. The phylum-directed primer, designated THAUM-494, showed >90% coverage for Thaumarchaeota and <1% tolerance to non-target taxa, indicating high specificity. To validate this result experimentally, PCRs were performed with THAUM-494 in combination with a universal archaeal primer (ARC917R or 1017FAR) and DNAs from five environmental samples to construct clone libraries. THAUM-494 showed a satisfactory specificity in empirical studies, as expected from the in silico results. Phylogenetic analysis of 859 cloned sequences obtained from 10 clone libraries revealed that >95% of the amplified sequences belonged to Thaumarchaeota. The most frequently sampled thaumarchaeotal subgroups in our samples were SCG, MG-I, and SAGMCG-I. To our knowledge, THAUM-494 is the first phylum-level primer for Thaumarchaeota. Furthermore, the high coverage and low tolerance of THAUM-494 will make it a potentially valuable tool in understanding the phylogenetic diversity and ecological niche of Thaumarchaeota.

Crenarchaeal heterotrophy in salt marsh sediments

Mesophilic Crenarchaeota (also known as Thaumarchaeota) are ubiquitous and abundant in marine habitats. However, very little is known about their metabolic function in situ. In this study, salt marsh sediments from New Jersey were screened via stable isotope probing (SIP) for heterotrophy by amending with a single (13)C-labeled compound (acetate, glycine or urea) or a complex (13)C-biopolymer (lipids, proteins or growth medium (ISOGRO)). SIP incubations were done at two substrate concentrations (30-150 μM; 2-10 mg ml(-1)), and (13)C-labeled DNA was analyzed by terminal restriction fragment length polymorphism (TRFLP) analysis of 16S rRNA genes. To test for autotrophy, an amendment with (13)C-bicarbonate was also performed. Our SIP analyses indicate salt marsh crenarchaea are heterotrophic, double within 2-3 days and often compete with heterotrophic bacteria for the same organic substrates. A clone library of (13)C-amplicons was screened to find matches to the (13)C-TRFLP peaks, with seven members of the Miscellaneous Crenarchaeal Group and seven members from the Marine Group 1.a Crenarchaeota being discerned. Some of these crenarchaea displayed a preference for particular carbon sources, whereas others incorporated nearly every (13)C-substrate provided. The data suggest salt marshes may be an excellent model system for studying crenarchaeal metabolic capabilities and can provide information on the competition between crenarchaea and other microbial groups to improve our understanding of microbial ecology.

Evolutionary analysis of a streamlined lineage of surface ocean Roseobacters

The vast majority of surface ocean bacteria are uncultivated. Compared with their cultured relatives, they frequently exhibit a streamlined genome, reduced G+C content and distinct gene repertoire. These genomic traits are relevant to environmental adaptation, and have generally been thought to become fixed in marine bacterial populations through selection. Using single-cell genomics, we sequenced four uncultivated cells affiliated with the ecologically relevant Roseobacter clade and used a composition-heterogeneous Bayesian phylogenomic model to resolve these single-cell genomes into a new clade. This lineage has no representatives in culture, yet accounts for ∼35% of Roseobacters in some surface ocean waters. Analyses of multiple genomic traits, including genome size, G+C content and percentage of noncoding DNA, suggest that these single cells are representative of oceanic Roseobacters but divergent from isolates. Population genetic analyses showed that substitution of physicochemically dissimilar amino acids and replacement of G+C-rich to G+C-poor codons are accelerated in the uncultivated clade, processes that are explained equally well by genetic drift as by the more frequently invoked explanation of natural selection. The relative importance of drift vs selection in this clade, and perhaps in other marine bacterial clades with streamlined G+C-poor genomes, remains unresolved until more evidence is accumulated.

Environmental distribution of two widespread uncultured freshwater Euryarchaeota clades unveiled by specific primers and quantitative PCR

Quantitative environmental distribution of two widely distributed uncultured freshwater Euryarchaeota with unknown functional role was explored by newly designed quantitative PCR primers targeting the 16S rRNA gene of clades Miscellaneous Euryarchaeota Group (MEG, containing the groups pMC2A384 and VALII/Eury4) and Deep-Sea Euryarchaeotal Groups (DSEG, targeting the cluster named VALIII containing the DHVE-3/DSEG, BC07-2A-27/DSEG-3 and DSEG-2 groups), respectively. The summer surface plankton of 28 lakes was analysed, and one additional dimictic deep alpine lake, Lake Redon, was temporally and vertically surveyed covering seasonal limnological variability. A trophic range between 0.2 and 5.2 μg l(-1) Chl a, and pH span from 3.8 to 9.5 was explored at altitudes between 632 and 2590 m above sea level. The primers showed to be highly selective with c. 85% coverage and 100% specificity. Only pH significantly explained the changes observed in gene abundances and environment. In Lake Redon, DSEG bloomed in deep stratified waters both in summer and early spring, and MEG at intermediate depths during the ice-cover period. Overall, MEG and DSEG showed a differential ecological distribution although correlational analyses indicated lack of coupling of both Euryarchaeota with phytoplankton (chlorophyll a). However, an intriguing positive and significant relationship was found between DSEG and putative ammonia oxidizing thaumarchaeota.

Archaea in biogeochemical cycles

Archaea constitute a considerable fraction of the microbial biomass on Earth. Like Bacteria they have evolved a variety of energy metabolisms using organic and/or inorganic electron donors and acceptors, and many of them are able to fix carbon from inorganic sources. Archaea thus play crucial roles in the Earth's global geochemical cycles and influence greenhouse gas emissions. Methanogenesis and anaerobic methane oxidation are important steps in the carbon cycle; both are performed exclusively by anaerobic archaea. Oxidation of ammonia to nitrite is performed by Thaumarchaeota. They represent the only archaeal group that resides in large numbers in the global aerobic terrestrial and marine environments on Earth. Sulfur-dependent archaea are confined mostly to hot environments, but metal leaching by acidophiles and reduction of sulfate by anaerobic, nonthermophilic methane oxidizers have a potential impact on the environment. The metabolisms of a large number of archaea, in particular those dominating the subsurface, remain to be explored.

Community transcriptomic assembly reveals microbes that contribute to deep-sea carbon and nitrogen cycling

The deep ocean is an important component of global biogeochemical cycles because it contains one of the largest pools of reactive carbon and nitrogen on earth. However, the microbial communities that drive deep-sea geochemistry are vastly unexplored. Metatranscriptomics offers new windows into these communities, but it has been hampered by reliance on genome databases for interpretation. We reconstructed the transcriptomes of microbial populations from Guaymas Basin, in the deep Gulf of California, through shotgun sequencing and de novo assembly of total community RNA. Many of the resulting messenger RNA (mRNA) contiguous sequences contain multiple genes, reflecting co-transcription of operons, including those from dominant members. Also prevalent were transcripts with only limited representation (2.8 times coverage) in a corresponding metagenome, including a considerable portion (1.2 Mb total assembled mRNA sequence) with similarity (96%) to a marine heterotroph, Alteromonas macleodii. This Alteromonas and euryarchaeal marine group II populations displayed abundant transcripts from amino-acid transporters, suggesting recycling of organic carbon and nitrogen from amino acids. Also among the most abundant mRNAs were catalytic subunits of the nitrite oxidoreductase complex and electron transfer components involved in nitrite oxidation. These and other novel genes are related to novel Nitrospirae and have limited representation in accompanying metagenomic data. High throughput sequencing of 16S ribosomal RNA (rRNA) genes and rRNA read counts confirmed that Nitrospirae are minor yet widespread members of deep-sea communities. These results implicate a novel bacterial group in deep-sea nitrite oxidation, the second step of nitrification. This study highlights metatranscriptomic assembly as a valuable approach to study microbial communities.

Archaeal amoA gene diversity points to distinct biogeography of ammonia-oxidizing Crenarchaeota in the ocean

Mesophilic ammonia-oxidizing Archaea (AOA) are abundant in a diverse range of marine environments, including the deep ocean, as revealed by the quantification of the archaeal amoA gene encoding the alpha-subunit of the ammonia monooxygenase. Using two different amoA primer sets, two distinct ecotypes of marine Crenarchaeota Group I (MCGI) were detected in the waters of the tropical Atlantic and the coastal Arctic. The HAC-AOA ecotype (high ammonia concentration AOA) was ≈ 8000 times and 15 times more abundant in the coastal Arctic and the top 300 m layer of the open equatorial Atlantic, respectively, than the LAC-AOA (low ammonia concentration AOA) ecotype. In contrast, the LAC-AOA ecotype dominated the lower meso- and bathypelagic waters of the tropical Atlantic (≈ 50 times more abundant than the HAC-AOA) where ammonia concentrations are well below the detection limit using conventional spectrophotometric or fluorometric methods. Cluster analysis of the sequences from the clone libraries obtained by the two amoA primer sets revealed two phylogenetically distinct clusters. Taken together, our results suggest the presence of two ecotypes of archaeal ammonia oxidizers corresponding to the medium (1.24 µM on average in the coastal Arctic) and low ammonia concentration (< 0.01 µM) in the shallow and the deep waters respectively.

Life in the dark: metagenomic evidence that a microbial slime community is driven by inorganic nitrogen metabolism

Beneath Australia's large, dry Nullarbor Plain lies an extensive underwater cave system, where dense microbial communities known as 'slime curtains' are found. These communities exist in isolation from photosynthetically derived carbon and are presumed to be chemoautotrophic. Earlier work found high levels of nitrite and nitrate in the cave waters and a high relative abundance of Nitrospirae in bacterial 16S rRNA clone libraries. This suggested that these communities may be supported by nitrite oxidation, however, details of the inorganic nitrogen cycling in these communities remained unclear. Here we report analysis of 16S rRNA amplicon and metagenomic sequence data from the Weebubbie cave slime curtain community. The microbial community is comprised of a diverse assortment of bacterial and archaeal genera, including an abundant population of Thaumarchaeota. Sufficient thaumarchaeotal sequence was recovered to enable a partial genome sequence to be assembled, which showed considerable synteny with the corresponding regions in the genome of the autotrophic ammonia oxidiser Nitrosopumilus maritimus SCM1. This partial genome sequence, contained regions with high sequence identity to the ammonia mono-oxygenase operon and carbon fixing 3-hydroxypropionate/4-hydroxybutyrate cycle genes of N. maritimus SCM1. Additionally, the community, as a whole, included genes encoding key enzymes for inorganic nitrogen transformations, including nitrification and denitrification. We propose that the Weebubbie slime curtain community represents a distinctive microbial ecosystem, in which primary productivity is due to the combined activity of archaeal ammonia-oxidisers and bacterial nitrite oxidisers.

Ocean currents shape the microbiome of Arctic marine sediments

Prokaryote communities were investigated on the seasonally stratified Alaska Beaufort Shelf (ABS). Water and sediment directly underlying water with origin in the Arctic, Pacific or Atlantic oceans were analyzed by pyrosequencing and length heterogeneity-PCR in conjunction with physicochemical and geographic distance data to determine what features structure ABS microbiomes. Distinct bacterial communities were evident in all water masses. Alphaproteobacteria explained similarity in Arctic surface water and Pacific derived water. Deltaproteobacteria were abundant in Atlantic origin water and drove similarity among samples. Most archaeal sequences in water were related to unclassified marine Euryarchaeota. Sediment communities influenced by Pacific and Atlantic water were distinct from each other and pelagic communities. Firmicutes and Chloroflexi were abundant in sediment, although their distribution varied in Atlantic and Pacific influenced sites. Thermoprotei dominated archaea in Pacific influenced sediments and Methanomicrobia dominated in methane-containing Atlantic influenced sediments. Length heterogeneity-PCR data from this study were analyzed with data from methane-containing sediments in other regions. Pacific influenced ABS sediments clustered with Pacific sites from New Zealand and Chilean coastal margins. Atlantic influenced ABS sediments formed another distinct cluster. Density and salinity were significant structuring features on pelagic communities. Porosity co-varied with benthic community structure across sites and methane did not. This study indicates that the origin of water overlying sediments shapes benthic communities locally and globally and that hydrography exerts greater influence on microbial community structure than the availability of methane.

Comparative genomic analysis and benzene, toluene, ethylbenzene, and o-, m-, and p-xylene (BTEX) degradation pathways of Pseudoxanthomonas spadix BD-a59

Pseudoxanthomonas spadix BD-a59, isolated from gasoline-contaminated soil, has the ability to degrade all six BTEX (benzene, toluene, ethylbenzene, and o-, m-, and p-xylene) compounds. The genomic features of strain BD-a59 were analyzed bioinformatically and compared with those of another fully sequenced Pseudoxanthomonas strain, P. suwonensis 11-1, which was isolated from cotton waste compost. The genome of strain BD-a59 differed from that of strain 11-1 in many characteristics, including the number of rRNA operons, dioxygenases, monooxygenases, genomic islands (GIs), and heavy metal resistance genes. A high abundance of phage integrases and GIs and the patterns in several other genetic measures (e.g., GC content, GC skew, Karlin signature, and clustered regularly interspaced short palindromic repeat [CRISPR] gene homology) indicated that strain BD-a59's genomic architecture may have been altered through horizontal gene transfers (HGT), phage attack, and genetic reshuffling during its evolutionary history. The genes for benzene/toluene, ethylbenzene, and xylene degradations were encoded on GI-9, -13, and -21, respectively, which suggests that they may have been acquired by HGT. We used bioinformatics to predict the biodegradation pathways of the six BTEX compounds, and these pathways were proved experimentally through the analysis of the intermediates of each BTEX compound using a gas chromatograph and mass spectrometry (GC-MS). The elevated abundances of dioxygenases, monooxygenases, and rRNA operons in strain BD-a59 (relative to strain 11-1), as well as other genomic characteristics, likely confer traits that enhance ecological fitness by enabling strain BD-a59 to degrade hydrocarbons in the soil environment.

Untangling genomes from metagenomes: revealing an uncultured class of marine Euryarchaeota

Ecosystems are shaped by complex communities of mostly unculturable microbes. Metagenomes provide a fragmented view of such communities, but the ecosystem functions of major groups of organisms remain mysterious. To better characterize members of these communities, we developed methods to reconstruct genomes directly from mate-paired short-read metagenomes. We closed a genome representing the as-yet uncultured marine group II Euryarchaeota, assembled de novo from 1.7% of a metagenome sequenced from surface seawater. The genome describes a motile, photo-heterotrophic cell focused on degradation of protein and lipids and clarifies the origin of proteorhodopsin. It also demonstrates that high-coverage mate-paired sequence can overcome assembly difficulties caused by interstrain variation in complex microbial communities, enabling inference of ecosystem functions for uncultured members.

Diversity, physiology, and niche differentiation of ammonia-oxidizing archaea

Nitrification, the aerobic oxidation of ammonia to nitrate via nitrite, has been suggested to have been a central part of the global biogeochemical nitrogen cycle since the oxygenation of Earth. The cultivation of several ammonia-oxidizing archaea (AOA) as well as the discovery that archaeal ammonia monooxygenase (amo)-like gene sequences are nearly ubiquitously distributed in the environment and outnumber their bacterial counterparts in many habitats fundamentally revised our understanding of nitrification. Surprising insights into the physiological distinctiveness of AOA are mirrored by the recognition of the phylogenetic uniqueness of these microbes, which fall within a novel archaeal phylum now known as Thaumarchaeota. The relative importance of AOA in nitrification, compared to ammonia-oxidizing bacteria (AOB), is still under debate. This minireview provides a synopsis of our current knowledge of the diversity and physiology of AOA, the factors controlling their ecology, and their role in carbon cycling as well as their potential involvement in the production of the greenhouse gas nitrous oxide. It emphasizes the importance of activity-based analyses in AOA studies and formulates priorities for future research.

Metagenomic analysis of hadopelagic microbial assemblages thriving at the deepest part of Mediterranean Sea, Matapan-Vavilov Deep

The marine pelagic zone situated > 200 m below the sea level (bls) is the largest marine subsystem, comprising more than two-thirds of the oceanic volume. At the same time, it is one of the least explored ecosystems on Earth. Few large-scale environmental genomics studies have been undertaken to examine the phylogenetic diversity and functional gene repertoire of planktonic microbes present in mesopelagic and bathypelagic environments. Here, we present the description of the deep-sea microbial community thriving at > 4900 m depth in Matapan-Vavilov Deep (MVD). This canyon is the deepest site of Mediterranean Sea, with a deepest point located at approximately 5270 m, 56 km SW of city Pylos (Greece) in the Ionian Sea (36°34.00N, 21°07.44E). Comparative analysis of whole-metagenomic data revealed that unlike other deep-sea metagenomes, the prokaryotic diversity in MVD was extremely poor. The decline in the dark primary production rates, measured at 4908 m depth, was coincident with overwhelming dominance of copiotrophic Alteromonas macleodii'deep-ecotype' AltDE at the expense of other prokaryotes including those potentially involved in both autotrophic and anaplerotic CO(2) fixation. We also demonstrate the occurrence in deep-sea metagenomes of several clustered regularly interspaced short palindromic repeats systems.

The metatranscriptome of a deep-sea hydrothermal plume is dominated by water column methanotrophs and lithotrophs

Microorganisms mediate geochemical processes in deep-sea hydrothermal vent plumes, which are a conduit for transfer of elements and energy from the subsurface to the oceans. Despite this important microbial influence on marine geochemistry, the ecology and activity of microbial communities in hydrothermal plumes is largely unexplored. Here, we use a coordinated metagenomic and metatranscriptomic approach to compare microbial communities in Guaymas Basin hydrothermal plumes to background waters above the plume and in the adjacent Carmen Basin. Despite marked increases in plume total RNA concentrations (3-4 times) and microbially mediated manganese oxidation rates (15-125 times), plume and background metatranscriptomes were dominated by the same groups of methanotrophs and chemolithoautotrophs. Abundant community members of Guaymas Basin seafloor environments (hydrothermal sediments and chimneys) were not prevalent in the plume metatranscriptome. De novo metagenomic assembly was used to reconstruct genomes of abundant populations, including Marine Group I archaea, Methylococcaceae, SAR324 Deltaproteobacteria and SUP05 Gammaproteobacteria. Mapping transcripts to these genomes revealed abundant expression of genes involved in the chemolithotrophic oxidation of ammonia (amo), methane (pmo) and sulfur (sox). Whereas amo and pmo gene transcripts were abundant in both plume and background, transcripts of sox genes for sulfur oxidation from SUP05 groups displayed a 10-20-fold increase in plumes. We conclude that the biogeochemistry of Guaymas Basin hydrothermal plumes is mediated by microorganisms that are derived from seawater rather than from seafloor hydrothermal environments such as chimneys or sediments, and that hydrothermal inputs serve as important electron donors for primary production in the deep Gulf of California.

Metatranscriptomics of the marine sponge Geodia barretti: tackling phylogeny and function of its microbial community

Geodia barretti is a marine cold-water sponge harbouring high numbers of microorganisms. Significant rates of nitrification have been observed in this sponge, indicating a substantial contribution to nitrogen turnover in marine environments with high sponge cover. In order to get closer insights into the phylogeny and function of the active microbial community and the interaction with its host G. barretti, a metatranscriptomic approach was employed, using the simultaneous analysis of rRNA and mRNA. Of the 262 298 RNA-tags obtained by pyrosequencing, 92% were assigned to ribosomal RNA (ribo-tags). A total of 109 325 SSU rRNA ribo-tags revealed a detailed picture of the community, dominated by group SAR202 of Chloroflexi, candidate phylum Poribacteria and Acidobacteria, which was different in its composition from that obtained in clone libraries prepared form the same samples. Optimized assembly strategies allowed the reconstruction of full-length rRNA sequences from the short ribo-tags for more detailed phylogenetic studies of the dominant taxa. Cells of several phyla were visualized by FISH analyses for confirmation. Of the remaining 21 325 RNA-tags, 10 023 were assigned to mRNA-tags, based on similarities to genes in the databases. A wide range of putative functional gene transcripts from over 10 different phyla were identified among the bacterial mRNA-tags. The most abundant mRNAs were those encoding key metabolic enzymes of nitrification from ammonia-oxidizing archaea as well as candidate genes involved in related processes. Our analysis demonstrates the potential and limits of using a combined rRNA and mRNA approach to explore the microbial community profile, phylogenetic assignments and metabolic activities of a complex, but little explored microbial community.

Polyclonality of concurrent natural populations of Alteromonas macleodii

We have analyzed a natural population of the marine bacterium, Alteromonas macleodii, from a single sample of seawater to evaluate the genomic diversity present. We performed full genome sequencing of four isolates and 161 metagenomic fosmid clones, all of which were assigned to A. macleodii by sequence similarity. Out of the four strain genomes, A. macleodii deep ecotype (AltDE1) represented a different genome, whereas AltDE2 and AltDE3 were identical to the previously described AltDE. Although the core genome (~80%) had an average nucleotide identity of 98.51%, both AltDE and AltDE1 contained flexible genomic islands (fGIs), that is, genomic islands present in both genomes in the same genomic context but having different gene content. Some of the fGIs encode cell surface receptors known to be phage recognition targets, such as the O-chain of the lipopolysaccharide, whereas others have genes involved in physiological traits (e.g., nutrient transport, degradation, and metal resistance) denoting microniche specialization. The presence in metagenomic fosmids of genomic fragments differing from the sequenced strain genomes, together with the presence of new fGIs, indicates that there are at least two more A. macleodii clones present. The availability of three or more sequences overlapping the same genomic region also allowed us to estimate the frequency and distribution of recombination events among these different clones, indicating that these clustered near the genomic islands. The results indicate that this natural A. macleodii population has multiple clones with a potential for different phage susceptibility and exploitation of resources, within a seemingly unstructured habitat.