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

Expanded phylogeny of extremely halophilic archaea shows multiple independent adaptations to hypersaline environments

Extremely halophilic archaea (Haloarchaea, Nanohaloarchaeota, Methanonatronarchaeia and Halarchaeoplasmatales) thrive in saturating salt concentrations where they must maintain osmotic equilibrium with their environment. The evolutionary history of adaptations enabling salt tolerance remains poorly understood, in particular because the phylogeny of several lineages is conflicting. Here we present a resolved phylogeny of extremely halophilic archaea obtained using improved taxon sampling and state-of-the-art phylogenetic approaches designed to cope with the strong compositional biases of their proteomes. We describe two uncultured lineages, Afararchaeaceae and Asbonarchaeaceae, which break the long branches at the base of Haloarchaea and Nanohaloarchaeota, respectively. We obtained 13 metagenome-assembled genomes (MAGs) of these archaea from metagenomes of hypersaline aquatic systems of the Danakil Depression (Ethiopia). Our phylogenomic analyses including these taxa show that at least four independent adaptations to extreme halophily occurred during archaeal evolution. Gene-tree/species-tree reconciliation suggests that gene duplication and horizontal gene transfer played an important role in this process, for example, by spreading key genes (such as those encoding potassium transporters) across extremely halophilic lineages.

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

IMPORTANCE

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

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

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

Pan-genome analysis of the Burkholderia gladioli PV. Cocovenenans reveal the extent of variation in the toxigenic gene cluster

Burkholderia gladioli has been reported as the pathogen responsible for cases of foodborne illness in many countries. The poisonous bongkrekic acid (BA) produced by B. gladioli was linked to a gene cluster absent in non-pathogenic strains. The whole genome sequence of eight bacteria strains, which were screened from the collected 175 raw food and environmental samples, were assembled and analyzed to detect a significant association of 19 protein-coding genes with the pathogenic status. Except for the common BA synthesis-related gene, several other genes, including the toxin-antitoxin genes, were also absent in the non-pathogenic strains. The bacteria strains with the BA gene cluster were found to form a single cluster in the analysis of all B. gladioli genome assemblies for the variants in the gene cluster. Divergence of this cluster was detected in the analysis for both the flanking sequences and those of the whole genome level, which indicates its complex origin. Genome recombination was found to cause a precise sequence deletion in the gene cluster region, which was found to be predominant in the non-pathogenic strains indicating the possible effect of horizontal gene transfer. Our study provided new information and resources for understanding the evolution and divergence of the B. gladioli species.

Metagenome-based metabolic modelling predicts unique microbial interactions in deep-sea hydrothermal plume microbiomes

Deep-sea hydrothermal vents are abundant on the ocean floor and play important roles in ocean biogeochemistry. In vent ecosystems such as hydrothermal plumes, microorganisms rely on reduced chemicals and gases in hydrothermal fluids to fuel primary production and form diverse and complex microbial communities. However, microbial interactions that drive these complex microbiomes remain poorly understood. Here, we use microbiomes from the Guaymas Basin hydrothermal system in the Pacific Ocean to shed more light on the key species in these communities and their interactions. We built metabolic models from metagenomically assembled genomes (MAGs) and infer possible metabolic exchanges and horizontal gene transfer (HGT) events within the community. We highlight possible archaea-archaea and archaea-bacteria interactions and their contributions to the robustness of the community. Cellobiose, D-Mannose 1-phosphate, O2, CO2, and H2S were among the most exchanged metabolites. These interactions enhanced the metabolic capabilities of the community by exchange of metabolites that cannot be produced by any other community member. Archaea from the DPANN group stood out as key microbes, benefiting significantly as acceptors in the community. Overall, our study provides key insights into the microbial interactions that drive community structure and organisation in complex hydrothermal plume microbiomes.

Bacterial origins of thymidylate metabolism in Asgard archaea and Eukarya

Asgard archaea include the closest known archaeal relatives of eukaryotes. Here, we investigate the evolution and function of Asgard thymidylate synthases and other folate-dependent enzymes required for the biosynthesis of DNA, RNA, amino acids and vitamins, as well as syntrophic amino acid utilization. Phylogenies of Asgard folate-dependent enzymes are consistent with their horizontal transmission from various bacterial groups. We experimentally validate the functionality of thymidylate synthase ThyX of the cultured 'Candidatus Prometheoarchaeum syntrophicum'. The enzyme efficiently uses bacterial-like folates and is inhibited by mycobacterial ThyX inhibitors, even though the majority of experimentally tested archaea are known to use carbon carriers distinct from bacterial folates. Our phylogenetic analyses suggest that the eukaryotic thymidylate synthase, required for de novo DNA synthesis, is not closely related to archaeal enzymes and might have been transferred from bacteria to protoeukaryotes during eukaryogenesis. Altogether, our study suggests that the capacity of eukaryotic cells to duplicate their genetic material is a sum of archaeal (replisome) and bacterial (thymidylate synthase) characteristics. We also propose that recent prevalent lateral gene transfer from bacteria has markedly shaped the metabolism of Asgard archaea.

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

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

Recovery of Lutacidiplasmatales archaeal order genomes suggests convergent evolution in Thermoplasmatota

The Terrestrial Miscellaneous Euryarchaeota Group has been identified in various environments, and the single genome investigated thus far suggests that these archaea are anaerobic sulfite reducers. We assemble 35 new genomes from this group that, based on genome analysis, appear to possess aerobic and facultative anaerobic lifestyles and may oxidise rather than reduce sulfite. We propose naming this order (representing 16 genera) "Lutacidiplasmatales" due to their occurrence in various acidic environments and placement within the phylum Thermoplasmatota. Phylum-level analysis reveals that Thermoplasmatota evolution had been punctuated by several periods of high levels of novel gene family acquisition. Several essential metabolisms, such as aerobic respiration and acid tolerance, were likely acquired independently by divergent lineages through convergent evolution rather than inherited from a common ancestor. Ultimately, this study describes the terrestrially prevalent Lutacidiciplasmatales and highlights convergent evolution as an important driving force in the evolution of archaeal lineages.

mOTUpan: a robust Bayesian approach to leverage metagenome-assembled genomes for core-genome estimation

