Growth temperature and genome size in bacteria are negatively correlated, suggesting genomic streamlining during thermal adaptation

Diversification, niche adaptation, and evolution of a candidate phylum thriving in the deep Critical Zone

The deep subsurface soil microbiome encompasses a vast amount of understudied phylogenetic diversity and metabolic novelty, and the metabolic capabilities and ecological roles of these communities remain largely unknown. We observed a widespread and relatively abundant bacterial phylum (CSP1-3) in deep soils and evaluated its phylogeny, ecology, metabolism, and evolutionary history. Genome analysis indicated that members of CSP1-3 were actively replicating in situ and were widely involved in the carbon, nitrogen, and sulfur cycles. We identified potential adaptive traits of CSP1-3 members for the oligotrophic deep soil environments, including a mixotrophic lifestyle, flexible energy metabolisms, and conservation pathways. The ancestor of CSP1-3 likely originated in an aquatic environment, subsequently colonizing topsoil and, later, deep soil environments, with major CSP1-3 clades adapted to each of these distinct niches. The transition into the terrestrial environment was associated with genome expansion, including the horizontal acquisition of a range of genes for carbohydrate and energy metabolism and, in one lineage, high-affinity terminal oxidases to support a microaerophilic lifestyle. Our results highlight the ecology and genome evolution of microbes in the deep Critical Zone.

Intraspecific divergence within Microcystis aeruginosa mediates the dynamics of freshwater harmful algal blooms under climate warming scenarios

Intraspecific biodiversity can have ecosystem-level consequences and may affect the accuracy of ecological forecasting. For example, rare genetic variants may have traits that prove beneficial under future environmental conditions. The cyanobacterium responsible for most freshwater harmful algal blooms worldwide, Microcystis aeruginosa, occurs in at least three types. While the dominant type occurs in eutrophic environments and is adapted to thrive in nutrient-rich conditions, two additional types have recently been discovered that inhabit oligotrophic and eutrophic environments and have genomic adaptations for survival under nutrient limitation. Here, we show that these oligotrophic types are widespread throughout the Eastern USA. By pairing an experimental warming study with gene expression analyses, we found that the eutrophic type may be most susceptible to climate warming. In comparison, oligotrophic types maintained their growth better and persisted longer under warming. As a mechanistic explanation for these patterns, we found that oligotrophic types responded to warming by widespread elevated expression of heat shock protein genes. Reduction of nutrient loading has been a historically effective mitigation strategy for controlling harmful algal blooms. Our results suggest that climate warming may benefit oligotrophic types of M. aeruginosa, potentially reducing the effectiveness of such mitigation efforts. In-depth study of intraspecific variation may therefore improve forecasting for understanding future whole ecosystem dynamics.

Genome evolution following an ecological shift in nectar-dwelling Acinetobacter

The bacterial genus Acinetobacter includes species found in environmental habitats like soil and water, as well as taxa adapted to be host-associated or pathogenic. High genetic diversity may allow for this habitat flexibility, but the specific genes underlying switches between habitats are poorly understood. One lineage of Acinetobacter has undergone a substantial habitat change by evolving from a presumed soil-dwelling ancestral state to thrive in floral nectar. Here, we compared the genomes of floral-dwelling and pollinator-associated Acinetobacter, including newly described species, with genomes from relatives found in other environments to determine the genomic changes associated with this ecological shift. Following one evolutionary origin of floral nectar adaptation, nectar-dwelling Acinetobacter taxa have undergone reduction in genome size compared with relatives and have experienced dynamic gene gains and losses as they diversified. Gene content changes suggest a shift to metabolism of monosaccharides rather than diverse carbohydrates, and scavenging of nitrogen sources, which we predict to be beneficial in nectar environments. Gene gains appear to result from duplication events, evolutionary divergence, and horizontal gene transfer. Most notably, nectar-dwelling Acinetobacter acquired the ability to degrade pectin from plant pathogens, and the genes underlying this ability have duplicated and are under selection within the clade. We hypothesize that this ability was a key trait for adaptation to floral nectar, as it could improve access to nutrients in the nutritionally unbalanced habitat of nectar. These results identify the genomic changes and traits coinciding with a dramatic habitat switch from soil to floral nectar.

IMPORTANCE

Many bacteria, including the genus Acinetobacter, commonly evolve to exploit new habitats. However, the genetic changes that underlie habitat switches are often unknown. Floral nectar is home to specialized microbes that can grow in this nutritionally unbalanced habitat. Several specialized Acinetobacter species evolved from soil-dwelling relatives to become common and abundant in floral nectar. Here, we investigate the genomic adaptations required to successfully colonize a novel habitat like floral nectar. We performed comparative genomics analyses between nectar-dwelling Acinetobacter and Acinetobacter species from other environments, like soil and water. We find that although gene loss coincided with the switch to living in nectar, gains of specific genes from other bacteria may have been particularly important for this ecological change. Acinetobacter living in nectar gained genes for degrading pectin, a plant polysaccharide, which may improve access to nutrients in their environment. These findings shed light on how evolutionary novelty evolves in bacteria.

Bacterial Growth Temperature as a Horizontally Acquired Polygenic Trait

Evolutionary events leading to organismal preference for a specific growth temperature, as well as genes whose products are needed for a proper function at that temperature, are poorly understood. Using 64 bacteria from phylum Thermotogota as a model system, we examined how optimal growth temperature changed throughout Thermotogota history. We inferred that Thermotogota's last common ancestor was a thermophile and that some Thermotogota evolved the mesophilic and hyperthermophilic lifestyles secondarily. By modeling gain and loss of genes throughout Thermotogota history and by reconstructing their phylogenies, we demonstrated that adaptations to lower and higher growth temperature involve both the acquisition of necessary genes and loss of unnecessary genes. Via a pangenome-wide association study, we correlated presence/absence of 68 genes with specific optimal growth temperature intervals. While some of these genes are poorly characterized, most are involved in metabolism of amino acids, nucleotides, carbohydrates, and lipids, as well as in signal transduction and regulation of transcription. Most of the 68 genes have a history of horizontal gene transfer with other bacteria and archaea that often grow at similar temperatures, suggesting that parallel acquisitions of genes likely promote independent adaptations of different Thermotogota species to specific growth temperatures.

Genome-streamlined SAR202 bacteria are widely present and active in the euphotic ocean

SAR202 bacteria are a diverse group of bacteria in the ocean. The SAR202 lineages dominate the bacterial community and evolve specialized metabolisms for oxidizing recalcitrant organic compounds in the dark ocean. SAR202 bacteria are also present in the euphotic oceans; however, their ecological roles and metabolic potential remain poorly understood. In this study, we collected 392 non-redundant metagenome-assembled genomes from different oceans, with 18% of these SAR202 genomes characterized by small genome sizes (<2 Mbp), low GC content (<40%), and high gene density. The 70 genome-streamlined SAR202 bacteria constitute more than an average of 90% of SAR202 in the euphotic zone and exhibit streamlined metabolic features compared to the dark ocean SAR202. Genome-streamlined SAR202 is distributed in many major SAR202 lineages (i.e. I, II, III, and VI). Phylogenomic analysis shows that the genome-streamlined SAR202 clades diverged from the non-genome-streamlined SAR202 lineages and evolved independently within the same clades. Certain genes are enriched in genome-streamlined SAR202, such as proteorhodopsin genes and the coding genes of major facilitator superfamily transporters, nucleoside transporters, and deoxyribodipyrimidine photo-lyase, indicating their adaptation to sunlit oligotrophic water. A detailed comparison between genome-streamlined SAR202 and non-genome-streamlined SAR202 was made to illustrate their distinct niche distribution and metabolic buildup. In addition, the metatranscriptomic analysis supports that genome-streamlined SAR202 bacteria are active in the upper ocean. This study represents a systematic study of streamlined SAR202 bacteria that occupy the euphotic ocean and provides a comprehensive view of the ecological roles of SAR202 bacteria in the ocean.

Adaptation of Archaeal Communities to Summer Hypoxia in the Sediment of Bohai Sea

