Acquisition of 1,000 eubacterial genes physiologically transformed a methanogen at the origin of Haloarchaea

Navigating the archaeal frontier: insights and projections from bioinformatic pipelines

Archaea continues to be one of the least investigated domains of life, and in recent years, the advent of metagenomics has led to the discovery of many new lineages at the phylum level. For the majority, only automatic genomic annotations can provide information regarding their metabolic potential and role in the environment. Here, genomic data from 2,978 archaeal genomes was used to perform automatic annotations using bioinformatics tools, alongside synteny analysis. These automatic classifications were done to assess how good these different tools perform in relation to archaeal data. Our study revealed that even with lowered cutoffs, several functional models do not capture the recently discovered archaeal diversity. Moreover, our investigation revealed that a significant portion of archaeal genomes, approximately 42%, remain uncharacterized. In comparison, within 3,235 bacterial genomes, a diverse range of unclassified proteins is obtained, with well-studied organisms like Escherichia coli having a substantially lower proportion of uncharacterized regions, ranging from <5 to 25%, and less studied lineages being comparable to archaea with the range of 35-40% of unclassified regions. Leveraging this analysis, we were able to identify metabolic protein markers, thereby providing insights into the metabolism of the archaea in our dataset. Our findings underscore a substantial gap between automatic classification tools and the comprehensive mapping of archaeal metabolism. Despite advances in computational approaches, a significant portion of archaeal genomes remains unexplored, highlighting the need for extensive experimental validation in this domain, as well as more refined annotation methods. This study contributes to a better understanding of archaeal metabolism and underscores the importance of further research in elucidating the functional potential of archaeal genomes.

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.

Identification and characterization of a novel type of ketohexokinase from the haloarchaeon Haloferax volcanii

Ketohexokinase (KHK) catalyzes the ATP-dependent phosphorylation of fructose, forming fructose-1-phosphate and ADP. The enzyme is well studied in Eukarya, in particular in humans and other vertebrates, but homologs have not been identified in Bacteria and Archaea. Here we report the identification of a novel type of KHK from the haloarchaeon Haloferax volcanii (HvKHK). The encoding gene khk was identified as HVO_1812. The gene was expressed as a 90-kDa homodimeric protein, catalyzing the phosphorylation of fructose with a Vmax value of 59 U/mg and apparent KM values for ATP and fructose of 0.47 and 1.29 mM, respectively. Homologs of HvKHK were only identified in a few haloarchaea and halophilic Bacteria. The protein showed low sequence identity to characterized KHKs from Eukarya and phylogenetic analyses indicate that haloarchaeal KHKs are largely separated from eukaryal KHKs. This is the first report of the identification of KHKs in prokaryotes that form a novel cluster of sugar kinases within the ribokinase/pfkB superfamily.

Members of the class Candidatus Ordosarchaeia imply an alternative evolutionary scenario from methanogens to haloarchaea

The origin of methanogenesis can be traced to the common ancestor of non-DPANN archaea, whereas haloarchaea (or Halobacteria) are believed to have evolved from a methanogenic ancestor through multiple evolutionary events. However, due to the accelerated evolution and compositional bias of proteins adapting to hypersaline habitats, Halobacteria exhibit substantial evolutionary divergence from methanogens, and the identification of the closest methanogen (either Methanonatronarchaeia or other taxa) to Halobacteria remains a subject of debate. Here, we obtained five metagenome-assembled genomes with high completeness from soda-saline lakes on the Ordos Plateau in Inner Mongolia, China, and we proposed the name Candidatus Ordosarchaeia for this novel class. Phylogenetic analyses revealed that Ca. Ordosarchaeia is firmly positioned near the median position between the Methanonatronarchaeia and Halobacteria-Hikarchaeia lineages. Functional predictions supported the transitional status of Ca. Ordosarchaeia with the metabolic potential of nonmethanogenic and aerobic chemoheterotrophy, as did remnants of the gene sequences of methylamine/dimethylamine/trimethylamine metabolism and coenzyme M biosynthesis. Based on the similarity of the methyl-coenzyme M reductase genes mcrBGADC in Methanonatronarchaeia with the phylogenetically distant methanogens, an alternative evolutionary scenario is proposed, in which Methanonatronarchaeia, Ca. Ordosarchaeia, Ca. Hikarchaeia, and Halobacteria share a common ancestor that initially lost mcr genes. However, certain members of Methanonatronarchaeia subsequently acquired mcr genes through horizontal gene transfer from distantly related methanogens. This hypothesis is supported by amalgamated likelihood estimation, phylogenetic analysis, and gene arrangement patterns. Altogether, Ca. Ordosarchaeia genomes clarify the sisterhood of Methanonatronarchaeia with Halobacteria and provide new insights into the evolution from methanogens to haloarchaea.

Environmental selection and evolutionary process jointly shape genomic and functional profiles of mangrove rhizosphere microbiomes

Mangrove reforestation with introduced species has been an important strategy to restore mangrove ecosystem functioning. However, how such activities affect microbially driven methane (CH4), nitrogen (N), and sulfur (S) cycling of rhizosphere microbiomes remains unclear. To understand the effect of environmental selection and the evolutionary process on microbially driven biogeochemical cycles in native and introduced mangrove rhizospheres, we analyzed key genomic and functional profiles of rhizosphere microbiomes from native and introduced mangrove species by metagenome sequencing technologies. Compared with the native mangrove (Kandelia obovata, KO), the introduced mangrove (Sonneratia apetala, SA) rhizosphere microbiome had significantly (p < 0.05) higher average genome size (AGS) (5.8 vs. 5.5 Mb), average 16S ribosomal RNA gene copy number (3.5 vs. 3.1), relative abundances of mobile genetic elements, and functional diversity in terms of the Shannon index (7.88 vs. 7.84) but lower functional potentials involved in CH4 cycling (e.g., mcrABCDG and pmoABC), N2 fixation (nifHDK), and inorganic S cycling (dsrAB, dsrC, dsrMKJOP, soxB, sqr, and fccAB). Similar results were also observed from the recovered Proteobacterial metagenome-assembled genomes with a higher AGS and distinct functions in the introduced mangrove rhizosphere. Additionally, salinity and ammonium were identified as the main environmental drivers of functional profiles of mangrove rhizosphere microbiomes through deterministic processes. This study advances our understanding of microbially mediated biogeochemical cycling of CH4, N, and S in the mangrove rhizosphere and provides novel insights into the influence of environmental selection and evolutionary processes on ecosystem functions, which has important implications for future mangrove reforestation.

Dissecting the HGT network of carbon metabolic genes in soil-borne microbiota