Recent advances in sequencing and bioinformatics have expanded the tree of life by providing genomes for uncultured environmentally relevant clades, either through metagenome-assembled genomes or through single-cell genomes. While this expanded diversity can provide novel insights into microbial population structure, most tools available for core-genome estimation are sensitive to genome completeness. Consequently, a major portion of the huge phylogenetic diversity uncovered by environmental genomic approaches remains excluded from such analyses. We present mOTUpan, a novel iterative Bayesian method for computing the core genome for sets of genomes of highly diverse completeness range. The likelihood for each gene cluster to belong to core or accessory genome is estimated by computing the probability of its presence/absence pattern in the target genome set. The core-genome prediction is computationally efficient and can be scaled up to thousands of genomes. It has shown comparable estimates to state-of-the-art tools Roary and PPanGGOLiN for high-quality genomes and is capable of using genomes at lower completeness thresholds. mOTUpan wraps a bootstrapping procedure to estimate the quality of a specific core-genome prediction, as the accuracy of each run will depend on the specific completeness distribution and the number of genomes in the dataset under scrutiny. mOTUpan is implemented in the mOTUlizer software package, and available at github.com/moritzbuck/mOTUlizer, under GPL 3.0 license.

Horizontal Gene Transfer in Archaea-From Mechanisms to Genome Evolution

Archaea remains the least-studied and least-characterized domain of life despite its significance not just to the ecology of our planet but also to the evolution of eukaryotes. It is therefore unsurprising that research into horizontal gene transfer (HGT) in archaea has lagged behind that of bacteria. Indeed, several archaeal lineages may owe their very existence to large-scale HGT events, and thus understanding both the molecular mechanisms and the evolutionary impact of HGT in archaea is highly important. Furthermore, some mechanisms of gene exchange, such as plasmids that transmit themselves via membrane vesicles and the formation of cytoplasmic bridges that allows transfer of both chromosomal and plasmid DNA, may be archaea specific. This review summarizes what we know about HGT in archaea, and the barriers that restrict it, highlighting exciting recent discoveries and pointing out opportunities for future research. Expected final online publication date for the Annual Review of Microbiology, Volume 76 is September 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

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.

Genome- and community-level interaction insights into the ecological role of archaea in rare earth element mine drainage in South China

Microbial communities play crucial roles in mine drainage generation and remediation. Despite the wide distribution of archaea in the mine ecosystem, their diversity and ecological roles remain less understood than bacteria. Here, we retrieved 56 archaeal metagenome-assembled genomes from a river impacted by rare earth element (REE) mining activities in South China. Genomic analysis showed that archaea represented four distinct lineages, including phyla of Thaumarchaeota, Micrarchaeota, Nanoarchaeota and Thermoplasmata. These archaea represented a considerable fraction (up to 40%) of the total prokaryote community, which might contribute to nitrogen and sulfur cycling in the REE mine drainage. Reconstructed metabolic potential among diverse archaea taxa revealed that archaea were involved in the network of ammonia oxidation, denitrification, sulfate redox reaction, and required substrates supplied by other community members. As the dominant driver of ammonia oxidation, Thaumarchaeota might provide substrates to support the survival of two nano-sized archaea belonging to Micrarchaeota and Nanoarchaeota. Despite the absence of biosynthesis pathways for amino acids and nucleotides, the potential capacity for nitrite reduction (nirD) was observed in Micrarchaeota, indicating that these nano-sized archaea encompassed diverse metabolisms. Moreover, Thermoplasmata, as keystone taxa in community, might be the main genetic donor for the other three archaeal phyla, transferring many environmental resistance related genes (e.g., V/A-type ATPase and Vitamin B12-transporting ATPase). The genetic interactions within archaeal community through horizontal gene transfer might be the key to the formation of archaeal resistance and functional partitioning. This study provides putative metabolic and genetic insights into the diverse archaea taxa from community-level perspectives, and highlights the ecological roles of archaea in REE contaminated aquatic environment.

Reconciling Asgardarchaeota Phylogenetic Proximity to Eukaryotes and Planctomycetes Cellular Features in the Evolution of Life

The relationship between the three domains of life—Archaea, Bacteria, and Eukarya—is one of Biology’s greatest mysteries. Current favored models imply two ancestral domains, Bacteria and Archaea, with eukaryotes originating within Archaea. This type of models has been supported by the recent description of the Asgardarchaeota, the closest prokaryotic relatives of eukaryotes. However, there are many problems associated with any scenarios implying that eukaryotes originated from within the Archaea, including genome mosaicism, phylogenies, the cellular organization of the Archaea, and their ancestral character. By contrast, all models of eukaryogenesis fail to consider two relevant discoveries: the detection of membrane coat proteins, and of phagocytosis-related processes in Planctomycetes, which are among the bacteria with the most developed endomembrane system.

Consideration of these often overlooked features and others found in Planctomycetes and related bacteria suggest an evolutionary model based on a single ancestral domain. In this model, the proximity of Asgard and eukaryotes is not rejected but instead, Asgard are considered as diverging away from a common ancestor instead of on the way toward the eukaryotic ancestor. This model based on a single ancestral domain solves most of the ambiguities associated with the ones based on two ancestral domains. The single-domain model is better suited to explain the origin and evolution of all three domains of life, blurring the distinctions between them. Support for this model as well as the opportunities that it presents not only for reinterpreting previous results, but also for planning future experiments, are explored.

(S)-3-Hydroxybutyryl-CoA Dehydrogenase From the Autotrophic 3-Hydroxypropionate/4-Hydroxybutyrate Cycle in Nitrosopumilus maritimus

