Teaching old dogs new tricks: genetic engineering methanogens.

Methanogenic archaea, which are integral to global carbon and nitrogen cycling, currently face challenges in genetic manipulation due to unique physiology and limited genetic tools. This review provides a survey of current and past developments in the genetic engineering of methanogens, including selection and counterselection markers, reporter systems, shuttle vectors, mutagenesis methods, markerless genetic exchange, and gene expression control. This review discusses genetic tools and emphasizes challenges tied to tool scarcity for specific methanogenic species. Mutagenesis techniques for methanogens, including physicochemical, transposon-mediated, liposome-mediated mutagenesis, and natural transformation, are outlined, along with achievements and challenges. Markerless genetic exchange strategies, such as homologous recombination and CRISPR/Cas-mediated genome editing, are also detailed. Finally, the review concludes by examining the control of gene expression in methanogens. The information presented underscores the urgent need for refined genetic tools in archaeal research. Despite historical challenges, recent advancements, notably CRISPR-based systems, hold promise for overcoming obstacles, with implications for global health, agriculture, climate change, and environmental engineering. This comprehensive review aims to bridge existing gaps in the literature, guiding future research in the expanding field of archaeal genetic engineering.

Temporal dynamics of volatile fatty acids profile, methane production, and prokaryotic community in an in vitro rumen fermentation system fed with maize silage.

Anaerobic in vitro fermentation is widely used to simulate rumen kinetics and study the microbiome and metabolite profiling in a controlled lab environment. However, a better understanding of the interplay between the temporal dynamics of fermentation kinetics, metabolic profiles, and microbial composition in in vitro rumen fermentation batch systems is required. To fill that knowledge gap, we conducted three in vitro rumen fermentations with maize silage as the substrate, monitoring total gas production (TGP), dry matter degradability (dDM), and methane (CH4) concentration at 6, 12, 24, 36, and 48 h in each fermentation. At each time point, we collected rumen fluid samples for microbiome analysis and volatile fatty acid (VFA) analysis. Amplicon sequencing of 16S rRNA genes (V4 region) was used to profile the prokaryotic community structure in the rumen during the fermentation process. As the fermentation time increased, dDM, TGP, VFA concentrations, CH4 concentration, and yield (mL CH4 per g DM at standard temperature and pressure (STP)) significantly increased. For the dependent variables, CH4 concentration and yield, as well as the independent variables TGP and dDM, polynomial equations were fitted. These equations explained over 85% of the data variability (R2 > 0.85) and suggest that TGP and dDM can be used as predictors to estimate CH4 production in rumen fermentation systems. Microbiome analysis revealed a dominance of Bacteroidota, Cyanobacteria, Desulfobacterota, Euryarchaeota, Fibrobacterota, Firmicutes, Patescibacteria, Proteobacteria, Spirochaetota, and Verrucomicrobiota. Significant temporal variations in Bacteroidota, Campylobacterota, Firmicutes, Proteobacteria, and Spirochaetota were detected. Estimates of alpha diversity based on species richness and the Shannon index showed no variation between fermentation time points. This study demonstrated that the in vitro fermentation characteristics of a given feed type (e.g., maize silage) can be predicted from a few parameters (CH4 concentration and yield, tVFA, acetic acid, and propionic acid) without running the actual in vitro trial if the rumen fluid is collected from similar donor cows. Although the dynamics of the rumen prokaryotes changed remarkably over time and in accordance with the fermentation kinetics, more time points between 0 and 24 h are required to provide more details about the microbial temporal dynamics at the onset of the fermentation.

Draft genome sequence of Thermococcus waiotapuensis WT1T, a thermophilic sulfur-dependent archaeon from the order Thermococcales.

Thermococcus waiotapuensis WT1T is a thermophilic, peptide, and amino acid-fermenting archaeon from the order Thermococcales. It was isolated from Waiotapu, Aotearoa-New Zealand, and has a genome size of 1.80 Mbp. The genome contains 2,000 total genes, of which 1,913 encode proteins and 46 encode tRNA.

Influence of aphotic haloclines and euxinia on organic biomarkers and microbial communities in a thalassohaline and alkaline volcanic crater lake.

