HGTector: an automated method facilitating genome-wide discovery of putative horizontal gene transfers

GutMetaNet: an integrated database for exploring horizontal gene transfer and functional redundancy in the human gut microbiome

Metagenomic studies have revealed the critical roles of complex microbial interactions, including horizontal gene transfer (HGT) and functional redundancy (FR), in shaping the gut microbiome's functional capacity and resilience. However, the lack of comprehensive data integration and systematic analysis approaches has limited the in-depth exploration of HGT and FR dynamics across large-scale gut microbiome datasets. To address this gap, we present GutMetaNet (https://gutmetanet.deepomics.org/), a first-of-its-kind database integrating extensive human gut microbiome data with comprehensive HGT and FR analyses. GutMetaNet contains 21 567 human gut metagenome samples with whole-genome shotgun sequencing data related to various health conditions. Through systematic analysis, we have characterized the taxonomic profiles and FR profiles, and identified 14 636 HGT events using a shared reference genome database across the collected samples. These HGT events have been curated into 8049 clusters, which are annotated with categorized mobile genetic elements, including transposons, prophages, integrative mobilizable elements, genomic islands, integrative conjugative elements and group II introns. Additionally, GutMetaNet incorporates automated analyses and visualizations for the HGT events and FR, serving as an efficient platform for in-depth exploration of the interactions among gut microbiome taxa and their implications for human health.

Reticulate evolution: Detection and utility in the phylogenomics era

Phylogenomics has enriched our understanding that the Tree of Life can have network-like or reticulate structures among some taxa and genes. Two non-vertical modes of evolution - hybridization/introgression and horizontal gene transfer - deviate from a strictly bifurcating tree model, causing non-treelike patterns. However, these reticulate processes can produce similar patterns to incomplete lineage sorting or recombination, potentially leading to ambiguity. Here, we present a brief overview of a phylogenomic workflow for inferring organismal histories and compare methods for distinguishing modes of reticulate evolution. We discuss how the timing of coalescent events can help disentangle introgression from incomplete lineage sorting and how horizontal gene transfer events can help determine the relative timing of speciation events. In doing so, we identify pitfalls of certain methods and discuss how to extend their utility across the Tree of Life. Workflows, methods, and future directions discussed herein underscore the need to embrace reticulate evolutionary patterns for understanding the timing and rates of evolutionary events, providing a clearer view of life's history.

Comprehensive genomic analysis of Brevibacillus brevis BF19 reveals its biocontrol potential against bitter gourd wilt

Bitter gourd wilt, a severe vascular disease triggered by the soilborne pathogen Fusarium oxysporum f. sp. momordicae (FOM), markedly constrains bitter gourd yield. In this study, a novel strain BF19 of Brevibacillus brevis was isolated and identified, exhibiting strong antimicrobial activity against FOM through in vivo and in vitro experiments. To comprehensively assess the biocontrol potential of strain BF19, we conducted phenotypic, phylogenetic, and comparative genomics analyses. Phenotypic analysis revealed that BF19 exhibited 53.33% biocontrol efficacy and significantly increased the average plant height, root fresh weight, and dry weight. Whole-genome sequencing and comparative genomic analysis revealed numerous potential genes associated with biocontrol mechanisms in BF19. Importantly, the integration of metabolic cluster prediction with liquid chromatography‒tandem mass spectrometry (LC‒MS/MS) revealed the presence of a macrobrevin antibiotic, a product of polyketide synthases (PKSs), predominantly in BF19 fermentation products. The effectiveness of the Br. brevis strain BF19 and its crude extract against bitter gourd wilt has also been confirmed. This study provides a genetic framework for future investigations on PKSs and establishes a scientific basis for optimizing field applications of microbial biopesticides derived from Br. brevis BF19.

Unveiling intraspecific diversity and evolutionary dynamics of the foodborne pathogen Bacillus paranthracis through high-quality pan-genome analysis

Understanding the evolutionary dynamics of foodborne pathogens throughout host-associated habitats is of utmost importance. Bacterial pan-genomes, as dynamic entities, are strongly influenced by ecological lifestyles. As a phenotypically diverse species in the Bacillus cereus group, Bacillus paranthracis is recognized as an emerging foodborne pathogen and a probiotic simultaneously. This poorly understood species is a suitable study model for adaptive pan-genome evolution. In this study, we determined the biogeographic distribution, abundance, genetic diversity, and genotypic profiles of key genetic elements of B. paranthracis. Metagenomic read recruitment analyses demonstrated that B. paranthracis members are globally distributed and abundant in host-associated habitats. A high-quality pan-genome of B. paranthracis was subsequently constructed to analyze the evolutionary dynamics involved in ecological adaptation comprehensively. The open pan-genome indicated a flexible gene repertoire with extensive genetic diversity. Significant divergences in the phylogenetic relationships, functional enrichment, and degree of selective pressure between the different components demonstrated different evolutionary dynamics between the core and accessory genomes driven by ecological forces. Purifying selection and gene loss are the main signatures of evolutionary dynamics in B. paranthracis pan-genome. The plasticity of the accessory genome is characterized by horizontal gene transfer (HGT), massive gene losses, and weak purifying or positive selection, which might contribute to niche-specific adaptation. In contrast, although the core genome dominantly undergoes purifying selection, its association with HGT and positively selected mutations indicates its potential role in ecological diversification. Furthermore, host fitness-related dynamics are characterized by the loss of secondary metabolite biosynthesis gene clusters (BGCs) and CAZyme-encoding genes and the acquisition of antimicrobial resistance (AMR) and virulence genes via HGT. This study offers a case study of pan-genome evolution to investigate the ecological adaptations reflected by biogeographical characteristics, thereby advancing the understanding of intraspecific diversity and evolutionary dynamics of foodborne pathogens.

Distinct microbial nitrogen cycling processes in the deepest part of the ocean

The Mariana Trench (MT) is the deepest part of the ocean on Earth. Previous studies have described the microbial community structures and functional potential in the seawater and surface sediment of MT. Still, the metabolic features and adaptation strategies of the microorganisms involved in nitrogen cycling processes are poorly understood. In this study, comparative metagenomic approaches were used to study microbial nitrogen cycling in three MT habitats, including hadal seawater [9,600-10,500 m below sea level (mbsl)], surface sediments [0-46 cm below seafloor (cmbsf) at a water depth between 7,143 and 8,638 mbsl], and deep sediments (200-306 cmbsf at a water depth of 8,300 mbsl). We identified five new nitrite-oxidizing bacteria (NOB) lineages that had adapted to the oligotrophic MT slope sediment, via their CO2 fixation capability through the reductive tricarboxylic acid (rTCA) or Calvin-Benson-Bassham (CBB) cycle; an anammox bacterium might perform aerobic respiration and utilize sedimentary carbohydrates for energy generation because it contains genes encoding type A cytochrome c oxidase and complete glycolysis pathway. In seawater, abundant alkane-oxidizing Ketobacter species can fix inert N2 released from other denitrifying and/or anammox bacteria. This study further expands our understanding of microbial life in the largely unexplored deepest part of the ocean.

IMPORTANCE

The metabolic features and adaptation strategies of the nitrogen cycling microorganisms in the deepest part of the ocean are largely unknown. This study revealed that anammox bacteria might perform aerobic respiration in response to nutrient limitation or O2 fluctuations in the Mariana Trench sediments. Meanwhile, an abundant alkane-oxidizing Ketobacter species could fix N2 in hadal seawater. This study provides new insights into the roles of hadal microorganisms in global nitrogen biogeochemical cycles. It substantially expands our understanding of the microbial life in the largely unexplored deepest part of the ocean.

Multi-omics signatures reveal genomic and functional heterogeneity of Cutibacterium acnes in normal and diseased skin

Cutibacterium acnes is the most abundant bacterium of the human skin microbiome since adolescence, participating in both skin homeostasis and diseases. Here, we demonstrate individual and niche heterogeneity of C. acnes from 1,234 isolate genomes. Skin disease (atopic dermatitis and acne) and body site shape genomic differences of C. acnes, stemming from horizontal gene transfer and selection pressure. C. acnes harbors characteristic metabolic functions, fewer antibiotic resistance genes and virulence factors, and a more stable genome compared with Staphylococcus epidermidis. Integrated genome, transcriptome, and metabolome analysis at the strain level unveils the functional characteristics of C. acnes. Consistent with the transcriptome signature, C. acnes in a sebum-rich environment induces toxic and pro-inflammatory effects on keratinocytes. L-carnosine, an anti-oxidative stress metabolite, is up-regulated in the C. acnes metabolome from atopic dermatitis and attenuates skin inflammation. Collectively, our study reveals the joint impact of genes and the microenvironment on C. acnes function.

Pan-genome insights into adaptive evolution of bacterial symbionts in mixed host-microbe symbioses represented by human gut microbiota Bacteroides cellulosilyticus