Ammonia-oxidizing archaea of the phylum Thaumarchaeota are among the most abundant organisms that exert primary control of oceanic and soil nitrification and are responsible for a large part of dark ocean primary production. They assimilate inorganic carbon via an energetically efficient version of the 3-hydroxypropionate/4-hydroxybutyrate cycle. In this cycle, acetyl-CoA is carboxylated to succinyl-CoA, which is then converted to two acetyl-CoA molecules with 4-hydroxybutyrate as the key intermediate. This conversion includes the (S)-3-hydroxybutyryl-CoA dehydrogenase reaction. Here, we heterologously produced the protein Nmar_1028 catalyzing this reaction in thaumarchaeon Nitrosopumilus maritimus, characterized it biochemically and performed its phylogenetic analysis. This NAD-dependent dehydrogenase is highly active with its substrate, (S)-3-hydroxybutyryl-CoA, and its low K_m value suggests that the protein is adapted to the functioning in the 3-hydroxypropionate/4-hydroxybutyrate cycle. Nmar_1028 is homologous to the dehydrogenase domain of crotonyl-CoA hydratase/(S)-3-hydroxybutyryl-CoA dehydrogenase that is present in many Archaea. Apparently, the loss of the dehydratase domain of the fusion protein in the course of evolution was accompanied by lateral gene transfer of 3-hydroxypropionyl-CoA dehydratase/crotonyl-CoA hydratase from Bacteria. Although (S)-3-hydroxybutyryl-CoA dehydrogenase studied here is neither unique nor characteristic for the HP/HB cycle, Nmar_1028 appears to be the only (S)-3-hydroxybutyryl-CoA dehydrogenase in N. maritimus and is thus essential for the functioning of the 3-hydroxypropionate/4-hydroxybutyrate cycle and for the biology of this important marine archaeon.

Rhizomal Reclassification of Living Organisms

Living organisms interact with each other during their lifetime, leading to genomes rearrangement and sequences transfer. These well-known phenomena give these organisms mosaic genomes, which challenge their classification. Moreover, many findings occurred between the IXXth and XXIst century, especially the discovery of giant viruses and candidate phyla radiation (CPR). Here, we tried to provide an updated classification, which integrates 216 representative genomes of the current described organisms. The reclassification was expressed through a genetic network based on the total genomic content, not on a single gene to represent the tree of life. This rhizomal exploration represents, more accurately, the evolutionary relationships among the studied species. Our analyses show a separated branch named fifth TRUC (Things Resisting Uncompleted Classifications). This taxon groups CPRs together, independently from Bacteria, Archaea (which regrouped also Nanoarchaeota and Asgard members), Eukarya, and the giant viruses (recognized recently as fourth TRUC). Finally, the broadening of analysis methods will lead to the discovery of new organisms, which justify the importance of updating the classification at every opportunity. In this perspective, our pragmatic representation could be adjusted along with the progress of evolutionary studies.

Comparative Genomic Analysis of the Foodborne Pathogen Burkholderia gladioli pv. cocovenenans Harboring a Bongkrekic Acid Biosynthesis Gene Cluster

The environmental bacterium Burkholderia gladioli pv. cocovenenans (B. cocovenenans) has been linked to fatal food poisoning cases in Asia and Africa. Bongkrekic acid (BA), a mitochondrial toxin produced by B. cocovenenans, is thought to be responsible for these outbreaks. While there are over 80 species in the Burkholderia genus, B. cocovenenans is the only pathovar capable of producing BA and causing human death. However, the genomic features of B. gladioli and the evolution of the BA biosynthesis gene cluster, bon, in B. cocovenenans remain elusive. In this study, 239 whole genome sequences (WGSs) of B. gladioli, isolated from 12 countries collected over 100 years, were used to analyze the intra-species genomic diversity and phylogenetic relationships of B. gladioli and to explore the origin and evolution of the bon gene cluster. Our results showed that the genome-wide average nucleotide identity (ANI) values were above 97.29% for pairs of B. gladioli genomes. Thirty-six of the 239 (15.06%) B. gladioli genomes, isolated from corn, rice, fruits, soil, and patients from Asia, Europe, North America, and South America, contained the bon gene cluster and formed three clades within the phylogenetic tree. Pan- and core-genome analysis suggested that the BA biosynthesis genes were recently acquired. Comparative genome analysis of the bon gene cluster showed that complex recombination events contributed to this toxin biosynthesis gene cluster's evolution and formation. This study suggests that a better understanding of the genomic diversity and evolution of this lethal foodborne pathovar will potentially contribute to B. cocovenenans food poisoning outbreak prevention.

The evolution pathway of ammonia-oxidizing archaea shaped by major geological events

Primordial nitrification processes have been studied extensively using geochemical approaches, but the biological origination of nitrification remains unclear. Ammonia-oxidizing archaea (AOA) are widely distributed nitrifiers and implement the rate-limiting step in nitrification. They are hypothesized to have been important players in the global nitrogen cycle in Earth’s early history. We performed systematic phylogenomic and marker gene analyses to elucidate the diversification timeline of AOA evolution. Our results suggested that the AOA ancestor experienced terrestrial geothermal environments at ∼1,165 Ma (1,928–880 Ma), and gradually evolved into mesophilic soil at ∼652 Ma (767–554 Ma) before diversifying into marine settings at ∼509 Ma (629–412 Ma) and later into shallow and deep oceans, respectively. Corroborated by geochemical evidence and modeling, the timing of key diversification nodes can be linked to the global magmatism and glaciation associated with the assembly and breakup of the supercontinent Rodinia, and the later oxygenation of the deep ocean. Results of this integrated study shed light on the geological forces that may have shaped the evolutionary pathways of the AOA, which played an important role in the ancient global nitrogen cycle.

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.

Integrating pan-genome with metagenome for microbial community profiling

Advances in sequencing technology have led to the increased availability of genomes and metagenomes, which has greatly facilitated microbial pan-genome and metagenome analysis in the community. In line with this trend, studies on microbial genomes and phenotypes have gradually shifted from individuals to environmental communities.

Pan-genomics and metagenomics are powerful strategies for in-depth profiling study of microbial communities. Pan-genomics focuses on genetic diversity, dynamics, and phylogeny at the multi-genome level, while metagenomics profiles the distribution and function of culture-free microbial communities in special environments. Combining pan-genome and metagenome analysis can reveal the microbial complicated connections from an individual complete genome to a mixture of genomes, thereby extending the catalog of traditional individual genomic profile to community microbial profile. Therefore, the combination of pan-genome and metagenome approaches has become a promising method to track the sources of various microbes and decipher the population-level evolution and ecosystem functions.

This review summarized the pan-genome and metagenome approaches, the combined strategies of pan-genome and metagenome, and applications of these combined strategies in studies of microbial dynamics, evolution, and function in communities. We discussed emerging strategies for the study of microbial communities that integrate information in both pan-genome and metagenome. We emphasized studies in which the integrating pan-genome with metagenome approach improved the understanding of models of microbial community profiles, both structural and functional. Finally, we illustrated future perspectives of microbial community profile: more advanced analytical techniques, including big-data based artificial intelligence, will lead to an even better understanding of the patterns of microbial communities.