The microbiota inhabiting soil plays a significant role in essential life-supporting element cycles. Here, we investigated the occurrence of horizontal gene transfer (HGT) and established the HGT network of carbon metabolic genes in 764 soil-borne microbiota genomes. Our study sheds light on the crucial role of HGT components in microbiological diversification that could have far-reaching implications in understanding how these microbial communities adapt to changing environments, ultimately impacting agricultural practices. In the overall HGT network of carbon metabolic genes in soil-borne microbiota, a total of 6,770 nodes and 3,812 edges are present. Among these nodes, phyla Proteobacteria, Actinobacteriota, Bacteroidota, and Firmicutes are predominant. Regarding specific classes, Actinobacteria, Gammaproteobacteria, Alphaproteobacteria, Bacteroidia, Actinomycetia, Betaproteobacteria, and Clostridia are dominant. The Kyoto Encyclopedia of Genes and Genomes (KEGG) functional assignments of glycosyltransferase (18.5%), glycolysis/gluconeogenesis (8.8%), carbohydrate-related transporter (7.9%), fatty acid biosynthesis (6.5%), benzoate degradation (3.1%) and butanoate metabolism (3.0%) are primarily identified. Glycosyltransferase involved in cell wall biosynthesis, glycosylation, and primary/secondary metabolism (with 363 HGT entries), ranks first overwhelmingly in the list of most frequently identified carbon metabolic HGT enzymes, followed by pimeloyl-ACP methyl ester carboxylesterase, alcohol dehydrogenase, and 3-oxoacyl-ACP reductase. Such HGT events mainly occur in the peripheral functions of the carbon metabolic pathway instead of the core section. The inter-microbe HGT genetic traits in soil-borne microbiota genetic sequences that we recognized, as well as their involvement in the metabolism and regulation processes of carbon organic, suggest a pervasive and substantial effect of HGT on the evolution of microbes.

A Transiently Hypersaline Microbial Mat Harbors a Diverse and Stable Archaeal Community in the Cuatro Cienegas Basin, Mexico

Microbial mats are biologically diverse communities that are analogs to some of the earliest ecosystems on Earth. In this study, we describe a unique transiently hypersaline microbial mat uncovered in a shallow pond within the Cuatro Cienegas Basin (CCB) in northern México. The CCB is an endemism-rich site that harbors living stromatolites that have been studied to understand the conditions of the Precambrian Earth. These microbial mats form elastic domes filled with biogenic gas, and the mats have a relatively large and stable subpopulation of archaea. For this reason, this site has been termed archaean domes (AD). The AD microbial community was analyzed by metagenomics over three seasons. The mat exhibited a highly diverse prokaryotic community dominated by bacteria. Bacterial sequences are represented in 37 phyla, mainly Proteobacteria, Firmicutes, and Actinobacteria, that together comprised >50% of the sequences from the mat. Archaea represented up to 5% of the retrieved sequences, with up to 230 different archaeal species that belong to 5 phyla (Euryarchaeota, Crenarchaeota, Thaumarchaeota, Korarchaeota, and Nanoarchaeota). The archaeal taxa showed low variation despite fluctuations in water and nutrient availability. In addition, predicted functions highlight stress responses to extreme conditions present in the AD, including salinity, pH, and water/drought fluctuation. The observed complexity of the AD mat thriving in high pH and fluctuating water and salt conditions within the CCB provides an extant model of great value for evolutionary studies, as well as a suitable analog to the early Earth and Mars.

Antimony efflux underpins phosphorus cycling and resistance of phosphate-solubilizing bacteria in mining soils

Microorganisms play crucial roles in phosphorus (P) turnover and P bioavailability increases in heavy metal-contaminated soils. However, microbially driven P-cycling processes and mechanisms of their resistance to heavy metal contaminants remain poorly understood. Here, we examined the possible survival strategies of P-cycling microorganisms in horizontal and vertical soil samples from the world's largest antimony (Sb) mining site, which is located in Xikuangshan, China. We found that total soil Sb and pH were the primary factors affecting bacterial community diversity, structure and P-cycling traits. Bacteria with the gcd gene, encoding an enzyme responsible for gluconic acid production, largely correlated with inorganic phosphate (Pi) solubilization and significantly enhanced soil P bioavailability. Among the 106 nearly complete bacterial metagenome-assembled genomes (MAGs) recovered, 60.4% carried the gcd gene. Pi transportation systems encoded by pit or pstSCAB were widely present in gcd-harboring bacteria, and 43.8% of the gcd-harboring bacteria also carried the acr3 gene encoding an Sb efflux pump. Phylogenetic and potential horizontal gene transfer (HGT) analyses of acr3 indicated that Sb efflux could be a dominant resistance mechanism, and two gcd-harboring MAGs appeared to acquire acr3 through HGT. The results indicated that Sb efflux could enhance P cycling and heavy metal resistance in Pi-solubilizing bacteria in mining soils. This study provides novel strategies for managing and remediating heavy metal-contaminated ecosystems.

Modular structure of complex II: An evolutionary perspective

Succinate dehydrogenases (SDHs) and fumarate reductases (FRDs) catalyse the interconversion of succinate and fumarate, a reaction highly conserved in all domains of life. The current classification of SDH/FRDs is based on the structure of the membrane anchor subunits and their cofactors. It is, however, unknown whether this classification would hold in the context of evolution. In this work, a large-scale comparative genomic analysis of complex II addresses the questions of its taxonomic distribution and phylogeny. Our findings report that for types C, D, and F, structural classification and phylogeny go hand in hand, while for types A, B and E the situation is more complex, highlighting the possibility for their classification into subgroups. Based on these findings, we proposed a revised version of the evolutionary scenario for these enzymes in which a primordial soluble module, corresponding to the cytoplasmatic subunits, would give rise to the current diversity via several independent membrane anchor attachment events.

Endosymbiotic selective pressure at the origin of eukaryotic cell biology

The dichotomy that separates prokaryotic from eukaryotic cells runs deep. The transition from pro- to eukaryote evolution is poorly understood due to a lack of reliable intermediate forms and definitions regarding the nature of the first host that could no longer be considered a prokaryote, the first eukaryotic common ancestor, FECA. The last eukaryotic common ancestor, LECA, was a complex cell that united all traits characterising eukaryotic biology including a mitochondrion. The role of the endosymbiotic organelle in this radical transition towards complex life forms is, however, sometimes questioned. In particular the discovery of the asgard archaea has stimulated discussions regarding the pre-endosymbiotic complexity of FECA. Here we review differences and similarities among models that view eukaryotic traits as isolated coincidental events in asgard archaeal evolution or, on the contrary, as a result of and in response to endosymbiosis. Inspecting eukaryotic traits from the perspective of the endosymbiont uncovers that eukaryotic cell biology can be explained as having evolved as a solution to housing a semi-autonomous organelle and why the addition of another endosymbiont, the plastid, added no extra compartments. Mitochondria provided the selective pressures for the origin (and continued maintenance) of eukaryotic cell complexity. Moreover, they also provided the energetic benefit throughout eukaryogenesis for evolving thousands of gene families unique to eukaryotes. Hence, a synthesis of the current data lets us conclude that traits such as the Golgi apparatus, the nucleus, autophagosomes, and meiosis and sex evolved as a response to the selective pressures an endosymbiont imposes.

Conserved and lineage-specific hypothetical proteins may have played a central role in the rise and diversification of major archaeal groups

BACKGROUND

Archaea play fundamental roles in the environment, for example by methane production and consumption, ammonia oxidation, protein degradation, carbon compound turnover, and sulfur compound transformations. Recent genomic analyses have profoundly reshaped our understanding of the distribution and functionalities of Archaea and their roles in eukaryotic evolution.

RESULTS

Here, 1179 representative genomes were selected from 3197 archaeal genomes. The representative genomes clustered based on the content of 10,866 newly defined archaeal protein families (that will serve as a community resource) recapitulates archaeal phylogeny. We identified the co-occurring proteins that distinguish the major lineages. Those with metabolic roles were consistent with experimental data. However, two families specific to Asgard were determined to be new eukaryotic signature proteins. Overall, the blocks of lineage-specific families are dominated by proteins that lack functional predictions.