Understanding the adaptation of archaea to hypoxia is essential for deciphering the functions and mechanisms of microbes when suffering environmental changes. However, the dynamics and responses of archaea to the sedimentary hypoxia in Bohai Sea are still unclear. In this study, the diversity, composition, and distribution of archaeal community in sediment along an inshore-offshore transect across the oxygen-depleted area in the Bohai Sea were investigated in June, July, and August of 2021 by employing high-throughput sequencing of 16S rRNA gene. Results indicated that the archaeal communities were dominated by Thermoproteota (80.61%), Asgardarchaeota (8.70%), and Thermoplasmatota (5.27%). Dissolved oxygen (DO) and NO3 - were the two key factors shaping the distribution of archaeal communities, accounting for 49.5% and 38.3% of the total variabilities (p < 0.05), respectively. With the intensity of oxygen depletion, the diversity of archaeal communities increased significantly. Microbial networks revealed that Bathyarchaeia played a key role in interacting with both bacteria and other archaeal groups. Furthermore, adaptions to hypoxia of archaea were also displayed by variation in relative abundance of the predicted ecological functions and the metabolic pathways. The enrichment of specific nitrogen transformation enzymes showed the potential for nitrogen fixation and removal, which might contribute to the balance of N budget and thus facilitate the ecological restoration under eutrophication in Bohai Sea. Our results provided a new picture on ecological and metabolic adaptions to hypoxia by archaea, which will be beneficial to further investigations in extreme environments both theoretically and practically.

Bigger genomes provide environment-dependent growth benefits in grasses

Increasing genome size (GS) has been associated with slower rates of DNA replication and greater cellular nitrogen (N) and phosphorus demands. Despite most plant species having small genomes, the existence of larger GS species suggests that such costs may be negligible or represent benefits under certain conditions. Focussing on the widespread and diverse grass family (Poaceae), we used data on species' climatic niches and growth rates under different environmental conditions to test for growth costs or benefits associated with GS. The influence of photosynthetic pathway, life history and evolutionary history on grass GS was also explored. We found that evolutionary history, photosynthetic pathway and life history all influence the distribution of grass species' GS. Genomes were smaller in annual and C4 species, the latter allowing for small cells necessary for C4 leaf anatomy. We found larger GS were associated with high N availability and, for perennial species, low growth-season temperature. Our findings reveal that GS is a globally important predictor of grass performance dependent on environmental conditions. The benefits for species with larger GS are likely due to associated larger cell sizes, allowing rapid biomass production where soil fertility meets N demands and/or when growth occurs via temperature-independent cell expansion.

Microbe-aided thermophilic composting accelerates manure fermentation

Aerobic composting is a key strategy to the sustainable use of livestock manure, which is however constrained by the slow kinetics. Microbe-aided thermophilic composting provides an attractive solution to this problem. In this study, we identified key thermophilic bacteria capable of accelerating manure composting based on the deciphering of manure bacterial community evolution in a thermophilic system. High-throughput sequencing showed a significant evolution of manure bacterial community structure with the increasing heating temperature. Firmicutes were substantially enriched by the heating, particularly some known thermotolerant bacterial species, such as Novibacillus thermophiles, Bacillus thermolactis, and Ammoniibacillus agariperforans. Correspondingly, through function prediction, we found bacterial taxa with cellulolytic and xylanolytic activities were significantly higher in the thermophilic process relative to the initial stage. Subsequently, a total of 47 bacteria were isolated in situ and their phylogenetic affiliation and degradation capacity were determined. Three isolates were back inoculated to the manure, resulting in shortened composting process from 5 to 3 days with Germination Index increased up to 134%, and improved compost quality particularly in wheat growth promoting. Comparing to the mesophilic and thermophilic Bacillus, the genomes of the three isolates manifested some features similar to the thermophiles, including smaller genome size and mutation of specific genes that enhance heat tolerance. This study provide robust evidence that microbe-aided thermophilic composting is capable to accelerate manure composting and improve the quality of compost, which represents a new hope to the sustainable use of manure from the meat industry.

High-quality genome of a novel Thermosynechococcaceae species from Namibia and characterization of its protein expression patterns at elevated temperatures

Thermophilic cyanobacteria thrive in extreme environments, making their thermoresistant enzymes valuable for industrial applications. Common habitats include hot springs, which act as evolutionary accelerators for speciation due to geographical isolation. The family Thermosynechococcaceae comprises thermophilic cyanobacteria known for their ability to thrive in high-temperature environments. These bacteria are notable for their photosynthetic capabilities, significantly contributing to primary production in extreme habitats. Members of Thermosynechococcaceae exhibit unique adaptations that allow them to perform photosynthesis efficiently at elevated temperatures, making them subjects of interest for studies on microbial ecology, evolution, and potential biotechnological applications. In this study, the genome of a thermophilic cyanobacterium, isolated from a hot spring near Okahandja in Namibia, was sequenced using a PacBio Sequel IIe long-read platform. Cultivations were performed at elevated temperatures of 40, 50, and 55°C, followed by proteome analyses based on the annotated genome. Phylogenetic investigations, informed by the 16S rRNA gene and aligned nucleotide identity (ANI), suggest that the novel cyanobacterium is a member of the family Thermosynechococcaceae. Furthermore, the new species was assigned to a separate branch, potentially representing a novel genus. Whole-genome alignments supported this finding, revealing few conserved regions and multiple genetic rearrangement events. Additionally, 129 proteins were identified as differentially expressed in a temperature-dependent manner. The results of this study broaden our understanding of cyanobacterial adaptation to extreme environments, providing a novel high-quality genome of Thermosynechococcaceae cyanobacterium sp. Okahandja and several promising candidate proteins for expression and characterization studies.

Metabolic interdependencies in thermophilic communities are revealed using co-occurrence and complementarity networks

Microbial communities exhibit intricate interactions underpinned by metabolic dependencies. To elucidate these dependencies, we present a workflow utilizing random matrix theory on metagenome-assembled genomes to construct co-occurrence and metabolic complementarity networks. We apply this approach to a temperature gradient hot spring, unraveling the interplay between thermal stress and metabolic cooperation. Our analysis reveals an increase in the frequency of metabolic interactions with rising temperatures. Amino acids, coenzyme A derivatives, and carbohydrates emerge as key exchange metabolites, forming the foundation for syntrophic dependencies, in which commensalistic interactions take a greater proportion than mutualistic ones. These metabolic exchanges are most prevalent between phylogenetically distant species, especially archaea-bacteria collaborations, as a crucial adaptation to harsh environments. Furthermore, we identify a significant positive correlation between basal metabolite exchange and genome size disparity, potentially signifying a means for streamlined genomes to leverage cooperation with metabolically richer partners. This phenomenon is also confirmed by another composting system which has a similar wide range of temperature fluctuations. Our workflow provides a feasible way to decipher the metabolic complementarity mechanisms underlying microbial interactions, and our findings suggested environmental stress regulates the cooperative strategies of thermophiles, while these dependencies have been potentially hardwired into their genomes during co-evolutions.

Comparative genomic analysis of thermophilic fungi reveals convergent evolutionary adaptations and gene losses

Thermophily is a trait scattered across the fungal tree of life, with its highest prevalence within three fungal families (Chaetomiaceae, Thermoascaceae, and Trichocomaceae), as well as some members of the phylum Mucoromycota. We examined 37 thermophilic and thermotolerant species and 42 mesophilic species for this study and identified thermophily as the ancestral state of all three prominent families of thermophilic fungi. Thermophilic fungal genomes were found to encode various thermostable enzymes, including carbohydrate-active enzymes such as endoxylanases, which are useful for many industrial applications. At the same time, the overall gene counts, especially in gene families responsible for microbial defense such as secondary metabolism, are reduced in thermophiles compared to mesophiles. We also found a reduction in the core genome size of thermophiles in both the Chaetomiaceae family and the Eurotiomycetes class. The Gene Ontology terms lost in thermophilic fungi include primary metabolism, transporters, UV response, and O-methyltransferases. Comparative genomics analysis also revealed higher GC content in the third base of codons (GC3) and a lower effective number of codons in fungal thermophiles than in both thermotolerant and mesophilic fungi. Furthermore, using the Support Vector Machine classifier, we identified several Pfam domains capable of discriminating between genomes of thermophiles and mesophiles with 94% accuracy. Using AlphaFold2 to predict protein structures of endoxylanases (GH10), we built a similarity network based on the structures. We found that the number of disulfide bonds appears important for protein structure, and the network clusters based on protein structures correlate with the optimal activity temperature. Thus, comparative genomics offers new insights into the biology, adaptation, and evolutionary history of thermophilic fungi while providing a parts list for bioengineering applications.

Lessons from Extremophiles: Functional Adaptations and Genomic Innovations across the Eukaryotic Tree of Life