Ancestral Reconstructions Decipher Major Adaptations of Ammonia-Oxidizing Archaea upon Radiation into Moderate Terrestrial and Marine Environments

Unlike all other archaeal lineages, ammonia-oxidizing archaea (AOA) of the phylum Thaumarchaeota are widespread and abundant in all moderate and oxic environments on Earth. The evolutionary adaptations that led to such unprecedented ecological success of a microbial clade characterized by highly conserved energy and carbon metabolisms have, however, remained underexplored. Here, we reconstructed the genomic content and growth temperature of the ancestor of all AOA, as well as the ancestors of the marine and soil lineages, based on 39 available complete or nearly complete genomes of AOA. Our evolutionary scenario depicts an extremely thermophilic, autotrophic, aerobic ancestor from which three independent lineages of a marine and two terrestrial groups radiated into moderate environments. Their emergence was paralleled by (i) a continuous acquisition of an extensive collection of stress tolerance genes mostly involved in redox maintenance and oxygen detoxification, (ii) an expansion of regulatory capacities in transcription and central metabolic functions, and (iii) an extended repertoire of cell appendages and modifications related to adherence and interactions with the environment. Our analysis provides insights into the evolutionary transitions and key processes that enabled the conquest of the diverse environments in which contemporary AOA are found.

Lysogenic host-virus interactions in SAR11 marine bacteria

Host-virus interactions structure microbial communities, drive biogeochemical cycles and enhance genetic diversity in nature1,2. Hypotheses proposed to explain the range of interactions that mediate these processes often invoke lysogeny3-6, a latent infection strategy used by temperate bacterial viruses to replicate in host cells until an induction event triggers the production and lytic release of free viruses. Most cultured bacteria harbour temperate viruses in their genomes (prophage)7. The absence of prophages in cultures of the dominant lineages of marine bacteria has contributed to an ongoing debate over the ecological significance of lysogeny and other viral life strategies in nature6,8-15. Here, we report the discovery of prophages in cultured SAR11, the ocean's most abundant clade of heterotrophic bacteria16,17. We show the concurrent production of cells and viruses, with enhanced virus production under carbon-limiting growth conditions. Evidence that related prophages are broadly distributed in the oceans suggests that similar interactions have contributed to the evolutionary success of SAR11 in nutrient-limited systems.

Origin and Evolution of Polycyclic Triterpene Synthesis

Polycyclic triterpenes are members of the terpene family produced by the cyclization of squalene. The most representative polycyclic triterpenes are hopanoids and sterols, the former are mostly found in bacteria, whereas the latter are largely limited to eukaryotes, albeit with a growing number of bacterial exceptions. Given their important role and omnipresence in most eukaryotes, contrasting with their scant representation in bacteria, sterol biosynthesis was long thought to be a eukaryotic innovation. Thus, their presence in some bacteria was deemed to be the result of lateral gene transfer from eukaryotes. Elucidating the origin and evolution of the polycyclic triterpene synthetic pathways is important to understand the role of these compounds in eukaryogenesis and their geobiological value as biomarkers in fossil records. Here, we have revisited the phylogenies of the main enzymes involved in triterpene synthesis, performing gene neighborhood analysis and phylogenetic profiling. Squalene can be biosynthesized by two different pathways containing the HpnCDE or Sqs proteins. Our results suggest that the HpnCDE enzymes are derived from carotenoid biosynthesis ones and that they assembled in an ancestral squalene pathway in bacteria, while remaining metabolically versatile. Conversely, the Sqs enzyme is prone to be involved in lateral gene transfer, and its emergence is possibly related to the specialization of squalene biosynthesis. The biosynthesis of hopanoids seems to be ancestral in the Bacteria domain. Moreover, no triterpene cyclases are found in Archaea, invoking a potential scenario in which eukaryotic genes for sterol biosynthesis assembled from ancestral bacterial contributions in early eukaryotic lineages.

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.

Prokaryotic Response to Phytodetritus-Derived Organic Material in Epi- and Mesopelagic Antarctic Waters

Particulate organic matter (POM) export represents the underlying principle of the biological carbon pump, driving the carbon flux from the sunlit to the dark ocean. The efficiency of this process is tightly linked to the prokaryotic community, as >70% of POM respiration is carried out by particle-associated prokaryotes. In the Ross Sea, one of the most productive areas of the Southern Ocean, up to 50% of the surface primary production is exported to the mesopelagic ocean as POM. Recent evidence suggests that a significant fraction of the POM in this area is composed of intact phytoplankton cells. During austral summer 2017, we set up bottle enrichment experiments in which we amended free-living surface and deep prokaryotic communities with organic matter pools generated from native microplankton, mimicking the particle export that may derive from mild (1 μg of Chlorophyll a L-1) and intense (10 μg of Chlorophyll a L-1) phytoplankton bloom. Over a course of 4 days, we followed free-living and particle-attached prokaryotes' abundance, the degradation rates of polysaccharides, proteins and lipids, heterotrophic production as well as inorganic carbon utilization and prokaryotic community structure dynamics. Our results showed that several rare or undetected taxa in the initial community became dominant during the time course of the incubations and that different phytodetritus-derived organic matter sources induced specific changes in microbial communities, selecting for peculiar degradation and utilization processes spectra. Moreover, the features of the supplied detritus (in terms of microplankton taxa composition) determined different colonization dynamics and organic matter processing modes. Our study provides insights into the mechanisms underlying the prokaryotic utilization of phytodetritus, a significant pool of organic matter in the dark ocean.

Inorganic and Organic Carbon Uptake Processes and Their Connection to Microbial Diversity in Meso- and Bathypelagic Arctic Waters (Eastern Fram Strait)