Studies on microbial communities, and their associated organic biomarkers, that are found thriving in the aphotic euxinic waters in modern stratified ecosystems are scarce compared to those undertaken in euxinic photic zones. The Dziani Dzaha (Mayotte, Indian Ocean) is a tropical, saline, alkaline crater lake that has recently been presented as a modern analog of Proterozoic Oceans due to its thalassohaline classification (having water of marine origin) and specific biogeochemical characteristics. Continuous intense photosynthetic production and microbial mineralization keep most of the water column permanently aphotic and anoxic preventing the development of a euxinic (sulfidic and anoxic) photic zone despite a high sulfide/sulfate ratio and the presence of permanent or seasonal haloclines. In this study, the molecular composition of the organic matter in Lake Dziani Dzaha was investigated and compared to the microbial diversity evaluated through 16S rRNA gene amplicon sequencing, over two contrasting seasons (rainy vs. dry) that influence water column stratification. Depth profiles of organic biomarker concentrations (chlorophyll-a and lipid biomarkers) and bacterial and archaeal OTU abundances appeared to be strongly dependent on the presence of aphotic haloclines and euxinia. OTU abundances revealed the importance of specific haloalkaliphilic bacterial and archaeal assemblages in phytoplanktonic biomass recycling and the biogeochemical functioning of the lake, suggesting new haloalkaline non-phototrophic anaerobic microbial precursors for some of the lipid biomarkers. Uncultured Firmicutes from the family Syntrophomonadaceae (Clostridiales), and Bacteroidetes from the ML635J-40 aquatic group, emerged as abundant chemotrophic bacterial members in the anoxic or euxinic waters and were probably responsible for the production of short-chain n-alkenes, wax esters, diplopterol, and tetrahymanol. Halocline-dependent euxinia also had a strong impact on the archaeal community which was dominated by Woesearchaeota in the sulfide-free waters. In the euxinic waters, methanogenic Euryarchaeota from the Methanomicrobia, Thermoplasmata, and WSA2 classes dominated and were likely at the origin of common hydrocarbon biomarkers of methanogens (phytane, pentamethyl-eicosenes, and partially hydrogenated squalene).

Influence of Environmental Factors on the Variability of Archaeal Communities in a Karst Wetland.

Archaea are ubiquitous and play an important role in elemental cycles in Earth's biosphere; but little is known about their diversity, distribution, abundance, and impact in karst environments. The present study investigated the effect of environmental factors on the variability of archaeal communities in the sediment of the Huixian karst wetland, the largest karst wetland in South China. Sediment cores were obtained from four sampling sites with different water depths and macrophyte inhabitants in both the winter of 2016 and the summer of 2018. The community analysis was based on PacBio sequencing and quantitative PCR of the archaeal 16S rRNA gene. The results showed that Euryarchaeota (57.4%) and Bathyarchaeota (38.7%) were dominant in all the samples. Methanogenic Methanosarcinales (25.1%) and Methanomicrobiales (13.7%), and methanotrophic archaea ANME-2d (9.0%) were the dominant Euryarchaeota; MCG-11 (16.5%), MCG-6 (9.1%), and MCG-5b (5.5%) were the dominant Bathyarchaeota. The community composition remained stable between summer and winter, and the vertical distributions of the archaeal phyla conformed to two patterns among the four sampling sites. In the winter samples, the archaeal 16S rRNA gene abundance was approximately 1.0E+10 copies/g of wet sediment and the Shannon index was 7.3±5, which were significantly higher than in the summer samples and in other karst environments. A correlation analysis showed that the moisture content and pH were the factors that mostly affected the archaeal communities. The prevalence of nitrate in the summer may be a key factor causing a significant decrease in archaeal abundance and diversity. Two features specific to karst environments, calcium-richness and weak alkalescence of the water supplies, may benefit the prevalence of bathyarchaeotal subgroups MCG-11, MCG-5b, and MCG-6. These results suggest that in karst wetlands, most of the archaea belong to clades that have significant roles in carbon turnover; their composition remains stable, but their abundance and diversity vary significantly from season to season.

Crop rotation reduces the frequency of anaerobic soil bacteria in Red Latosol of Brazil.

Crop diversity affects the processes of soil physical structuring and most likely provokes changes in the frequencies of soil microbial communities. The study was conducted for soil prokaryotic diversity sequencing 16S rDNA genes from a 25-year no-tillage experiment comprised of two crop systems: crop succession (Triticum aestivum-Glycine max) and rotation (Vicia sativa-Zea mays-Avena sativa-Glycine max-Triticum aestivum-Glycine max). The hypothesis was that a crop system with higher crop diversification (rotation) would affect the frequencies of prokaryotic taxa against a less diverse crop system (succession) altering the major soil functions guided by bacterial diversity. Soils in both crop systems were dominated by Proteobacteria (31%), Acidobacteria (23%), Actinobacteria (10%), and Gemmatimonadetes (7.2%), among other common copiotrophic soil bacteria. Crop systems did not affect the richness and diversity indexes of soil bacteria and soil archaea. However, the crop rotation system reduced only the frequencies of anaerobic metabolism bacteria Chloroacidobacteria, Holophagae, Spirochaetes, Euryarchaeota, and Crenarchaeota. It can be concluded that crop succession, a system that is poorer in root diversity over time, may have conditioned the soil to lower oxygen diffusion and built up ecological niches that suitable for anaerobic bacteria tolerating lower levels of oxygen. On the other hand, it appeared that crop rotation has restructured the soil over the years while enabling copiotrophic aerobic bacteria to dominate the soil ecosystem. The changes prompted by crop succession have implications for efficient soil organic matter decomposition, reduced greenhouse gas emissions, higher root activity, and overall soil productivity, which compromise to agriculture sustainability.