From hydrothermal vents, to glaciers, to deserts, research in extreme environments has reshaped our understanding of how and where life can persist. Contained within the genomes of extremophilic organisms are the blueprints for a toolkit to tackle the multitude of challenges of survival in inhospitable environments. As new sequencing technologies have rapidly developed, so too has our understanding of the molecular and genomic mechanisms that have facilitated the success of extremophiles. Although eukaryotic extremophiles remain relatively understudied compared to bacteria and archaea, an increasing number of studies have begun to leverage 'omics tools to shed light on eukaryotic life in harsh conditions. In this perspective paper, we highlight a diverse breadth of research on extremophilic lineages across the eukaryotic tree of life, from microbes to macrobes, that are collectively reshaping our understanding of molecular innovations at life's extremes. These studies are not only advancing our understanding of evolution and biological processes but are also offering a valuable roadmap on how emerging technologies can be applied to identify cellular mechanisms of adaptation to cope with life in stressful conditions, including high and low temperatures, limited water availability, and heavy metal habitats. We shed light on patterns of molecular and organismal adaptation across the eukaryotic tree of life and discuss a few promising research directions, including investigations into the role of horizontal gene transfer in eukaryotic extremophiles and the importance of increasing phylogenetic diversity of model systems.

Eco-evolutionary strategies for relieving carbon limitation under salt stress differ across microbial clades

With the continuous expansion of saline soils under climate change, understanding the eco-evolutionary tradeoff between the microbial mitigation of carbon limitation and the maintenance of functional traits in saline soils represents a significant knowledge gap in predicting future soil health and ecological function. Through shotgun metagenomic sequencing of coastal soils along a salinity gradient, we show contrasting eco-evolutionary directions of soil bacteria and archaea that manifest in changes to genome size and the functional potential of the soil microbiome. In salt environments with high carbon requirements, bacteria exhibit reduced genome sizes associated with a depletion of metabolic genes, while archaea display larger genomes and enrichment of salt-resistance, metabolic, and carbon-acquisition genes. This suggests that bacteria conserve energy through genome streamlining when facing salt stress, while archaea invest in carbon-acquisition pathways to broaden their resource usage. These findings suggest divergent directions in eco-evolutionary adaptations to soil saline stress amongst microbial clades and serve as a foundation for understanding the response of soil microbiomes to escalating climate change.

Genome streamlining in Parcubacteria transitioning from soil to groundwater

BACKGROUND

To better understand the influence of habitat on the genetic content of bacteria, with a focus on members of Candidate Phyla Radiation (CPR) bacteria, we studied the effects of transitioning from soil via seepage waters to groundwater on genomic composition of ultra-small Parcubacteria, the dominating CPR class in seepage waters, using genome resolved metagenomics.

RESULTS

Bacterial metagenome-assembled genomes (MAGs), (318 total, 32 of Parcubacteria) were generated from seepage waters and compared directly to groundwater counterparts. The estimated average genome sizes of members of major phyla Proteobacteria, Bacteroidota and Cand. Patescibacteria (Candidate Phyla Radiation - CPR bacteria) were significantly higher in soil-seepage water as compared to their groundwater counterparts. Seepage water Parcubacteria (Paceibacteria) exhibited 1.18-fold greater mean genome size and 2-fold lower mean proportion of pseudogenes than those in groundwater. Bacteroidota and Proteobacteria also showed a similar trend of reduced genomes in groundwater compared to seepage. While exploring gene loss and adaptive gains in closely related CPR lineages in groundwater, we identified a membrane protein, and a lipoglycopeptide resistance gene unique to a seepage Parcubacterium genome. A nitrite reductase gene was also identified and was unique to the groundwater Parcubacteria genomes, likely acquired from other planktonic microbes via horizontal gene transfer.

CONCLUSIONS

Overall, our data suggest that bacteria in seepage waters, including ultra-small Parcubacteria, have significantly larger genomes and higher metabolic enrichment than their groundwater counterparts, highlighting possible genome streamlining of the latter in response to habitat selection in an oligotrophic environment.

Hot springs viruses at Yellowstone National Park have ancient origins and are adapted to thermophilic hosts

Geothermal springs house unicellular red algae in the class Cyanidiophyceae that dominate the microbial biomass at these sites. Little is known about host-virus interactions in these environments. We analyzed the virus community associated with red algal mats in three neighboring habitats (creek, endolithic, soil) at Lemonade Creek, Yellowstone National Park (YNP), USA. We find that despite proximity, each habitat houses a unique collection of viruses, with the giant viruses, Megaviricetes, dominant in all three. The early branching phylogenetic position of genes encoded on metagenome assembled virus genomes (vMAGs) suggests that the YNP lineages are of ancient origin and not due to multiple invasions from mesophilic habitats. The existence of genomic footprints of adaptation to thermophily in the vMAGs is consistent with this idea. The Cyanidiophyceae at geothermal sites originated ca. 1.5 Bya and are therefore relevant to understanding biotic interactions on the early Earth.

Adaptation in Unstable Environments and Global Gene Losses: Small but Stable Gene Networks by the May-Wigner Theory

Although gene loss is common in evolution, it remains unclear whether it is an adaptive process. In a survey of seven major mangrove clades that are woody plants in the intertidal zones of daily environmental perturbations, we noticed that they generally evolved reduced gene numbers. We then focused on the largest clade of Rhizophoreae and observed the continual gene set reduction in each of the eight species. A great majority of gene losses are concentrated on environmental interaction processes, presumably to cope with the constant fluctuations in the tidal environments. Genes of the general processes for woody plants are largely retained. In particular, fewer gene losses are found in physiological traits such as viviparous seeds, high salinity, and high tannin content. Given the broad and continual genome reductions, we propose the May-Wigner theory (MWT) of system stability as a possible mechanism. In MWT, the most effective solution for buffering continual perturbations is to reduce the size of the system (or to weaken the total genic interactions). Mangroves are unique as immovable inhabitants of the compound environments in the land-sea interface, where environmental gradients (such as salinity) fluctuate constantly, often drastically. Extending MWT to gene regulatory network (GRN), computer simulations and transcriptome analyses support the stabilizing effects of smaller gene sets in mangroves vis-à-vis inland plants. In summary, we show the adaptive significance of gene losses in mangrove plants, including the specific role of promoting phenotype innovation and a general role in stabilizing GRN in unstable environments as predicted by MWT.

The global distribution of angiosperm genome size is shaped by climate

Angiosperms, which inhabit diverse environments across all continents, exhibit significant variation in genome sizes, making them an excellent model system for examining hypotheses about the global distribution of genome size. These include the previously proposed large genome constraint, mutational hazard, polyploidy-mediated, and climate-mediated hypotheses. We compiled the largest genome size dataset to date, encompassing 16 017 (> 5% of known) angiosperm species, and analyzed genome size distribution using a comprehensive geographic distribution dataset for all angiosperms. We observed that angiosperms with large range sizes generally had small genomes, supporting the large genome constraint hypothesis. Climate was shown to exert a strong influence on genome size distribution along the global latitudinal gradient, while the frequency of polyploidy and the type of growth form had negligible effects. In contrast to the unimodal patterns along the global latitudinal gradient shown by plant size traits and polyploid proportions, the increase in angiosperm genome size from the equator to 40-50°N/S is probably mediated by different (mostly climatic) mechanisms than the decrease in genome sizes observed from 40 to 50°N northward. Our analysis suggests that the global distribution of genome sizes in angiosperms is mainly shaped by climatically mediated purifying selection, genetic drift, relaxed selection, and environmental filtering.

Idiosyncratic genome evolution of the thermophilic cyanobacterium Synechococcus at the limits of phototrophy

Thermophilic microorganisms are expected to have smaller cells and genomes compared with mesophiles, a higher proportion of horizontally acquired genes, and distinct nucleotide and amino acid composition signatures. Here, we took an integrative approach to investigate these apparent correlates of thermophily for Synechococcus A/B cyanobacteria, which include the most heat-tolerant phototrophs on the planet. Phylogenomics confirmed a unique origin of different thermotolerance ecotypes, with low levels of continued gene flow between ecologically divergent but overlapping populations, which has shaped the distribution of phenotypic traits along these geothermal gradients. More thermotolerant strains do have smaller genomes, but genome reduction is associated with a decrease in community richness and metabolic diversity, rather than with cell size. Horizontal gene transfer played only a limited role during Synechococcus evolution, but, the most thermotolerant strains have acquired a Thermus tRNA modification enzyme that may stabilize translation at high temperatures. Although nucleotide base composition was not associated with thermotolerance, we found a general replacement of aspartate with glutamate, as well as a dramatic remodeling of amino acid composition at the highest temperatures that substantially differed from previous predictions. We conclude that Synechococcus A/B genome diversification largely does not conform to the standard view of temperature adaptation. In addition, carbon fixation was more thermolabile than photosynthetic oxygen evolution for the most thermotolerant strains compared with less tolerant lineages. This suggests that increased flow of reducing power generated during the light reactions to an electron sink(s) beyond carbon dioxide has emerged during temperature adaptation of these bacteria.

How Do Thermophiles Organize Their Genomes?