The deep Arctic Ocean is increasingly vulnerable to climate change effects, yet our understanding of its microbial processes is limited. We collected samples from shelf waters, mesopelagic Atlantic Waters (AW) and bathypelagic Norwegian Sea Deep Waters (NSDW) in the eastern Fram Strait, along coast-to-offshore transects off Svalbard during boreal summer. We measured community respiration, heterotrophic carbon production (HCP), and dissolved inorganic carbon utilization (DICu) together with prokaryotic abundance, diversity, and metagenomic predictions. In deep samples, HCP was significantly faster in AW than in NSDW, while we observed no differences in DICu rates. Organic carbon uptake was higher than its inorganic counterpart, suggesting a major reliance of deep microbial Arctic communities on heterotrophic metabolism. Community structure and spatial distribution followed the hydrography of water masses. Distinct from other oceans, the most abundant OTU in our deep samples was represented by the archaeal MG-II. To address the potential biogeochemical role of each water mass-specific microbial community, as well as their link with the measured rates, PICRUSt-based predicted metagenomes were built. The results showed that pathways of auto- and heterotrophic carbon utilization differed between the deep water masses, although this was not reflected in measured DICu rates. Our findings provide new insights to understand microbial processes and diversity in the dark Arctic Ocean and to progress toward a better comprehension of the biogeochemical cycles and their trends in light of climate changes.

Division of labour in a matrix, rather than phagocytosis or endosymbiosis, as a route for the origin of eukaryotic cells

Abstract

Two apparently irreconcilable models dominate research into the origin of eukaryotes. In one model, amitochondrial proto-eukaryotes emerged autogenously from the last universal common ancestor of all cells. Proto-eukaryotes subsequently acquired mitochondrial progenitors by the phagocytic capture of bacteria. In the second model, two prokaryotes, probably an archaeon and a bacterial cell, engaged in prokaryotic endosymbiosis, with the species resident within the host becoming the mitochondrial progenitor. Both models have limitations. A search was therefore undertaken for alternative routes towards the origin of eukaryotic cells. The question was addressed by considering classes of potential pathways from prokaryotic to eukaryotic cells based on considerations of cellular topology. Among the solutions identified, one, called here the “third-space model”, has not been widely explored. A version is presented in which an extracellular space (the third-space), serves as a proxy cytoplasm for mixed populations of archaea and bacteria to “merge” as a transitionary complex without obligatory endosymbiosis or phagocytosis and to form a precursor cell. Incipient nuclei and mitochondria diverge by division of labour. The third-space model can accommodate the reorganization of prokaryote-like genomes to a more eukaryote-like genome structure. Nuclei with multiple chromosomes and mitosis emerge as a natural feature of the model. The model is compatible with the loss of archaeal lipid biochemistry while retaining archaeal genes and provides a route for the development of membranous organelles such as the Golgi apparatus and endoplasmic reticulum. Advantages, limitations and variations of the “third-space” models are discussed.

Reviewers

This article was reviewed by Damien Devos, Buzz Baum and Michael Gray.

Development of a Vanillate Biosensor for the Vanillin Biosynthesis Pathway in E. coli

The engineered de novo vanillin biosynthesis pathway constructed in Escherichia coli is industrially relevant but limited by the reaction catalyzed by catechol O-methyltransferase, which is intended to catalyze the conversion of protocatechuate to vanillate. To identify alternative O-methyltransferases, we constructed a vanillate sensor based on the Caulobacter crescentus VanR-VanO system. Using an E. coli promoter library, we achieved greater than 14-fold dynamic range in our best rationally constructed sensor. We found that this construct and an evolved variant demonstrate remarkable substrate selectivity, exhibiting no detectable response to the regioisomer byproduct isovanillate and minimal response to structurally similar pathway intermediates. We then harnessed the evolved biosensor to conduct rapid bioprospecting of natural catechol O-methyltransferases and identified three previously uncharacterized but active O-methyltransferases. Collectively, these efforts enrich our knowledge of how biosensing can aid metabolic engineering and constitute the foundation for future improvements in vanillin pathway productivity.

Depth distributions of nitrite reductase (nirK) gene variants reveal spatial dynamics of thaumarchaeal ecotype populations in coastal Monterey Bay

Ammonia-oxidizing archaea (AOA) of the phylum Thaumarchaeota are key players in nutrient cycling, yet large gaps remain in our understanding of their ecology and metabolism. Despite multiple lines of evidence pointing to a central role for copper-containing nitrite reductase (NirK) in AOA metabolism, the thaumarchaeal nirK gene is rarely studied in the environment. In this study, we examine the diversity of nirK in the marine pelagic environment, in light of previously described ecological patterns of pelagic thaumarchaeal populations. Phylogenetic analyses show that nirK better resolves diversification patterns of marine Thaumarchaeota, compared to the conventionally used marker gene amoA. Specifically, we demonstrate that the three major phylogenetic clusters of marine nirK correspond to the three 'ecotype' populations of pelagic Thaumarchaeota. In this context, we further examine the relative distributions of the three variant groups in metagenomes and metatranscriptomes representing two depth profiles in coastal Monterey Bay. Our results reveal that nirK effectively tracks the dynamics of thaumarchaeal ecotype populations, particularly finer-scale diversification patterns within major lineages. We also find evidence for multiple copies of nirK per genome in a fraction of thaumarchaeal cells in the water column, which must be taken into account when using it as a molecular marker.

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.

Phylogenomics suggests oxygen availability as a driving force in Thaumarchaeota evolution

Ammonia-oxidizing archaea (AOA) of the phylum Thaumarchaeota are widespread in marine and terrestrial habitats, playing a major role in the global nitrogen cycle. However, their evolutionary history remains unexplored, which limits our understanding of their adaptation mechanisms. Here, our comprehensive phylogenomic tree of Thaumarchaeota supports three sequential events: origin of AOA from terrestrial non-AOA ancestors, colonization of the shallow ocean, and expansion to the deep ocean. Careful molecular dating suggests that these events coincided with the Great Oxygenation Event around 2300 million years ago (Mya), and oxygenation of the shallow and deep ocean around 800 and 635-560 Mya, respectively. The first transition was likely enabled by the gain of an aerobic pathway for energy production by ammonia oxidation and biosynthetic pathways for cobalamin and biotin that act as cofactors in aerobic metabolism. The first transition was also accompanied by the loss of dissimilatory nitrate and sulfate reduction, loss of oxygen-sensitive pyruvate oxidoreductase, which reduces pyruvate to acetyl-CoA, and loss of the Wood-Ljungdahl pathway for anaerobic carbon fixation. The second transition involved gain of a K⁺ transporter and of the biosynthetic pathway for ectoine, which may function as an osmoprotectant. The third transition was accompanied by the loss of the uvr system for repairing ultraviolet light-induced DNA lesions. We conclude that oxygen availability drove the terrestrial origin of AOA and their expansion to the photic and dark oceans, and that the stressors encountered during these events were partially overcome by gene acquisitions from Euryarchaeota and Bacteria, among other sources.

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.

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

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

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.