Animal hosts harbor diverse assemblages of microbial symbionts that play crucial roles in the host's lifestyle. The link between microbial symbiosis and host development remains poorly understood. In particular, little is known about the adaptive evolution of gut bacteria in host-microbe symbioses. Recently, symbiotic relationships have been categorized as open, closed, or mixed, reflecting their modes of inter-host transmission and resulting in distinct genomic features. Members of the genus Bacteroides are the most abundant human gut microbiota and possess both probiotic and pathogenic potential, providing an excellent model for studying pan-genome evolution in symbiotic systems. Here, we determined the complete genome of an novel clinical strain PL2022, which was isolated from a blood sample and performed pan-genome analyses on a representative set of Bacteroides cellulosilyticus strains to quantify the influence of the symbiotic relationship on the evolutionary dynamics. B. cellulosilyticus exhibited correlated genomic features with both open and closed symbioses, suggesting a mixed symbiosis. An open pan-genome is characterized by abundant accessory gene families, potential horizontal gene transfer (HGT), and diverse mobile genetic elements (MGEs), indicating an innovative gene pool, mainly associated with genomic islands and plasmids. However, massive parallel gene loss, weak purifying selection, and accumulation of positively selected mutations were the main drivers of genome reduction in B. cellulosilyticus. Metagenomic read recruitment analyses showed that B. cellulosilyticus members are globally distributed and active in human gut habitats, in line with predominant vertical transmission in the human gut. However, existence and/or high abundance were also detected in non-intestinal tissues, other animal hosts, and non-host environments, indicating occasional horizontal transmission to new niches, thereby creating arenas for the acquisition of novel genes. This case study of adaptive evolution under a mixed host-microbe symbiosis advances our understanding of symbiotic pan-genome evolution. Our results highlight the complexity of genetic evolution in this unusual intestinal symbiont.

Comprehensive analysis of genomic variation, pan-genome and biosynthetic potential of Corynebacterium glutamicum strains

Corynebacterium glutamicum is a non-pathogenic species of the Corynebacteriaceae family. It has been broadly used in industrial biotechnology for the production of valuable products. Though it is widely accepted at the industrial level, knowledge about the genomic diversity of the strains is limited. Here, we investigated the comparative genomic features of the strains and pan-genomic characteristics. We also observed phylogenetic relationships among the strains based on average nucleotide identity (ANI). We found diversity between strains at the genomic and pan-genomic levels. Less than one-third of the C. glutamicum pan-genome consists of core genes and soft-core genes. Whereas, a large number of strain-specific genes covered about half of the total pan-genome. Besides, C. glutamicum pan-genome is open and expanding, which indicates the possible addition of new gene families to the pan-genome. We also investigated the distribution of biosynthetic gene clusters (BGCs) among the strains. We discovered slight variations of BGCs at the strain level. Several BGCs with the potential to express novel bioactive secondary metabolites have been identified. Therefore, by utilizing the characteristic advantages of C. glutamicum, different strains can be potential applicants for natural drug discovery.

Speciation Features of Ferdinandcohnia quinoae sp. nov to Adapt to the Plant Host

The bacterial strain SECRCQ15T was isolated from seeds of Chenopodium quinoa in Spain. Phylogenetic, chemotaxonomic, and phenotypic analyses, as well as genome similarity indices, support the classification of the strain into a novel species of the genus Ferdinandcohnia, for which we propose the name Ferdinandcohnia quinoae sp. nov. To dig deep into the speciation features of the strain SECRCQ15T, we performed a comparative genomic analysis of the genome of this strain and those of the type strains of species from the genus Ferdinandcohnia. We found several genes related with plant growth-promoting mechanisms within the SECRCQ15T genome. We also found that singletons of F. quinoae SECRCQ15T are mainly related to the use of carbohydrates, which is a common trait of plant-associated bacteria. To further reveal speciation events in this strain, we revealed genes undergoing diversifying selection (e.g., genes encoding ribosomal proteins) and functions likely lost due to pseudogenization. Also, we found that this novel species contains 138 plant-associated gene-cluster functions that are unique within the genus Ferdinandcohnia. These features may explain both the ecological and taxonomical differentiation of this new taxon.

The global speciation continuum of the cyanobacterium Microcoleus

Speciation is a continuous process driven by genetic, geographic, and ecological barriers to gene flow. It is widely investigated in multicellular eukaryotes, yet we are only beginning to comprehend the relative importance of mechanisms driving the emergence of barriers to gene flow in microbial populations. Here, we explored the diversification of the nearly ubiquitous soil cyanobacterium Microcoleus. Our dataset consisted of 291 genomes, of which 202 strains and eight herbarium specimens were sequenced for this study. We found that Microcoleus represents a global speciation continuum of at least 12 lineages, which radiated during Eocene/Oligocene aridification and exhibit varying degrees of divergence and gene flow. The lineage divergence has been driven by selection, geographical distance, and the environment. Evidence of genetic divergence and selection was widespread across the genome, but we identified regions of exceptional differentiation containing candidate genes associated with stress response and biosynthesis of secondary metabolites.

Pangenome analysis of Shewanella xiamenensis revealed important genetic traits concerning genetic diversity, pathogenicity and antibiotic resistance

BACKGROUND

Shewanella xiamenensis, widely distributed in natural environments, has long been considered as opportunistic pathogen. Recently, significant changes in the resistance spectrum have been observed in S. xiamenensis, due to acquired antibiotic resistance genes. Therefore, a pan-genome analysis was conducted to illuminate the genomic changes in S. xiamenensis.

RESULTS

Phylogenetic analysis revealed three major clusters and three singletons, among which close relationship between several strains was discovered, regardless of their host and niches. The "open" genomes with diversity of accessory and strain-specific genomes took advantage towards diversity environments. The purifying selection pressure was the main force on genome evolution, especially in conservative genes. Only 53 gene families were under positive selection pressure. Phenotypic resistance analysis revealed 21 strains were classified as multi-drug resistance (MDR). Ten types of antibiotic resistance genes and two heavy metal resistance operons were discovered in S. xiamenensis. Mobile genetic elements and horizontal gene transfer increased genome diversity and were closely related to MDR strains. S. xiamenensis carried a variety of virulence genes and macromolecular secretion systems, indicating their important roles in pathogenicity and adaptability. Type IV secretion system was discovered in 15 genomes with various sequence structures, indicating it was originated from different donors through horizontal gene transfer.

CONCLUSIONS

This study provided with a detailed insight into the changes in the pan-genome of S. xiamenensis, highlighting its capability to acquire new mobile genetic elements and resistance genes for its adaptation to environment and pathogenicity to human and animals.

Genomic analysis of Ralstonia pickettii reveals the genetic features for potential pathogenicity and adaptive evolution in drinking water

Ralstonia pickettii, the most critical clinical pathogen of the genus Ralstonia, has been identified as a causative agent of numerous harmful infections. Additionally, Ralstonia pickettii demonstrates adaptability to extreme environmental conditions, such as those found in drinking water. In this study, we conducted a comprehensive genomic analysis to investigate the genomic characteristics related to potential pathogenicity and adaptive evolution in drinking water environments of Ralstonia pickettii. Through phylogenetic analysis and population genetic analysis, we divided Ralstonia pickettii into five Groups, two of which were associated with drinking water environments. The open pan-genome with a large and flexible gene repertoire indicated a high genetic plasticity. Significant differences in functional enrichment were observed between the core- and pan-genome of different groups. Diverse mobile genetic elements (MGEs), extensive genomic rearrangements, and horizontal gene transfer (HGT) events played a crucial role in generating genetic diversity. In drinking water environments, Ralstonia pickettii exhibited strong adaptability, and the acquisition of specific adaptive genes was potentially facilitated by genomic islands (GIs) and HGT. Furthermore, environmental pressures drove the adaptive evolution of Ralstonia pickettii, leading to the accumulation of unique mutations in key genes. These mutations may have a significant impact on various physiological functions, particularly carbon metabolism and energy metabolism. The presence of virulence-related elements associated with macromolecular secretion systems, virulence factors, and antimicrobial resistance indicated the potential pathogenicity of Ralstonia pickettii, making it capable of causing multiple nosocomial infections. This study provides comprehensive insights into the potential pathogenicity and adaptive evolution of Ralstonia pickettii in drinking water environments from a genomic perspective.

Detecting horizontal gene transfer among microbiota: an innovative pipeline for identifying co-shared genes within the mobilome through advanced comparative analysis

Horizontal gene transfer (HGT) is a key driver in the evolution of bacterial genomes. The acquisition of genes mediated by HGT may enable bacteria to adapt to ever-changing environmental conditions. Long-term application of antibiotics in intensive agriculture is associated with the dissemination of antibiotic resistance genes among bacteria with the consequences causing public health concern. Commensal farm-animal-associated gut microbiota are considered the reservoir of the resistance genes. Therefore, in this study, we identified known and not-yet characterized mobilized genes originating from chicken and porcine fecal samples using our innovative pipeline followed by network analysis to provide appropriate visualization to support proper interpretation.

Comamonas resistens sp. nov. and Pseudomonas triclosanedens sp. nov., two members of the phylum Pseudomonadota isolated from the wastewater treatment system of a pharmaceutical factory

