Considerations for Detecting Organic Indicators of Metabolism on Enceladus.

Enceladus is of interest to astrobiology and the search for life since it is thought to host active hydrothermal activity and habitable conditions. It is also possible that the organics detected on Enceladus may indicate an active prebiotic or biotic system; in particular, the conditions on Enceladus may favor mineral-driven protometabolic reactions. When including metabolism-related biosignatures in Enceladus mission concepts, it is necessary to base these in a clearer understanding of how these signatures could also be produced prebiotically. In addition, postulating which biological metabolisms to look for on Enceladus requires a non-Earth-centric approach since the details of biological metabolic pathways are heavily shaped by adaptation to geochemical conditions over the planet's history. Creating metabolism-related organic detection objectives for Enceladus missions, therefore, requires consideration of how metabolic systems may operate differently on another world, while basing these speculations on observed Earth-specific microbial processes. In addition, advances in origin-of-life research can play a critical role in distinguishing between interpretations of any future organic detections on Enceladus, and the discovery of an extant prebiotic system would be a transformative astrobiological event in its own right.

Assessing the Influence of HGT on the Evolution of Stress Responses in Microbial Communities from Shark Bay, Western Australia.

Pustular microbial mats in Shark Bay, Western Australia, are modern analogs of microbial systems that colonized peritidal environments before the evolution of complex life. To understand how these microbial communities evolved to grow and metabolize in the presence of various environmental stresses, the horizontal gene transfer (HGT) detection tool, MetaCHIP, was used to identify the horizontal transfer of genes related to stress response in 83 metagenome-assembled genomes from a Shark Bay pustular mat. Subsequently, maximum-likelihood phylogenies were constructed using these genes and their most closely related homologs from other environments in order to determine the likelihood of these HGT events occurring within the pustular mat. Phylogenies of several stress-related genes-including those involved in response to osmotic stress, oxidative stress and arsenic toxicity-indicate a potentially long history of HGT events and are consistent with these transfers occurring outside of modern pustular mats. The phylogeny of a particular osmoprotectant transport gene reveals relatively recent adaptations and suggests interactions between Planctomycetota and Myxococcota within these pustular mats. Overall, HGT phylogenies support a potentially broad distribution in the relative timing of the HGT events of stress-related genes and demonstrate ongoing microbial adaptations and evolution in these pustular mat communities.

Fast-Evolving Alignment Sites Are Highly Informative for Reconstructions of Deep Tree of Life Phylogenies.

The trimming of fast-evolving sites, often known as "slow-fast" analysis, is broadly used in microbial phylogenetic reconstruction under the assumption that fast-evolving sites do not retain an accurate phylogenetic signal due to substitution saturation. Therefore, removing sites that have experienced multiple substitutions would improve the signal-to-noise ratio in phylogenetic analyses, with the remaining slower-evolving sites preserving a more reliable record of evolutionary relationships. Here, we show that, contrary to this assumption, even the fastest-evolving sites present in the conserved proteins often used in Tree of Life studies contain reliable and valuable phylogenetic information, and that the trimming of such sites can negatively impact the accuracy of phylogenetic reconstruction. Simulated alignments modeled after ribosomal protein datasets used in Tree of Life studies consistently show that slow-evolving sites are less likely to recover true bipartitions than even the fastest-evolving sites. Furthermore, site-specific substitution rates are positively correlated with the frequency of accurately recovered short-branched bipartitions, as slowly evolving sites are less likely to have experienced substitutions along these intervals. Using published Tree of Life sequence alignment datasets, we also show that both slow- and fast-evolving sites contain similarly inconsistent phylogenetic signals, and that, for fast-evolving sites, this inconsistency can be attributed to poor alignment quality. Furthermore, trimming fast sites, slow sites, or both is shown to have a substantial impact on phylogenetic reconstruction across multiple evolutionary models. This is perhaps most evident in the resulting placements of the Eukarya and Asgardarchaeota groups, which are especially sensitive to the implementation of different trimming schemes.

Late acquisition of the rTCA carbon fixation pathway by Chlorobi.