A phylogenomic and ecological analysis of the globally abundant Marine Group II archaea (Ca. Poseidoniales ord. nov.)

Marine Group II (MGII) archaea represent the most abundant planktonic archaeal group in ocean surface waters, but our understanding of the group has been limited by a lack of cultured representatives and few sequenced genomes. Here, we conducted a comparative phylogenomic analysis of 270 recently available MGII metagenome-assembled genomes (MAGs) to investigate their evolution and ecology. Based on a rank-normalised genome phylogeny, we propose that MGII is an order-level lineage for which we propose the name Candidatus Poseidoniales (after Gr. n. Poseidon, God of the sea), comprising the families Candidatus Poseidonaceae fam. nov. (formerly subgroup MGIIa) and Candidatus Thalassarchaeaceae fam. nov. (formerly subgroup MGIIb). Within these families, 21 genera could be resolved, many of which had distinct biogeographic ranges and inferred nutrient preferences. Phylogenetic analyses of key metabolic functions suggest that the ancestor of Ca. Poseidoniales was a surface water-dwelling photoheterotroph that evolved to occupy multiple related ecological niches based primarily on spectral tuning of proteorhodopsin genes. Interestingly, this adaptation appears to involve an overwrite mechanism whereby an existing single copy of the proteorhodopsin gene is replaced by a horizontally transferred copy, which in many instances should allow an abrupt change in light absorption capacity. Phototrophy was lost entirely from five Ca. Poseidoniales genera coinciding with their adaptation to deeper aphotic waters. We also report the first instances of nitrate reductase in two genera acquired via horizontal gene transfer (HGT), which was a potential adaptation to oxygen limitation. Additional metabolic traits differentiating families and genera include flagellar-based adhesion, transporters, and sugar, amino acid, and peptide degradation. Our results suggest that HGT has shaped the evolution of Ca. Poseidoniales to occupy a variety of ecological niches and to become the most successful archaeal lineage in ocean surface waters.

Convergent evolution of unusual complex I homologs with increased proton pumping capacity: energetic and ecological implications

Respiratory complex I is part of a large family of homologous enzymes that carry out the transfer of electrons between soluble cytoplasmic electron carriers and membrane-bound electron carriers. These complexes are vital bioenergetic enzymes that serve as the entry points into electron transport chains for a wide variety of microbial metabolisms, and electron transfer is coupled to proton translocation. The core complex of this enzyme is made up of 11 protein subunits, with three major proton pumping subunits. Here, we document a large number of modified complex I gene cassettes found in genome sequences from diverse cultured bacteria, shotgun metagenomics, and environmentally derived archaeal fosmids all of which encode a fourth proton pumping subunit. The incorporation of this extra subunit into a functional protein complex is supported by large amino acid insertions in the amphipathic helix that runs the length of the protein complex. Phylogenetic analyses reveal that these modified complexes appear to have arisen independently multiple times in a remarkable case of convergent molecular evolution. From an energetic perspective, we hypothesize that this modification on the canonical complex I architecture allows for the translocation of a fifth proton per reaction cycle-the physiological utility of this modified complex is discussed.

The Rhizome of Lokiarchaeota Illustrates the Mosaicity of Archaeal Genomes

Genome remodeling and exchange of sequences are widespread in the prokaryotic world and mosaic genomes challenge the classification of prokaryotes, which cannot be properly achieved in terms of a single gene or group of genes. Here, we studied individually the gene collection of the archaic microorganism Lokiarchaeum sp., suggested as an archaeal host close to the emergence of the eukaryotes. The network or rhizome of all Lokiarchaeum sp. genes revealed that the genomic repertoire is mainly composed of genes from archaeal (∼36%) and bacterial origin (∼28%), distantly followed by components of eukaryotic origin (∼2%). Thirty-three percent of genes were unique to this species (ORFans). The mosaicity of archaea was also supported by studying Methanomassiliicoccus luminyensis, an archaea from the gut, in which 67% of the genomic repertoire arised from archaea and 22% from bacteria. Our results illustrate the intricate evolutionary relationships of the archaeal genome repertoire and highlight the rhizome-like processes of evolution in archaea, their mosaicity, and chimeric origin composed of different domains of life, questioning the reality of a tree of life.

Does the Promoter Constitute a Barrier in the Horizontal Transposon Transfer Process? Insight from Bari Transposons

The contribution of the transposons’ promoter in the horizontal transfer process is quite overlooked in the scientific literature. To shed light on this aspect we have mimicked the horizontal transfer process in laboratory and assayed in a wide range of hosts (fly, human, yeast and bacteria) the promoter activity of the 5′ terminal sequences in Bari1 and Bari3, two Drosophila transposons belonging to the Tc1-mariner superfamily. These sequences are able to drive the transcription of a reporter gene even in distantly related organisms at least at the episomal level. By combining bioinformatics and experimental approaches, we define two distinct promoter sequences for each terminal sequence analyzed, which allow transcriptional activity in prokaryotes and eukaryotes, respectively. We propose that the Bari family of transposons, and possibly other members of the Tc1-mariner superfamily, might have evolved “blurry promoters,” which have facilitated their diffusion in many living organisms through horizontal transfer.

The enigmatic SAR202 cluster up close: shedding light on a globally distributed dark ocean lineage involved in sulfur cycling

The dark ocean microbiota represents the unknown majority in the global ocean waters. The SAR202 cluster belonging to the phylum Chloroflexi was the first microbial lineage discovered to specifically inhabit the aphotic realm, where they are abundant and globally distributed. The absence of SAR202 cultured representatives is a significant bottleneck towards understanding their metabolic capacities and role in the marine environment. In this work, we use a combination of metagenome-assembled genomes from deep-sea datasets and publicly available single-cell genomes to construct a genomic perspective of SAR202 phylogeny, metabolism and biogeography. Our results suggest that SAR202 cluster members are medium sized, free-living cells with a heterotrophic lifestyle, broadly divided into two distinct clades. We present the first evidence of vertical stratification of these microbes along the meso- and bathypelagic ocean layers. Remarkably, two distinct species of SAR202 cluster are highly abundant in nearly all deep bathypelagic metagenomic datasets available so far. SAR202 members metabolize multiple organosulfur compounds, many appear to be sulfite-oxidizers and are predicted to play a major role in sulfur turnover in the dark water column. This concomitantly suggests an unsuspected availability of these nutrient sources to allow for the high abundance of these microbes in the deep sea.

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.