Five strains of two novel species were isolated from the wastewater treatment systems of a pharmaceutical factory located in Zhejiang province, PR China. Strains ZM22T and Y6 were identified as belonging to a potential novel species of the genus Comamonas, whereas strains ZM23T, ZM24 and ZM25 were identified as belonging to a novel species of the genus Pseudomonas. These strains were characterized by polyphasic approaches including 16S rRNA gene analysis, multi-locus sequence analysis, average nucleotide identity (ANI), in silico DNA-DNA hybridization (isDDH), physiological and biochemical tests, as well as chemotaxonomic analysis. Genome-based phylogenetic analysis further confirmed that strains ZM22T and Y6 form a distinct clade closely related to Comamonas testosteroni ATCC 11996T and Comamonas thiooxydans DSM 17888T. Strains ZM23T, ZM24 and ZM25 were grouped as a separate clade closely related to Pseudomonas nitroreducens DSM 14399T and Pseudomonas nicosulfuronedens LAM1902T. The orthoANI and isDDH results indicated that strains ZM22T and Y6 belong to the same species. In addition, genomic DNA fingerprinting demonstrated that these strains do not originate from a single clone. The same results were observed for strains ZM23T, ZM24 and ZM25. Strains ZM22T and Y6 were resistant to multiple antibiotics, whereas strains ZM23T, ZM24 and ZM25 were able to degrade an emerging pollutant, triclosan. The phylogenetic, physiological and biochemical characteristics, as well as chemotaxonomy, allowed these strains to be distinguished from their genus, and we therefore propose the names Comamonas resistens sp. nov. (type strain ZM22=MCCC 1K08496T=KCTC 82561T) and Pseudomonas triclosanedens sp. nov. (type strain ZM23T=MCCC 1K08497T=JCM 36056T), respectively.

Prediction and characterization of prophages of Stenotrophomonas maltophilia reveals a remarkable phylogenetic diversity of prophages

Prophages, which enables bacterial hosts to acquire novel traits, and increase genetic variation and evolutionary innovation, are considered to be one of the greatest drivers of bacterial diversity and evolution. Stenotrophomonas maltophilia is widely distributed and one of the most important multidrug resistant bacteria in hospitals. However, the distribution and genetic diversity of S. maltophilia prophages have not been elucidated. In this study, putative prophages were predicted in S. maltophilia genomes by using virus prediction tools, and the genetic diversity and phylogeny of S. maltophilia and the prophages they harbor were further analyzed. A total of 356 prophage regions were predicted from 88 S. maltophilia genomes. Among them, 144 were intact prophages, but 77.09% of the intact prophages did not match any known phage sequences in the public database. The number of prophage carried by S. maltophilia is related to its host habitat and is an important factor affecting the size of the host genome, but it is not related to the genetic diversity of the prophage. The prediction of auxiliary genes encoded by prophage showed that antibiotic resistance genes was not predicted for any of the prophages except for one questionable prophage, while 53 virulence genes and 169 carbohydrate active enzymes were predicted from 11.24 and 44.1% prophages, respectively. Most of the prophages (72.29%) mediated horizontal gene transfer of S. maltophilia genome, but only involved in 6.25% of the horizontal gene transfer events. In addition, CRISPR prediction indicated 97.75% S. maltophilia strains contained the CRISPR-Cas system containing 818 spacer sequences. However, these spacer sequences did not match any known S. maltophilia phages, and only a few S. maltophilia prophages. Comparative genomic analysis revealed a highly conserved and syntenic organization with genomic rearrangement between the prophages and the known related S. maltophilia phages. Our results indicate a high prevalence and genetic diversity of prophages in the genome of S. maltophilia, as well as the presence of a large number of uncharacterized phages. It provides an important complement to understanding the diversity and biological characteristics of phages, as well as the interactions and evolution between bacteria and phages.

Comparative genomic analysis of Chryseobacterium species: deep insights into plant-growth-promoting and halotolerant capacities

Members of the genus Chryseobacterium have attracted great interest as beneficial bacteria that can promote plant growth and biocontrol. Given the recent risks of climate change, it is important to develop tolerance strategies for efficient applications of plant-beneficial bacteria in saline environments. However, the genetic determinants of plant-growth-promoting and halotolerance effects in Chryseobacterium have not yet been investigated at the genomic level. Here, a comparative genomic analysis was conducted with seven Chryseobacterium species. Phylogenetic and phylogenomic analyses revealed niche-specific evolutionary distances between soil and freshwater Chryseobacterium species, consistent with differences in genomic statistics, indicating that the freshwater bacteria have smaller genome sizes and fewer genes than the soil bacteria. Phosphorus- and zinc-cycling genes (required for nutrient acquisition in plants) were universally present in all species, whereas nitrification and sulphite reduction genes (required for nitrogen- and sulphur-cycling, respectively) were distributed only in soil bacteria. A pan-genome containing 6842 gene clusters was constructed, which reflected the general features of the core, accessory and unique genomes. Halotolerant species with an accessory genome shared a Kdp potassium transporter and biosynthetic pathways for branched-chain amino acids and the carotenoid lycopene, which are associated with countermeasures against salt stress. Protein-protein interaction network analysis was used to define the genetic determinants of Chryseobacterium salivictor NBC122 that reduce salt damage in bacteria and plants. Sixteen hub genes comprised the aromatic compound degradation and Por secretion systems, which are required to cope with complex stresses associated with saline environments. Horizontal gene transfer and CRISPR-Cas analyses indicated that C. salivictor NBC122 underwent more evolutionary events when interacting with different environments. These findings provide deep insights into genomic adaptation to dynamic interactions between plant-growth-promoting Chryseobacterium and salt stress.

Insights into early evolutionary adaptations of the Akkermansia genus to the vertebrate gut

Akkermansia, a relevant mucin degrader from the vertebrate gut microbiota, is a member of the deeply branched Verrucomicrobiota, as well as the only known member of this phylum to be described as inhabitants of the gut. Only a few Akkermansia species have been officially described so far, although there is genomic evidence addressing the existence of more species-level variants for this genus. This niche specialization makes Akkermansia an interesting model for studying the evolution of microorganisms to their adaptation to the gastrointestinal tract environment, including which kind of functions were gained when the Akkermansia genus originated or how the evolutionary pressure functions over those genes. In order to gain more insight into Akkermansia adaptations to the gastrointestinal tract niche, we performed a phylogenomic analysis of 367 high-quality Akkermansia isolates and metagenome-assembled genomes, in addition to other members of Verrucomicrobiota. This work was focused on three aspects: the definition of Akkermansia genomic species clusters and the calculation and functional characterization of the pangenome for the most represented species; the evolutionary relationship between Akkermansia and their closest relatives from Verrucomicrobiota, defining the gene families which were gained or lost during the emergence of the last Akkermansia common ancestor (LAkkCA) and; the evaluation of the evolutionary pressure metrics for each relevant gene family of main Akkermansia species. This analysis found 25 Akkermansia genomic species clusters distributed in two main clades, divergent from their non-Akkermansia relatives. Pangenome analyses suggest that Akkermansia species have open pangenomes, and the gene gain/loss model indicates that genes associated with mucin degradation (both glycoside hydrolases and peptidases), (micro)aerobic metabolism, surface interaction, and adhesion were part of LAkkCA. Specifically, mucin degradation is a very ancestral innovation involved in the origin of Akkermansia. Horizontal gene transfer detection suggests that Akkermansia could receive genes mostly from unknown sources or from other Gram-negative gut bacteria. Evolutionary metrics suggest that Akkemansia species evolved differently, and even some conserved genes suffered different evolutionary pressures among clades. These results suggest a complex evolutionary landscape of the genus and indicate that mucin degradation could be an essential feature in Akkermansia evolution as a symbiotic species.

Host shift promotes divergent evolution between closely related holoparasitic species

Distinct hosts have been hypothesized to possess the potential for affecting species differentiation and genome evolution of parasitic organisms. However, what host shift history is experienced by the closely related parasites and whether disparate evolution of their genomes occur remain largely unknown. Here, we screened horizontal gene transfer (HGT) events in a pair of sister species of holoparasitic Boschniakia (Orobanchaceae) having obligate hosts from distinct families to recall the former host-parasite associations and performed a comparative analysis to investigate the difference of their organelle genomes. Except those from the current hosts (Ericaceae and Betulaceae), we identified a number of HGTs from Rosaceae supporting the occurrence of unexpected ancient host shifts. Different hosts transfer functional genes which changed nuclear genomes of this sister species. Likewise, different donors transferred sequences to their mitogenomes, which vary in size due to foreign and repetitive elements rather than other factors found in other parasites. The plastomes are both severely reduced, and the degree of difference in reduction syndrome reaches the intergeneric level. Our findings provide new insights into the genome evolution of parasites adapting to different hosts and extend the mechanism of host shift promoting species differentiation to parasitic plant lineages.

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.

Comparative genome identification of accessory genes associated with strong biofilm formation in Vibrio parahaemolyticus

Vibrio parahaemolyticus biofilms on the seafood processing plant surfaces are a potential source of seafood contamination and subsequent food poisoning. Strains differ in their ability to form biofilm, but little is known about the genetic characteristics responsible for biofilm development. In this study, pangenome and comparative genome analysis of V. parahaemolyticus strains reveals genetic attributes and gene repertoire that contribute to robust biofilm formation. The study identified 136 accessory genes that were exclusively present in strong biofilm forming strains and these were functionally assigned to the Gene Ontology (GO) pathways of cellulose biosynthesis, rhamnose metabolic and catabolic processes, UDP-glucose processes and O antigen biosynthesis (p < 0.05). Strategies of CRISPR-Cas defence and MSHA pilus-led attachment were implicated via Kyoto Encyclopedia of Genes and Genomes (KEGG) annotation. Higher levels of horizontal gene transfer (HGT) were inferred to confer more putatively novel properties on biofilm-forming V. parahaemolyticus. Furthermore, cellulose biosynthesis, a neglected potential virulence factor, was identified as being acquired from within the order Vibrionales. The cellulose synthase operons in V. parahaemolyticus were examined for their prevalence (22/138, 15.94 %) and were found to consist of the genes bcsG, bcsE, bcsQ, bcsA, bcsB, bcsZ, bcsC. This study provides insights into robust biofilm formation of V. parahaemolyticus at the genomic level and facilitates: identification of key attributes for robust biofilm formation, elucidation of biofilm formation mechanisms and development of potential targets for novel control strategies of persistent V. parahaemolyticus.