The reverse tricarboxylic acid (rTCA) cycle is touted as a primordial mode of carbon fixation due to its autocatalytic propensity and oxygen intolerance. Despite this inferred antiquity, however, the earliest rock record affords scant supporting evidence. In fact, based on the chimeric inheritance of rTCA cycle steps within the Chlorobiaceae, even the use of the chemical fossil record of this group is now subject to question. While the 1.64-billion-year-old Barney Creek Formation contains chemical fossils of the earliest known putative Chlorobiaceae-derived carotenoids, interferences from the accompanying hydrocarbon matrix have hitherto precluded the carbon isotope measurements necessary to establish the physiology of the organisms that produced them. Overcoming this obstacle, here we report a suite of compound-specific carbon isotope measurements identifying a cyanobacterially dominated ecosystem featuring heterotrophic bacteria. We demonstrate chlorobactane is 13C-depleted when compared to contemporary equivalents, showing only slight 13C-enrichment over co-existing cyanobacterial carotenoids. The absence of this diagnostic isotopic fingerprint, in turn, confirms phylogenomic hypotheses that call for the late assembly of the rTCA cycle and, thus, the delayed acquisition of autotrophy within the Chlorobiaceae. We suggest that progressive oxygenation of the Earth System caused an increase in the marine sulfate inventory thereby providing the selective pressure to fuel the Neoproterozoic shift towards energy-efficient photoautotrophy within the Chlorobiaceae.

Metaproteogenomic Profile of a Mesopelagic Adenylylsulfate Reductase: Course-Based Discovery Using the Ocean Protein Portal.

Adenylylsulfate reductase (Apr) is a flavoprotein with a dissimilatory sulfate reductase function. Its ability to catalyze the reverse reaction in sulfur oxidizers has propelled a complex phylogenetic history of transfers with sulfate reducers and made this enzyme an important protein in ocean sulfur cycling. As part of a graduate course, we analyzed metaproteomic data from the Ocean Protein Portal and observed evidence of Apr alpha (AprA) and beta (AprB) subunits in the Central Pacific Ocean. The protein was originally taxonomically attributed toChlorobium tepidum TLS, a green sulfur bacterium. However, our phylogenomic and oceanographic contextual analysis contradicted this label, instead showing that this protein is consistent with the genomic material from the newly discovered Candidatus Lambdaproteobacteriaclass, implying that the ecological role of this lineage in oxygen minimum twilight zones is underappreciated. This study illustrates how metaproteogenomic analysis can contribute to more accurate metagenomic/proteomic annotations and comprehensive ocean biogeochemical processes conducive to course-based research experiences.

Chitin utilization by marine picocyanobacteria and the evolution of a planktonic lifestyle.

Marine picocyanobacteria Prochlorococcus and Synechococcus, the most abundant photosynthetic cells in the oceans, are generally thought to have a primarily single-celled and free-living lifestyle. However, while studying the ability of picocyanobacteria to supplement photosynthetic carbon fixation with the use of exogenous organic carbon, we found the widespread occurrence of genes for breaking down chitin, an abundant source of organic carbon that exists primarily as particles. We show that cells that encode a chitin degradation pathway display chitin degradation activity, attach to chitin particles, and show enhanced growth under low light conditions when exposed to chitosan, a partially deacetylated soluble form of chitin. Marine chitin is largely derived from arthropods, which underwent major diversifications 520 to 535 Mya, close to when marine picocyanobacteria are inferred to have appeared in the ocean. Phylogenetic analyses confirm that the chitin utilization trait was acquired at the root of marine picocyanobacteria. Together this leads us to postulate that attachment to chitin particles allowed benthic cyanobacteria to emulate their mat-based lifestyle in the water column, initiating their expansion into the open ocean, seeding the rise of modern marine ecosystems. Subsequently, transitioning to a constitutive planktonic life without chitin associations led to cellular and genomic streamlining along a major early branch within Prochlorococcus. Our work highlights how the emergence of associations between organisms from different trophic levels, and their coevolution, creates opportunities for colonizing new environments. In this view, the rise of ecological complexity and the expansion of the biosphere are deeply intertwined processes.

A Proterozoic microbial origin of extant cyanide-hydrolyzing enzyme diversity.

