Regulatory sequence-based discovery of anti-defense genes in archaeal viruses

In silico identification of viral anti-CRISPR proteins (Acrs) has relied largely on the guilt-by-association method using known Acrs or anti-CRISPR associated proteins (Acas) as the bait. However, the low number and limited spread of the characterized archaeal Acrs and Aca hinders our ability to identify Acrs using guilt-by-association. Here, based on the observation that the few characterized archaeal Acrs and Aca are transcribed immediately post viral infection, we hypothesize that these genes, and many other unidentified anti-defense genes (ADG), are under the control of conserved regulatory sequences including a strong promoter, which can be used to predict anti-defense genes in archaeal viruses. Using this consensus sequence based method, we identify 354 potential ADGs in 57 archaeal viruses and 6 metagenome-assembled genomes. Experimental validation identified a CRISPR subtype I-A inhibitor and the first virally encoded inhibitor of an archaeal toxin-antitoxin based immune system. We also identify regulatory proteins potentially akin to Acas that can facilitate further identification of ADGs combined with the guilt-by-association approach. These results demonstrate the potential of regulatory sequence analysis for extensive identification of ADGs in viruses of archaea and bacteria.

Mechanisms of extracellular electron transfer in anaerobic methanotrophic archaea

Anaerobic methanotrophic (ANME) archaea are environmentally important, uncultivated microorganisms that oxidize the potent greenhouse gas methane. During methane oxidation, ANME archaea engage in extracellular electron transfer (EET) with other microbes, metal oxides, and electrodes through unclear mechanisms. Here, we cultivate ANME-2d archaea ('Ca. Methanoperedens') in bioelectrochemical systems and observe strong methane-dependent current (91-93% of total current) associated with high enrichment of 'Ca. Methanoperedens' on the anode (up to 82% of the community), as determined by metagenomics and transmission electron microscopy. Electrochemical and metatranscriptomic analyses suggest that the EET mechanism is similar at various electrode potentials, with the possible involvement of an uncharacterized short-range electron transport protein complex and OmcZ nanowires.

Distribution and genomic variation of ammonia-oxidizing archaea in abyssal and hadal surface sediments

Ammonia-oxidizing archaea of the phylum Thaumarchaeota play a central role in the biogeochemical cycling of nitrogen in benthic sediments, at the interface between pelagic and subsurface ecosystems. However, our understanding of their niche separation and of the processes controlling their population structure in hadal and abyssal surface sediments is still limited. Here, we reconstructed 47 AOA metagenome-assembled genomes (MAGs) from surface sediments of the Atacama and Kermadec trench systems. They formed deep-sea-specific groups within the family Nitrosopumilaceae and were assigned to six amoA gene-based clades. MAGs from different clades had distinct distribution patterns along oxygen-ammonium counter gradients in surface sediments. At the species level, MAGs thus seemed to form different ecotypes and follow deterministic niche-based distributions. In contrast, intraspecific population structure, defined by patterns of Single Nucleotide Variants (SNV), seemed to reflect more complex contributions of both deterministic and stochastic processes. Firstly, the bathymetric range had a strong effect on population structure, with distinct populations in abyssal plains and hadal trenches. Then, hadal populations were clearly separated by trench system, suggesting a strong isolation-by-topography effect, whereas abyssal populations were rather controlled by sediment depth or geographic distances, depending on the clade considered. Interestingly, genetic variability between samples was lowest in sediment layers where the mean MAG coverage was highest, highlighting the importance of selective pressure linked with each AOA clade's ecological niche. Overall, our results show that deep-sea AOA genome distributions seem to follow both deterministic and stochastic processes, depending on the genomic variability scale considered.

Genomics of soil depth niche partitioning in the Thaumarchaeota family Gagatemarchaeaceae