HGTphyloDetect: facilitating the identification and phylogenetic analysis of horizontal gene transfer

BACKGROUND

Horizontal gene transfer (HGT) is an important driver in genome evolution, gain-of-function, and metabolic adaptation to environmental niches. Genome-wide identification of putative HGT events has become increasingly practical, given the rapid growth of genomic data. However, existing HGT analysis toolboxes are not widely used, limited by their inability to perform phylogenetic reconstruction to explore potential donors, and the detection of HGT from both evolutionarily distant and closely related species.

RESULTS

In this study, we have developed HGTphyloDetect, which is a versatile computational toolbox that combines high-throughput analysis with phylogenetic inference, to facilitate comprehensive investigation of HGT events. Two case studies with Saccharomyces cerevisiae and Candida versatilis demonstrate the ability of HGTphyloDetect to identify horizontally acquired genes with high accuracy. In addition, HGTphyloDetect enables phylogenetic analysis to illustrate a likely path of gene transmission among the evolutionarily distant or closely related species.

CONCLUSIONS

The HGTphyloDetect computational toolbox is designed for ease of use and can accurately find HGT events with a very low false discovery rate in a high-throughput manner. The HGTphyloDetect toolbox and its related user tutorial are freely available at https://github.com/SysBioChalmers/HGTphyloDetect.

Pan-Chromosome and Comparative Analysis of Agrobacterium fabrum Reveal Important Traits Concerning the Genetic Diversity, Evolutionary Dynamics, and Niche Adaptation of the Species

Agrobacterium fabrum has been critical for the development of plant genetic engineering and agricultural biotechnology due to its ability to transform eukaryotic cells. However, the gene composition, evolutionary dynamics, and niche adaptation of this species is still unknown. Therefore, we established a comparative genomic analysis based on a pan-chromosome data set to evaluate the genetic diversity of A. fabrum. Here, 25 A. fabrum genomes were selected for analysis by core genome phylogeny combined with the average nucleotide identity (ANI), amino acid identity (AAI), and in silico DNA-DNA hybridization (DDH) values. An open pan-genome of A. fabrum exhibits genetic diversity with variable accessorial genes as evidenced by a consensus pan-genome of 12 representative genomes. The genomic plasticity of A. fabrum is apparent in its putative sequences for mobile genetic elements (MGEs), limited horizontal gene transfer barriers, and potentially horizontally transferred genes. The evolutionary constraints and functional enrichment in the pan-chromosome were measured by the Clusters of Orthologous Groups (COG) categories using eggNOG-mapper software, and the nonsynonymous/synonymous rate ratio (dN/dS) was determined using HYPHY software. Comparative analysis revealed significant differences in the functional enrichment and the degree of purifying selection between the core genome and non-core genome. We demonstrate that the core gene families undergo stronger purifying selection but have a significant bias to contain one or more positively selected sites. Furthermore, although they shared similar genetic diversity, we observed significant differences between chromosome 1 (Chr I) and the chromid in their functional features and evolutionary constraints. We demonstrate that putative genetic elements responsible for plant infection, ecological adaptation, and speciation represent the core genome, highlighting their importance in the adaptation of A. fabrum to plant-related niches. Our pan-chromosome analysis of A. fabrum provides comprehensive insights into the genetic properties, evolutionary patterns, and niche adaptation of the species. IMPORTANCE Agrobacterium spp. live in diverse plant-associated niches such as soil, the rhizosphere, and vegetation, which are challenged by multiple stressors such as diverse energy sources, plant defenses, and microbial competition. They have evolved the ability to utilize diverse resources, escape plant defenses, and defeat competitors. However, the underlying genetic diversity and evolutionary dynamics of Agrobacterium spp. remain unexplored. We examined the phylogeny and pan-genome of A. fabrum to define intraspecies evolutionary relationships. Our results indicate an open pan-genome and numerous MGEs and horizontally transferred genes among A. fabrum genomes, reflecting the flexibility of the chromosomes and the potential for genetic exchange. Furthermore, we observed significant differences in the functional features and evolutionary constraints between the core and accessory genomes and between Chr I and the chromid, respectively.

Bacterial origins of thymidylate metabolism in Asgard archaea and Eukarya

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

Genomic characterization and molecular dating of the novel bacterium Permianibacter aggregans HW001T, which originated from Permian ground water

The Permian Basin is a unique ecosystem located in the southwest of the USA. An unanswered question is whether the bacteria in the Permian Basin adapted to the changing paleomarine environment and survived in the remnants of Permian groundwater. In our previous study, a novel bacterial strain, Permianibacter aggregans HW001T, was isolated from microalgae cultures incubated with Permian Basin waters, and was shown to originate from the Permian Ocean. In this study, strain HW001T was shown to be the representative strain of a novel family, classified as 'Permianibacteraceae'. The results of molecular dating suggested that the strain HW001T diverged ~ 447 million years ago (mya), which is the early Permian period (~ 250 mya). Genome analysis was used to access its potential energy utilization and biosynthesis capacity. A large number of transporters, carbohydrate-active enzymes and protein-degradation related genes have been annotated in the genome of strain HW001T. In addition, a series of important metabolic pathways, such as peptidoglycan biosynthesis, osmotic stress response system and multifunctional quorum sensing were annotated, which may confer the ability to adapt to various unfavorable environmental conditions. Finally, the evolutionary history of strain HW001T was reconstructed and the horizontal transfer of genes was predicted, indicating that the adaptation of P. aggregans to a changing marine environment depends on the evolution of their metabolic capabilities, especially in signal transmission. In conclusion, the results of this study provide genomic information for revealing the adaptive mechanism of strain HW001T to the changing ancient oceans.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-023-00164-3.

One to host them all: genomics of the diverse bacterial endosymbionts of the spider Oedothorax gibbosus

Bacterial endosymbionts of the groups Wolbachia, Cardinium and Rickettsiaceae are well known for their diverse effects on their arthropod hosts, ranging from mutualistic relationships to reproductive phenotypes. Here, we analysed a unique system in which the dwarf spider Oedothorax gibbosus is co-infected with up to five different endosymbionts affiliated with Wolbachia, 'Candidatus Tisiphia' (formerly Torix group Rickettsia), Cardinium and Rhabdochlamydia. Using short-read genome sequencing data, we show that the endosymbionts are heterogeneously distributed among O. gibbosus populations and are frequently found co-infecting spider individuals. To study this intricate host-endosymbiont system on a genome-resolved level, we used long-read sequencing to reconstruct closed genomes of the Wolbachia, 'Ca. Tisiphia' and Cardinium endosymbionts. We provide insights into the ecology and evolution of the endosymbionts and shed light on the interactions with their spider host. We detected high quantities of transposable elements in all endosymbiont genomes and provide evidence that ancestors of the Cardinium, 'Ca. Tisiphia' and Wolbachia endosymbionts have co-infected the same hosts in the past. Our findings contribute to broadening our knowledge about endosymbionts infecting one of the largest animal phyla on Earth and show the usefulness of transposable elements as an evolutionary 'contact-tracing' tool.

Leveraging comparative genomics to uncover alien genes in bacterial genomes

A significant challenge in bacterial genomics is to catalogue genes acquired through the evolutionary process of horizontal gene transfer (HGT). Both comparative genomics and sequence composition-based methods have often been invoked to quantify horizontally acquired genes in bacterial genomes. Comparative genomics methods rely on completely sequenced genomes and therefore the confidence in their predictions increases as the databases become more enriched in completely sequenced genomes. Recent developments including in microbial genome sequencing call for reassessment of alien genes based on information-rich resources currently available. We revisited the comparative genomics approach and developed a new algorithm for alien gene detection. Our algorithm compared favourably with the existing comparative genomics-based methods and is capable of detecting both recent and ancient transfers. It can be used as a standalone tool or in concert with other complementary algorithms for comprehensively cataloguing alien genes in bacterial genomes.

Insights into Genomic Evolution and the Potential Genetic Basis of Klebsiella variicola subsp. variicola ZH07 Reveal Its Potential for Plant Growth Promotion and Autotoxin Degradation