All cells must maintain the structural and functional integrity of the genome under a wide range of environments. High temperatures pose a formidable challenge to cells by denaturing the DNA double helix, causing chemical damage to DNA, and increasing the random thermal motion of chromosomes. Thermophiles, predominantly classified as bacteria or archaea, exhibit an exceptional capacity to mitigate these detrimental effects and prosper under extreme thermal conditions, with some species tolerating temperatures higher than 100°C. Their genomes are mainly characterized by the presence of reverse gyrase, a unique topoisomerase that introduces positive supercoils into DNA. This enzyme has been suggested to maintain the genome integrity of thermophiles by limiting DNA melting and mediating DNA repair. Previous studies provided significant insights into the mechanisms by which NAPs, histones, SMC superfamily proteins, and polyamines affect the 3D genomes of thermophiles across different scales. Here, I discuss current knowledge of the genome organization in thermophiles and pertinent research questions for future investigations.

Investigation of lysing excess sludge slurry using hydrolase secreting thermophilic bacterial communities

Sludge reduction is a critical challenge in biological wastewater treatment. Combining excess sludge slurry lysis technology with traditional activated sludge processes is a promising approach for in-situ sludge reduction. Here, a strategy for excess sludge slurry lysis based on thermophilic bacterial communities (LTBC) was proposed. This investigation focused on the process of excess sludge slurry lysis dominated by thermophilic bacterial communities domesticated at different temperatures (55-75 °C). The evolution of sludge lysate was analyzed, and the mechanism of excess sludge slurry lysis under the action of thermophilic bacterial communities was elucidated through amplicon sequencing analysis. The results demonstrated that the aerobic thermophilic bacterial communities adapted to 75 °C exhibit the highest efficiency in sludge slurry lysis. During LTBC process, the removal efficiency of volatile suspended solids reached 53.9 ± 1.8% within 2 h, and 97.0 ± 1.0% of the protein and 96.0 ± 1.0% of the polysaccharide in the extracellular polymers was solubilized, and bacterial cell walls in sludge were disrupted. Fourier transform infrared spectroscopy and excitation-emission matrix spectroscopy of the sludge lysate demonstrated that the LTBC process was accompanied by humification process. The accumulation of humic acid primarily occurred at 55 °C and 65 °C, while fulvic acid occurred at 75 °C. The thermophilic bacterial communities adapted to 75 °C were dominated by Thermus and Thermaerobacter. Phylogenetic studies showed that the LTBC hydrolase system comprises enzymes related to protein hydrolysis, carbohydrate hydrolysis, and peptidoglycan hydrolysis, including metalopeptidase MepB, neutral α-glucosidase C, N-acetyl Muramyl-L-alanine amidase, and others enzymes. These results provide a theoretical basis for the application of LTBC technology in the reduction of sludge which generated in traditional waste water activated sludge processes.

Metagenomic profiles of archaea and bacteria within thermal and geochemical gradients of the Guaymas Basin deep subsurface

Previous studies of microbial communities in subseafloor sediments reported that microbial abundance and diversity decrease with sediment depth and age, and microbes dominating at depth tend to be a subset of the local seafloor community. However, the existence of geographically widespread, subsurface-adapted specialists is also possible. Here, we use metagenomic and metatranscriptomic analyses of the hydrothermally heated, sediment layers of Guaymas Basin (Gulf of California, Mexico) to examine the distribution and activity patterns of bacteria and archaea along thermal, geochemical and cell count gradients. We find that the composition and distribution of metagenome-assembled genomes (MAGs), dominated by numerous lineages of Chloroflexota and Thermoproteota, correlate with biogeochemical parameters as long as temperatures remain moderate, but downcore increasing temperatures beyond ca. 45 ºC override other factors. Consistently, MAG size and diversity decrease with increasing temperature, indicating a downcore winnowing of the subsurface biosphere. By contrast, specific archaeal MAGs within the Thermoproteota and Hadarchaeota increase in relative abundance and in recruitment of transcriptome reads towards deeper, hotter sediments, marking the transition towards a specialized deep, hot biosphere.

Hydrogenimonas cancrithermarum sp. nov., a hydrogen- and thiosulfate-oxidizing mesophilic chemolithoautotroph isolated from diffuse-flow fluids on the East Pacific Rise, and an emended description of the genus Hydrogenimonas

A novel mesophilic, hydrogen- and thiosulfate-oxidizing bacterium, strain ISO32T, was isolated from diffuse-flow hydrothermal fluids from the Crab Spa vent on the East Pacific Rise. Cells of ISO32T were rods, being motile by means of a single polar flagellum. The isolate grew at a temperature range between 30 and 55 °C (optimum, 43 °C), at a pH range between 5.3 and 7.6 (optimum, pH 5.8) and in the presence of 2.0-4.0 % NaCl (optimum, 2.5 %). The isolate was able to grow chemolithoautotrophically with molecular hydrogen, thiosulfate or elemental sulfur as the sole electron donor. Thiosulfate, elemental sulfur, nitrate and molecular oxygen were each used as a sole electron acceptor. Phylogenetic analysis of 16S rRNA gene sequences placed ISO32T in the genus Hydrogenimonas of the class Epsilonproteobacteria, with Hydrogenimonas thermophila EP1-55-1 %T as its closest relative (95.95 % similarity). On the basis of the phylogenetic, physiological and genomic characteristics, it is proposed that the organism represents a novel species within the genus Hydrogenimonas, Hydrogenimonas cancrithermarum sp. nov. The type strain is ISO32T (=JCM 39185T =KCTC 25252T). Furthermore, the genomic properties of members of the genus Hydrogenimonas are distinguished from those of members of other thermophilic genera in the orders Campylobacterales (Nitratiruptor and Nitrosophilus) and Nautiliales (Caminibacter, Nautilia and Lebetimonas), with larger genome sizes and lower 16S rRNA G+C content values. Comprehensive metabolic comparisons based on genomes revealed that genes responsible for the Pta-AckA pathway were observed exclusively in members of mesophilic genera in the order Campylobacterales and of the genus Hydrogenimonas. Our results indicate that the genus Hydrogenimonas contributes to elucidating the evolutionary history of Epsilonproteobacteria in terms of metabolism and transition from a thermophilic to a mesophilic lifestyle.

Bacterial and archaeal community distributions and cosmopolitanism across physicochemically diverse hot springs

Terrestrial hot springs harbor diverse microbial communities whose compositions are shaped by the wide-ranging physico-chemistries of individual springs. The effect of enormous physico-chemical differences on bacterial and archaeal distributions and population structures is little understood. We therefore analysed the prevalence and relative abundance of bacteria and archaea in the sediments (n = 76) of hot spring features, in the Taupō Volcanic Zone (New Zealand), spanning large differences in major anion water chemistry, pH (2.0-7.5), and temperature (17.5-92.9 °C). Community composition, based on 16S rRNA amplicon sequence variants (ASVs) was strongly influenced by both temperature and pH. However, certain lineages characterized diverse hot springs. At the domain level, bacteria and archaea shared broadly equivalent community abundances across physico-chemically diverse springs, despite slightly lower bacteria-to-archaea ratios and microbial 16S rRNA gene concentrations at higher temperatures. Communities were almost exclusively dominated by Proteobacteria, Euryarchaeota or Crenarchaeota. Eight archaeal and bacterial ASVs from Thermoplasmatales, Desulfurellaceae, Mesoaciditogaceae and Acidithiobacillaceae were unusually prevalent (present in 57.9-84.2% of samples) and abundant (1.7-12.0% sample relative abundance), and together comprised 44% of overall community abundance. Metagenomic analyses generated multiple populations associated with dominant ASVs, and showed characteristic traits of each lineage for sulfur, nitrogen and hydrogen metabolism. Differences in metabolic gene composition and genome-specific metabolism delineated populations from relatives. Genome coverage calculations showed that populations associated with each lineage were distributed across a physicochemically broad range of hot springs. Results imply that certain bacterial and archaeal lineages harbor different population structures and metabolic potentials for colonizing diverse hot spring environments.

Hot spring distribution and survival mechanisms of thermophilic comammox Nitrospira