CONCLUSIONS

Given that these hypothetical proteins are near ubiquitous within major archaeal groups, we propose that they were important in the origin of most of the major archaeal lineages. Interestingly, although there were clearly phylum-specific co-occurring proteins, no such blocks of protein families were shared across superphyla, suggesting a burst-like origin of new lineages early in archaeal evolution.

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.

Microbial diversity in extreme environments

A wide array of microorganisms, including many novel, phylogenetically deeply rooted taxa, survive and thrive in extreme environments. These unique and reduced-complexity ecosystems offer a tremendous opportunity for studying the structure, function and evolution of natural microbial communities. Marker gene surveys have resolved patterns and ecological drivers of these extremophile assemblages, revealing a vast uncultured microbial diversity and the often predominance of archaea in the most extreme conditions. New omics studies have uncovered linkages between community function and environmental variables, and have enabled discovery and genomic characterization of major new lineages that substantially expand microbial diversity and change the structure of the tree of life. These efforts have significantly advanced our understanding of the diversity, ecology and evolution of microorganisms populating Earth's extreme environments, and have facilitated the exploration of microbiota and processes in more complex ecosystems.

Distinct Expansion of Group II Introns During Evolution of Prokaryotes and Possible Factors Involved in Its Regulation

Group II introns (G2Is) are ribozymes that have retroelement characteristics in prokaryotes. Although G2Is are suggested to have been an important evolutionary factor in the prokaryote-to-eukaryote transition, comprehensive analyses of these introns among the tens of thousands of prokaryotic genomes currently available are still limited. Here, we developed a bioinformatic pipeline that systematically collects G2Is and applied it to prokaryotic genomes. We found that in bacteria, 25% (447 of 1,790) of the total representative genomes had an average of 5.3 G2Is, and in archaea, 9% (28 of 296) of the total representative genomes had an average of 3.0 G2Is. The greatest number of G2Is per genome was 101 in Arthrospira platensis (phylum Cyanobacteriota). A comprehensive sequence analysis of the intron-encoded protein (IEP) in each G2I sequence was conducted and resulted in the addition of three new IEP classes (U1-U3) to the previous classification. This analysis suggested that about 30% of all IEPs are non-canonical IEPs. The number of G2Is per genome was defined almost at the phylum level, and at least in the following two phyla, Firmicutes, and Cyanobacteriota, the type of IEP was largely associated as a factor in the G2I increase, i.e., there was an explosive increase in G2Is with bacterial C-type IEPs, mainly in the phylum Firmicutes, and in G2Is with CL-type IEPs, mainly in the phylum Cyanobacteriota. We also systematically analyzed the relationship between genomic signatures and the mechanism of these increases in G2Is. This is the first study to systematically characterize G2Is in the prokaryotic phylogenies.

Archaeal Communities: The Microbial Phylogenomic Frontier

Archaea are a unique system for investigating the diversity of life. There are the most diverse group of organisms with the longest evolutionary history of life on Earth. Phylogenomic investigations reveal the complex evolutionary history of Archaea, overturning longstanding views of the history of life. They exist in the harshest environments and benign conditions, providing a system to investigate the basis for living in extreme environments. They are frequently members of microbial communities, albeit generally rare. Archaea were central in the evolution of Eukaryotes and can be used as a proxy for studying life on other planets. Future advances will depend not only upon phylogenomic studies but also on a better understanding of isolation and cultivation techniques.

What makes a megaplasmid?

Naturally occurring plasmids come in different sizes. The smallest are less than a kilobase of DNA, while the largest can be over three orders of magnitude larger. Historically, research has tended to focus on smaller plasmids that are usually easier to isolate, manipulate and sequence, but with improved genome assemblies made possible by long-read sequencing, there is increased appreciation that very large plasmids—known as megaplasmids—are widespread, diverse, complex, and often encode key traits in the biology of their host microorganisms. Why are megaplasmids so big? What other features come with large plasmid size that could affect bacterial ecology and evolution? Are megaplasmids 'just' big plasmids, or do they have distinct characteristics? In this perspective, we reflect on the distribution, diversity, biology, and gene content of megaplasmids, providing an overview to these large, yet often overlooked, mobile genetic elements.

This article is part of the theme issue ‘The secret lives of microbial mobile genetic elements’.

Large-scale k-mer-based analysis of the informational properties of genomes, comparative genomics and taxonomy

Information theoretic approaches are ubiquitous and effective in a wide variety of bioinformatics applications. In comparative genomics, alignment-free methods, based on short DNA words, or k-mers, are particularly powerful. We evaluated the utility of varying k-mer lengths for genome comparisons by analyzing their sequence space coverage of 5805 genomes in the KEGG GENOME database. In subsequent analyses on four k-mer lengths spanning the relevant range (11, 21, 31, 41), hierarchical clustering of 1634 genus-level representative genomes using pairwise 21- and 31-mer Jaccard similarities best recapitulated a phylogenetic/taxonomic tree of life with clear boundaries for superkingdom domains and high subtree similarity for named taxons at lower levels (family through phylum). By analyzing ~14.2M prokaryotic genome comparisons by their lowest-common-ancestor taxon levels, we detected many potential misclassification errors in a curated database, further demonstrating the need for wide-scale adoption of quantitative taxonomic classifications based on whole-genome similarity.

Deciphering Symbiotic Interactions of "Candidatus Aenigmarchaeota" with Inferred Horizontal Gene Transfers and Co-occurrence Networks

"Candidatus Aenigmarchaeota" ("Ca. Aenigmarchaeota") represents one of the earliest proposed evolutionary branches within the Diapherotrites, Parvarchaeota, Aenigmarchaeota, Nanoarchaeota, and Nanohaloarchaeota (DPANN) superphylum. However, their ecological roles and potential host-symbiont interactions are still poorly understood. Here, eight metagenome-assembled genomes (MAGs) were reconstructed from hot spring ecosystems, and further in-depth comparative and evolutionary genomic analyses were conducted on these MAGs and other genomes downloaded from public databases. Although with limited metabolic capacities, we reported that "Ca. Aenigmarchaeota" in thermal environments harbor more genes related to carbohydrate metabolism than "Ca. Aenigmarchaeota" in nonthermal environments. Evolutionary analyses suggested that members from the Thaumarchaeota, Aigarchaeota, Crenarchaeota, and Korarchaeota (TACK) superphylum and Euryarchaeota contribute substantially to the niche expansion of "Ca. Aenigmarchaeota" via horizontal gene transfer (HGT), especially genes related to virus defense and stress responses. Based on co-occurrence network results and recent genetic exchanges among community members, we conjectured that "Ca. Aenigmarchaeota" may be symbionts associated with one MAG affiliated with the genus Pyrobaculum, though host specificity might be wide and variable across different "Ca. Aenigmarchaeota" organisms. This study provides significant insight into possible DPANN-host interactions and ecological roles of "Ca. Aenigmarchaeota." IMPORTANCE Recent advances in sequencing technology promoted the blowout discovery of super tiny microbes in the Diapherotrites, Parvarchaeota, Aenigmarchaeota, Nanoarchaeota, and Nanohaloarchaeota (DPANN) superphylum. However, the unculturable properties of the majority of microbes impeded our investigation of their behavior and symbiotic lifestyle in the corresponding community. By integrating horizontal gene transfer (HGT) detection and co-occurrence network analysis on "Candidatus Aenigmarchaeota" ("Ca. Aenigmarchaeota"), we made one of the first attempts to infer their putative interaction partners and further decipher the potential functional and genetic interactions between the symbionts. We revealed that HGTs contributed by members from the Thaumarchaeota, Aigarchaeota, Crenarchaeota, and Korarchaeota (TACK) superphylum and Euryarchaeota conferred "Ca. Aenigmarchaeota" with the ability to survive under different environmental stresses, such as virus infection, high temperature, and oxidative stress. This study demonstrates that the interaction partners might be inferable by applying informatics analyses on metagenomic sequencing data.