The accumulation of autotoxins in soil causes continuous cropping obstacle stress in crops, and the bioremediation of autotoxins by microorganisms is an efficient process. In this study, strain ZH07 was isolated from the peanut rhizosphere and was found to be utilizing multiple autotoxins as its carbon sources. Based on its genomic characteristics and a phylogenetic analysis, ZH07 represents a member of Klebsiella variicola subsp. variicola. A comparative genomic analysis exhibited evolutionary dynamics exhibited by mobile genetic elements (MGEs), strain-specific genes, potential horizontal genes, and evolutionary constraints driven by purifying selection, which facilitated its genomic adaptation to rhizosphere soil. Genome mining revealed the potential genomic properties associated with plant growth promotion, such as nitrogen fixation, indole acetic acid synthesis, phosphonate solubilization and assimilation, siderophore production, and secondary metabolite synthesis. Moreover, abundant genes putatively responsible for the biodegradation of aromatic xenobiotics, including benzoic acid, cinnamic acid, vanillic acid, protocatechuic acid, phenylacetic acid, and p-hydroxybenzoic acid were also observed in the ZH07 genome. Compared to autotoxin stress alone, the combination of ZH07 application promoted peanut germination and seedling growth. Our analysis revealed the genetic adaptation of ZH07 to the rhizosphere environment and the potential genetic basis and effectiveness of the isolate to serve as a plant growth stimulator. IMPORTANCE Continuous cropping obstacles reduce the production and quality of agricultural products, and the application of rhizosphere beneficial microbes is an important strategy. Strain ZH07 showed autotoxin-degrading and plant growth-promoting capacities. The objectives of this study were to characterize its genomic evolution and the potential genetic basis of the autotoxin degradation and plant growth promotion. ZH07 represents a member of Klebsiella variicola subsp. variicola, based on genomic and phylogenetic analyses. Its genomic components have undergone different degrees of purifying selection, and the disparity in the evolutionary rate may be associated with its niche adaptation. A systematic analysis of the ZH07 genome identified the potential genetic basis that contributes to plant growth promotion and to aromatic xenobiotic biodegradation. This study demonstrates that plant growth-promoting rhizobacteria (PGPR) play important roles in autotoxin biodegradation and can be used as biofertilizers to enhance the growth of peanuts in response to continuous cropping obstacle stress.

Comparative Genomics of the Genus Pseudomonas Reveals Host- and Environment-Specific Evolution

Each Earth ecosystem has unique microbial communities. Pseudomonas bacteria have evolved to occupy a plethora of different ecological niches, including living hosts, such as animals and plants. Many genes necessary for the Pseudomonas-niche interaction and their encoded functions remain unknown. Here, we describe a comparative genomic study of 3,274 genomes with 19,056,667 protein-coding sequences from Pseudomonas strains isolated from diverse environments. We detected functional divergence of Pseudomonas that depends on the niche. Each group of strains from a certain environment harbored a distinctive set of metabolic pathways or functions. The horizontal transfer of genes, which mainly proceeded between closely related taxa, was dependent on the isolation source. Finally, we detected thousands of undescribed proteins and functions associated with each Pseudomonas lifestyle. This research represents an effort to reveal the mechanisms underlying the ecology, pathogenicity, and evolution of Pseudomonas, and it will enable clinical, ecological, and biotechnological advances. IMPORTANCE Microbes play important roles in the health of living beings and in the environment. The knowledge of these functions may be useful for the development of new clinical and biotechnological applications and the restoration and preservation of natural ecosystems. However, most mechanisms implicated in the interaction of microbes with the environment remain poorly understood; thus, this field of research is very important. Here, we try to understand the mechanisms that facilitate the differential adaptation of Pseudomonas-a large and ubiquitous bacterial genus-to the environment. We analyzed more than 3,000 Pseudomonas genomes and searched for genetic patterns that can be related with their coevolution with different hosts (animals, plants, or fungi) and environments. Our results revealed that thousands of genes and genetic features are associated with each niche. Our data may be useful to develop new technical and theoretical advances in the fields of ecology, health, and industry.

Comparative Genomic Insights into the Evolution of Halobacteria-Associated "Candidatus Nanohaloarchaeota"

Members of the phylum "Candidatus Nanohaloarchaeota," a representative lineage within the DPANN superphylum, are characterized by their nanosized cells and symbiotic lifestyle with Halobacteria. However, the development of the symbiosis remains unclear. Here, we propose two novel families, "Candidatus Nanoanaerosalinaceae" and "Candidatus Nanohalalkaliarchaeaceae" in "Ca. Nanohaloarchaeota," represented by five dereplicated metagenome-assembled genomes obtained from hypersaline sediments or related enrichment cultures of soda-saline lakes. Phylogenetic analyses reveal that the two novel families are placed at the root of the family "Candidatus Nanosalinaceae," including the cultivated taxa. The two novel families prefer hypersaline sediments, and the acid shift of predicted proteomes indicates a "salt-in" strategy for hypersaline adaptation. They contain a lower proportion of putative horizontal gene transfers from Halobacteria than "Ca. Nanosalinaceae," suggesting a weaker association with Halobacteria. Functional prediction and historical events reconstruction disclose that they exhibit divergent potentials in carbohydrate and organic acid metabolism and environmental responses. Globally, comparative genomic analyses based on the new families enrich the taxonomic and functional diversity of "Ca. Nanohaloarchaeota" and provide insights into the evolutionary process of "Ca. Nanohaloarchaeota" and their symbiotic relationship with Halobacteria. IMPORTANCE The DPANN superphylum is a group of archaea widely distributed in various habitats. They generally have small cells and have a symbiotic lifestyle with other archaea. The archaeal symbiotic interaction is vital to understanding microbial communities. However, the formation and evolution of the symbiosis between the DPANN lineages and other diverse archaea remain unclear. Based on phylogeny, habitat distribution, hypersaline adaptation, host prediction, functional potentials, and historical events of "Ca. Nanohaloarchaeota," a representative phylum within the DPANN superphylum, we report two novel families representing intermediate stages, and we infer the evolutionary process of "Ca. Nanohaloarchaeota" and their Halobacteria-associated symbiosis. Altogether, this research helps in understanding the evolution of symbiosis in "Ca. Nanohaloarchaeota" and provides a model for the evolution of other DPANN lineages.

Expanding the genomic encyclopedia of Actinobacteria with 824 isolate reference genomes

The phylum Actinobacteria includes important human pathogens like Mycobacterium tuberculosis and Corynebacterium diphtheriae and renowned producers of secondary metabolites of commercial interest, yet only a small part of its diversity is represented by sequenced genomes. Here, we present 824 actinobacterial isolate genomes in the context of a phylum-wide analysis of 6,700 genomes including public isolates and metagenome-assembled genomes (MAGs). We estimate that only 30%-50% of projected actinobacterial phylogenetic diversity possesses genomic representation via isolates and MAGs. A comparison of gene functions reveals novel determinants of host-microbe interaction as well as environment-specific adaptations such as potential antimicrobial peptides. We identify plasmids and prophages across isolates and uncover extensive prophage diversity structured mainly by host taxonomy. Analysis of >80,000 biosynthetic gene clusters reveals that horizontal gene transfer and gene loss shape secondary metabolite repertoire across taxa. Our observations illustrate the essential role of and need for high-quality isolate genome sequences.

Reconstructing horizontal gene flow network to understand prokaryotic evolution

Horizontal gene transfer (HGT) is a major source of phenotypic innovation and a mechanism of niche adaptation in prokaryotes. Quantification of HGT is critical to decipher its myriad roles in microbial evolution and adaptation. Advances in genome sequencing and bioinformatics have augmented our ability to understand the microbial world, particularly the direct or indirect influence of HGT on diverse life forms. Methods for detecting HGT can be classified into phylogenetic-based and parametric or composition-based approaches. Here, we exploited the complementary strengths of both the approaches to construct a high confidence horizontal gene flow network. Our network is unique in its ability to detect the transfer of native genes of a genome to genomes from other taxa, thus establishing donor and recipient organisms (taxa), rather than through a post hoc analysis as is the practice with several other approaches. The scale-free horizontal gene flow network presented here provides new insights into modes of transfer for the exchange of genetic information and also illuminates differential gene flow across phyla.

Reconstructing horizontal gene flow network to understand prokaryotic evolution

Horizontal gene transfer (HGT) is a major source of phenotypic innovation and a mechanism of niche adaptation in prokaryotes. Quantification of HGT is critical to decipher its myriad roles in microbial evolution and adaptation. Advances in genome sequencing and bioinformatics have augmented our ability to understand the microbial world, particularly the direct or indirect influence of HGT on diverse life forms. Methods for detecting HGT can be classified into phylogenetic-based and parametric or composition-based approaches. Here, we exploited the complementary strengths of both the approaches to construct a high confidence horizontal gene flow network. Our network is unique in its ability to detect the transfer of native genes of a genome to genomes from other taxa, thus establishing donor and recipient organisms (taxa), rather than through a post hoc analysis as is the practice with several other approaches. The scale-free horizontal gene flow network presented here provides new insights into modes of transfer for the exchange of genetic information and also illuminates differential gene flow across phyla.

Next Generation Sequencing and Comparative Genomic Analysis Reveal Extreme Plasticity of Two Burkholderia glumae Strains HN1 and HN2

Burkholderia glumae is an important rice pathogen, thus the genomic and evolutionary history may be helpful to control this notorious pathogen. Here, we present two complete genomes of the B. glumae strains HN1 and HN2, which were isolated from diseased rice seed in China. Average nucleotide identity (ANI) analysis shows greater than 99% similarity of the strains HN1 and HN2 with other published B. glumae genomes. Genomic annotation revealed that the genome of strain HN1 consists of five replicons (6,680,415 bp) with an overall G + C content of 68.06%, whereas the genome of strain HN2 comprises of three replicons (6,560,085 bp) with an overall G + C content of 68.34%. The genome of HN1 contains 5434 protein-coding genes, 351 pseudogenes, and 1 CRISPR, whereas the genome of HN2 encodes 5278 protein-coding genes, 357 pseudogenes, and 2 CRISPR. Both strains encode many pathogenic-associated genes (143 genes in HN1 vs. 141 genes in HN2). Moreover, comparative genomic analysis shows the extreme plasticity of B. glumae, which may contribute to its pathogenicity. In total, 259 single-copy genes were affected by positive selection. These genes may contribute to the adaption to different environments. Notably, six genes were characterized as virulence factors which may be an additional way to assist the pathogenicity of B. glumae.