Knowledge of deeply-rooted non-ammonia oxidising Thaumarchaeota lineages from terrestrial environments is scarce, despite their abundance in acidic soils. Here, 15 new deeply-rooted thaumarchaeotal genomes were assembled from acidic topsoils (0-15 cm) and subsoils (30-60 cm), corresponding to two genera of terrestrially prevalent Gagatemarchaeaceae (previously known as thaumarchaeotal Group I.1c) and to a novel genus of heterotrophic terrestrial Thaumarchaeota. Unlike previous predictions, metabolic annotations suggest Gagatemarchaeaceae perform aerobic respiration and use various organic carbon sources. Evolutionary divergence between topsoil and subsoil lineages happened early in Gagatemarchaeaceae history, with significant metabolic and genomic trait differences. Reconstruction of the evolutionary mechanisms showed that the genome expansion in topsoil Gagatemarchaeaceae resulted from extensive early lateral gene acquisition, followed by progressive gene duplication throughout evolutionary history. Ancestral trait reconstruction using the expanded genomic diversity also did not support the previous hypothesis of a thermophilic last common ancestor of the ammonia-oxidising archaea. Ultimately, this study provides a good model for studying mechanisms driving niche partitioning between spatially related ecosystems.

Stepwise pathway for early evolutionary assembly of dissimilatory sulfite and sulfate reduction

Microbial dissimilatory sulfur metabolism utilizing dissimilatory sulfite reductases (Dsr) influenced the biochemical sulfur cycle during Earth's history and the Dsr pathway is thought to be an ancient metabolic process. Here we performed comparative genomics, phylogenetic, and synteny analyses of several Dsr proteins involved in or associated with the Dsr pathway across over 195,000 prokaryotic metagenomes. The results point to an archaeal origin of the minimal DsrABCMK(N) protein set, having as primordial function sulfite reduction. The acquisition of additional Dsr proteins (DsrJOPT) increased the Dsr pathway complexity. Archaeoglobus would originally possess the archaeal-type Dsr pathway and the archaeal DsrAB proteins were replaced with the bacterial reductive-type version, possibly at the same time as the acquisition of the QmoABC and DsrD proteins. Further inventions of two Qmo complex types, which are more spread than previously thought, allowed microorganisms to use sulfate as electron acceptor. The ability to use the Dsr pathway for sulfur oxidation evolved at least twice, with Chlorobi and Proteobacteria being extant descendants of these two independent adaptations.

Mcr-dependent methanogenesis in Archaeoglobaceae enriched from a terrestrial hot spring

The preeminent source of biological methane on Earth is methyl coenzyme M reductase (Mcr)-dependent archaeal methanogenesis. A growing body of evidence suggests a diversity of archaea possess Mcr, although experimental validation of hypothesized methane metabolisms has been missing. Here, we provide evidence of a functional Mcr-based methanogenesis pathway in a novel member of the family Archaeoglobaceae, designated Methanoglobus nevadensis, which we enriched from a terrestrial hot spring on the polysaccharide xyloglucan. Our incubation assays demonstrate methane production that is highly sensitive to the Mcr inhibitor bromoethanesulfonate, stimulated by xyloglucan and xyloglucan-derived sugars, concomitant with the consumption of molecular hydrogen, and causing a deuterium fractionation in methane characteristic of hydrogenotrophic and methylotrophic methanogens. Combined with the recovery and analysis of a high-quality M. nevadensis metagenome-assembled genome encoding a divergent Mcr and diverse potential electron and carbon transfer pathways, our observations suggest methanogenesis in M. nevadensis occurs via Mcr and is fueled by the consumption of cross-fed byproducts of xyloglucan fermentation mediated by other community members. Phylogenetic analysis shows close affiliation of the M. nevadensis Mcr with those from Korarchaeota, Nezhaarchaeota, Verstraetearchaeota, and other Archaeoglobales that are divergent from well-characterized Mcr. We propose these archaea likely also use functional Mcr complexes to generate methane on the basis of our experimental validation in M. nevadensis. Thus, divergent Mcr-encoding archaea may be underestimated sources of biological methane in terrestrial and marine hydrothermal environments.