Comparative Genomics Provides Insights into the Genetic Diversity and Evolution of the DPANN Superphylum

DPANN is known as highly diverse, globally widespread, and mostly ectosymbiotic archaeal superphylum. However, this group of archaea was overlooked for a long time, and there were limited in-depth studies reported. In this investigation, 41 metagenome-assembled genomes (MAGs) belonging to the DPANN superphylum were recovered (18 MAGs had average nucleotide identity [ANI] values of <95% and a percentage of conserved proteins [POCP] of >50%, while 14 MAGs showed a POCP of <50%), which were analyzed comparatively with 515 other published DPANN genomes. Mismatches to known 16S rRNA gene primers were identified among 16S rRNA genes of DPANN archaea. Numbers of gene families lost (mostly related to energy and amino acid metabolism) were over three times greater than those gained in the evolution of DPANN archaea. Lateral gene transfer (LGT; ∼45.5% was cross-domain) had facilitated niche adaption of the DPANN archaea, ensuring a delicate equilibrium of streamlined genomes with efficient niche-adaptive strategies. For instance, LGT-derived cytochrome bd ubiquinol oxidase and arginine deiminase in the genomes of “Candidatus Micrarchaeota” could help them better adapt to aerobic acidic mine drainage habitats. In addition, most DPANN archaea acquired enzymes for biosynthesis of extracellular polymeric substances (EPS) and transketolase/transaldolase for the pentose phosphate pathway from Bacteria.

IMPORTANCE The domain Archaea is a key research model for gaining insights into the origin and evolution of life, as well as the relevant biogeochemical processes. The discovery of nanosized DPANN archaea has overthrown many aspects of microbiology. However, the DPANN superphylum still contains a vast genetic novelty and diversity that need to be explored. Comprehensively comparative genomic analysis on the DPANN superphylum was performed in this study, with an attempt to illuminate its metabolic potential, ecological distribution and evolutionary history. Many interphylum differences within the DPANN superphylum were found. For example, Altiarchaeota had the biggest genome among DPANN phyla, possessing many pathways missing in other phyla, such as formaldehyde assimilation and the Wood-Ljungdahl pathway. In addition, LGT acted as an important force to provide DPANN archaeal genetic flexibility that permitted the occupation of diverse niches. This study has advanced our understanding of the diversity and genome evolution of archaea.

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.

Respiratory Heme A-Containing Oxidases Originated in the Ancestors of Iron-Oxidizing Bacteria

Respiration is a major trait shaping the biology of many environments. Cytochrome oxidase containing heme A (COX) is a common terminal oxidase in aerobic bacteria and is the only one in mammalian mitochondria. The synthesis of heme A is catalyzed by heme A synthase (CtaA/Cox15), an enzyme that most likely coevolved with COX. The evolutionary origin of COX in bacteria has remained unknown. Using extensive sequence and phylogenetic analysis, we show that the ancestral type of heme A synthases is present in iron-oxidizing Proteobacteria such as Acidithiobacillus spp. These bacteria also contain a deep branching form of the major COX subunit (COX1) and an ancestral variant of CtaG, a protein that is specifically required for COX biogenesis. Our work thus suggests that the ancestors of extant iron-oxidizers were the first to evolve COX. Consistent with this conclusion, acidophilic iron-oxidizing prokaryotes lived on emerged land around the time for which there is the earliest geochemical evidence of aerobic respiration on earth. Hence, ecological niches of iron oxidation have apparently promoted the evolution of aerobic respiration.

Gene Duplications Trace Mitochondria to the Onset of Eukaryote Complexity

The last eukaryote common ancestor (LECA) possessed mitochondria and all key traits that make eukaryotic cells more complex than their prokaryotic ancestors, yet the timing of mitochondrial acquisition and the role of mitochondria in the origin of eukaryote complexity remain debated. Here, we report evidence from gene duplications in LECA indicating an early origin of mitochondria. Among 163,545 duplications in 24,571 gene trees spanning 150 sequenced eukaryotic genomes, we identify 713 gene duplication events that occurred in LECA. LECA's bacterial-derived genes include numerous mitochondrial functions and were duplicated significantly more often than archaeal-derived and eukaryote-specific genes. The surplus of bacterial-derived duplications in LECA most likely reflects the serial copying of genes from the mitochondrial endosymbiont to the archaeal host's chromosomes. Clustering, phylogenies and likelihood ratio tests for 22.4 million genes from 5,655 prokaryotic and 150 eukaryotic genomes reveal no evidence for lineage-specific gene acquisitions in eukaryotes, except from the plastid in the plant lineage. That finding, and the functions of bacterial genes duplicated in LECA, suggests that the bacterial genes in eukaryotes are acquisitions from the mitochondrion, followed by vertical gene evolution and differential loss across eukaryotic lineages, flanked by concomitant lateral gene transfer among prokaryotes. Overall, the data indicate that recurrent gene transfer via the copying of genes from a resident mitochondrial endosymbiont to archaeal host chromosomes preceded the onset of eukaryotic cellular complexity, favoring mitochondria-early over mitochondria-late hypotheses for eukaryote origin.

Glucose Metabolism and Acetate Switch in Archaea: the Enzymes in Haloferax volcanii

The halophilic archaeon Haloferax volcanii has been proposed to degrade glucose via the semiphosphorylative Entner-Doudoroff (spED) pathway. Following our previous studies on key enzymes of this pathway, we now focus on the characterization of enzymes involved in 3-phosphoglycerate conversion to pyruvate, in anaplerosis, and in acetyl coenzyme A (acetyl-CoA) formation from pyruvate. These enzymes include phosphoglycerate mutase, enolase, pyruvate kinase, phosphoenolpyruvate carboxylase, and pyruvate-ferredoxin oxidoreductase. The essential function of these enzymes were shown by transcript analyses and growth experiments with respective deletion mutants. Furthermore, we show that H. volcanii-during aerobic growth on glucose-excreted significant amounts of acetate, which was consumed in the stationary phase (acetate switch). The enzyme catalyzing the conversion of acetyl-CoA to acetate as part of the acetate overflow mechanism, an ADP-forming acetyl-CoA synthetase (ACD), was characterized. The functional involvement of ACD in acetate formation and of AMP-forming acetyl-CoA synthetases (ACSs) in activation of excreted acetate was proven by using respective deletion mutants. Together, the data provide a comprehensive analysis of enzymes of the spED pathway and of anaplerosis and report the first genetic evidence of the functional involvement of enzymes of the acetate switch in archaea.IMPORTANCE In this work, we provide a comprehensive analysis of glucose degradation via the semiphosphorylative Entner-Doudoroff pathway in the haloarchaeal model organism Haloferax volcanii The study includes transcriptional analyses, growth experiments with deletion mutants. and characterization of all enzymes involved in the conversion of 3-phosphoglycerate to acetyl coenzyme A (acetyl-CoA) and in anaplerosis. Phylogenetic analyses of several enzymes indicate various lateral gene transfer events from bacteria to haloarchaea. Furthermore, we analyzed the key players involved in the acetate switch, i.e., in the formation (overflow) and subsequent consumption of acetate during aerobic growth on glucose. Together, the data provide novel aspects of glucose degradation, anaplerosis, and acetate switch in H. volcanii and thus expand our understanding of the unusual sugar metabolism in archaea.