Insights into Evolutionary, Genomic, and Biogeographic Characterizations of Chryseobacterium nepalense Represented by a Polyvinyl Alcohol-Degrading Bacterium, AC3

Chryseobacterium spp. are Gram-negative rods found ubiquitously in the environment, with certain species being reported as having unusual degrading properties. Polyvinyl alcohol (PVA) is used widely in industry but causes serious global environmental pollution. Here, we report the complete genome sequence of a novel bacterium, AC3, that efficiently degrades PVA. As the representative genome of Chryseobacterium nepalense, key genomic characteristics (e.g., mobile genetic elements, horizontal genes, genome-scale metabolic network, secondary metabolite biosynthesis gene clusters, and carbohydrate-active enzymes) were comprehensively investigated to reveal the potential genetic features of this species. Core genome phylogenetic analysis in combination with average nucleotide identity, average amino acid identity, and in silico DNA-DNA hybridization values provided an accurate taxonomic position of C. nepalense in the genus Chryseobacterium. Comparative genomic analysis of AC3 with closely related species suggested evolutionary dynamics characterized by a species-specific genetic repertoire, dramatic rearrangements, and evolutionary constraints driven by selective pressure, which facilitated the speciation and adaptative evolution of C. nepalense. Biogeographic characterization indicated that this species is ubiquitously distributed not only in soil habitats but also in a variety of other source niches. Bioinformatic analysis revealed the potential genetic basis of PVA degradation in AC3, which included six putative genes associated with the synthesis of PVA dehydrogenase, cytochrome c, oxidized PVA hydrolase, and secondary alcohol dehydrogenase. Our study reports the first complete genome of C. nepalense with PVA-degrading properties, providing comprehensive insights into the genomic characteristics of this species and increasing our understanding of the microbial degradation of PVA. IMPORTANCE Although PVA is a biodegradable polymer, the widespread use of PVA in global industrialization has resulted in serious environmental problems. To date, knowledge of effective and applicable PVA-degrading bacteria is limited, and thus, the discovery of novel PVA biodegraders is pertinent. Here, we isolated a novel bacterial strain, AC3, which efficiently degraded PVA. The complete genome of AC3 was sequenced as the first genome sequence of the species C. nepalense. Comparative genomic analysis was performed to comprehensively investigate the phylogenetic relationships, genome-scale metabolic network, key genomic characteristics associated with genomic evolution, evolutionary dynamics between AC3 and its close relatives, and biogeographic characterization of C. nepalense, particularly regarding the potential genetic basis of PVA degradation. These findings could advance our understanding of the genomic characteristics of C. nepalense and PVA bioremediation.

Pan-genome evolution and its association with divergence of metabolic functions in Bifidobacterium genus

Previous studies were mainly focused on genomic evolution and diversity of type species of Bifidobacterium genus due to their health-promoting effect on host. However, those studies were mainly based on species-level taxonomic resolution, adaptation, and characterization of carbohydrate metabolic features of the bifidobacterial species. Here, a comprehensive analysis of the type strain genome unveils the association of pan-genome evolution with the divergence of metabolic function of the Bifidobacterium genus. This study has also demonstrated that horizontal gene transfer, as well as genome expansion and reduction events, leads to the divergence of metabolic functions in Bifidobacterium genus. Furthermore, the genome-based search of probiotic traits among all the available bifidobacterial type strains gives hints on type species, that could confer health benefits to nutrient-deficient individuals. Altogether, the present study provides insight into the developments of genomic evolution, functional divergence, and potential probiotic type species of the Bifidobacterium genus.

An Ancient Respiratory System in the Widespread Sedimentary Archaea Thermoprofundales

Thermoprofundales, formerly Marine Benthic Group D (MBG-D), is a ubiquitous archaeal lineage found in sedimentary environments worldwide. However, its taxonomic classification, metabolic pathways, and evolutionary history are largely unexplored because of its uncultivability and limited number of sequenced genomes. In this study, phylogenomic analysis and average amino acid identity values of a collection of 146 Thermoprofundales genomes revealed five Thermoprofundales subgroups (A-E) with distinct habitat preferences. Most of the microorganisms from Subgroups B and D were thermophiles inhabiting hydrothermal vents and hot spring sediments, whereas those from Subgroup E were adapted to surface environments where sunlight is available. H2 production may be featured in Thermoprofundales as evidenced by a gene cluster encoding the ancient membrane-bound hydrogenase (MBH) complex. Interestingly, a unique structure separating the MBH gene cluster into two modular units was observed exclusively in the genomes of Subgroup E, which included a peripheral arm encoding the [NiFe] hydrogenase domain and a membrane arm encoding the Na+/H+ antiporter domain. These two modular structures were confirmed to function independently by detecting the H2-evolving activity in vitro and salt tolerance to 0.2 M NaCl in vivo, respectively. The peripheral arm of Subgroup E resembles the proposed common ancestral respiratory complex of modern respiratory systems, which plays a key role in the early evolution of life. In addition, molecular dating analysis revealed that Thermoprofundales is an early emerging archaeal lineage among the extant MBH-containing microorganisms, indicating new insights into the evolution of this ubiquitous archaea lineage.

Genome-guided prediction of acid resistance mechanisms in acidophilic methanotrophs of phylogenetically deep-rooted Verrucomicrobia isolated from geothermal environments

Verrucomicrobia are a group of microorganisms that have been proposed to be deeply rooted in the Tree of Life. Some are methanotrophs that oxidize the potent greenhouse gas methane and are thus important in decreasing atmospheric concentrations of the gas, potentially ameliorating climate change. They are widespread in various environments including soil and fresh or marine waters. Recently, a clade of extremely acidophilic Verrucomicrobia, flourishing at pH < 3, were described from high-temperature geothermal ecosystems. This novel group could be of interest for studies about the emergence of life on Earth and to astrobiologists as homologs for possible extraterrestrial life. In this paper, we describe predicted mechanisms for survival of this clade at low pH and suggest its possible evolutionary trajectory from an inferred neutrophilic ancestor. Extreme acidophiles are defined as organisms that thrive in extremely low pH environments (≤ pH 3). Many are polyextremophiles facing high temperatures and high salt as well as low pH. They are important to study for both providing fundamental insights into biological mechanisms of survival and evolution in such extreme environments and for understanding their roles in biotechnological applications such as industrial mineral recovery (bioleaching) and mitigation of acid mine drainage. They are also, potentially, a rich source of novel genes and pathways for the genetic engineering of microbial strains. Acidophiles of the Verrucomicrobia phylum are unique as they are the only known aerobic methanotrophs that can grow optimally under acidic (pH 2–3) and moderately thermophilic conditions (50–60°C). Three moderately thermophilic genera, namely Methylacidiphilum, Methylacidimicrobium, and Ca. Methylacidithermus, have been described in geothermal environments. Most of the investigations of these organisms have focused on their methane oxidizing capabilities (methanotrophy) and use of lanthanides as a protein cofactor, with no extensive study that sheds light on the mechanisms that they use to flourish at extremely low pH. In this paper, we extend the phylogenetic description of this group of acidophiles using whole genome information and we identify several mechanisms, potentially involved in acid resistance, including “first line of defense” mechanisms that impede the entry of protons into the cell. These include the presence of membrane-associated hopanoids, multiple copies of the outer membrane protein (Slp), and inner membrane potassium channels (kup, kdp) that generate a reversed membrane potential repelling the intrusion of protons. Acidophilic Verrucomicrobia also display a wide array of proteins potentially involved in the “second line of defense” where protons that evaded the first line of defense and entered the cell are expelled or neutralized, such as the glutamate decarboxylation (gadAB) and phosphate-uptake systems. An exclusive N-type ATPase F₀-F₁ was identified only in acidophiles of Verrucomicrobia and is predicted to be a specific adaptation in these organisms. Phylogenetic analyses suggest that many predicted mechanisms are evolutionarily conserved and most likely entered the acidophilic lineage of Verrucomicrobia by vertical descent from a common ancestor. However, it is likely that some defense mechanisms such as gadA and kup entered the acidophilic Verrucomicrobia lineage by horizontal gene transfer.

Genomic insights into versatile lifestyle of three new bacterial candidate phyla

Metagenomic explorations of the Earth's biosphere enable the discovery of previously unknown bacterial lineages of phylogenetic and ecological significance. Here, we retrieved 11 metagenomic-assembled genomes (MAGs) affiliated to three new monophyletic bacterial lineages from the seawater of the Yap Trench. Phylogenomic analysis revealed that each lineage is a new bacterial candidate phylum, subsequently named Candidatus Qinglongiota, Candidatus Heilongiota, and Candidatus Canglongiota. Metabolic reconstruction of genomes from the three phyla suggested that they adopt a versatile lifestyle, with the potential to utilize various types of sugars, proteins, and/or short-chain fatty acids through anaerobic pathways. This was further confirmed by a global distribution map of the three phyla, indicating a preference for oxygen-limited or particle-attached niches, such as anoxic sedimentary environments. Of note, Candidatus Canglongiota genomes harbor genes for the complete Wood- Ljungdahl pathway and sulfate reduction that are similar to those identified in some sulfate-reducing bacteria. Evolutionary analysis indicated that gene gain and loss events, and horizontal gene transfer (HGT) play important roles in shaping the genomic and metabolic features of the three new phyla. This study presents the genomic insight into the ecology, metabolism, and evolution of three new phyla, which broadens the phylum-level diversity within the domain Bacteria.