In addition to its role as a toxic environmental contaminant, cyanide has been hypothesized to play a key role in prebiotic chemistry and early biogeochemical evolution. While cyanide-hydrolyzing enzymes have been studied and engineered for bioremediation, the extant diversity of these enzymes remains underexplored. Additionally, the age and evolution of microbial cyanide metabolisms is poorly constrained. Here we provide comprehensive phylogenetic and molecular clock analyses of the distribution and evolution of the Class I nitrilases, thiocyanate hydrolases, and nitrile hydratases. Molecular clock analyses indicate that bacterial cyanide-reducing nitrilases were present by the Paleo- to Mesoproterozoic, and were subsequently horizontally transferred into eukaryotes. These results present a broad diversity of microbial enzymes that could be optimized for cyanide bioremediation.

DaTeR: error-correcting phylogenetic chronograms using relative time constraints.

MOTIVATION: A chronogram is a dated phylogenetic tree whose branch lengths have been scaled to represent time. Such chronograms are computed based on available date estimates (e.g. from dated fossils), which provide absolute time constraints for one or more nodes of an input undated phylogeny, coupled with an appropriate underlying model for evolutionary rates variation along the branches of the phylogeny. However, traditional methods for phylogenetic dating cannot take into account relative time constraints, such as those provided by inferred horizontal transfer events. In many cases, chronograms computed using only absolute time constraints are inconsistent with known relative time constraints. RESULTS: In this work, we introduce a new approach, Dating Trees using Relative constraints (DaTeR), for phylogenetic dating that can take into account both absolute and relative time constraints. The key idea is to use existing Bayesian approaches for phylogenetic dating to sample posterior chronograms satisfying desired absolute time constraints, minimally adjust or 'error-correct' these sampled chronograms to satisfy all given relative time constraints, and aggregate across all error-corrected chronograms. DaTeR uses a constrained optimization framework for the error-correction step, finding minimal deviations from previously assigned dates or branch lengths. We applied DaTeR to a biological dataset of 170 Cyanobacterial taxa and a reliable set of 24 transfer-based relative constraints, under six different molecular dating models. Our extensive analysis of this dataset demonstrates that DaTeR is both highly effective and scalable and that its application can significantly improve estimated chronograms. AVAILABILITY AND IMPLEMENTATION: Freely available from https://compbio.engr.uconn.edu/software/dater/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Chimeric inheritance and crown-group acquisitions of carbon fixation genes within Chlorobiales: Origins of autotrophy in Chlorobiales and implication for geological biomarkers.

The geological record of microbial metabolisms and ecologies primarily consists of stable isotope fractionations and the diagenetic products of biogenic lipids. Carotenoid lipid biomarkers are particularly useful proxies for reconstructing this record, providing information on microbial phototroph primary productivity, redox couples, and oxygenation. The biomarkers okenane, chlorobactane, and isorenieratene are generally considered to be evidence of anoxygenic phototrophs, and provide a record that extends to 1.64 Ga. The utility of the carotenoid biomarker record may be enhanced by examining the carbon isotopic ratios in these products, which are diagnostic for specific pathways of biological carbon fixation found today within different microbial groups. However, this joint inference assumes that microbes have conserved these pathways across the duration of the preserved biomarker record. Testing this hypothesis, we performed phylogenetic analyses of the enzymes constituting the reductive tricarboxylic acid (rTCA) cycle in Chlorobiales, the group of anoxygenic phototrophic bacteria usually implicated in the deposition of chlorobactane and isorenieretane. We find phylogenetically incongruent patterns of inheritance across all enzymes, indicative of horizontal gene transfers to both stem and crown Chlorobiales from multiple potential donor lineages. This indicates that a complete rTCA cycle was independently acquired at least twice within Chlorobiales and was not present in the last common ancestor. When combined with recent molecular clock analyses, these results predict that the Mesoproterzoic lipid biomarker record diagnostic for Chlorobiales should not preserve isotopic fractionations indicative of a full rTCA cycle. Furthermore, we conclude that coupling isotopic and biomarker records is insufficient for reliably reconstructing microbial paleoecologies in the absence of a complementary and consistent phylogenomic narrative.

Early Nitrogenase Ancestors Encompassed Novel Active Site Diversity.