Crystal structure of the 6-phosphogluconate dehydrogenase from Gluconobacter oxydans reveals tetrameric 6PGDHs as the crucial intermediate in the evolution of structure and cofactor preference in the 6PGDH family

Background: The enzyme 6-phosphogluconate dehydrogenase (6PGDH) is the central enzyme of the oxidative pentose phosphate pathway. Members of the 6PGDH family belong to different classes: either homodimeric enzymes assembled from long-chain subunits or homotetrameric ones assembled from short-chain subunits. Dimeric 6PGDHs bear an internal duplication absent in tetrameric 6PGDHs and distant homologues of the β-hydroxyacid dehydrogenase (βHADH) superfamily. Methods: We use X-ray crystallography to determine the structure of the apo form of the 6PGDH from Gluconobacter oxydans (Go6PGDH). We carried out a structural and phylogenetic analysis of short and long-chain 6PGDHs. We put forward an evolutionary hypothesis explaining the differences seen in oligomeric state vs. dinucleotide preference of the 6PGDH family. We determined the cofactor preference of Go6PGDH at different 6-phosphogluconate concentrations, characterizing the wild-type enzyme and three-point mutants of residues in the cofactor binding site of Go6PGDH. Results: The structural comparison suggests that the 6PG binding site initially evolved by exchanging C-terminal α-helices between subunits. An internal duplication event changed the quaternary structure of the enzyme from a tetrameric to a dimeric arrangement. The phylogenetic analysis suggests that 6PGDHs have spread from Bacteria to Archaea and Eukarya on multiple occasions by lateral gene transfer. Sequence motifs consistent with NAD+- and NADP+-specificity are found in the β2-α2 loop of dimeric and tetrameric 6PGDHs. Site-directed mutagenesis of Go6PGDH inspired by this analysis fully reverses dinucleotide preference. One of the mutants we engineered has the highest efficiency and specificity for NAD+ so far described for a 6PGDH. Conclusions: The family 6PGDH comprises dimeric and tetrameric members whose active sites are conformed by a C-terminal α-helix contributed from adjacent subunits. Dimeric 6PGDHs have evolved from the duplication-fusion of the tetrameric C-terminal domain before independent transitions of cofactor specificity. Changes in the conserved β2-α2 loop are crucial to modulate the cofactor specificity in Go6PGDH.

Metagenomic insights into the metabolism and evolution of a new Thermoplasmata order (Candidatus Gimiplasmatales)

Thermoplasmata is a widely distributed and ecologically important archaeal class in the phylum Euryarchaeota. Because few cultures and genomes are available, uncharacterized Thermoplasmata metabolisms remain unexplored. In this study, we obtained four medium- to high-quality archaeal metagenome-assembled genomes (MAGs) from the filamentous fragments of black-odorous aquatic sediments (Foshan, Guangdong, China). Based on their 16S rRNA gene and ribosomal protein phylogenies, the four MAGs belong to the previously unnamed Thermoplasmata UBA10834 clade. We propose that this clade (five reference genomes from the Genome Taxonomy Database (GTDB) and four MAGs from this study) be considered a new order, Candidatus Gimiplasmatales. Metabolic pathway reconstructions indicated that the Ca. Gimiplasmatales MAGs can biosynthesize isoprenoids and nucleotides de novo. Additionally, some taxa have genes for formaldehyde and acetate assimilation, and the Wood-Ljungdahl CO₂ -fixation pathway, indicating a mixotrophic lifestyle. Sulfur reduction, hydrogen metabolism, and arsenic detoxification pathways were predicted, indicating sulfur-, hydrogen-, and arsenic-transformation potentials. Comparative genomics indicated that the H₄ F Wood-Ljungdahl pathway of both Ca. Gimiplasmatales and Methanomassiliicoccales was likely obtained by the interdomain lateral gene transfer from the Firmicutes. Collectively, this study elucidates the taxonomic and potential metabolic diversity of the new order Ca. Gimiplasmatales and the evolution of this subgroup and its sister lineage Methanomassiliicoccales.

Pan-genome analysis and ancestral state reconstruction of class halobacteria: probability of a new super-order

Halobacteria, a class of Euryarchaeota are extremely halophilic archaea that can adapt to a wide range of salt concentration generally from 10% NaCl to saturated salt concentration of 32% NaCl. It consists of the orders: Halobacteriales, Haloferaciales and Natriabales. Pan-genome analysis of class Halobacteria was done to explore the core (300) and variable components (Softcore: 998, Cloud:36531, Shell:11784). The core component revealed genes of replication, transcription, translation and repair, whereas the variable component had a major portion of environmental information processing. The pan-gene matrix was mapped onto the core-gene tree to find the ancestral (44.8%) and derived genes (55.1%) of the Last Common Ancestor of Halobacteria. A High percentage of derived genes along with presence of transformation and conjugation genes indicate the occurrence of horizontal gene transfer during the evolution of Halobacteria. A Core and pan-gene tree were also constructed to infer a phylogeny which implicated on the new super-order comprising of Natrialbales and Halobacteriales.

Anomalous Phylogenetic Behavior of Ribosomal Proteins in Metagenome-Assembled Asgard Archaea

Metagenomic studies permit the exploration of microbial diversity in a defined habitat, and binning procedures enable phylogenomic analyses, taxon description, and even phenotypic characterizations in the absence of morphological evidence. Such lineages include asgard archaea, which were initially reported to represent archaea with eukaryotic cell complexity, although the first images of such an archaeon show simple cells with prokaryotic characteristics. However, these metagenome-assembled genomes (MAGs) might suffer from data quality problems not encountered in sequences from cultured organisms due to two common analytical procedures of bioinformatics: assembly of metagenomic sequences and binning of assembled sequences on the basis of innate sequence properties and abundance across samples. Consequently, genomic sequences of distantly related taxa, or domains, can in principle be assigned to the same MAG and result in chimeric sequences. The impacts of low-quality or chimeric MAGs on phylogenomic and metabolic prediction remain unknown. Debates that asgard archaeal data are contaminated with eukaryotic sequences are overshadowed by the lack of evidence indicating that individual asgard MAGs stem from the same chromosome. Here, we show that universal proteins including ribosomal proteins of asgard archaeal MAGs fail to meet the basic phylogenetic criterion fulfilled by genome sequences of cultured archaea investigated to date: These proteins do not share common evolutionary histories to the same extent as pure culture genomes do, pointing to a chimeric nature of asgard archaeal MAGs. Our analysis suggests that some asgard archaeal MAGs represent unnatural constructs, genome-like patchworks of genes resulting from assembly and/or the binning process.

A spectrum of verticality across genes