Living Organisms Author Their Read-Write Genomes in Evolution

Evolutionary variations generating phenotypic adaptations and novel taxa resulted from complex cellular activities altering genome content and expression: (i) Symbiogenetic cell mergers producing the mitochondrion-bearing ancestor of eukaryotes and chloroplast-bearing ancestors of photosynthetic eukaryotes; (ii) interspecific hybridizations and genome doublings generating new species and adaptive radiations of higher plants and animals; and, (iii) interspecific horizontal DNA transfer encoding virtually all of the cellular functions between organisms and their viruses in all domains of life. Consequently, assuming that evolutionary processes occur in isolated genomes of individual species has become an unrealistic abstraction. Adaptive variations also involved natural genetic engineering of mobile DNA elements to rewire regulatory networks. In the most highly evolved organisms, biological complexity scales with "non-coding" DNA content more closely than with protein-coding capacity. Coincidentally, we have learned how so-called "non-coding" RNAs that are rich in repetitive mobile DNA sequences are key regulators of complex phenotypes. Both biotic and abiotic ecological challenges serve as triggers for episodes of elevated genome change. The intersections of cell activities, biosphere interactions, horizontal DNA transfers, and non-random Read-Write genome modifications by natural genetic engineering provide a rich molecular and biological foundation for understanding how ecological disruptions can stimulate productive, often abrupt, evolutionary transformations.

Proposed Role for KaiC-Like ATPases as Major Signal Transduction Hubs in Archaea

All organisms must adapt to ever-changing environmental conditions and accordingly have evolved diverse signal transduction systems. In bacteria, the most abundant networks are built around the two-component signal transduction systems that include histidine kinases and receiver domains. In contrast, eukaryotic signal transduction is dominated by serine/threonine/tyrosine protein kinases. Both of these systems are also found in archaea, but they are not as common and diversified as their bacterial and eukaryotic counterparts, suggesting the possibility that archaea have evolved other, still uncharacterized signal transduction networks. Here we propose a role for KaiC family ATPases, known to be key components of the circadian clock in cyanobacteria, in archaeal signal transduction. The KaiC family is notably expanded in most archaeal genomes, and although most of these ATPases remain poorly characterized, members of the KaiC family have been shown to control archaellum assembly and have been found to be a stable component of the gas vesicle system in Halobacteria Computational analyses described here suggest that KaiC-like ATPases and their homologues with inactivated ATPase domains are involved in many other archaeal signal transduction pathways and comprise major hubs of complex regulatory networks. We predict numerous input and output domains that are linked to KaiC-like proteins, including putative homologues of eukaryotic DEATH domains that could function as adapters in archaeal signaling networks. We further address the relationships of the archaeal family of KaiC homologues to the bona fide KaiC of cyanobacteria and implications for the existence of a KaiC-based circadian clock apparatus in archaea.IMPORTANCE Little is currently known about signal transduction pathways in Archaea Recent studies indicate that KaiC-like ATPases, known as key components of the circadian clock apparatus in cyanobacteria, are involved in the regulation of archaellum assembly and, likely, type IV pili and the gas vesicle system in Archaea We performed comprehensive comparative genomic analyses of the KaiC family. A vast protein interaction network was revealed, with KaiC family proteins as hubs for numerous input and output components, many of which are shared with two-component signal transduction systems. Putative KaiC-based signal transduction systems are predicted to regulate the activities of membrane-associated complexes and individual proteins, such as signal recognition particle and membrane transporters, and also could be important for oxidative stress response regulation. KaiC-centered signal transduction networks are predicted to play major roles in archaeal physiology, and this work is expected to stimulate their experimental characterization.

Integrative modeling of gene and genome evolution roots the archaeal tree of life

A root for the archaeal tree is essential for reconstructing the metabolism and ecology of early cells and for testing hypotheses that propose that the eukaryotic nuclear lineage originated from within the Archaea; however, published studies based on outgroup rooting disagree regarding the position of the archaeal root. Here we constructed a consensus unrooted archaeal topology using protein concatenation and a multigene supertree method based on 3,242 single gene trees, and then rooted this tree using a recently developed model of genome evolution. This model uses evidence from gene duplications, horizontal transfers, and gene losses contained in 31,236 archaeal gene families to identify the most likely root for the tree. Our analyses support the monophyly of DPANN (Diapherotrites, Parvarchaeota, Aenigmarchaeota, Nanoarchaeota, Nanohaloarchaea), a recently discovered cosmopolitan and genetically diverse lineage, and, in contrast to previous work, place the tree root between DPANN and all other Archaea. The sister group to DPANN comprises the Euryarchaeota and the TACK Archaea, including Lokiarchaeum, which our analyses suggest are monophyletic sister lineages. Metabolic reconstructions on the rooted tree suggest that early Archaea were anaerobes that may have had the ability to reduce CO2 to acetate via the Wood-Ljungdahl pathway. In contrast to proposals suggesting that genome reduction has been the predominant mode of archaeal evolution, our analyses infer a relatively small-genomed archaeal ancestor that subsequently increased in complexity via gene duplication and horizontal gene transfer.

Genome and epigenome of a novel marine Thaumarchaeota strain suggest viral infection, phosphorothioation DNA modification and multiple restriction systems