Corrigendum: High Representation of Archaea Across All Depths in Oxic and Low-pH Sediment Layers Underlying an Acidic Stream.

[This corrects the article DOI: 10.3389/fmicb.2020.576520.].

ADP-dependent glucose/glucosamine kinase from Thermococcus kodakarensis: cloning and characterization.

The genome sequence of Thermococcus kodakarensis contains an open reading frame, TK1110, annotated as ADP-dependent glucokinase. The encoding gene was expressed in Escherichia coli and the gene product, TK-GLK, was produced in soluble and active form. The recombinant enzyme was extremely thermostable. Thermostability was increased significantly in the presence of ammonium sulfate. ADP was the preferred co-factor for TK-GLK, which could be replaced with CDP but with a 60% activity. TK-GLK was a metal ion-dependent enzyme which exhibited glucokinase, glucosamine kinase and glucose 6-phosphatase activities. It catalyzed the phosphorylation of both glucose and glucosamine with nearly the same rate and affinity. The apparent Km values for glucose and glucosamine were 0.48 ± 0.03 and 0.47 ± 0.09 mM, respectively. The catalytic efficiency (kcat/Km) values against these two substrates were 6.2 × 105 ± 0.25 and 5.8 × 105 ± 0.75 M-1 s-1. The apparent Km value for dephosphorylation of glucose 6-phosphate was ~14-fold higher than that of glucose phosphorylation. Similarly, catalytic efficiency (kcat/Km) for phosphatase reaction was ~19-fold lower than that for the kinase reaction. To the best of our knowledge, this is the first report that describes the reversible nature of a euryarchaeal ADP-dependent glucokinase.

Euryarchaeal genomes are folded into SMC-dependent loops and domains, but lack transcription-mediated compartmentalization.

Hi-C has become a routine method for probing the 3D organization of genomes. However, when applied to prokaryotes and archaea, the current protocols are expensive and limited in their resolution. We develop a cost-effective Hi-C protocol to explore chromosome conformations of these two kingdoms at the gene or operon level. We first validate it on E. coli and V. cholera, generating sub-kilobase-resolution contact maps, and then apply it to the euryarchaeota H. volcanii, Hbt. salinarum, and T. kodakaraensis. With a resolution of up to 1 kb, we explore the diversity of chromosome folding in this phylum. In contrast to crenarchaeota, these euryarchaeota lack (active/inactive) compartment-like structures. Instead, their genomes are composed of self-interacting domains and chromatin loops. In H. volcanii, these structures are regulated by transcription and the archaeal structural maintenance of chromosomes (SMC) protein, further supporting the ubiquitous role of these processes in shaping the higher-order organization of genomes.

High Representation of Archaea Across All Depths in Oxic and Low-pH Sediment Layers Underlying an Acidic Stream.

Parys Mountain or Mynydd Parys (Isle of Anglesey, United Kingdom) is a mine-impacted environment, which accommodates a variety of acidophilic organisms. Our previous research of water and sediments from one of the surface acidic streams showed a high proportion of archaea in the total microbial community. To understand the spatial distribution of archaea, we sampled cores (0-20 cm) of sediment and conducted chemical analyses and taxonomic profiling of microbiomes using 16S rRNA gene amplicon sequencing in different core layers. The taxonomic affiliation of sequencing reads indicated that archaea represented between 6.2 and 54% of the microbial community at all sediment depths. Majority of archaea were associated with the order Thermoplasmatales, with the most abundant group of sequences being clustered closely with the phylotype B_DKE, followed by "E-plasma," "A-plasma," other yet uncultured Thermoplasmatales with Ferroplasma and Cuniculiplasma spp. represented in minor proportions. Thermoplasmatales were found at all depths and in the whole range of chemical conditions with their abundance correlating with sediment Fe, As, Cr, and Mn contents. The bacterial microbiome component was largely composed in all layers of sediment by members of the phyla Proteobacteria, Actinobacteria, Nitrospirae, Firmicutes, uncultured Chloroflexi (AD3 group), and Acidobacteria. This study has revealed a high abundance of Thermoplasmatales in acid mine drainage-affected sediment layers and pointed at these organisms being the main contributors to carbon, and probably to iron and sulfur cycles in this ecosystem.