The recent discovery of Nitrospira species capable of complete ammonia oxidation (comammox) in non-marine natural and engineered ecosystems under mesothermal conditions has changed our understanding of microbial nitrification. However, little is known about the occurrence of comammox bacteria or their ability to survive in moderately thermal and/or hyperthermal habitats. Here, we report the wide distribution of comammox Nitrospira in five terrestrial hot springs at temperatures ranging from 36 to 80°C and provide metagenome-assembled genomes of 11 new comammox strains. Interestingly, the identification of dissimilatory nitrate reduction to ammonium (DNRA) in thermophilic comammox Nitrospira lineages suggests that they have versatile ecological functions as both sinks and sources of ammonia, in contrast to the described mesophilic comammox lineages, which lack the DNRA pathway. Furthermore, the in situ expression of key genes associated with nitrogen metabolism, thermal adaptation, and oxidative stress confirmed their ability to survive in the studied hot springs and their contribution to nitrification in these environments. Additionally, the smaller genome size and higher GC content, less polar and more charged amino acids in usage profiles, and the expression of a large number of heat shock proteins compared to mesophilic comammox strains presumably confer tolerance to thermal stress. These novel insights into the occurrence, metabolic activity, and adaptation of comammox Nitrospira in thermal habitats further expand our understanding of the global distribution of comammox Nitrospira and have significant implications for how these unique microorganisms have evolved thermal tolerance strategies.

High speciation rate of niche specialists in hot springs

Ecological and evolutionary processes simultaneously regulate microbial diversity, but the evolutionary processes and their driving forces remain largely unexplored. Here we investigated the ecological and evolutionary characteristics of microbiota in hot springs spanning a broad temperature range (54.8-80 °C) by sequencing the 16S rRNA genes. Our results demonstrated that niche specialists and niche generalists are embedded in a complex interaction of ecological and evolutionary dynamics. On the thermal tolerance niche axis, thermal (T) sensitive (at a specific temperature) versus T-resistant (at least in five temperatures) species were characterized by different niche breadth, community abundance and dispersal potential, consequently differing in potential evolutionary trajectory. The niche-specialized T-sensitive species experienced strong temperature barriers, leading to completely species shift and high fitness but low abundant communities at each temperature ("home niche"), and such trade-offs thus reinforced peak performance, as evidenced by high speciation across temperatures and increasing diversification potential with temperature. In contrast, T-resistant species are advantageous of niche expansion but with poor local performance, as shown by wide niche breadth with high extinction, indicating these niche generalists are "jack-of-all-trades, master-of-none". Despite of such differences, the T-sensitive and T-resistant species are evolutionarily interacted. Specifically, the continuous transition from T-sensitive to T-resistant species insured the exclusion probability of T-resistant species at a relatively constant level across temperatures. The co-evolution and co-adaptation of T-sensitive and T-resistant species were in line with the red queen theory. Collectively, our findings demonstrate that high speciation of niche specialists could alleviate the environmental-filtering-induced negative effect on diversity.

Inference and reconstruction of the heimdallarchaeial ancestry of eukaryotes

In the ongoing debates about eukaryogenesis-the series of evolutionary events leading to the emergence of the eukaryotic cell from prokaryotic ancestors-members of the Asgard archaea play a key part as the closest archaeal relatives of eukaryotes1. However, the nature and phylogenetic identity of the last common ancestor of Asgard archaea and eukaryotes remain unresolved2-4. Here we analyse distinct phylogenetic marker datasets of an expanded genomic sampling of Asgard archaea and evaluate competing evolutionary scenarios using state-of-the-art phylogenomic approaches. We find that eukaryotes are placed, with high confidence, as a well-nested clade within Asgard archaea and as a sister lineage to Hodarchaeales, a newly proposed order within Heimdallarchaeia. Using sophisticated gene tree and species tree reconciliation approaches, we show that analogous to the evolution of eukaryotic genomes, genome evolution in Asgard archaea involved significantly more gene duplication and fewer gene loss events compared with other archaea. Finally, we infer that the last common ancestor of Asgard archaea was probably a thermophilic chemolithotroph and that the lineage from which eukaryotes evolved adapted to mesophilic conditions and acquired the genetic potential to support a heterotrophic lifestyle. Our work provides key insights into the prokaryote-to-eukaryote transition and a platform for better understanding the emergence of cellular complexity in eukaryotic cells.

Thermocaproicibacter melissae gen. nov., sp. nov., a thermophilic chain-elongating bacterium, producing n-caproate from polymeric carbohydrates

Strain MDTJ8^T is a chain-elongating thermophilic bacterium isolated from a thermophilic acidogenic anaerobic digestor treating human waste while producing the high commodity chemical n-caproate. The strain grows and produces formate, acetate, n-butyrate, n-caproate and lactate from mono-, di- and polymeric saccharides at 37-60 °C (optimum, 50-55 °C) and at pH 5.0-7.0 (optimum, pH 6.5). The organism is an obligate anaerobe, is motile and its cells form rods (0.3-0.5×1.0-3.0 µm) that stain Gram-positive and occur primarily as chains. Phylogenetic analysis of both the 16S rRNA gene and full genome sequence shows that strain MDTJ8^T belongs to a group that consists of mesophylic chain-elongating bacteria within the family Oscillospiraceae, being nearest to Caproicibacter fermentans EA1^T (94.8 %) and Caproiciproducens galactitolivorans BS-1^T (93.7 %). Its genome (1.96 Mbp) with a G+C content of 49.6 mol% is remarkably smaller than those of other chain-elongating bacteria of the family Oscillospiraceae. Pairwise average nucleotide identity and DNA-DNA hybridization values between strain MDJT8^T and its mesophilic family members are less than 70 and 35 %, respectively, while pairwise average amino acid identity values are less than 68 %. In addition, strain MDJT8^T uses far less carbohydrate and non-carbohydrate substrates compared to its nearest family members. The predominant cellular fatty acids of strain MDTJ8^T are C_14 : 0, C_14 : 0 DMA (dimethyl acetal) and C_16 : 0, while its polar lipid profile shows three unidentified glycophospholipids, 11 glycolipids, 13 phospholipids and six unidentified lipids. No respiratory quinones and polyamines are detected. Based on its phylogenetic, genotypic, morphological, physiological, biochemical and chemotaxonomic characteristics, strain MDTJ8^T represents a novel species and novel genus of the family Oscillospiraceae and Thermocaproicibacter melissae gen. nov., sp. nov. is proposed as its name. The type strain is MDTJ8^T (=DSM 114174^T=LMG 32615^T=NCCB 100883^T).

Metagenomic Discovery of "Candidatus Parvarchaeales"-Related Lineages Sheds Light on Adaptation and Diversification from Neutral-Thermal to Acidic-Mesothermal Environments

"Candidatus Parvarchaeales" microbes, representing a DPANN archaeal group with limited metabolic potential and reliance on hosts for their growth, were initially found in acid mine drainage (AMD). Due to the lack of representatives, however, their ecological roles and adaptation to extreme habitats such as AMD as well as how they diverge across the lineage remain largely unexplored. By applying genome-resolved metagenomics, 28 Parvarchaeales-associated metagenome-assembled genomes (MAGs) representing two orders and five genera were recovered. Among them, we identified three new genera and proposed the names "Candidatus Jingweiarchaeum," "Candidatus Haiyanarchaeum," and "Candidatus Rehaiarchaeum," with the former two belonging to a new order, "Candidatus Jingweiarchaeales." Further analyses of the metabolic potentials revealed substantial niche differentiation between Jingweiarchaeales and Parvarchaeales. Jingweiarchaeales may rely on fermentation, salvage pathways, partial glycolysis, and the pentose phosphate pathway (PPP) for energy conservation reservation, while the metabolic potentials of Parvarchaeales might be more versatile. Comparative genomic analyses suggested that Jingweiarchaeales favor habitats with higher temperatures and that Parvarchaeales are better adapted to acidic environments. We further revealed that the thermal adaptation of these lineages, especially Haiyanarchaeum, might rely on genomic features such as the usage of specific amino acids, genome streamlining, and hyperthermophile featured genes such as rgy. Notably, the adaptation of Parvarchaeales to acidic environments was possibly driven by horizontal gene transfer (HGT). The reconstruction of ancestral states demonstrated that both may have originated from thermal and neutral environments and later spread to mesothermal and acidic environments. These evolutionary processes may also be accompanied by adaptation to oxygen-rich environments via HGT. IMPORTANCE "Candidatus Parvarchaeales" microbes may represent a lineage uniquely distributed in extreme environments such as AMD and hot springs. However, little is known about the strategies and processes of how they adapted to these extreme environments. By the discovery of potential new order-level lineages, "Ca. Jingweiarchaeales," and in-depth comparative genomic analysis, we unveiled the functional differentiation of these lineages. Furthermore, we show that the adaptation of these lineages to high-temperature and acidic environments was driven by different strategies, with the former relying more on genomic characteristics such as genome streamlining and amino acid compositions and the latter relying more on the acquisition of genes associated with acid tolerance. Finally, by the reconstruction of the ancestral states of the optimal growth temperature (OGT) and isoelectric point (pI), we showed the potential evolutionary process of Parvarchaeales-related lineages with regard to the shift from the high-temperature environment of their common ancestors to low-temperature (potentially acidic) environments.

