Unexpectedly high mutation rate of a deep-sea hyperthermophilic anaerobic archaeon

Challenges in estimating effective population sizes from metagenome-assembled genomes

Effective population size (Ne) plays a critical role in shaping the relative efficiency between natural selection and genetic drift, thereby serving as a cornerstone for understanding microbial ecological dynamics. Direct Ne estimation relies on neutral genetic diversity within closely related genomes, which is, however, often constrained by the culturing difficulties for the vast majority of prokaryotic lineages. Metagenome-assembled genomes (MAGs) offer a high-throughput alternative for genomic data acquisition, yet their accuracy in Ne estimation has not been fully verified. This study examines the Thermococcus genus, comprising 66 isolated strains and 29 MAGs, to evaluate the reliability of MAGs in Ne estimation. Despite the even distribution across the Thermococcus phylogeny and the comparable internal average nucleotide identity (ANI) between isolate populations and MAG populations, our results reveal consistently lower Ne estimates from MAG populations. This trend of underestimation is also observed in various MAG populations across three other bacterial genera. The underrepresentation of genetic variation in MAGs, including loss of allele frequency data and variable genomic segments, likely contributes to the underestimation of Ne. Our findings underscore the necessity for caution when employing MAGs for evolutionary studies, which often depend on high-quality genome assemblies and nucleotide-level diversity.

Similar mutation rates but different mutation spectra in moderate and extremely halophilic archaea

Archaea are a major part of Earth's microbiota and extremely diverse. Yet, we know very little about the process of mutation that drives such diversification. To expand beyond previous work with the moderate halophilic archaeal species Haloferax volcanii, we performed a mutation-accumulation experiment followed by whole-genome sequencing in the extremely halophilic archaeon Halobacterium salinarum. Although Hfx. volcanii and Hbt. salinarum have different salt requirements, both species have highly polyploid genomes and similar GC content. We accumulated mutations for an average of 1250 generations in 67 mutation accumulation lines of Hbt. salinarum, and revealed 84 single-base substitutions and 10 insertion-deletion mutations. The estimated base-substitution mutation rate of 3.99 × 10-10 per site per generation or 1.0 × 10-3 per genome per generation in Hbt. salinarum is similar to that reported for Hfx. volcanii (1.2 × 10-3 per genome per generation), but the genome-wide insertion-deletion rate and spectrum of mutations are somewhat dissimilar in these archaeal species. The spectra of spontaneous mutations were AT biased in both archaea, but they differed in significant ways that may be related to differences in the fidelity of DNA replication/repair mechanisms or a simple result of the different salt concentrations.

Proteome-wide 3D structure prediction provides insights into the ancestral metabolism of ancient archaea and bacteria

Ancestral metabolism has remained controversial due to a lack of evidence beyond sequence-based reconstructions. Although prebiotic chemists have provided hints that metabolism might originate from non-enzymatic protometabolic pathways, gaps between ancestral reconstruction and prebiotic processes mean there is much that is still unknown. Here, we apply proteome-wide 3D structure predictions and comparisons to investigate ancestorial metabolism of ancient bacteria and archaea, to provide information beyond sequence as a bridge to the prebiotic processes. We compare representative bacterial and archaeal strains, which reveal surprisingly similar physiological and metabolic characteristics via microbiological and biophysical experiments. Pairwise comparison of protein structures identify the conserved metabolic modules in bacteria and archaea, despite interference from overly variable sequences. The conserved modules (for example, middle of glycolysis, partial TCA, proton/sulfur respiration, building block biosynthesis) constitute the basic functions that possibly existed in the archaeal-bacterial common ancestor, which are remarkably consistent with the experimentally confirmed protometabolic pathways. These structure-based findings provide a new perspective to reconstructing the ancestral metabolism and understanding its origin, which suggests high-throughput protein 3D structure prediction is a promising approach, deserving broader application in future ancestral exploration.

Comparison of prokaryotes between Mount Everest and the Mariana Trench

BACKGROUND

Mount Everest and the Mariana Trench represent the highest and deepest places on Earth, respectively. They are geographically separated, with distinct extreme environmental parameters that provide unique habitats for prokaryotes. Comparison of prokaryotes between Mount Everest and the Mariana Trench will provide a unique perspective to understanding the composition and distribution of environmental microbiomes on Earth.

RESULTS

Here, we compared prokaryotic communities between Mount Everest and the Mariana Trench based on shotgun metagenomic analysis. Analyzing 25 metagenomes and 1176 metagenome-assembled genomes showed distinct taxonomic compositions between Mount Everest and the Mariana Trench, with little taxa overlap, and significant differences in genome size, GC content, and predicted optimal growth temperature. However, community metabolic capabilities exhibited striking commonality, with > 90% of metabolic modules overlapping among samples of Mount Everest and the Mariana Trench, with the only exception for CO₂ fixations (photoautotrophy in Mount Everest but chemoautotrophy in the Mariana Trench). Most metabolic pathways were common but performed by distinct taxa in the two extreme habitats, even including some specialized metabolic pathways, such as the versatile degradation of various refractory organic matters, heavy metal metabolism (e.g., As and Se), stress resistance, and antioxidation. The metabolic commonality indicated the overall consistent roles of prokaryotes in elemental cycling and common adaptation strategies to overcome the distinct stress conditions despite the intuitively huge differences in Mount Everest and the Mariana Trench.

CONCLUSION

Our results, the first comparison between prokaryotes in the highest and the deepest habitats on Earth, may highlight the principles of prokaryotic diversity: although taxa are habitat-specific, primary metabolic functions could be always conserved. Video abstract.

Strategy of micro-environmental adaptation to cold seep among different brittle stars’ colonization