Unveiling the Biodiversity of Hyperthermophilic Archaea in Jharia Coal Mines: Potential Threat to Methanogenesis?

AIM: To examine the biodiversity of archaeal sulfate reducers and methanogens present in the underground coal mines of Jharia using metagenomics and pyrosequencing. OBJECTIVES: 1) Bioinformatical analysis of the metagenomic data related to a taxonomic analysis obtained from the coal to investigate complete archaeal taxonomic features of the coal bed methane (CBM) microbiome. 2) Bioinformatical analysis of the metagenomic data related to a functional analysis obtained from the coal to investigate functional features relating to taxonomic diversity of the CBM microbiome. 3) The functional attributes have been examined specifically for ORFs related to sulfite reduction and methanogenesis.The taxonomic and functional biodiversity related to euryarchaeota will help in a better understanding of the obstacles associated with methane production imposed by the sulfate reducers. BACKGROUND: The microbial methanogenesis in the coal microbiome is a resultant of substrate utilization by primarily fermentative bacteria and methanogens. The present work reveals the biodiversity of archaeal sulfate reducers and methanogens present in the underground coal mines of Jharia using metagenomics and pyrosequencing. METHODOLOGY: Bioinformatical analysis for structural and functional attributes was accomplished using MG-RAST. The structural analysis was accomplished using RefSeq database, whereas the functional analysis was done via CoG database with a cut off value, a sequence percent identity, and sequence alignment length cut off of 1e-5, 60% and 45, respectively. RESULTS: Attained communities revealed the dominance of hyperthermophilic archaea Pyrococcus furiosus along with Thermococcus kodakarensis in the coal metagenome.The obtained results also suggest the presence of dissimilatory sulfite reductase and formylmethanofuran dehydrogenase, formylmethanofuran: tetrahydromethanopterin formyltransferase involved in sulfite reduction and methanogenesis, respectively, in the microbiome. CONCLUSION: This report is the first attempt to showcase the existence of specific euryarchaeal diversity and their related functional attributes from Jharia coal mines through high throughput sequencing. The study helps in developing a better understanding of the presence of indigenous microbes (archaea) and their functions in the coal microbiome, which can be utilized further to resolve the energy crisis.

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.

Evolutionary considerations of the oligosaccharyltransferase AglB and other aspects of N-glycosylation across Archaea.

Various biological markers in members of the TACK and Asgard archaeal super-phyla show Eukarya-like traits. These include the oligosaccharyltransferase, responsible for transferring glycans from the lipid carrier upon which they are assembled onto selected asparagine residues of target proteins during N-glycosylation. In Archaea, oligosaccharyltransferase activity is catalyzed by AglB. To gain deeper insight into AglB and N-glycosylation across archaeal phylogeny, bioinformatics approaches were employed to address variability in AglB sequence motifs involved in enzyme activity, construct a phylogenetic tree based on AglB sequences, search for archaeal homologues of non-catalytic subunits of the multimeric eukaryal oligosaccharyltransferase complex and predict the presence of aglB-based clusters of glycosylation-related genes in the Euryarchaeota and the DPANN, TACK and Asgard super-phyla. In addition, site-directed mutagenesis and mass spectrometry were employed to study the natural variability in the WWDXG motif central to oligosaccharyltransferase activity seen in archaeal AglB. The results clearly distinguish AglB from members of the DPANN super-phylum and the Euryarchaeota from the same enzyme in members of the TACK and Asgard super-phyla, which showed considerable similarity to its eukaryal homologue Stt3. The results thus support the evolutionary proximity of Eukarya and the TACK and Asgard archaea.

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

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

In situ absorbance measurements: a new means to study respiratory electron transfer in chemolithotrophic microorganisms.