Abiotic selection of microbial genome size in the global ocean

Strong purifying selection is considered a major evolutionary force behind small microbial genomes in the resource-poor photic ocean. However, very little is currently known about how the size of prokaryotic genomes evolves in the global ocean and whether patterns reflect shifts in resource availability in the epipelagic and relatively stable deep-sea environmental conditions. Using 364 marine microbial metagenomes, we investigate how the average genome size of uncultured planktonic prokaryotes varies across the tropical and polar oceans to the hadal realm. We find that genome size is highest in the perennially cold polar ocean, reflecting elongation of coding genes and gene dosage effects due to duplications in the interior ocean microbiome. Moreover, the rate of change in genome size due to temperature is 16-fold higher than with depth up to 200 m. Our results demonstrate how environmental factors can influence marine microbial genome size selection and ecological strategies of the microbiome.

Warmer and drier ecosystems select for smaller bacterial genomes in global soils

Bacterial genome size reflects bacterial evolutionary processes and metabolic lifestyles, with implications for microbial community assembly and ecosystem functions. However, to understand the extent of genome-mediated microbial responses to environmental selections, we require studies that observe genome size distributions along environmental gradients representing different conditions that soil bacteria normally encounter. In this study, we used surface soils collected from 237 sites across the globe and analyzed how environmental conditions (e.g., soil carbon and nutrients, aridity, pH, and temperature) affect soil bacterial occurrences and genome size at the community level using bacterial community profiling. We used a joint species distribution model to quantify the effects of environments on species occurrences and found that aridity was a major regulator of genome size with warmer and drier environments selecting bacteria with smaller genomes. Drought-induced physiological constraints on bacterial growth (e.g., water scarcity for cell metabolisms) may have led to these correlations. This finding suggests that increasing cover by warmer and drier ecosystems may result in bacterial genome simplifications by a reduction of genome size.

Evolutionary diversification of methanotrophic ANME-1 archaea and their expansive virome

'Candidatus Methanophagales' (ANME-1) is an order-level clade of archaea responsible for anaerobic methane oxidation in deep-sea sediments. The diversity, ecology and evolution of ANME-1 remain poorly understood. In this study, we use metagenomics on deep-sea hydrothermal samples to expand ANME-1 diversity and uncover the effect of virus-host dynamics. Phylogenetic analyses reveal a deep-branching, thermophilic family, 'Candidatus Methanospirareceae', closely related to short-chain alkane oxidizers. Global phylogeny and near-complete genomes show that hydrogen metabolism within ANME-1 is an ancient trait that was vertically inherited but differentially lost during lineage diversification. Metagenomics also uncovered 16 undescribed virus families so far exclusively targeting ANME-1 archaea, showing unique structural and replicative signatures. The expansive ANME-1 virome contains a metabolic gene repertoire that can influence host ecology and evolution through virus-mediated gene displacement. Our results suggest an evolutionary continuum between anaerobic methane and short-chain alkane oxidizers and underscore the effects of viruses on the dynamics and evolution of methane-driven ecosystems.

Isolation and characterization of a thermophilic chain elongating bacterium that produces the high commodity chemical n-caproate from polymeric carbohydrates

A thermophilic chain elongating bacterium, strain MDTJ8, was isolated from a thermophilic acidogenic anaerobic digestor producing n-caproate from human waste, growing optimally at 50-55 °C and pH 6.5. 16S rRNA gene analysis suggests that MDTJ8 represents a new species/genus within a group currently composed of mesophilic chain elongators of the Oscillospiraceae family. Genome analysis showed that strain MDTJ8 contains homologues of genes encoding for chain elongation and energy conservation but also indicated n-caproate production from carbohydrates including polymeric substances. This was confirmed by culturing experiments in which MDTJ8 converted, at pH 6.5 and 55 °C, mono-, di- and polymeric carbohydrates (starch and hemicellulose) to n-caproate reaching concentrations up to 283 mg/L and accounting for up to 10 % of the measured fermentation products. MDTJ8 is the first axenic organism that thermophilically performs chain elongation, opening doors to understand and intensify thermophilic bioprocesses targeting anaerobic digestion towards the production of the value-added chemical n-caproate.

Growth rate-associated transcriptome reorganization in response to genomic, environmental, and evolutionary interruptions

The genomic, environmental, and evolutionary interruptions caused the changes in bacterial growth, which were stringently associated with changes in gene expression. The growth and gene expression changes remained unclear in response to these interruptions that occurred combinative. As a pilot study, whether and how bacterial growth was affected by the individual and dual interruptions of genome reduction, environmental stress, and adaptive evolution were investigated. Growth assay showed that the presence of the environmental stressors, i.e., threonine and chloramphenicol, significantly decreased the growth rate of the wild-type Escherichia coli, whereas not that of the reduced genome. It indicated a canceling effect in bacterial growth due to the dual interruption of the genomic and environmental changes. Experimental evolution of the reduced genome released the canceling effect by improving growth fitness. Intriguingly, the transcriptome architecture maintained a homeostatic chromosomal periodicity regardless of the genomic, environmental, and evolutionary interruptions. Negative epistasis in transcriptome reorganization was commonly observed in response to the dual interruptions, which might contribute to the canceling effect. It was supported by the changes in the numbers of differentially expressed genes (DEGs) and the enriched regulons and functions. Gene network analysis newly constructed 11 gene modules, one out of which was correlated to the growth rate. Enrichment of DEGs in these modules successfully categorized them into three types, i.e., conserved, responsive, and epistatic. Taken together, homeostasis in transcriptome architecture was essential to being alive, and it might be attributed to the negative epistasis in transcriptome reorganization and the functional differentiation in gene modules. The present study directly connected bacterial growth fitness with transcriptome reorganization and provided a global view of how microorganisms responded to genomic, environmental, and evolutionary interruptions for survival from wild nature.

Compensatory Transcriptional Response of Fischerella thermalis to Thermal Damage of the Photosynthetic Electron Transfer Chain

Key organisms in the environment, such as oxygenic photosynthetic primary producers (photosynthetic eukaryotes and cyanobacteria), are responsible for fixing most of the carbon globally. However, they are affected by environmental conditions, such as temperature, which in turn affect their distribution. Globally, the cyanobacterium Fischerella thermalis is one of the main primary producers in terrestrial hot springs with thermal gradients up to 60 °C, but the mechanisms by which F. thermalis maintains its photosynthetic activity at these high temperatures are not known. In this study, we used molecular approaches and bioinformatics, in addition to photophysiological analyses, to determine the genetic activity associated with the energy metabolism of F. thermalis both in situ and in high-temperature (40 °C to 65 °C) cultures. Our results show that photosynthesis of F. thermalis decays with temperature, while increased transcriptional activity of genes encoding photosystem II reaction center proteins, such as PsbA (D1), could help overcome thermal damage at up to 60 °C. We observed that F. thermalis tends to lose copies of the standard G4 D1 isoform while maintaining the recently described D1INT isoform, suggesting a preference for photoresistant isoforms in response to the thermal gradient. The transcriptional activity and metabolic characteristics of F. thermalis, as measured by metatranscriptomics, further suggest that carbon metabolism occurs in parallel with photosynthesis, thereby assisting in energy acquisition under high temperatures at which other photosynthetic organisms cannot survive. This study reveals that, to cope with the harsh conditions of hot springs, F. thermalis has several compensatory adaptations, and provides emerging evidence for mixotrophic metabolism as being potentially relevant to the thermotolerance of this species. Ultimately, this work increases our knowledge about thermal adaptation strategies of cyanobacteria.

An insight into the mechanisms of homeostasis in extremophiles

The homeostasis of extremophiles is one that is a diamond hidden in the rough. The way extremophiles adapt to their extreme environments gives a clue into the true extent of what is possible when it comes to life. The discovery of new extremophiles is ever-expanding and an explosion of knowledge surrounding their successful existence in extreme environments is obviously perceived in scientific literature. The present review paper aims to provide a comprehensive view on the different mechanisms governing the extreme adaptations of extremophiles, along with insights and discussions on what the limits of life can possibly be. The membrane adaptations that are vital for survival are discussed in detail. It was found that there are many alterations in the genetic makeup of such extremophiles when compared to their mesophilic counterparts. Apart from the several proteins involved, the significance of chaperones, efflux systems, DNA repair proteins and a host of other enzymes that adapt to maintain functionality, are enlisted, and explained. A deeper understanding of the underlying mechanisms could have a plethora of applications in the industry. There are cases when certain microbes can withstand extreme doses of antibiotics. Such microbes accumulate numerous genetic elements (or plasmids) that possess genes for multiple drug resistance (MDR). A deeper understanding of such mechanisms helps in the development of potential approaches and therapeutic schemes for treating pathogen-mediated outbreaks. An in-depth analysis of the parameters - radiation, pressure, temperature, pH value and metal resistance - are discussed in this review, and the key to survival in these precarious niches is described.