Ancestral sequence reconstruction (ASR) infers predicted ancestral states for sites within sequences and can constrain the functions and properties of ancestors of extant protein families. Here, we compare the likely sequences of inferred nitrogenase ancestors to extant nitrogenase sequence diversity. We show that the most-likely combinations of ancestral states for key substrate channel residues are not represented in extant sequence space, and rarely found within a more broadly defined physiochemical space-supporting that the earliest ancestors of extant nitrogenases likely had alternative substrate channel composition. These differences may indicate differing environmental selection pressures acting on nitrogenase substrate specificity in ancient environments. These results highlight ASR's potential as an in silico tool for developing hypotheses about ancestral enzyme functions, as well as improving hypothesis testing through more targeted in vitro and in vivo experiments.

Standard Candles for Dating Microbial Lineages.

Molecular clock analyses are challenging for microbial phylogenies, due to a lack of fossil calibrations that can reliably provide absolute time constraints. An alternative source of temporal constraints for microbial groups is provided by the inheritance of proteins that are specific for the utilization of eukaryote-derived substrates, which have often been dispersed across the Tree of Life via horizontal gene transfer. In particular, animal, algal, and plant-derived substrates are often produced by groups with more precisely known divergence times, providing an older-bound on their availability within microbial environments. Therefore, these ages can serve as "standard candles" for dating microbial groups across the Tree of Life, expanding the reach of informative molecular clock investigations. Here, we formally develop the concept of substrate standard candles and describe how they can be propagated and applied using both microbial species trees and individual gene family phylogenies. We also provide detailed evaluations of several candidate standard candles and discuss their suitability in light of their often complex evolutionary and metabolic histories.

Inferred ancestry of scytonemin biosynthesis proteins in cyanobacteria indicates a response to Paleoproterozoic oxygenation.

Protection from radiation damage is an important adaptation for phototrophic microbes. Living in surface, shallow water, and peritidal environments, cyanobacteria are especially exposed to long-wavelength ultraviolet (UVA) radiation. Several groups of cyanobacteria within these environments are protected from UVA damage by the production of the pigment scytonemin. Paleontological evidence of cyanobacteria in UVA-exposed environments from the Proterozoic, and possibly as early as the Archaean, suggests a long evolutionary history of radiation protection within this group. We show that phylogenetic analyses of enzymes in the scytonemin biosynthesis pathway support this hypothesis and reveal a deep history of vertical inheritance of this pathway within extant cyanobacterial diversity. Referencing this phylogeny to cyanobacterial molecular clocks suggests that scytonemin production likely appeared during the early Proterozoic, soon after the Great Oxygenation Event. This timing is consistent with an adaptive scenario for the evolution of scytonemin production, wherein the threat of UVA-generated reactive oxygen species becomes significantly greater once molecular oxygen is more pervasive across photosynthetic environments.

Oxidative metabolisms catalyzed Earth's oxygenation.

The burial of organic carbon, which prevents its remineralization via oxygen-consuming processes, is considered one of the causes of Earth's oxygenation. Yet, higher levels of oxygen are thought to inhibit burial. Here we propose a resolution of this conundrum, wherein Earth's initial oxygenation is favored by oxidative metabolisms generating partially oxidized organic matter (POOM), increasing burial via interaction with minerals in sediments. First, we introduce the POOM hypothesis via a mathematical argument. Second, we reconstruct the evolutionary history of one key enzyme family, flavin-dependent Baeyer-Villiger monooxygenases, that generates POOM, and show the temporal consistency of its diversification with the Proterozoic and Phanerozoic atmospheric oxygenation. Finally, we propose that the expansion of oxidative metabolisms instigated a positive feedback, which was amplified by the chemical changes to minerals on Earth's surface. Collectively, these results suggest that Earth's oxygenation is an autocatalytic transition induced by a combination of biological innovations and geological changes.

Novel nitrite reductase domain structure suggests a chimeric denitrification repertoire in the phylum Chloroflexi.