Phylogenies of the 16S rRNA gene and its hypervariable regions lack concordance with core genome phylogenies

BACKGROUND

The 16S rRNA gene is used extensively in bacterial phylogenetics, in species delineation, and now widely in microbiome studies. However, the gene suffers from intragenomic heterogeneity, and reports of recombination and an unreliable phylogenetic signal are accumulating. Here, we compare core gene phylogenies to phylogenies constructed using core gene concatenations to estimate the strength of signal for the 16S rRNA gene, its hypervariable regions, and all core genes at the intra- and inter-genus levels. Specifically, we perform four intra-genus analyses (Clostridium, n = 65; Legionella, n = 47; Staphylococcus, n = 36; and Campylobacter, n = 17) and one inter-genus analysis [41 core genera of the human gut microbiome (31 families, 17 orders, and 12 classes), n = 82].

RESULTS

At both taxonomic levels, the 16S rRNA gene was recombinant and subject to horizontal gene transfer. At the intra-genus level, the gene showed one of the lowest levels of concordance with the core genome phylogeny (50.7% average). Concordance for hypervariable regions was lower still, with entropy masking providing little to no benefit. A major factor influencing concordance was SNP count, which showed a positive logarithmic association. Using this relationship, we determined that 690 ± 110 SNPs were required for 80% concordance (average 16S rRNA gene SNP count was 254). We also found a wide range in 16S-23S-5S rRNA operon copy number among genomes (1-27). At the inter-genus level, concordance for the whole 16S rRNA gene was markedly higher (73.8% - 10th out of 49 loci); however, the most concordant hypervariable regions (V4, V3-V4, and V1-V2) ranked in the third quartile (62.5 to 60.0%).

CONCLUSIONS

Ramifications of a poor phylogenetic performance for the 16S rRNA gene are far reaching. For example, in addition to incorrect species/strain delineation and phylogenetic inference, it has the potential to confound community diversity metrics if phylogenetic information is incorporated - for example, with popular approaches such as Faith's phylogenetic diversity and UniFrac. Our results highlight the problematic nature of these approaches and their use (along with entropy masking) is discouraged. Lastly, the wide range in 16S rRNA gene copy number among genomes also has a strong potential to confound diversity metrics. Video Abstract.

Genomic Insights Into the Interspecific Diversity and Evolution of Mobiluncus, a Pathogen Associated With Bacterial Vaginosis

Bacterial vaginosis (BV) is a common vaginal infection and has been associated with increased risk for a wide array of health issues. BV is linked with a variety of heterogeneous pathogenic anaerobic bacteria, among which Mobiluncus is strongly associated with BV diagnosis. However, their genetic features, pathogenicity, interspecific diversity, and evolutionary characters have not been illustrated at genomic level. The current study performed phylogenomic and comparative genomic analyses of Mobiluncus. Phylogenomic analyses revealed remarkable phylogenetic distinctions among different species. Compared with M. curtisii, M. mulieris had a larger genome and pangenome size with more insertion sequences but less CRISPR-Cas systems. In addition, these two species were diverse in profile of virulence factors, but harbored similar antibiotic resistance genes. Statistically different functional genome profiles between strains from the two species were determined, as well as correlations of some functional genes/pathways with putative pathogenicity. We also showed that high levels of horizontal gene transfer might be an important strategy for species diversification and pathogenicity. Collectively, this study provides the first genome sequence level description of Mobiluncus, and may shed light on its virulence/pathogenicity, functional diversification, and evolutionary dynamics. Our study could facilitate the further investigations of this important pathogen, and might improve the future treatment of BV.

Comparative Genomics and Pan-Genome Driven Prediction of a Reduced Genome of Akkermansia muciniphila

Akkermanisia muciniphila imparts important health benefits and is considered a next-generation probiotic. It is imperative to understand the genomic diversity and metabolic potential of the species for safer applications as probiotics. As it resides with both health-promoting and pathogenic bacteria, understanding the evolutionary patterns are crucial, but this area remains largely unexplored. Moreover, pan-genome has previously been established based on only a limited number of strains and without careful strain selection. The pan-genomics have become very important for understanding species diversity and evolution. In the current study, a systematic approach was used to find a refined pan-genome profile of A. muciniphila by excluding too-diverse strains based on average nucleotide identity-based species demarcation. The strains were divided into four phylogroups using a variety of clustering techniques. Horizontal gene transfer and recombination patterns were also elucidated. Evolutionary patterns revealed that different phylogroups were expanding differently. Furthermore, a comparative evaluation of the metabolic potential of the pan-genome and its subsections was performed. Lastly, the study combines functional annotation, persistent genome, and essential genes to devise an approach to determine a minimal genome that can systematically remove unwanted genes, including virulent factors. The selection of one strain to be used as a chassis for the prediction of a reduced genome was very carefully performed by analyzing several genomic parameters, including the number of unique genes and the resistance and pathogenic potential of the strains. The strategy could be applied to other microbes, including human-associated microbiota, towards a common goal of predicting a minimal or a reduced genome.

Evolutionary Divergence of Marinobacter Strains in Cryopeg Brines as Revealed by Pangenomics

Marinobacter spp. are cosmopolitan in saline environments, displaying a diverse set of metabolisms that allow them to competitively occupy these environments, some of which can be extreme in both salinity and temperature. Here, we introduce a distinct cluster of Marinobacter genomes, composed of novel isolates and in silico assembled genomes obtained from subzero, hypersaline cryopeg brines, relic seawater-derived liquid habitats within permafrost sampled near Utqiaġvik, Alaska. Using these new genomes and 45 representative publicly available genomes of Marinobacter spp. from other settings, we assembled a pangenome to examine how the new extremophile members fit evolutionarily and ecologically, based on genetic potential and environmental source. This first genus-wide genomic analysis revealed that Marinobacter spp. in general encode metabolic pathways that are thermodynamically favored at low temperature, cover a broad range of organic compounds, and optimize protein usage, e.g., the Entner–Doudoroff pathway, the glyoxylate shunt, and amino acid metabolism. The new isolates contributed to a distinct clade of subzero brine-dwelling Marinobacter spp. that diverged genotypically and phylogenetically from all other Marinobacter members. The subzero brine clade displays genomic characteristics that may explain competitive adaptations to the extreme environments they inhabit, including more abundant membrane transport systems (e.g., for organic substrates, compatible solutes, and ions) and stress-induced transcriptional regulatory mechanisms (e.g., for cold and salt stress) than in the other Marinobacter clades. We also identified more abundant signatures of potential horizontal transfer of genes involved in transcription, the mobilome, and a variety of metabolite exchange systems, which led to considering the importance of this evolutionary mechanism in an extreme environment where adaptation via vertical evolution is physiologically rate limited. Assessing these new extremophile genomes in a pangenomic context has provided a unique view into the ecological and evolutionary history of the genus Marinobacter, particularly with regard to its remarkable diversity and its opportunism in extremely cold and saline environments.

Current Methods for Recombination Detection in Bacteria

The role of genetic exchanges, i.e., homologous recombination (HR) and horizontal gene transfer (HGT), in bacteria cannot be overestimated for it is a pivotal mechanism leading to their evolution and adaptation, thus, tracking the signs of recombination and HGT events is importance both for fundamental and applied science. To date, dozens of bioinformatics tools for revealing recombination signals are available, however, their pros and cons as well as the spectra of solvable tasks have not yet been systematically reviewed. Moreover, there are two major groups of software. One aims to infer evidence of HR, while the other only deals with horizontal gene transfer (HGT). However, despite seemingly different goals, all the methods use similar algorithmic approaches, and the processes are interconnected in terms of genomic evolution influencing each other. In this review, we propose a classification of novel instruments for both HR and HGT detection based on the genomic consequences of recombination. In this context, we summarize available methodologies paying particular attention to the type of traceable events for which a certain program has been designed.

Metagenomic analysis further extends the role of Chloroflexi in fundamental biogeochemical cycles

Chloroflexi members are ubiquitous and have been extensively studied; however, the evolution and metabolic pathways of Chloroflexi members have long been debated. In the present study, the evolution and the metabolic potentials of 17 newly obtained Chloroflexi metagenome-assembled genomes (MAGs) were evaluated using genome and horizontal gene transfer (HGT) analysis. Taxonomic analysis suggests that the MAGs of the present study might be novel. One MAG encodes genes for anoxygenic phototrophy. The HGT analysis suggest that genes responsible for anoxygenic phototrophy in the MAG might have been transferred from Proteobacteria/Chlorobi. The evolution of anaerobic photosynthesis, which has long been questioned, has now been shown to be the result of HGT events. An incomplete Wood-Ljungdahl pathway (with missing genes metF, acsE, fdh, and acsA) was reported in Dehalococcoidetes members. In the present study, MAGs that were not the Dehalococcoidetes members encode genes acsA, acsB, metF and acsE. The genes responsible for sulfate reduction (sat, cysC and sir), dissimilatory sulfite reductase (dsrA and dsrB), and aerobic and anaerobic carbon monoxide oxidation (coxSML and cooSF) were detected in the present study MAGs. The present study expands our knowledge of the possible metabolic potentials of the phylum Chloroflexi and clarifies the evolution of anaerobic photosynthesis.

Comprehensive Analysis Reveals the Genetic and Pathogenic Diversity of Ralstonia solanacearum Species Complex and Benefits Its Taxonomic Classification