Genomic attributes of thermophilic and hyperthermophilic bacteria and archaea

Thermophiles and hyperthermophiles are immensely useful in understanding the evolution of life, besides their utility in environmental and industrial biotechnology. Advancements in sequencing technologies have revolutionized the field of microbial genomics. The massive generation of data enhances the sequencing coverage multi-fold and allows to analyse the entire genomic features of microbes efficiently and accurately. The mandate of a pure isolate can also be bypassed where whole metagenome-assembled genomes and single cell-based sequencing have fulfilled the majority of the criteria to decode various attributes of microbial genomes. A boom has, therefore, been seen in analysing the extremophilic bacteria and archaea using sequence-based approaches. Due to extensive sequence analysis, it becomes easier to understand the gene flow and their evolution among the members of bacteria and archaea. For instance, sequencing unveiled that Thermotoga maritima shares around 24% of genes of archaeal origin. Comparative and functional genomics provide an analytical view to understanding the microbial diversity of thermophilic bacteria and archaea, their interactions with other microbes, their adaptations, gene flow, and evolution over time. In this review, the genomic features of thermophilic bacteria and archaea are dealt with comprehensively.

The microbiome of cryospheric ecosystems

The melting of the cryosphere is among the most conspicuous consequences of climate change, with impacts on microbial life and related biogeochemistry. However, we are missing a systematic understanding of microbiome structure and function across cryospheric ecosystems. Here, we present a global inventory of the microbiome from snow, ice, permafrost soils, and both coastal and freshwater ecosystems under glacier influence. Combining phylogenetic and taxonomic approaches, we find that these cryospheric ecosystems, despite their particularities, share a microbiome with representatives across the bacterial tree of life and apparent signatures of early and constrained radiation. In addition, we use metagenomic analyses to define the genetic repertoire of cryospheric bacteria. Our work provides a reference resource for future studies on climate change microbiology.

The cryosphere includes those parts of Earth where water or soil is frozen, such as snow, ice, glaciers and permafrost soils. Here, the authors present a global inventory of cryospheric microbial communities and their genetic repertoires.

Comparative genomic analysis of Thermus provides insights into the evolutionary history of an incomplete denitrification pathway

Biological denitrification is a crucial process in the nitrogen biogeochemical cycle, and Thermus has been reported to be a significant heterotrophic denitrifier in terrestrial geothermal environments. However, neither the denitrification potential nor the evolutionary history of denitrification genes in the genus Thermus or phylum Deinococcota is well understood. Here, we performed a comparative analysis of 23 Thermus genomes and identified denitrification genes in 15 Thermus strains. We confirmed that Thermus harbors an incomplete denitrification pathway as none of the strains contain the nosZ gene. Ancestral character state reconstructions and phylogenetic analyses showed that narG, nirS, and norB genes were acquired by the last common ancestor of Thermales and were inherited vertically. In contrast, nirK of Thermales was acquired via two distinct horizontal gene transfers from Proteobacteria to the genus Caldithermus and from an unknown donor to the common ancestor of all known Thermus species except Thermus filiformis. This study expands our understanding of the genomic potential for incomplete denitrification in Thermus, revealing a largely vertical evolutionary history of the denitrification pathway in the Thermaceae, and supporting the important role for Thermus as an important heterotrophic denitrifier in geothermal environments.

Translation stalling proline motifs are enriched in slow-growing, thermophilic, and multicellular bacteria

Rapid bacterial growth depends on the speed at which ribosomes can translate mRNA into proteins. mRNAs that encode successive stretches of proline can cause ribosomes to stall, substantially reducing translation speed. Such stalling is especially detrimental for species that must grow and divide rapidly. Here, we focus on di-prolyl motifs (XXPPX) and ask whether their prevalence varies with growth rate. To find out we conducted a broad survey of such motifs in >3000 bacterial genomes across 35 phyla. Indeed, fast-growing species encode fewer motifs than slow-growing species, especially in highly expressed proteins. We also found many di-prolyl motifs within thermophiles, where prolines can help maintain proteome stability. Moreover, bacteria with complex, multicellular lifecycles also encode many di-prolyl motifs. This is especially evident in the slow-growing phylum Myxococcota. Bacteria in this phylum encode many serine-threonine kinases, and many di-prolyl motifs at potential phosphorylation sites within these kinases. Serine-threonine kinases are involved in cell signaling and help regulate developmental processes linked to multicellularity in the Myxococcota. Altogether, our observations suggest that weakened selection on translational rate, whether due to slow or thermophilic growth, may allow di-prolyl motifs to take on new roles in biological processes that are unrelated to translational rate.

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.

Discovery of the non-cosmopolitan lineages in Candidatus Thermoprofundales

The recently proposed order Candidatus Thermoprofundales, currently containing only one family-level lineage Marine Benthic Group-D (MBG-D), is distributed in global subsurface ecosystems and ecologically important, but its diversity, evolution and metabolism remain largely unknown. Here we described two novel family-level specialized lineages in Ca. Thermoprofundales, JdFR-43 and HyVt, which are restricted to specific biotopes (primarily in marine hydrothermal vents and occasionally in oil reservoirs and hot springs) in contrast to the cosmopolitan lineage MBG-D. The comparative genomics revealed that the specialized lineages have streamlined genomes, higher GC contents, enriched genes associated with nucleotide biosynthesis, ribosome biogenesis and DNA repair and additional thermostable aminopeptidases, enabling them to adapt to high-temperature habitats such as marine hydrothermal vents, deep subsurface oil reservoirs and hot springs. On the contrary, the unique metabolic traits of the cosmopolitan MBG-D, motility, glycolysis, butanoate metabolism, secondary metabolites production and additional genes for specific peptides and carbohydrates degradation potentially enhance its response to environmental change. Substrate preference is found for most MAGs across all lineages with the ability to utilize both polysaccharides (chitin and starch) and proteinaceous substances, whereas JdFR-43 members from oil reservoirs can only utilize proteins. These results expand the diversity of Ca. Thermoprofundales significantly and further improve our understandings of the adaptations of Ca. Thermoprofundales to various environments. This article is protected by copyright. All rights reserved.

Environmental stress leads to genome streamlining in a widely distributed species of soil bacteria

Bacteria have highly flexible pangenomes, which are thought to facilitate evolutionary responses to environmental change, but the impacts of environmental stress on pangenome evolution remain unclear. Using a landscape pangenomics approach, I demonstrate that environmental stress leads to consistent, continuous reduction in genome content along four environmental stress gradients (acidity, aridity, heat, salinity) in naturally occurring populations of Bradyrhizobium diazoefficiens (widespread soil-dwelling plant mutualists). Using gene-level network and duplication functional traits to predict accessory gene distributions across environments, genes predicted to be superfluous are more likely lost in high stress, while genes with multi-functional roles are more likely retained. Genes with higher probabilities of being lost with stress contain significantly higher proportions of codons under strong purifying and positive selection. Gene loss is widespread across the entire genome, with high gene-retention hotspots in close spatial proximity to core genes, suggesting Bradyrhizobium has evolved to cluster essential-function genes (accessory genes with multifunctional roles and core genes) in discrete genomic regions, which may stabilise viability during genomic decay. In conclusion, pangenome evolution through genome streamlining are important evolutionary responses to environmental change. This raises questions about impacts of genome streamlining on the adaptive capacity of bacterial populations facing rapid environmental change.

Overlapping genes in natural and engineered genomes

Modern genome-scale methods that identify new genes, such as proteogenomics and ribosome profiling, have revealed, to the surprise of many, that overlap in genes, open reading frames and even coding sequences is widespread and functionally integrated into prokaryotic, eukaryotic and viral genomes. In parallel, the constraints that overlapping regions place on genome sequences and their evolution can be harnessed in bioengineering to build more robust synthetic strains and constructs. With a focus on overlapping protein-coding and RNA-coding genes, this Review examines their discovery, topology and biogenesis in the context of their genome biology. We highlight exciting new uses for sequence overlap to control translation, compress synthetic genetic constructs, and protect against mutation.

Variation in genomic traits of microbial communities among ecosystems