Denitrification plays a central role in the global nitrogen cycle, reducing and removing nitrogen from marine and terrestrial ecosystems. The flux of nitrogen species through this pathway has a widespread impact, affecting ecological carrying capacity, agriculture, and climate. Nitrite reductase (Nir) and nitric oxide reductase (NOR) are the two central enzymes in this pathway. Here we present a previously unreported Nir domain architecture in members of phylum Chloroflexi. Phylogenetic analyses of protein domains within Nir indicate that an ancestral horizontal transfer and fusion event produced this chimeric domain architecture. We also identify an expanded genomic diversity of a rarely reported NOR subtype, eNOR. Together, these results suggest a greater diversity of denitrification enzyme arrangements exist than have been previously reported.

An Efficient, Nonphylogenetic Method for Detecting Genes Sharing Evolutionary Signals in Phylogenomic Data Sets.

Assessing the compatibility between gene family phylogenies is a crucial and often computationally demanding step in many phylogenomic analyses. Here, we describe the Evolutionary Similarity Index (IES), a means to assess shared evolution between gene families using a weighted orthogonal distance regression model applied to sequence distances. The utilization of pairwise distance matrices circumvents comparisons between gene tree topologies, which are inherently uncertain and sensitive to evolutionary model choice, phylogenetic reconstruction artifacts, and other sources of error. Furthermore, IES enables the many-to-many pairing of multiple copies between similarly evolving gene families. This is done by selecting non-overlapping pairs of copies, one from each assessed family, and yielding the least sum of squared residuals. Analyses of simulated gene family data sets show that IES's accuracy is on par with popular tree-based methods while also less susceptible to noise introduced by sequence alignment and evolutionary model fitting. Applying IES to an empirical data set of 1,322 genes from 42 archaeal genomes identified eight major clusters of gene families with compatible evolutionary trends. The most cohesive cluster consisted of 62 genes with compatible evolutionary signal, which occur as both single-copy and multiple homologs per genome; phylogenetic analysis of concatenated alignments from this cluster produced a tree closely matching previously published species trees for Archaea. Four other clusters are mainly composed of accessory genes with limited distribution among Archaea and enriched toward specific metabolic functions. Pairwise evolutionary distances obtained from these accessory gene clusters suggest patterns of interphyla horizontal gene transfer. An IES implementation is available at https://github.com/lthiberiol/evolSimIndex.

Assessing the accuracy of phylogenetic rooting methods on prokaryotic gene families.

Almost all standard phylogenetic methods for reconstructing gene trees result in unrooted trees; yet, many of the most useful applications of gene trees require that the gene trees be correctly rooted. As a result, several computational methods have been developed for inferring the root of unrooted gene trees. However, the accuracy of such methods has never been systematically evaluated on prokaryotic gene families, where horizontal gene transfer is often one of the dominant evolutionary events driving gene family evolution. In this work, we address this gap by conducting a thorough comparative evaluation of five different rooting methods using large collections of both simulated and empirical prokaryotic gene trees. Our simulation study is based on 6000 true and reconstructed gene trees on 100 species and characterizes the rooting accuracy of the four methods under 36 different evolutionary conditions and 3 levels of gene tree reconstruction error. The empirical study is based on a large, carefully designed data set of 3098 gene trees from 504 bacterial species (406 Alphaproteobacteria and 98 Cyanobacteria) and reveals insights that supplement those gleaned from the simulation study. Overall, this work provides several valuable insights into the accuracy of the considered methods that will help inform the choice of rooting methods to use when studying microbial gene family evolution. Among other findings, this study identifies parsimonious Duplication-Transfer-Loss (DTL) rooting and Minimal Ancestor Deviation (MAD) rooting as two of the most accurate gene tree rooting methods for prokaryotes and specifies the evolutionary conditions under which these methods are most accurate, demonstrates that DTL rooting is highly sensitive to high evolutionary rates and gene tree error, and that rooting methods based on branch-lengths are generally robust to gene tree reconstruction error.

An Expanded Ribosomal Phylogeny of Cyanobacteria Supports a Deep Placement of Plastids.