Diffusing fluid from methane seepage in cold seep field creates zones with physicochemical gradients and divergent ecosystems like the mussel beds and clam beds. Three species of brittle stars (Ophiuroidea) were discovered in the Haima cold seep fields, of which Ophiophthalmus serratus and Histampica haimaensis were found on top of or within mussel beds and clam beds, whereas Amphiura sp. was only collected from muds in the clam bed assemblage. Here, we evaluated the genetic signatures of micro-environmental adaptation of brittle stars to cold seep through the comparison of mitogenomes. This study provided two complete mitogenome sequences of O. serratus and Amphiura sp. and compared with those of H. haimaensis and other non-seep species. We found that the split events of the seep and non-seep species were as ancient as the Cretaceous period (∼148–98 Mya). O. serratus and H. haimaensis display rapid residue mutation and mitogenome rearrangements compared to their shallow or deep-sea relatives, in contrast, Amphiura sp. only show medium, regardless of nucleotide mutation rate or mitogenome rearrangement, which may correlate with their adaptation to one or two micro-ecosystems. Furthermore, we identified 10 positively selected residues in ND4 in the Amphiura sp. lineage, suggesting important roles of the dehydrogenase complex in Amphiura sp. adaptive to the cold seep environment. Our results shed light on the different evolutionary strategies during colonization in different micro-environments.

A Large-Scale Genome-Based Survey of Acidophilic Bacteria Suggests That Genome Streamlining Is an Adaption for Life at Low pH

The genome streamlining theory suggests that reduction of microbial genome size optimizes energy utilization in stressful environments. Although this hypothesis has been explored in several cases of low-nutrient (oligotrophic) and high-temperature environments, little work has been carried out on microorganisms from low-pH environments, and what has been reported is inconclusive. In this study, we performed a large-scale comparative genomics investigation of more than 260 bacterial high-quality genome sequences of acidophiles, together with genomes of their closest phylogenetic relatives that live at circum-neutral pH. A statistically supported correlation is reported between reduction of genome size and decreasing pH that we demonstrate is due to gene loss and reduced gene sizes. This trend is independent from other genome size constraints such as temperature and G + C content. Genome streamlining in the evolution of acidophilic bacteria is thus supported by our results. The analyses of predicted Clusters of Orthologous Genes (COG) categories and subcellular location predictions indicate that acidophiles have a lower representation of genes encoding extracellular proteins, signal transduction mechanisms, and proteins with unknown function but are enriched in inner membrane proteins, chaperones, basic metabolism, and core cellular functions. Contrary to other reports for genome streamlining, there was no significant change in paralog frequencies across pH. However, a detailed analysis of COG categories revealed a higher proportion of genes in acidophiles in the following categories: "replication and repair," "amino acid transport," and "intracellular trafficking". This study brings increasing clarity regarding the genomic adaptations of acidophiles to life at low pH while putting elements, such as the reduction of average gene size, under the spotlight of streamlining theory.

Prochlorococcus have low global mutation rate and small effective population size

Prochlorococcus are the most abundant free-living photosynthetic carbon-fixing organisms in the ocean. Prochlorococcus show small genome sizes, low genomic G+C content, reduced DNA repair gene pool and fast evolutionary rates, which are typical features of endosymbiotic bacteria. Nevertheless, their evolutionary mechanisms are believed to be different. Evolution of endosymbiotic bacteria is dominated by genetic drift owing to repeated population bottlenecks, whereas Prochlorococcus are postulated to have extremely large effective population sizes (N_e) and thus drift has rarely been considered. However, accurately extrapolating N_e requires measuring an unbiased global mutation rate through mutation accumulation, which is challenging for Prochlorococcus. Here, we managed this experiment over 1,065 days using Prochlorococcus marinus AS9601, sequenced genomes of 141 mutant lines and determined its mutation rate to be 3.50 × 10^-10 per site per generation. Extrapolating N_e additionally requires identifying population boundaries, which we defined using PopCOGenT and over 400 genomes related to AS9601. Accordingly, we calculated its N_e to be 1.68 × 10⁷, which is only reasonably greater than that of endosymbiotic bacteria but surprisingly smaller than that of many free-living bacteria extrapolated using the same approach. Our results therefore suggest that genetic drift is a key driver of Prochlorococcus evolution.

Mechanisms driving genome reduction of a novel Roseobacter lineage

Members of the marine Roseobacter group are key players in the global carbon and sulfur cycles. While over 300 species have been described, only 2% possess reduced genomes (mostly 3-3.5 Mbp) compared to an average roseobacter (>4 Mbp). These taxonomic minorities are phylogenetically diverse but form a Pelagic Roseobacter Cluster (PRC) at the genome content level. Here, we cultivated eight isolates constituting a novel Roseobacter lineage which we named 'CHUG'. Metagenomic and metatranscriptomic read recruitment analyses showed that CHUG members are globally distributed and active in marine pelagic environments. CHUG members possess some of the smallest genomes (~2.6 Mb) among all known roseobacters, but they do not exhibit canonical features of typical bacterioplankton lineages theorized to have undergone genome streamlining processes, like higher coding density, fewer paralogues and rarer pseudogenes. While CHUG members form a genome content cluster with traditional PRC members, they show important differences. Unlike other PRC members, neither the relative abundances of CHUG members nor their relative gene expression levels are correlated with chlorophyll a concentration across the global samples. CHUG members cannot utilize most phytoplankton-derived metabolites or synthesize vitamin B₁₂, a key metabolite mediating the roseobacter-phytoplankton interactions. This combination of features is evidence for the hypothesis that CHUG members may have evolved a free-living lifestyle decoupled from phytoplankton. This ecological transition was accompanied by the loss of signature genes involved in roseobacter-phytoplankton symbiosis, suggesting that relaxation of purifying selection owing to lifestyle shift is likely an important driver of genome reduction in CHUG.