Lateral gene transfer (LGT) has impacted prokaryotic genome evolution, yet the extent to which LGT compromises vertical evolution across individual genes and individual phyla is unknown, as are the factors that govern LGT frequency across genes. Estimating LGT frequency from tree comparisons is problematic when thousands of genomes are compared, because LGT becomes difficult to distinguish from phylogenetic artefacts. Here we report quantitative estimates for verticality across all genes and genomes, leveraging a well-known property of phylogenetic inference: phylogeny works best at the tips of trees. From terminal (tip) phylum level relationships, we calculate the verticality for 19,050,992 genes from 101,422 clusters in 5,655 prokaryotic genomes and rank them by their verticality. Among functional classes, translation, followed by nucleotide and cofactor biosynthesis, and DNA replication and repair are the most vertical. The most vertically evolving lineages are those rich in ecological specialists such as Acidithiobacilli, Chlamydiae, Chlorobi and Methanococcales. Lineages most affected by LGT are the α-, β-, γ-, and δ- classes of Proteobacteria and the Firmicutes. The 2,587 eukaryotic clusters in our sample having prokaryotic homologues fail to reject eukaryotic monophyly using the likelihood ratio test. The low verticality of α-proteobacterial and cyanobacterial genomes requires only three partners-an archaeal host, a mitochondrial symbiont, and a plastid ancestor-each with mosaic chromosomes, to directly account for the prokaryotic origin of eukaryotic genes. In terms of phylogeny, the 100 most vertically evolving prokaryotic genes are neither representative nor predictive for the remaining 97% of an average genome. In search of factors that govern LGT frequency, we find a simple but natural principle: Verticality correlates strongly with gene distribution density, LGT being least likely for intruding genes that must replace a preexisting homologue in recipient chromosomes. LGT is most likely for novel genetic material, intruding genes that encounter no competing copy.

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.

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.

The controversy on the ancestral arsenite oxidizing enzyme; deducing evolutionary histories with phylogeny and thermodynamics

The three presently known enzymes responsible for arsenic-using bioenergetic processes are arsenite oxidase (Aio), arsenate reductase (Arr) and alternative arsenite oxidase (Arx), all of which are molybdoenzymes from the vast group referred to as the Mo/W-bisPGD enzyme superfamily. Since arsenite is present in substantial amounts in hydrothermal environments, frequently considered as vestiges of primordial biochemistry, arsenite-based bioenergetics has long been predicted to be ancient. Conflicting scenarios, however, have been put forward proposing either Arr/Arx or Aio as operating in the ancestral metabolism. Phylogenetic data argue in favor of Aio whereas biochemical and physiological data led several authors to propose Arx/Arr as the most ancient anaerobic arsenite metabolizing enzymes. Here we combine phylogenetic approaches with physiological and biochemical experiments to demonstrate that the Arx/Arr enzymes could not have been functional in the Archaean geological eon. We propose that Arr reacts with menaquinones to reduce arsenate whereas Arx reacts with ubiquinone to oxidize arsenite, in line with thermodynamic considerations. The distribution of the quinone biosynthesis pathways, however, clearly indicates that the ubiquinone pathway is recent. An updated phylogeny of Arx furthermore reinforces the hypothesis of a recent emergence of this enzyme. We therefore conclude that anaerobic arsenite redox conversion in the Archaean must have been performed in a metabolism involving Aio.

On the evolution of cytochrome oxidases consuming oxygen

This review examines the current state of the art on the evolution of the families of Heme Copper Oxygen reductases (HCO) that oxidize cytochrome c and reduce oxygen to water, chiefly cytochrome oxidase, COX. COX is present in many bacterial and most eukaryotic lineages, but its origin has remained elusive. After examining previous proposals for COX evolution, the review summarizes recent insights suggesting that COX enzymes might have evolved in soil dwelling, probably iron-oxidizing bacteria which lived on emerged land over two billion years ago. These bacteria were the likely ancestors of extant acidophilic iron-oxidizers such as Acidithiobacillus spp., which belong to basal lineages of the phylum Proteobacteria. Proteobacteria may thus be considered the originators of COX, which was then laterally transferred to other prokaryotes. The taxonomy of bacteria is presented in relation to the current distribution of COX and C family oxidases, from which COX may have evolved.

A New Phylogenomic Approach For Quantifying Horizontal Gene Transfer Trends in Prokaryotes

It is well established nowadays that among prokaryotes, various families of orthologous genes exhibit conflicting evolutionary history. A prime factor for this conflict is horizontal gene transfer (HGT) - the transfer of genetic material not via vertical descent. Thus, the prevalence of HGT is challenging the meaningfulness of the classical Tree of Life concept. Here we present a comprehensive study of HGT representing the entire prokaryotic world. We mainly rely on a novel analytic approach for analyzing an aggregate of gene histories, by means of the quartet plurality distribution (QPD) that we develop. Through the analysis of real and simulated data, QPD is used to reveal evidence of a barrier against HGT, separating the archaea from the bacteria and making HGT between the two domains, in general, quite rare. In contrast, bacteria's confined HGT is substantially more frequent than archaea's. Our approach also reveals that despite intensive HGT, a strong tree-like signal can be extracted, corroborating several previous works. Thus, QPD, which enables one to analytically combine information from an aggregate of gene trees, can be used for understanding patterns and rates of HGT in prokaryotes, as well as for validating or refuting models of horizontal genetic transfers and evolution in general.

Horizontal Gene Transfer as a Source of Conflict and Cooperation in Prokaryotes

Horizontal gene transfer (HGT) is one of the most important processes in prokaryote evolution. The sharing of DNA can spread neutral or beneficial genes, as well as genetic parasites across populations and communities, creating a large proportion of the variability acted on by natural selection. Here, we highlight the role of HGT in enhancing the opportunities for conflict and cooperation within and between prokaryote genomes. We discuss how horizontally acquired genes can cooperate or conflict both with each other and with a recipient genome, resulting in signature patterns of gene co-occurrence, avoidance, and dependence. We then describe how interactions involving horizontally transferred genes may influence cooperation and conflict at higher levels (populations, communities, and symbioses). Finally, we consider the benefits and drawbacks of HGT for prokaryotes and its fundamental role in understanding conflict and cooperation from the gene-gene to the microbiome level.

Descent of Bacteria and Eukarya From an Archaeal Root of Life

The 3 biological domains delineated based on small subunit ribosomal RNAs (SSU rRNAs) are confronted by uncertainties regarding the relationship between Archaea and Bacteria, and the origin of Eukarya. The similarities between the paralogous valyl-tRNA and isoleucyl-tRNA synthetases in 5398 species estimated by BLASTP, which decreased from Archaea to Bacteria and further to Eukarya, were consistent with vertical gene transmission from an archaeal root of life close to Methanopyrus kandleri through a Primitive Archaea Cluster to an Ancestral Bacteria Cluster, and to Eukarya. The predominant similarities of the ribosomal proteins (rProts) of eukaryotes toward archaeal rProts relative to bacterial rProts established that an archaeal parent rather than a bacterial parent underwent genome merger with bacteria to generate eukaryotes with mitochondria. Eukaryogenesis benefited from the predominantly archaeal accelerated gene adoption (AGA) phenotype pertaining to horizontally transferred genes from other prokaryotes and expedited genome evolution via both gene-content mutations and nucleotidyl mutations. Archaeons endowed with substantial AGA activity were accordingly favored as candidate archaeal parents. Based on the top similarity bitscores displayed by their proteomes toward the eukaryotic proteomes of Giardia and Trichomonas, and high AGA activity, the Aciduliprofundum archaea were identified as leading candidates of the archaeal parent. The Asgard archaeons and a number of bacterial species were among the foremost potential contributors of eukaryotic-like proteins to Eukarya.