The phylum Cyanobacteria includes free-living bacteria and plastids, the descendants of cyanobacteria that were engulfed by the ancestral lineage of the major photosynthetic eukaryotic group Archaeplastida. Endosymbiotic events that followed this primary endosymbiosis spread plastids across diverse eukaryotic groups. The remnants of the ancestral cyanobacterial genome present in all modern plastids, enable the placement of plastids within Cyanobacteria using sequence-based phylogenetic analyses. To date, such phylogenetic studies have produced conflicting results and two competing hypotheses: (1) plastids diverge relatively recently in cyanobacterial evolution and are most closely related to nitrogen-fixing cyanobacteria, or (2) plastids diverge early in the evolutionary history of cyanobacteria, before the divergence of most cyanobacterial lineages. Here, we use phylogenetic analysis of ribosomal proteins from an expanded data set of cyanobacterial and representative plastid genomes to infer a deep placement for the divergence of the plastid ancestor lineage. We recover plastids as sister to Gloeomargarita and show that the group diverges from other cyanobacterial groups before Pseudanabaena, a previously unreported placement. The tree topologies and phylogenetic distances in our study have implications for future molecular clock studies that aim to model accurate divergence times, especially with respect to groups containing fossil calibrations. The newly sequenced cyanobacterial groups included here will also enable the use of novel cyanobacterial microfossil calibrations.

Paleozoic diversification of terrestrial chitin-degrading bacterial lineages.

BACKGROUND: Establishing the divergence times of groups of organisms is a major goal of evolutionary biology. This is especially challenging for microbial lineages due to the near-absence of preserved physical evidence (diagnostic body fossils or geochemical biomarkers). Horizontal gene transfer (HGT) can serve as a temporal scaffold between microbial groups and other fossil-calibrated clades, potentially improving these estimates. Specifically, HGT to or from organisms with fossil-calibrated age estimates can propagate these constraints to additional groups that lack fossils. While HGT is common between lineages, only a small subset of HGT events are potentially informative for dating microbial groups. RESULTS: Constrained by published fossil-calibrated studies of fungal evolution, molecular clock analyses show that multiple clades of Bacteria likely acquired chitinase homologs via HGT during the very late Neoproterozoic into the early Paleozoic. These results also show that, following these HGT events, recipient terrestrial bacterial clades likely diversified ~ 300-500 million years ago, consistent with established timescales of arthropod and plant terrestrialization. CONCLUSIONS: We conclude that these age estimates are broadly consistent with the dispersal of chitinase genes throughout the microbial world in direct response to the evolution and ecological expansion of detrital-chitin producing groups. The convergence of multiple lines of evidence demonstrates the utility of HGT-based dating methods in microbial evolution. The pattern of inheritance of chitinase genes in multiple terrestrial bacterial lineages via HGT processes suggests that these genes, and possibly other genes encoding substrate-specific enzymes, can serve as a "standard candle" for dating microbial lineages across the Tree of Life.

A Briefly Argued Case That Asgard Archaea Are Part of the Eukaryote Tree.

The recent discovery of the Lokiarchaeota and other members of the Asgard superphylum suggests that closer analysis of the cell biology and evolution of these groups may help shed light on the origin of the eukaryote cell. Asgard lineages often appear in molecular phylogenies as closely related to eukaryotes, and possess "Eukaryote Signature Proteins" coded by genes previously thought to be unique to eukaryotes. This phylogenetic affinity to eukaryotes has been widely interpreted as indicating that Asgard lineages are "eukaryote-like archaea," with eukaryotes evolving from within a paraphyletic Archaea. Guided by the established principles of systematics, we examine the potential implications of the monophyly of Asgard lineages and Eukarya. We show that a helpful parallel case is that of Synapsida, a group that includes modern mammals and their more "reptile-like" ancestors, united by shared derived characters that evolved in their common ancestor. While this group contains extinct members that share many similarities with modern reptiles and their extinct relatives, they are evolutionarily distinct from Sauropsida, the group which includes modern birds, reptiles, and all other amniotes. Similarly, Asgard lineages and eukaryotes are united by shared derived characters to the exclusion of all other groups. Consequently, the Asgard group is not only highly informative for our understanding of eukaryogenesis, but may be better understood as being early diverging members of a broader group including eukaryotes, for which we propose the name "Eukaryomorpha." Significantly, this means that the relationship between Eukarya and Asgard lineages cannot, on its own, resolve the debate over 2 vs. 3 Domains of life; instead, resolving this debate depends upon identifying the root of Archaea with respect to Bacteria.