A widespread group of large plasmids in methanotrophic Methanoperedens archaea

Anaerobic methanotrophic (ANME) archaea obtain energy from the breakdown of methane, yet their extrachromosomal genetic elements are little understood. Here we describe large plasmids associated with ANME archaea of the Methanoperedens genus in enrichment cultures and other natural anoxic environments. By manual curation we show that two of the plasmids are large (155,605 bp and 191,912 bp), circular, and may replicate bidirectionally. The plasmids occur in the same copy number as the main chromosome, and plasmid genes are actively transcribed. One of the plasmids encodes three tRNAs, ribosomal protein uL16 and elongation factor eEF2; these genes appear to be missing in the host Methanoperedens genome, suggesting an obligate interdependence between plasmid and host. Our work opens the way for the development of genetic vectors to shed light on the physiology and biochemistry of Methanoperedens, and potentially genetically edit them to enhance growth and accelerate methane oxidation rates.

Anaerobic carboxydotrophy in sulfur-respiring haloarchaea from hypersaline lakes

Anaerobic carboxydotrophy is a widespread catabolic trait in bacteria, with two dominant pathways: hydrogenogenic and acetogenic. The marginal mode by direct oxidation to CO₂ using an external e-acceptor has only a few examples. Use of sulfidic sediments from two types of hypersaline lakes in anaerobic enrichments with CO as an e-donor and elemental sulfur as an e-acceptor led to isolation of two pure cultures of anaerobic carboxydotrophs belonging to two genera of sulfur-reducing haloarchaea: Halanaeroarchaeum sp. HSR-CO from salt lakes and Halalkaliarchaeum sp. AArc-CO from soda lakes. Anaerobic growth of extremely halophilic archaea with CO was obligatory depended on the presence of elemental sulfur as the electron acceptor and yeast extract as the carbon source. CO served as a direct electron donor and H₂ was not generated from CO when cells were incubated with or without sulfur. The genomes of the isolates encode a catalytic Ni,Fe-CODH subunit CooS (distantly related to bacterial homologs) and its Ni-incorporating chaperone CooC (related to methanogenic homologs) within a single genomic locus. Similar loci were also present in a genome of the type species of Halalkaliarchaeum closely related to AArc-CO, and the ability for anaerobic sulfur-dependent carboxydotrophy was confirmed for three different strains of this genus. Moreover, similar proteins are encoded in three of the four genomes of recently described carbohydrate-utilizing sulfur-reducing haloarchaea belonging to the genus Halapricum and in two yet undescribed haloarchaeal species. Overall, this work demonstrated for the first time the potential for anaerobic sulfur-dependent carboxydotrophy in extremely halophilic archaea.

Comparative genomic analysis of Methanimicrococcus blatticola provides insights into host adaptation in archaea and the evolution of methanogenesis

Other than the Methanobacteriales and Methanomassiliicoccales, the characteristics of archaea that inhabit the animal microbiome are largely unknown. Methanimicrococcus blatticola, a member of the Methanosarcinales, currently reunites two unique features within this order: it is a colonizer of the animal digestive tract and can only reduce methyl compounds with H2 for methanogenesis, a increasingly recognized metabolism in the archaea and whose origin remains debated. To understand the origin of these characteristics, we have carried out a large-scale comparative genomic analysis. We infer the loss of more than a thousand genes in M. blatticola, by far the largest genome reduction across all Methanosarcinales. These include numerous elements for sensing the environment and adapting to more stable gut conditions, as well as a significant remodeling of the cell surface components likely involved in host and gut microbiota interactions. Several of these modifications parallel those previously observed in phylogenetically distant archaea and bacteria from the animal microbiome, suggesting large-scale convergent mechanisms of adaptation to the gut. Strikingly, M. blatticola has lost almost all genes coding for the H4MPT methyl branch of the Wood-Ljungdahl pathway (to the exception of mer), a phenomenon never reported before in any member of Class I or Class II methanogens. The loss of this pathway illustrates one of the evolutionary processes that may have led to the emergence of methyl-reducing hydrogenotrophic methanogens, possibly linked to the colonization of organic-rich environments (including the animal gut) where both methyl compounds and hydrogen are abundant.