Marine Thaumarchaeota are abundant ammonia-oxidizers but have few representative laboratory-cultured strains. We report the cultivation of Candidatus Nitrosomarinus catalina SPOT01, a novel strain that is less warm-temperature tolerant than other cultivated Thaumarchaeota. Using metagenomic recruitment, strain SPOT01 comprises a major portion of Thaumarchaeota (4-54%) in temperate Pacific waters. Its complete 1.36 Mbp genome possesses several distinguishing features: putative phosphorothioation (PT) DNA modification genes; a region containing probable viral genes; and putative urea utilization genes. The PT modification genes and an adjacent putative restriction enzyme (RE) operon likely form a restriction modification (RM) system for defence from foreign DNA. PacBio sequencing showed >98% methylation at two motifs, and inferred PT guanine modification of 19% of possible TGCA sites. Metagenomic recruitment also reveals the putative virus region and PT modification and RE genes are present in 18-26%, 9-14% and <1.5% of natural populations at 150 m with ≥85% identity to strain SPOT01. The presence of multiple probable RM systems in a highly streamlined genome suggests a surprising importance for defence from foreign DNA for dilute populations that infrequently encounter viruses or other cells. This new strain provides new insights into the ecology, including viral interactions, of this important group of marine microbes.

Symbiosis in eukaryotic evolution

Fifty years ago, Lynn Margulis, inspiring in early twentieth-century ideas that put forward a symbiotic origin for some eukaryotic organelles, proposed a unified theory for the origin of the eukaryotic cell based on symbiosis as evolutionary mechanism. Margulis was profoundly aware of the importance of symbiosis in the natural microbial world and anticipated the evolutionary significance that integrated cooperative interactions might have as mechanism to increase cellular complexity. Today, we have started fully appreciating the vast extent of microbial diversity and the importance of syntrophic metabolic cooperation in natural ecosystems, especially in sediments and microbial mats. Also, not only the symbiogenetic origin of mitochondria and chloroplasts has been clearly demonstrated, but improvement in phylogenomic methods combined with recent discoveries of archaeal lineages more closely related to eukaryotes further support the symbiogenetic origin of the eukaryotic cell. Margulis left us in legacy the idea of 'eukaryogenesis by symbiogenesis'. Although this has been largely verified, when, where, and specifically how eukaryotic cells evolved are yet unclear. Here, we shortly review current knowledge about symbiotic interactions in the microbial world and their evolutionary impact, the status of eukaryogenetic models and the current challenges and perspectives ahead to reconstruct the evolutionary path to eukaryotes.

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.

Horizontal gene transfer drives the evolution of Rh50 permeases in prokaryotes

Background

Rh50 proteins belong to the family of ammonia permeases together with their Amt/MEP homologs. Ammonia permeases increase the permeability of NH₃/NH₄⁺ across cell membranes and are believed to be involved in excretion of toxic ammonia and in the maintenance of pH homeostasis. RH50 genes are widespread in eukaryotes but absent in land plants and fungi, and remarkably rare in prokaryotes. The evolutionary history of RH50 genes in prokaryotes is just beginning to be unveiled.

Results

Here, a molecular phylogenetic approach suggests horizontal gene transfer (HGT) as a primary force driving the evolution and spread of RH50 among prokaryotes. In addition, the taxonomic distribution of the RH50 gene among prokaryotes turned out to be very narrow; a single-copy RH50 is present in the genome of only a small proportion of Bacteria, and, first evidence to date, in only three methanogens among Euryarchaea. The coexistence of RH50 and AMT in prokaryotes seems also a rare event. Finally, phylogenetic analyses were used to reconstruct the HGT network along which prokaryotic RH50 evolution has taken place.

Conclusions

The eukaryotic or bacterial “origin” of the RH50 gene remains unsolved. The RH50 prokaryotic HGT network suggests a preferential directionality of transfer from aerobic to anaerobic organisms. The observed HGT events between archaeal methanogens, anaerobic and aerobic ammonia-oxidizing bacteria suggest that syntrophic relationships play a major role in the structuring of the network, and point to oxygen minimum zones as an ecological niche that might be of crucial importance for HGT-driven evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0850-6) contains supplementary material, which is available to authorized users.

Identity of major sulfur-cycle prokaryotes in freshwater lake ecosystems revealed by a comprehensive phylogenetic study of the dissimilatory adenylylsulfate reductase

Adenylylsulfate reductase is a heterodimeric complex of two subunits, AprB and AprA, and is a key enzyme in dissimilatory sulfate reduction and sulfur oxidation. Common use of aprA as a functional marker gene has revealed the diversity of sulfur-cycle prokaryotes in diverse environments. In this study, we established a comprehensive sequence set of apr genes and employed it to reanalyze apr phylogeny, evaluate the coverage of a widely used primer set (AprA-1-FW/AprA-5-RV), and categorize environmental aprA sequences. Phylogenetic tree construction revealed new members of Apr lineage II and several previously unrecognized lateral gene transfer events. Using the established phylogenetic tree, we classified all previously reported aprA sequences amplified from freshwater lakes with the primer pair AprA-1-FW/AprA-5-RV in addition to the aprA sequences newly retrieved from freshwater lakes; the obtained results were complemented by 16S rRNA clone library analysis. Apr-based classifications of some of operational taxonomic units were supported by 16S rRNA-based analysis. This study updates our knowledge on the phylogeny of aprBA and shows the identities of several sulfur-cycle bacteria, which could not be classified to a known taxa until now. The established apr sequence set is publicly available and can be applied to assign environmental sequences to known lineages.

mdRNA-Seq analysis of marine microbial communities from the northern Red Sea

Metatranscriptomic differential RNA-Seq (mdRNA-Seq) identifies the suite of active transcriptional start sites at single-nucleotide resolution through enrichment of primary transcript 5′ ends. Here we analyzed the microbial community at 45 m depth at Station A in the northern Gulf of Aqaba, Red Sea, during 500 m deep mixing in February 2012 using mdRNA-Seq and a parallel classical RNA-Seq approach. We identified promoters active in situ for five different pico-planktonic genera (the SAR11 clade of Alphaproteobacteria, Synechococcus of Cyanobacteria, Euryarchaeota, Thaumarchaeota, and Micromonas as an example for picoeukaryotic algae), showing the applicability of this approach to highly diverse microbial communities. 16S rDNA quantification revealed that 24% of the analyzed community were group II marine Euryarchaeota in which we identified a highly abundant non-coding RNA, Tan1, and detected very high expression of genes encoding intrinsically disordered proteins, as well as enzymes for the synthesis of specific B vitamins, extracellular peptidases, carbohydrate-active enzymes, and transport systems. These results highlight previously unknown functions of Euryarchaeota with community-wide relevance. The complementation of metatranscriptomic studies with mdRNA-Seq provides substantial additional information regarding transcriptional start sites, promoter activities, and the identification of non-coding RNAs.