Older Than Genes: The Acetyl CoA Pathway and Origins

For decades, microbiologists have viewed the acetyl CoA pathway and organisms that use it for H2-dependent carbon and energy metabolism, acetogens and methanogens, as ancient. Classical evidence and newer evidence indicating the antiquity of the acetyl CoA pathway are summarized here. The acetyl CoA pathway requires approximately 10 enzymes, roughly as many organic cofactors, and more than 500 kDa of combined subunit molecular mass to catalyze the conversion of H2 and CO2 to formate, acetate, and pyruvate in acetogens and methanogens. However, a single hydrothermal vent alloy, awaruite (Ni3Fe), can convert H2 and CO2 to formate, acetate, and pyruvate under mild hydrothermal conditions on its own. The chemical reactions of H2 and CO2 to pyruvate thus have a natural tendency to occur without enzymes, given suitable inorganic catalysts. This suggests that the evolution of the enzymatic acetyl CoA pathway was preceded by-and patterned along-a route of naturally occurring exergonic reactions catalyzed by transition metal minerals that could activate H2 and CO2 by chemisorption. The principle of forward (autotrophic) pathway evolution from preexisting non-enzymatic reactions is generalized to the concept of patterned evolution of pathways. In acetogens, exergonic reduction of CO2 by H2 generates acyl phosphates by highly reactive carbonyl groups undergoing attack by inert inorganic phosphate. In that ancient reaction of biochemical energy conservation, the energy behind formation of the acyl phosphate bond resides in the carbonyl, not in phosphate. The antiquity of the acetyl CoA pathway is usually seen in light of CO2 fixation; its role in primordial energy coupling via acyl phosphates and substrate-level phosphorylation is emphasized here.

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.

Acid Experimental Evolution of the Haloarchaeon Halobacterium sp. NRC-1 Selects Mutations Affecting Arginine Transport and Catabolism

Halobacterium sp. NRC-1 (NRC-1) is an extremely halophilic archaeon that is adapted to multiple stressors such as UV, ionizing radiation and arsenic exposure; it is considered a model organism for the feasibility of microbial life in iron-rich brine on Mars. We conducted experimental evolution of NRC-1 under acid and iron stress. NRC-1 was serially cultured in CM+ medium modified by four conditions: optimal pH (pH 7.5), acid stress (pH 6.3), iron amendment (600 μM ferrous sulfate, pH 7.5), and acid plus iron (pH 6.3, with 600 μM ferrous sulfate). For each condition, four independent lineages of evolving populations were propagated. After 500 generations, 16 clones were isolated for phenotypic characterization and genomic sequencing. Genome sequences of all 16 clones revealed 378 mutations, of which 90% were haloarchaeal insertion sequences (ISH) and ISH-mediated large deletions. This proportion of ISH events in NRC-1 was five-fold greater than that reported for comparable evolution of Escherichia coli. One acid-evolved clone had increased fitness compared to the ancestral strain when cultured at low pH. Seven of eight acid-evolved clones had a mutation within or upstream of arcD, which encodes an arginine-ornithine antiporter; no non-acid adapted strains had arcD mutations. Mutations also affected the arcR regulator of arginine catabolism, which protects bacteria from acid stress by release of ammonia. Two acid-adapted strains shared a common mutation in bop, which encodes bacterio-opsin, apoprotein for the bacteriorhodopsin light-driven proton pump. Thus, in the haloarchaeon NRC-1, as in bacteria, pH adaptation was associated with genes involved in arginine catabolism and proton transport. Our study is among the first to report experimental evolution with multiple resequenced genomes of an archaeon. Haloarchaea are polyextremophiles capable of growth under environmental conditions such as concentrated NaCl and desiccation, but little is known about pH stress. Interesting parallels appear between the molecular basis of pH adaptation in NRC-1 and in bacteria, particularly the acid-responsive arginine-ornithine system found in oral streptococci.

On the Origin of Isoprenoid Biosynthesis

Isoprenoids and their derivatives represent the largest group of organic compounds in nature and are distributed universally in the three domains of life. Isoprenoids are biosynthesized from isoprenyl diphosphate units, generated by two distinctive biosynthetic pathways: mevalonate pathway and methylerthritol 4-phosphate pathway. Archaea and eukaryotes exclusively have the former pathway, while most bacteria have the latter. Some bacteria, however, are known to possess the mevalonate pathway genes. Understanding the evolutionary history of these two isoprenoid biosynthesis pathways in each domain of life is critical since isoprenoids are so interweaved in the architecture of life that they would have had indispensable roles in the early evolution of life. Our study provides a detailed phylogenetic analysis of enzymes involved in the mevalonate pathway and sheds new light on its evolutionary history. The results suggest that a potential mevalonate pathway is present in the recently discovered superphylum Candidate Phyla Radiation (CPR), and further suggest a strong evolutionary relationship exists between archaea and CPR. Interestingly, CPR harbors the characteristics of both the bacterial-type and archaeal-type mevalonate pathways and may retain signatures regarding the ancestral isoprenoid biosynthesis pathway in the last universal common ancestor. Our study supports the ancient origin of the mevalonate pathway in the three domains of life as previously inferred, but concludes that the evolution of the mevalonate pathway was more complex.

Networks Consolidate the Core Concepts of Evolution by Natural Selection

Microbiology has unraveled rich evidence of ongoing reticulate evolutionary processes and complex interactions both within and between cells. These phenomena feature real biological networks, which can logically be analyzed using network-based tools. It is thus not surprising that network sciences, a field independent from evolutionary biology and microbiology, have recently pervasively infused their methods into both fields. Importantly, network tools bring forward observations enhancing the understanding of three core evolutionary concepts: variation, fitness, and heredity. Consequently, our work shows how network sciences can enhance evolutionary theory by explaining the evolution by natural selection of a broad diversity of units of selection, while updating the popular figure of Darwin's tree of life with a comprehensive sketch of the networks of evolution.

Eukaryote Genes Are More Likely than Prokaryote Genes to Be Composites

The formation of new genes by combining parts of existing genes is an important evolutionary process. Remodelled genes, which we call composites, have been investigated in many species, however, their distribution across all of life is still unknown. We set out to examine the extent to which genomes from cells and mobile genetic elements contain composite genes. We identify composite genes as those that show partial homology to at least two unrelated component genes. In order to identify composite and component genes, we constructed sequence similarity networks (SSNs) of more than one million genes from all three domains of life, as well as viruses and plasmids. We identified non-transitive triplets of nodes in this network and explored the homology relationships in these triplets to see if the middle nodes were indeed composite genes. In total, we identified 221,043 (18.57%) composites genes, which were distributed across all genomic and functional categories. In particular, the presence of composite genes is statistically more likely in eukaryotes than prokaryotes.

The genomes of polyextremophilic cyanidiales contain 1% horizontally transferred genes with diverse adaptive functions