Comparative genomic analysis reveals metabolic flexibility of Woesearchaeota

The archaeal phylum Woesearchaeota, within the DPANN superphylum, includes phylogenetically diverse microorganisms that inhabit various environments. Their biology is poorly understood due to the lack of cultured isolates. Here, we analyze datasets of Woesearchaeota 16S rRNA gene sequences and metagenome-assembled genomes to infer global distribution patterns, ecological preferences and metabolic capabilities. Phylogenomic analyses indicate that the phylum can be classified into ten subgroups, termed A-J. While a symbiotic lifestyle is predicted for most, some members of subgroup J might be host-independent. The genomes of several Woesearchaeota, including subgroup J, encode putative [FeFe] hydrogenases (known to be important for fermentation in other organisms), suggesting that these archaea might be anaerobic fermentative heterotrophs.

Brockarchaeota, a novel archaeal phylum with unique and versatile carbon cycling pathways

Geothermal environments, such as hot springs and hydrothermal vents, are hotspots for carbon cycling and contain many poorly described microbial taxa. Here, we reconstructed 15 archaeal metagenome-assembled genomes (MAGs) from terrestrial hot spring sediments in China and deep-sea hydrothermal vent sediments in Guaymas Basin, Gulf of California. Phylogenetic analyses of these MAGs indicate that they form a distinct group within the TACK superphylum, and thus we propose their classification as a new phylum, 'Brockarchaeota', named after Thomas Brock for his seminal research in hot springs. Based on the MAG sequence information, we infer that some Brockarchaeota are uniquely capable of mediating non-methanogenic anaerobic methylotrophy, via the tetrahydrofolate methyl branch of the Wood-Ljungdahl pathway and reductive glycine pathway. The hydrothermal vent genotypes appear to be obligate fermenters of plant-derived polysaccharides that rely mostly on substrate-level phosphorylation, as they seem to lack most respiratory complexes. In contrast, hot spring lineages have alternate pathways to increase their ATP yield, including anaerobic methylotrophy of methanol and trimethylamine, and potentially use geothermally derived mercury, arsenic, or hydrogen. Their broad distribution and their apparent anaerobic metabolic versatility indicate that Brockarchaeota may occupy previously overlooked roles in anaerobic carbon cycling.

Genome wide transcriptomic analysis of the soil ammonia oxidizing archaeon Nitrososphaera viennensis upon exposure to copper limitation

Ammonia-oxidizing archaea (AOA) are widespread in nature and are involved in nitrification, an essential process in the global nitrogen cycle. The enzymes for ammonia oxidation and electron transport rely heavily on copper (Cu), which can be limited in nature. In this study the model soil archaeon Nitrososphaera viennensis was investigated via transcriptomic analysis to gain insight regarding possible Cu uptake mechanisms and compensation strategies when Cu becomes limiting. Upon Cu limitation, N. viennensis exhibited impaired nitrite production and thus growth, which was paralleled by downregulation of ammonia oxidation, electron transport, carbon fixation, nucleotide, and lipid biosynthesis pathway genes. Under Cu-limitation, 1547 out of 3180 detected genes were differentially expressed, with 784 genes upregulated and 763 downregulated. The most highly upregulated genes encoded proteins with a possible role in Cu binding and uptake, such as the Cu chelator and transporter CopC/D, disulfide bond oxidoreductase D (dsbD), and multicopper oxidases. While this response differs from the marine strain Nitrosopumilus maritimus, conserved sequence motifs in some of the Cu-responsive genes suggest conserved transcriptional regulation in terrestrial AOA. This study provides possible gene regulation and energy conservation mechanisms linked to Cu bioavailability and presents the first model for Cu uptake by a soil AOA.