Ralstonia solanacearum species complex (RSSC) is a diverse group of plant pathogens that attack a wide range of hosts and cause devastating losses worldwide. In this study, we conducted a comprehensive analysis of 131 RSSC strains to detect their genetic diversity, pathogenicity, and evolution dynamics. Average nucleotide identity analysis was performed to explore the genomic relatedness among these strains, and finally obtained an open pangenome with 32,961 gene families. To better understand the diverse evolution and pathogenicity, we also conducted a series of analyses of virulence factors (VFs) and horizontal gene transfer (HGT) in the pangenome and at the single genome level. The distribution of VFs and mobile genetic elements (MGEs) showed significant differences among different groups and strains, which were consistent with the new nomenclatures of the RSSC with three distinct species. Further functional analysis showed that most HGT events conferred from Burkholderiales and played a great role in shaping the genomic plasticity and genetic diversity of RSSC genomes. Our work provides insights into the genetic polymorphism, evolution dynamics, and pathogenetic variety of RSSC and provides strong supports for the new taxonomic classification, as well as abundant resources for studying host specificity and pathogen emergence.

Comparative Genomic Analysis of Vibrio cincinnatiensis Provides Insights into Genetic Diversity, Evolutionary Dynamics, and Pathogenic Traits of the Species

Vibrio cincinnatiensis is a poorly understood pathogenic Vibrio species, and the underlying mechanisms of its genetic diversity, genomic plasticity, evolutionary dynamics, and pathogenicity have not yet been comprehensively investigated. Here, a comparative genomic analysis of V. cincinnatiensis was constructed. The open pan-genome with a flexible gene repertoire exhibited genetic diversity. The genomic plasticity and stability were characterized by the determinations of diverse mobile genetic elements (MGEs) and barriers to horizontal gene transfer (HGT), respectively. Evolutionary divergences were exhibited by the difference in functional enrichment and selective pressure between the different components of the pan-genome. The evolution on the Chr I and Chr II core genomes was mainly driven by purifying selection. Predicted essential genes in V. cincinnatiensis were mainly found in the core gene families on Chr I and were subject to stronger evolutionary constraints. We identified diverse virulence-related elements, including the gene clusters involved in encoding flagella, secretion systems, several pili, and scattered virulence genes. Our results indicated the pathogenic potential of V. cincinnatiensis and highlighted that HGT events from other Vibrio species promoted pathogenicity. This pan-genome study provides comprehensive insights into this poorly understood species from the genomic perspective.

Comprehensive Analysis of the Nocardia cyriacigeorgica Complex Reveals Five Species-Level Clades with Different Evolutionary and Pathogenicity Characteristics

Nocardia cyriacigeorgica is a common etiological agent of nocardiosis that has increasingly been implicated in serious pulmonary infections, especially in immunocompromised individuals. However, the evolution, diversity, and pathogenesis of N. cyriacigeorgica have remained unclear. Here, we performed a comparative genomic analysis using 91 N. cyriacigeorgica strains, 45 of which were newly sequenced in this study. Phylogenetic and average nucleotide identity (ANI) analyses revealed that N. cyriacigeorgica contained five species-level clades (8.6 to 14.6% interclade genetic divergence), namely, the N. cyriacigeorgica complex (NCC). Further pan-genome analysis revealed extensive differences among the five clades in nine functional categories, such as energy production, lipid metabolism, secondary metabolites, and signal transduction mechanisms. All 2,935 single-copy core genes undergoing purifying selection were highly conserved across NCC. However, clades D and E exhibited reduced selective constraints, compared to clades A to C. Horizontal gene transfer (HGT) and mobile genetic elements contributed to genomic plasticity, and clades A and B had experienced a higher level of HGT events than other clades. A total of 129 virulence factors were ubiquitous across NCC, such as the mce operon, hemolysin, and type VII secretion system (T7SS). However, different distributions of three toxin-coding genes and two new types of mce operons were detected, which might contribute to pathogenicity differences among the members of the NCC. Overall, our study provides comprehensive insights into the evolution, genetic diversity, and pathogenicity of NCC, facilitating the prevention of infections.

IMPORTANCENocardia species are opportunistic bacterial pathogens that can affect all organ systems, primarily the skin, lungs, and brain. N. cyriacigeorgica is the most prevalent species within the genus, exhibits clinical significance, and can cause severe infections when disseminated throughout the body. However, the evolution, diversity, and pathogenicity of N. cyriacigeorgica remain unclear. Here, we have conducted a comparative genomic analysis of 91 N. cyriacigeorgica strains and revealed that N. cyriacigeorgica is not a single species but is composed of five closely related species. In addition, we discovered that these five species differ in many ways, involving selection pressure, horizontal gene transfer, functional capacity, pathogenicity, and antibiotic resistance. Overall, our work provides important clues in dissecting the evolution, genetic diversity, and pathogenicity of NCC, thereby advancing prevention measures against these infections.

Prediction and Inferred Evolution of Acid Tolerance Genes in the Biotechnologically Important Acidihalobacter Genus

Acidihalobacter is a genus of acidophilic, gram-negative bacteria known for its ability to oxidize pyrite minerals in the presence of elevated chloride ions, a capability rare in other iron-sulfur oxidizing acidophiles. Previous research involving Acidihalobacter spp. has focused on their applicability in saline biomining operations and their genetic arsenal that allows them to cope with chloride, metal and oxidative stress. However, an understanding of the molecular adaptations that enable Acidihalobacter spp. to thrive under both acid and chloride stress is needed to provide a more comprehensive understanding of how this genus can thrive in such extreme biomining conditions. Currently, four genomes of the Acidihalobacter genus have been sequenced: Acidihalobacter prosperus DSM 5130^T, Acidihalobacter yilgarnensis DSM 105917^T, Acidihalobacter aeolianus DSM 14174^T, and Acidihalobacter ferrooxydans DSM 14175^T. Phylogenetic analysis shows that the Acidihalobacter genus roots to the Chromatiales class consisting of mostly halophilic microorganisms. In this study, we aim to advance our knowledge of the genetic repertoire of the Acidihalobacter genus that has enabled it to cope with acidic stress. We provide evidence of gene gain events that are hypothesized to help the Acidihalobacter genus cope with acid stress. Potential acid tolerance mechanisms that were found in the Acidihalobacter genomes include multiple potassium transporters, chloride/proton antiporters, glutamate decarboxylase system, arginine decarboxylase system, urease system, slp genes, squalene synthesis, and hopanoid synthesis. Some of these genes are hypothesized to have entered the Acidihalobacter via vertical decent from an inferred non-acidophilic ancestor, however, horizontal gene transfer (HGT) from other acidophilic lineages is probably responsible for the introduction of many acid resistance genes.

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.

Metabacillus dongyingensis sp. nov. Is Represented by the Plant Growth-Promoting Bacterium BY2G20 Isolated from Saline-Alkaline Soil and Enhances the Growth of Zea mays L. under Salt Stress

A novel plant growth-promoting rhizobacterium (PGPR), which was designated strain BY2G20, was isolated from saline-alkaline soil in Dongying, China. Strain BY2G20 can grow at a NaCl range from 0 to 7% and a pH range from 7 to 9 and can prevent the growth of the phytopathogen Ralstonia solanacearum. Based on its phenotypic and genomic characteristics and phylogenetic analysis, strain BY2G20 represents a novel species of the genus Metabacillus, for which the name Metabacillus dongyingensis sp. nov. is proposed. Comparative genomic analysis of strain BY2G20 with its closely related species exhibited a high level of evolutionary plasticity derived by horizontal gene transfer, which facilitated adaptative evolution. Different evolutionary constraints have operated on the diverse functions of BY2G20, with the gene adapted to saline-alkaline ecosystems experiencing functional constraints. We determined the genetic properties of saline-alkaline tolerance and plant growth promotion, such as cation-proton antiporters, cation transporters, osmoprotectant synthesis and transport, H⁺-transporting F₁F₀-ATPase, indole-3-acetic acid production, and secondary metabolite synthesis. We also evaluated the effects of strain BY2G20 on the growth of Zea mays L. (maize) under salt stress. The physiological parameters of maize such as plant height, stem diameter, dry biomass, and fresh biomass were significantly higher after inoculating strain BY2G20 under salt stress, indicating that inoculation with BY2G20 enhanced the growth of maize in saline areas. This study demonstrates that M. dongyingensis sp. nov. BY2G20 is a potential candidate for organic agriculture biofertilizers in saline-alkaline areas.

IMPORTANCE Plant growth and yield are adversely affected by soil salinity. PGPRs can promote plant growth and enhance plant tolerance to salt stress. In this study, a saline-alkaline tolerant PGPR strain BY2G20 was isolated from the rhizosphere of Ulmus pumila in Dongying, China. Strain BY2G20 represents a novel species within the genus Metabacillus based on phenotypic, genomic, and phylogenetic analysis. Genomic components have undergone different functional constraints, and the disparity in the evolutionary rate may be associated with the adaptation to a specific niche. Genomic analysis revealed numerous adaptive features of strain BY2G20 to a saline-alkaline environment and rhizosphere, especially genes related to salt tolerance, pH adaptability, and plant growth promotion. Our work also exhibited that inoculation of strain BY2G20 enhanced the growth of maize under salt stress. This study demonstrates that PGPRs play an important role in stimulating salt tolerance in plants and can be used as biofertilizers to enhance the growth of crops in saline-alkaline areas.