Absorbance measurements on intact chemolithotrophic microorganisms that respire aerobically on soluble iron are described that used a novel integrating cavity absorption meter to eliminate the effects of light scattering on the experimental results. Steady state kinetic measurements on ferric iron production by intact cells revealed that the Michaelis Menten equation described the initial rates of product formation for at least 8 different chemolithotrophic microorganisms in 6 phyla distributed equally among the archaea and the Gram negative and Gram positive eubacteria. Cell-monitored turnover measurements during aerobic respiration on soluble iron by the same 12 intact microorganisms revealed six different patterns of iron-dependent absorbance changes, suggesting that there may be at least six different sets of prosthetic groups and biomolecules that can accomplish aerobic respiration on soluble iron. Detailed kinetic studies revealed that the 3-component iron respiratory chain of Acidithiobacillus ferrooxidans functioned as an ensemble with a single macroscopic rate constant when the iron-reduced proteins were oxidized in the presence of excess molecular oxygen. The principal member of this 3-component system was a cupredoxin called rusticyanin that was present in the periplasm of At. ferrooxidans at an approximate concentration of 350 mg/mL, an observation that provides new insights into the crowded environments in the periplasms of Gram negative eubacteria that conduct electrons across their periplasm. The ability to conduct direct spectrophotometric measurements under noninvasive physiological conditions represents a new and powerful approach to examine the rates and extents of biological events in situ without disrupting the complexity of the live cellular environment.

Biochemical characterization of a highly active ADP-dependent phosphofructokinase from Thermococcus kodakarensis.

The genome sequence of Thermococcus kodakarensis contains an open reading frame, TK0376, annotated as ADP-dependent phosphofructokinase belonging to pfkC family. The encoding gene was expressed in Escherichia coli and the gene product was characterized. The recombinant protein was produced in soluble and active form. Phosphofructokinase activity of TK0376 was metal-ion dependent and the highest activity (5090 µmol min-1 mg-1) was found in the presence of Co2+ followed by Mg2+ (3280 µmol min-1 mg-1) at 90°C and pH 7.5. TK0376 preferred ADP as phosphoryl donor, however, it could be replaced by ATP but with a 5-fold lower activity. It catalyzed the phosphorylation of fructose 6-phosphate and dephosphorylation of fructose 1,6-bisphosphate. In addition, it was able to phosphorylate glucose and nucleosides but with a much lower rate compared to that of fructose 6-phosphate. The apparent kcat and Km values against fructose 6-phosphate were 4238 s-1 and 0.74 mM, respectively. The rate of dephosphorylation of fructose 1,6-bisphosphate was 3-times lower at 50°C than the phosphorylation of fructose 6-phosphate. Similarly, the rate of phosphorylation of glucose was 450-fold lower than that of fructose 6-phosphate. Phosphofructokinase activity was not allosterically regulated, but it was slightly enhanced by phosphoenol pyruvate, and inhibited by ATP and AMP in a competitive manner.

Microbial Communities of Polymetallic Deposits' Acidic Ecosystems of Continental Climatic Zone With High Temperature Contrasts.

Acid mine drainage (AMD) systems are globally widespread and are an important source of metal pollution in riverine and coastal systems. Microbial AMD communities have been extensively studied for their ability to thrive under extremely acidic conditions and for their immense contribution to the dissolution of metal ores. However, little is known on microbial inhabitants of AMD systems subjected to extremely contrasting continental seasonal temperature patterns as opposed to maritime climate zones, experiencing much weaker annual temperature variations. Here, we investigated three types of AMD sites in Eastern Transbaikalia (Russia). In this region, all surface water bodies undergo a deep and long (up to 6 months) freezing, with seasonal temperatures varying between -33 and +24°C, which starkly contrasts the common well-studied AMD environments. We sampled acidic pit lake (Sherlovaya Gora site) located in the area of a polymetallic deposit, acidic drainage water from Bugdaya gold-molybdenum-tungsten deposit and Ulan-Bulak natural acidic spring. These systems showed the abundance of bacteria-derived reads mostly affiliated with Actinobacteria, Acidobacteria, Alpha- and Gammaproteobacteria, chloroplasts, Chloroflexi, Bacteroidetes, and Firmicutes. Furthermore, candidate taxa "Ca. Saccharibacteria" (previously known as TM7), "Ca. Parcubacteria" (OD1) and WPS-2 were represented in substantial quantities (10-20%). Heterotrophy and iron redox cycling can be considered as central processes of carbon and energy flow for majority of detected bacterial taxa. Archaea were detected in low numbers, with Terrestrial Miscellaneous Euryarchaeal Group (TMEG), to be most abundant (3%) in acidic spring Ulan-Bulak. Composition of these communities was found to be typical in comparison to other AMD sites; however, certain groups (as Ignavibacteriae) could be specifically associated with this area. This study provides insight into the microbial diversity patterns in acidic ecosystems formed in areas of polymetallic deposits in extreme continental climate zone with contrasting temperature parameters.

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 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.