Gene duplication drives genome expansion in a major lineage of Thaumarchaeota

Ammonia-oxidising archaea of the phylum Thaumarchaeota are important organisms in the nitrogen cycle, but the mechanisms driving their radiation into diverse ecosystems remain underexplored. Here, existing thaumarchaeotal genomes are complemented with 12 genomes belonging to the previously under-sampled Nitrososphaerales to investigate the impact of lateral gene transfer (LGT), gene duplication and loss across thaumarchaeotal evolution. We reveal a major role for gene duplication in driving genome expansion subsequent to early LGT. In particular, two large LGT events are identified into Nitrososphaerales and the fate of these gene families is highly lineage-specific, being lost in some descendant lineages, but undergoing extensive duplication in others, suggesting niche-specific roles. Notably, some genes involved in carbohydrate transport or coenzyme metabolism were duplicated, likely facilitating niche specialisation in soils and sediments. Overall, our results suggest that LGT followed by gene duplication drives Nitrososphaerales evolution, highlighting a previously under-appreciated mechanism of genome expansion in archaea.

Hikarchaeia demonstrate an intermediate stage in the methanogen-to-halophile transition

Halobacteria (henceforth: Haloarchaea) are predominantly aerobic halophiles that are thought to have evolved from anaerobic methanogens. This remarkable transformation most likely involved an extensive influx of bacterial genes. Whether it entailed a single massive transfer event or a gradual stream of transfers remains a matter of debate. To address this, genomes that descend from methanogen-to-halophile intermediates are necessary. Here, we present five such near-complete genomes of Marine Group IV archaea (Hikarchaeia), the closest known relatives of Haloarchaea. Their inclusion in gene tree-aware ancestral reconstructions reveals an intermediate stage that had already lost a large number of genes, including nearly all of those involved in methanogenesis and the Wood-Ljungdahl pathway. In contrast, the last Haloarchaea common ancestor gained a large number of genes and expanded its aerobic respiration and salt/UV resistance gene repertoire. Our results suggest that complex and gradual patterns of gain and loss shaped the methanogen-to-halophile transition.

Culture of Methanogenic Archaea from Human Colostrum and Milk

Archaeal sequences have been detected in human colostrum and milk, but no studies have determined whether living archaea are present in either of these fluids. Methanogenic archaea are neglected since they are not detected by usual molecular and culture methods. By using improved DNA detection protocols and microbial culture techniques associated with antioxidants previously developed in our center, we investigated the presence of methanogenic archaea using culture and specific Methanobrevibacter smithii and Methanobrevibacter oralis real-time PCR in human colostrum and milk. M. smithii was isolated from 3 colostrum and 5 milk (day 10) samples. M. oralis was isolated from 1 milk sample. For 2 strains, the genome was sequenced, and the rhizome was similar to that of strains previously isolated from the human mouth and gut. M. smithii was detected in the colostrum or milk of 5/13 (38%) and 37/127 (29%) mothers by culture and qPCR, respectively. The different distribution of maternal body mass index according to the detection of M. smithii suggested an association with maternal metabolic phenotype. M. oralis was not detected by molecular methods. Our results suggest that breastfeeding may contribute to the vertical transmission of these microorganisms and may be essential to seed the infant's microbiota with these neglected critical commensals from the first hour of life.

Divergent methyl-coenzyme M reductase genes in a deep-subseafloor Archaeoglobi