Free-living bacteria in nutrient limited environments often exhibit traits which may reduce the cost of reproduction, such as smaller genome size, low GC content and fewer sigma (σ) factor and 16S rRNA gene copies. Despite the potential utility of these traits to detect relationships between microbial communities and ecosystem-scale properties, few studies have assessed these traits on a community-scale. Here, we analysed these traits from publicly available metagenomes derived from marine, soil, host-associated and thermophilic communities. In marine and thermophilic communities, genome size and GC content declined in parallel, consistent with genomic streamlining, with GC content in thermophilic communities generally higher than in marine systems. In contrast, soil communities averaging smaller genomes featured higher GC content and were often from low-carbon environments, suggesting unique selection pressures in soil bacteria. The abundance of specific σ-factors varied with average genome size and ecosystem type. In oceans, abundance of fliA, a σ-factor controlling flagella biosynthesis, was positively correlated with community average genome size—reflecting known trade-offs between nutrient conservation and chemotaxis. In soils, a high abundance of the stress response σ-factor gene rpoS was associated with smaller average genome size and often located in harsh and/or carbon-limited environments—a result which tracks features observed in culture and indicates an increased capacity for stress response in nutrient-poor soils. This work shows how ecosystem-specific constraints are associated with trade-offs which are embedded in the genomic features of bacteria in microbial communities, and which can be detected at the community level, highlighting the importance of genomic features in microbial community analysis.

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 Reveals Thermal Adaptation and a High Metabolic Diversity in "Candidatus Bathyarchaeia"

"Candidatus Bathyarchaeia" is a phylogenetically diverse and widely distributed lineage often in high abundance in anoxic submarine sediments; however, their evolution and ecological roles in terrestrial geothermal habitats are poorly understood. In the present study, 35 Ca. Bathyarchaeia metagenome-assembled genomes (MAGs) were recovered from hot spring sediments in Tibet and Yunnan, China. Phylogenetic analysis revealed all MAGs of Ca. Bathyarchaeia can be classified into 7 orders and 15 families. Among them, 4 families have been first discovered in the present study, significantly expanding the known diversity of Ca. Bathyarchaeia. Comparative genomics demonstrated Ca. Bathyarchaeia MAGs from thermal habitats to encode a large variety of genes related to carbohydrate degradation, which are likely a metabolic adaptation of these organisms to a lifestyle at high temperatures. At least two families are potential methanogens/alkanotrophs, indicating a potential for the catalysis of short-chain hydrocarbons. Three MAGs from Family-7.3 are identified as alkanotrophs due to the detection of an Mcr complex. Family-2 contains the largest number of genes relevant to alkyl-CoM transformation, indicating the potential for methylotrophic methanogenesis, although their evolutionary history suggests the ancestor of Ca. Bathyarchaeia was unable to metabolize alkanes. Subsequent lineages have acquired the ability via horizontal gene transfer. Overall, our study significantly expands our knowledge and understanding of the metabolic capabilities, habitat adaptations, and evolution of Ca. Bathyarchaeia in thermal environments. IMPORTANCE Ca. Bathyarchaeia MAGs from terrestrial hot spring habitats are poorly revealed, though they have been studied extensively in marine ecosystems. In this study, we uncovered the metabolic capabilities and ecological role of Ca. Bathyarchaeia in hot springs and give a comprehensive comparative analysis between thermal and nonthermal habitats to reveal the thermal adaptability of Ca. Bathyarchaeia. Also, we attempt to determine the evolutionary history of methane/alkane metabolism in Ca. Bathyarchaeia, since it appears to be the first archaea beyond Euryarchaeota which contains the mcrABG genes. The reclassification of Ca. Bathyarchaeia and significant genomic differences among different lineages largely expand our knowledge on these cosmopolitan archaea, which will be beneficial in guiding the future studies.

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.

Decomposition of Microbial Necromass Is Divergent at the Individual Taxonomic Level in Soil

The turnover of microbial biomass plays an important part in providing a significant source of carbon (C) to soil organic C. However, whether the decomposition of microbial necromass (non-living microbial biomass) in the soil varies at the individual taxa level remains largely unknown. To fill up these gaps, we compared the necromass decomposition of bacterial and archaeal taxa by separating live microbial biomass with 18O-stable isotope probing from dead microbial biomass in soil. Our results showed that most of the microbial necromass at the operational taxonomic unit level (88.51%), which mainly belong to Acidobacteria, Actinobacteria, Gemmatimonadetes, and Proteobacteria, decomposed significantly after 30 days. In addition, there were great variations in necromass decomposition within each phylum, such as the decomposition of operational taxonomic units in Proteobacteria that ranged from 51% (Beijerinckia) to 92% (Nitrosospira). More importantly, the necromass decomposition was not related to the chemical composition of the cell wall but might positively correlate with the guanine-cytosine content of DNA and negatively correlated with genome size. This study provided a new insight that the decomposition of microbial necromass in soil was divergent at the individual taxonomic level and could not be fully explained by previously proposed mechanisms.

Genomic adaptations enabling Acidithiobacillus distribution across wide-ranging hot spring temperatures and pHs

BACKGROUND

Terrestrial hot spring settings span a broad spectrum of physicochemistries. Physicochemical parameters, such as pH and temperature, are key factors influencing differences in microbial composition across diverse geothermal areas. Nonetheless, analysis of hot spring pools from the Taupo Volcanic Zone (TVZ), New Zealand, revealed that some members of the bacterial genus, Acidithiobacillus, are prevalent across wide ranges of hot spring pHs and temperatures. To determine the genomic attributes of Acidithiobacillus that inhabit such diverse conditions, we assembled the genomes of 19 uncultivated hot spring Acidithiobacillus strains from six geothermal areas and compared these to 37 publicly available Acidithiobacillus genomes from various habitats.

RESULTS

Analysis of 16S rRNA gene amplicons from 138 samples revealed that Acidithiobacillus comprised on average 11.4 ± 16.8% of hot spring prokaryotic communities, with three Acidithiobacillus amplicon sequence variants (ASVs) (TVZ_G1, TVZ_G2, TVZ_G3) accounting for > 90% of Acidithiobacillus in terms of relative abundance, and occurring in 126 out of 138 samples across wide ranges of temperature (17.5-92.9 °C) and pH (1.0-7.5). We recovered 19 environmental genomes belonging to each of these three ASVs, as well as a fourth related group (TVZ_G4). Based on genome average nucleotide identities, the four groups (TVZ_G1-TVZ_G4) constitute distinct species (ANI < 96.5%) of which three are novel Acidithiobacillus species (TVZ_G2-TVZ_G4) and one belongs to Acidithiobacillus caldus (TVZ_G1). All four TVZ Acidithiobacillus groups were found in hot springs with temperatures above the previously known limit for the genus (up to 40 °C higher), likely due to significantly higher proline and GC contents than other Acidithiobacillus species, which are known to increase thermostability. Results also indicate hot spring-associated Acidithiobacillus have undergone genome streamlining, likely due to thermal adaptation. Moreover, our data suggest that Acidithiobacillus prevalence across varied hot spring pHs is supported by distinct strategies, whereby TVZ_G2-TVZ_G4 regulate pH homeostasis mostly through Na⁺/H⁺ antiporters and proton-efflux ATPases, whereas TVZ_G1 mainly relies on amino acid decarboxylases.

CONCLUSIONS

This study provides insights into the distribution of Acidithiobacillus species across diverse hot spring physichochemistries and determines genomic features and adaptations that potentially enable Acidithiobacillus species to colonize a broad range of temperatures and pHs in geothermal environments. Video Abstract.

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

Deep-sea hydrothermal vents resemble the early Earth, and thus the dominant Thermococcaceae inhabitants, which occupy an evolutionarily basal position of the archaeal tree and take an obligate anaerobic hyperthermophilic free-living lifestyle, are likely excellent models to study the evolution of early life. Here, we determined that unbiased mutation rate of a representative species, Thermococcus eurythermalis, exceeded that of all known free-living prokaryotes by 1-2 orders of magnitude, and thus rejected the long-standing hypothesis that low mutation rates were selectively favored in hyperthermophiles. We further sequenced multiple and diverse isolates of this species and calculated that T. eurythermalis has a lower effective population size than other free-living prokaryotes by 1-2 orders of magnitude. These data collectively indicate that the high mutation rate of this species is not selectively favored but instead driven by random genetic drift. The availability of these unusual data also helps explore mechanisms underlying microbial genome size evolution. We showed that genome size is negatively correlated with mutation rate and positively correlated with effective population size across 30 bacterial and archaeal lineages, suggesting that increased mutation rate and random genetic drift are likely two important mechanisms driving microbial genome reduction. Future determinations of the unbiased mutation rate of more representative lineages with highly reduced genomes such as Prochlorococcus and Pelagibacterales that dominate marine microbial communities are essential to test these hypotheses.