The role and extent of horizontal gene transfer (HGT) in eukaryotes are hotly disputed topics that impact our understanding of the origin of metabolic processes and the role of organelles in cellular evolution. We addressed this issue by analyzing 10 novel Cyanidiales genomes and determined that 1% of their gene inventory is HGT-derived. Numerous HGT candidates share a close phylogenetic relationship with prokaryotes that live in similar habitats as the Cyanidiales and encode functions related to polyextremophily. HGT candidates differ from native genes in GC-content, number of splice sites, and gene expression. HGT candidates are more prone to loss, which may explain the absence of a eukaryotic pan-genome. Therefore, the lack of a pan-genome and cumulative effects fail to provide substantive arguments against our hypothesis of recurring HGT followed by differential loss in eukaryotes. The maintenance of 1% HGTs, even under selection for genome reduction, underlines the importance of non-endosymbiosis related foreign gene acquisition.

New views on an old enzyme: allosteric regulation and evolution of archaeal pyruvate kinases

Pyruvate kinases (PKs) synthesize ATP as the final step of glycolysis in the three domains of life. PKs from most bacteria and eukarya are allosteric enzymes that are activated by sugar phosphates; for example, the feed-forward regulator fructose-1,6-bisphosphate, or AMP as a sensor of energy charge. Archaea utilize unusual glycolytic pathways, but the allosteric properties of PKs from these species are largely unknown. Here, we present an analysis of 24 PKs from most archaeal clades with respect to allosteric properties, together with phylogenetic analyses constructed using a novel mode of rooting protein trees. We find that PKs from many Thermoproteales, an order of crenarchaeota, are allosterically activated by 3-phosphoglycerate (3PG). We also identify five conserved amino acids that form the binding pocket for 3PG. 3PG is generated via an irreversible reaction in the modified glycolytic pathway of these archaea and therefore functions as a feed-forward regulator. We also show that PKs from hyperthermophilic Methanococcales, an order of euryarchaeota, are activated by AMP. Phylogenetic analyses indicate that 3PG-activated PKs form an evolutionary lineage that is distinct from that of sugar-phosphate activated PKs, and that sugar phosphate-activated PKs originated as AMP-regulated PKs in hyperthermophilic Methanococcales. Since the phospho group of sugar phosphates and 3PG overlap in the allosteric site, our data indicate that the allostery in PKs first started from a progenitor phosphate-binding site that evolved in two spatially distinct directions: one direction generated the canonical site that responds to sugar phosphates and the other gave rise to the 3PG site present in Thermoproteales. Overall, our data suggest an intimate connection between the allosteric properties and evolution of PKs.

Why the Lipid Divide? Membrane Proteins as Drivers of the Split between the Lipids of the Three Domains of Life

Recent results from engineered and natural samples show that the starkly different lipids of archaea and bacteria can form stable hybrid membranes. But if the two types can mix, why don't they? That is, why do most bacteria and all eukaryotes have only typically bacterial lipids, and archaea archaeal lipids? It is suggested here that the reason may lie on the other main component of cellular membranes: membrane proteins, and their close adaptation to the lipids. Archaeal lipids in modern bacteria could suggest that the last universal common ancestor (LUCA) had both lipid types. However, this would imply a rather elaborate evolutionary scenario, while negating simpler alternatives. In light of widespread horizontal gene transfer across the prokaryotic domains, hybrid membranes reveal that the lipid divide did not just occur once at the divergence of archaea and bacteria from LUCA. Instead, it continues to occur actively to this day.

Oxygen Reductases in Alphaproteobacterial Genomes: Physiological Evolution From Low to High Oxygen Environments

Oxygen reducing terminal oxidases differ with respect to their subunit composition, heme groups, operon structure, and affinity for O₂. Six families of terminal oxidases are currently recognized, all of which occur in alphaproteobacterial genomes, two of which are also present in mitochondria. Many alphaproteobacteria encode several different terminal oxidases, likely reflecting ecological versatility with respect to oxygen levels. Terminal oxidase evolution likely started with the advent of O₂ roughly 2.4 billion years ago and terminal oxidases diversified in the Proterozoic, during which oxygen levels remained low, around the Pasteur point (ca. 2 μM O₂). Among the alphaproteobacterial genomes surveyed, those from members of the Rhodospirillaceae reveal the greatest diversity in oxygen reductases. Some harbor all six terminal oxidase types, in addition to many soluble enzymes typical of anaerobic fermentations in mitochondria and hydrogenosomes of eukaryotes. Recent data have it that O₂ levels increased to current values (21% v/v or ca. 250 μM) only about 430 million years ago. Ecological adaptation brought forth different lineages of alphaproteobacteria and different lineages of eukaryotes that have undergone evolutionary specialization to high oxygen, low oxygen, and anaerobic habitats. Some have remained facultative anaerobes that are able to generate ATP with or without the help of oxygen and represent physiological links to the ancient proteobacterial lineage at the origin of mitochondria and eukaryotes. Our analysis reveals that the genomes of alphaproteobacteria appear to retain signatures of ancient transitions in aerobic metabolism, findings that are relevant to mitochondrial evolution in eukaryotes as well.

Bipartite Network Analysis of Gene Sharings in the Microbial World

Extensive microbial gene flows affect how we understand virology, microbiology, medical sciences, genetic modification, and evolutionary biology. Phylogenies only provide a narrow view of these gene flows: plasmids and viruses, lacking core genes, cannot be attached to cellular life on phylogenetic trees. Yet viruses and plasmids have a major impact on cellular evolution, affecting both the gene content and the dynamics of microbial communities. Using bipartite graphs that connect up to 149,000 clusters of homologous genes with 8,217 related and unrelated genomes, we can in particular show patterns of gene sharing that do not map neatly with the organismal phylogeny. Homologous genes are recycled by lateral gene transfer, and multiple copies of homologous genes are carried by otherwise completely unrelated (and possibly nested) genomes, that is, viruses, plasmids and prokaryotes. When a homologous gene is present on at least one plasmid or virus and at least one chromosome, a process of "gene externalization," affected by a postprocessed selected functional bias, takes place, especially in Bacteria. Bipartite graphs give us a view of vertical and horizontal gene flow beyond classic taxonomy on a single very large, analytically tractable, graph that goes beyond the cellular Web of Life.

Bioinformatics Analysis of Metabolism Pathways of Archaeal Energy Reserves

Energy storage compounds play crucial roles in prokaryotic physiology. Five chemical compounds have been identified in prokaryotes as energy reserves: polyphosphate (polyP), polyhydroxyalkanoates (PHAs), glycogen, wax ester (WE) and triacylglycerol (TAG). Currently, no systematic study of archaeal energy storage metabolism exists. In this study, we collected 427 archaeal reference sequences from UniProt database. A thorough pathway screening of energy reserves led to an overview of distribution patterns of energy metabolism in archaea. We also explored how energy metabolism might have impact on archaeal extremophilic phenotypes. Based on the systematic analyses of archaeal proteomes, we confirmed that metabolism pathways of polyP, PHAs and glycogen are present in archaea, but TAG and WE are completely absent. It was also confirmed that PHAs are tightly related to halophilic archaea with larger proteome size and higher GC contents, while polyP is mainly present in methanogens. In sum, this study systematically investigates energy storage metabolism in archaea and provides a clear correlation between energy metabolism and the ability to survive in extreme environments. With more genomic editing tools developed for archaea and molecular mechanisms unravelled for energy storage metabolisms (ESMs), there will be a better understanding of the unique lifestyle of archaea in extreme environments.