The methyl-coenzyme M reductase (MCR) complex is a key enzyme in archaeal methane generation and has recently been proposed to also be involved in the oxidation of short-chain hydrocarbons including methane, butane, and potentially propane. The number of archaeal clades encoding the MCR continues to grow, suggesting that this complex was inherited from an ancient ancestor, or has undergone extensive horizontal gene transfer. Expanding the representation of MCR-encoding lineages through metagenomic approaches will help resolve the evolutionary history of this complex. Here, a near-complete Archaeoglobi metagenome-assembled genome (MAG; Ca. Polytropus marinifundus gen. nov. sp. nov.) was recovered from the deep subseafloor along the Juan de Fuca Ridge flank that encodes two divergent McrABG operons similar to those found in Ca. Bathyarchaeota and Ca. Syntrophoarchaeum MAGs. Ca. P. marinifundus is basal to members of the class Archaeoglobi, and encodes the genes for β-oxidation, potentially allowing an alkanotrophic metabolism similar to that proposed for Ca. Syntrophoarchaeum. Ca. P. marinifundus also encodes a respiratory electron transport chain that can potentially utilize nitrate, iron, and sulfur compounds as electron acceptors. Phylogenetic analysis suggests that the Ca. P. marinifundus MCR operons were horizontally transferred, changing our understanding of the evolution and distribution of this complex in the Archaea.

Novel Large Sulfur Bacteria in the Metagenomes of Groundwater-Fed Chemosynthetic Microbial Mats in the Lake Huron Basin

Little is known about large sulfur bacteria (LSB) that inhabit sulfidic groundwater seeps in large lakes. To examine how geochemically relevant microbial metabolisms are partitioned among community members, we conducted metagenomic analysis of a chemosynthetic microbial mat in the Isolated Sinkhole, which is in a deep, aphotic environment of Lake Huron. For comparison, we also analyzed a white mat in an artesian fountain that is fed by groundwater similar to Isolated Sinkhole, but that sits in shallow water and is exposed to sunlight. De novo assembly and binning of metagenomic data from these two communities yielded near complete genomes and revealed representatives of two families of LSB. The Isolated Sinkhole community was dominated by novel members of the Beggiatoaceae that are phylogenetically intermediate between known freshwater and marine groups. Several of these Beggiatoaceae had 16S rRNA genes that contained introns previously observed only in marine taxa. The Alpena fountain was dominated by populations closely related to Thiothrix lacustris and an SM1 euryarchaeon known to live symbiotically with Thiothrix spp. The SM1 genomic bin contained evidence of H₂-based lithoautotrophy. Genomic bins of both the Thiothrix and Beggiatoaceae contained genes for sulfur oxidation via the rDsr pathway, H₂ oxidation via Ni-Fe hydrogenases, and the use of O₂ and nitrate as electron acceptors. Mats at both sites also contained Deltaproteobacteria with genes for dissimilatory sulfate reduction (sat, apr, and dsr) and hydrogen oxidation (Ni-Fe hydrogenases). Overall, the microbial mats at the two sites held low-diversity microbial communities, displayed evidence of coupled sulfur cycling, and did not differ largely in their metabolic potentials, despite the environmental differences. These results show that groundwater-fed communities in an artesian fountain and in submerged sinkholes of Lake Huron are a rich source of novel LSB, associated heterotrophic and sulfate-reducing bacteria, and archaea.

Extensive Inter-Domain Lateral Gene Transfer in the Evolution of the Human Commensal Methanosphaera stadtmanae

Methanosphaera stadtmanae is a commensal methanogenic archaeon found in the human gut. As most of its niche-neighbors are bacteria, it is expected that lateral gene transfer (LGT) from bacteria might have contributed to the evolutionary history of this organism. We performed a phylogenomic survey of putative LGT events in M. stadtmanae, using a phylogenetic pipeline. Our analysis indicates that a substantial fraction of the proteins of M. stadtmanae are inferred to have been involved in inter-domain LGT. Laterally acquired genes have had a large contribution to surface functions, by providing novel glycosyltransferase functions. In addition, several ABC transporters seem to be of bacterial origin, including the molybdate transporter. Thus, bacterial genes contributed to the adaptation of M. stadtmanae to a host-dependent lifestyle by allowing a larger variation in surface structures and increasing transport efficiency in the gut niche which is diverse and competitive.