Phylogeny of bacterial and archaeal genomes using conserved genes: supertrees and supermatrices

Pathogenomic profile and clonal diversity of potential zoonotic MRSA-CC7-ST789-t091-SCCmecV from human skin and soft tissue infections

The whole genome sequence (WGS) of prevalent MRSA strains harboring mecA gene obtained from skin and soft tissue infections (SSTIs) in Nigerian hospitals were profiled for pathogenomic structure and evaluated for clonal diversity. The two MRSA strains identified among 66 isolated multi-drug resistant S. aureus from a collection of 256 clinical samples were phenotyped for antibiotic resistance and genotyped for mecA, SCCmec, and spa types. The mecA positive MRSA was analysed using whole-genome sequencing for resistomes, virulomes, phylogenomic profiles and clonal diversity. The identified MRSA-CC7-ST789-t091-SCCmecV strains from a female child (aged 1 year) with severe otorrhea and an adult male (aged 23) with purulent wound abscess showed high-level resistance to streptomycin, vancomycin, kanamycin, sulfamethoxazole and ciprofloxacin. Both strains harbored abundant resistomes, inherent plasmids, chromosomal replicons and typical seven housekeeping genes (arc3, aroE4, glpF1, gmk4, pta4, tpi6, yqiL3). The most abundant putative virulomes were pathogenesis-associated proteins (included hemolysin gamma, leucocidins, proteases, staphylococcal superantigen/enterotoxin-like genes (Set/Ssl), capsule- and biofilm-associated genes, and hyaluronate lyase). Comparative phylogenomic analysis revealed the relatedness of the two clonal strains with prevalent MRSA-CC7 pathotypes observed in Italy (2013 and 2014), Denmark (2014), Thailand (2015 and 2016), USA (2018), and Nigeria (2016 and 2020); and share high genetic similarities with livestock strains from cow milk and cattle. Identified MRSA-CC7-ST789-t091-SCCmecV pathotypes implicated in SSTIs from Nigeria harboring repertoires of antibiotic resistance and virulence genes, and genetic relatedness with livestock strains; show the possibility of gene transfer between animal and human. Adequate hospital MRSA infection control and geno-epidemiological surveillance for animal and human transfer is required.

Phylogenomic evaluation of Mangrovimicrobium sediminis gen. nov. sp. nov., the first nitrogen fixing member of the family Halieaceae adapted to mangrove habitat and reclassification of Halioglobus pacificus to Pseudohaliglobus pacificus comb. nov

The taxonomic position and genomic characteristics of a nitrogen fixing and polymer degrading marine bacterium, strain SAOS 164 isolated from a mangrove sediment sample was investigated. Sequence analysis based on 16S rRNA gene identified it as a member of family Halieaceae with closest similarity to Haliea salexigens DSM 19537T (96.3 %), H. alexandrii LZ-16-2T (96.2 %) and Parahaliea maris HSLHS9T (96.0 %) but was distantly related to the genera Haliea, Parahaliea and Halioglobus in phylogenetic trees. In order to ascertain the exact taxonomic position, phylogeny based on RpoBC proteins, whole genome, core and orthologous genes, and comparative analysis of metabolic potential retrieved the strain in an independent lineage clustering along with the genera Halioglobus, Pseudohalioglobus and Seongchinamella. Further, various genome based delimitation parameters represented by mol % GC content, percentage of conserved proteins (POCP), and amino acid identity (AAI) along with chemotaxonomic markers (i.e. fatty acids and polar lipids) supported the inferences of genome based phylogeny and indicated that the strain SAOS 164 belongs to a novel genus. The genome was mapped to 4.8 Mb in size with 65.1 % DNA mol% G + C content. In-silico genomic investigation and phenotyping revealed diverse metabolite genes/pathways related to polymer hydrolysis, nitrogen fixation, light induced growth, carbohydrate, sulfur, phosphorus and amino acid metabolism, virulence factors, defense mechanism, and stress-responsive elements facilitating survival in the mangrove habitat. Based on polyphasic taxonomic approach including genome analyses, a novel genus Mangrovimicrobium sediminis gen. nov. sp. nov. (=SAOS 164T = MTCC 12907T = KCTC 52755T = JCM 32136T) is proposed. Additionally, the reclassification of Halioglobus pacificus (=DSM 27932T = KCTC 23430T = S1-72T) to Pseudhalioglobus pacificus comb. nov. is also proposed.

Advancements in prokaryotic systematics and the role of Bergey's International Society for Microbial Systematicsin addressing challenges in the meta-data era

Prokaryotes are ubiquitous in the biosphere, important for human health and drive diverse biological and environmental processes. Systematics of prokaryotes, whose origins can be traced to the discovery of microorganisms in the 17th century, has transitioned from a phenotype-based classification to a more comprehensive polyphasic taxonomy and eventually to the current genome-based taxonomic approach. This transition aligns with a foundational shift from studies focused on phenotypic traits that have limited comparative value to those using genome sequences. In this context, Bergey's Manual of Systematics of Archaea and Bacteria (BMSAB) and Bergey's International Society for Microbial Systematics (BISMiS) play a pivotal role in guiding prokaryotic systematics. This review focuses on the historical development of prokaryotic systematics with a focus on the roles of BMSAB and BISMiS. We also explore significant contributions and achievements by microbiologists, highlight the latest progress in the field and anticipate challenges and opportunities within prokaryotic systematics. Additionally, we outline five focal points of BISMiS that are aimed at addressing these challenges. In conclusion, our collaborative effort seeks to enhance ongoing advancements in prokaryotic systematics, ensuring its continued relevance and innovative characters in the contemporary landscape of genomics and bioinformatics.

Metagenomic insight into the pathogenic-related characteristics and resistome profiles within microbiome residing on the Angkor sandstone monuments in Cambodia

To reveal the characteristics of indigenous microbiome including the pathogenic-related ones on Angkor monuments in Cambodia and the distribution pattern of resistome at different locations, several sites, namely Angkor Wat, Bayon of Angkor Thom, and Prasat Preah Vihear with different exposure levels to tourists were selected to conduct the metagenomic analysis in this study. The general characteristics of the microbiome on these monuments were revealed, and the association between the environmental geo-ecological feature and the indigenous microbiome was delineated. The most common microbial groups included 6 phyla, namely Acidobacteria, Actinobacteria, Gemmatimonadetes, Nitrospirae, Proteobacteria and Verrucomicrobia on the monuments, but Firmicutes and Chlamydiae were the most dominant phyla found in bats droppings. The taxonomic family of Chitinophagaceae could serve as a signature microbial group for Preah Vihear, the less visited site. More importantly, the pathogenic-related characteristics of the microbiome residing on Angkor monuments were uncovered. A set of specific antibiotic resistance genes (ARGs) with cross-niches dispersal capacity (between the environmental microbiome and the microbiome within warm blood fauna) was identified to be high by the source tracking analysis based on ARGs profile varies in this study. Among the 10 ARG-types detected in this study, 6 of them are confined to resistance mechanism of antibiotic efflux-pump. The findings of this study provide new a new direction on public health management and implication globally at archaeological sites for tourism.

A Guide to Phylogenomic Inference

Phylogenomics aims at reconstructing the evolutionary histories of organisms taking into account whole genomes or large fractions of genomes. Phylogenomics has significant applications in fields such as evolutionary biology, systematics, comparative genomics, and conservation genetics, providing valuable insights into the origins and relationships of species and contributing to our understanding of biological diversity and evolution. This chapter surveys phylogenetic concepts and methods aimed at both gene tree and species tree reconstruction while also addressing common pitfalls, providing references to relevant computer programs. A practical phylogenomic analysis example including bacterial genomes is presented at the end of the chapter.

Taxonomic proposal for a deep branching bacterial phylogenetic lineage: transfer of the family Thermodesulfobiaceae to Thermodesulfobiales ord. nov., Thermodesulfobiia classis nov. and Thermodesulfobiota phyl. nov

The family Thermodesulfobiaceae, comprising one genus Thermodesulfobium with two validly published species, is currently assigned to order Thermoanaerobacterales within the class Clostridia of the phylum Bacillota. At the same time, the very first 16S rRNA gene sequence-based phylogenetic studies of representatives of the genus pointed out great differences between Thermodesulfobium and other members of the phylum Bacillota. Subsequent studies of new Thermodesulfobium representatives supported deep phylogenetic branching of this lineage within bacterial tree, implying that it represents a novel phylum. The results of the phylogenomic analysis performed in the frames of the present work confirm previous findings and suggest that Thermodesulfobium represents a distinct phylum-level lineage. Thus, we propose the transfer of the family Thermodesulfobiaceae to the new order Thermodesulfobiales within the new class Thermodesulfobiia and the new phylum Thermodesulfobiota.

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.

Evolution of the Inhibitory and Non-Inhibitory ε, ζ, and IF1 Subunits of the F1FO-ATPase as Related to the Endosymbiotic Origin of Mitochondria

The F1FO-ATP synthase nanomotor synthesizes >90% of the cellular ATP of almost all living beings by rotating in the “forward” direction, but it can also consume the same ATP pools by rotating in “reverse.” To prevent futile F1FO-ATPase activity, several different inhibitory proteins or domains in bacteria (ε and ζ subunits), mitochondria (IF1), and chloroplasts (ε and γ disulfide) emerged to block the F1FO-ATPase activity selectively. In this study, we analyze how these F1FO-ATPase inhibitory proteins have evolved. The phylogeny of the α-proteobacterial ε showed that it diverged in its C-terminal side, thus losing both the inhibitory function and the ATP-binding/sensor motif that controls this inhibition. The losses of inhibitory function and the ATP-binding site correlate with an evolutionary divergence of non-inhibitory α-proteobacterial ε and mitochondrial δ subunits from inhibitory bacterial and chloroplastidic ε subunits. Here, we confirm the lack of inhibitory function of wild-type and C-terminal truncated ε subunits of P. denitrificans. Taken together, the data show that ζ evolved to replace ε as the primary inhibitor of the F1FO-ATPase of free-living α-proteobacteria. However, the ζ inhibitory function was also partially lost in some symbiotic α-proteobacteria and totally lost in some strictly parasitic α-proteobacteria such as the Rickettsiales order. Finally, we found that ζ and IF1 likely evolved independently via convergent evolution before and after the endosymbiotic origin mitochondria, respectively. This led us to propose the ε and ζ subunits as tracer genes of the pre-endosymbiont that evolved into the actual mitochondria.

New Insights into the Taxonomy of Bacteria in the Genomic Era and a Case Study with Rhizobia

Since early studies, the history of prokaryotes taxonomy has dealt with many changes driven by the development of new and more robust technologies. As a result, the number of new taxa descriptions is exponentially increasing, while an increasing number of others has been subject of reclassification, demanding from the taxonomists more effort to maintain an organized hierarchical system. However, expectations are that the taxonomy of prokaryotes will acquire a more stable status with the genomic era. Other analyses may continue to be necessary to determine microbial features, but the use of genomic data might be sufficient to provide reliable taxa delineation, helping taxonomy to reach the goal of correct classification and identification. Here we describe the evolution of prokaryotes' taxonomy until the genomic era, emphasizing bacteria and taking as an example the history of rhizobia taxonomy. This example was chosen because of the importance of the symbiotic nitrogen fixation of legumes with rhizobia to the nitrogen input to both natural ecosystems and agricultural crops. This case study reports the technological advances and the methodologies used to classify and identify bacterial species and indicates the actual rules required for an accurate description of new taxa.

Testing of candidate probiotics to prevent dental caries induced by Streptococcus mutans in a mouse model

AIMS

We evaluated two species of human oral commensal streptococci in protection against dental caries induced by Streptococcus mutans.

METHODS AND RESULTS

Candidate probiotics, Streptococcus sp. A12, Streptococcus sanguinis BCC23 and an arginine deiminase mutant of BCC23 (∆arcADS) were tested for their ability to reduce S. mutans-induced caries in an established mouse model. Mice were colonized with a probiotic, challenged with S. mutans, then intermittently reinoculated with a probiotic strain. Oral colonization of each strain and autochthonous bacteria was assessed by qPCR. Both BCC23 strains, but not A12, were associated with markedly reduced sulcal caries, persistently colonized mucosal and dental biofilms, and significantly lowered S. mutans counts. All three strains enhanced mucosal colonization of autochthonous bacteria. In a follow-up experiment, when S. mutans was established first, dental and mucosal colonization of S. mutans was unaltered by a subsequent challenge with either BCC23 strain. Results between BCC23 and BCC23 ∆arcADS were equivalent.

CONCLUSIONS

BCC23 is a potential probiotic to treat patients at high caries risk. Its effectiveness is independent of ADS activity, but initial dental cleaning to enhance establishment in dental biofilms may be required.

SIGNIFICANCE AND IMPACT OF THE STUDY

In vivo testing of candidate probiotics is highly informative, as effectiveness is not always reflected by genotype or in vitro behaviors.

Testing Phylogenetic Stability with Variable Taxon Sampling

Over the past three decades, computational capabilities have grown at such a rapid rate that they have given rise to many computationally heavy science fields such as phylogenomics. As increasingly more genomes are sequenced in the three domains of life, larger and more species-complete phylogenetic tree reconstructions are leading to a better understanding of the tree of life and the evolutionary histories in deep times. However, these large datasets pose unique challenges from a modeling and computational perspective: accurately describing the evolutionary process of thousands of species is still beyond the capability of current models, while the computational burden limits our ability to test multiple hypotheses. Thus, it is common practice to reduce the size of a dataset by selecting species to represent a clade (taxon sampling). Unfortunately, this process is subjective, and comparisons of large tree of life studies show that choice and number of species used in a dataset can alter the topology obtained. Thus, taxon sampling is, in itself, a process that needs to be fully investigated to determine its effect on phylogenetic stability. Here, we present the theory and practical application of an automated pipeline that can be easily implemented to explore the effect of taxon sampling on phylogenetic reconstructions. The application of this approach was recently discussed in a study of Terrabacteria and shows its power in investigating the accuracy of deep nodes of a phylogeny.

Evolutionary and Comparative Analysis of Bacterial Nonhomologous End Joining Repair

DNA double-strand breaks (DSBs) are a threat to genome stability. In all domains of life, DSBs are faithfully fixed via homologous recombination. Recombination requires the presence of an uncut copy of duplex DNA which is used as a template for repair. Alternatively, in the absence of a template, cells utilize error-prone nonhomologous end joining (NHEJ). Although ubiquitously found in eukaryotes, NHEJ is not universally present in bacteria. It is unclear as to why many prokaryotes lack this pathway. Toward understanding what could have led to the current distribution of bacterial NHEJ, we carried out comparative genomics and phylogenetic analysis across ∼6,000 genomes. Our results show that this pathway is sporadically distributed across the phylogeny. Ancestral reconstruction further suggests that NHEJ was absent in the eubacterial ancestor and can be acquired via specific routes. Integrating NHEJ occurrence data for archaea, we also find evidence for extensive horizontal exchange of NHEJ genes between the two kingdoms as well as across bacterial clades. The pattern of occurrence in bacteria is consistent with correlated evolution of NHEJ with key genome characteristics of genome size and growth rate; NHEJ presence is associated with large genome sizes and/or slow growth rates, with the former being the dominant correlate. Given the central role these traits play in determining the ability to carry out recombination, it is possible that the evolutionary history of bacterial NHEJ may have been shaped by requirement for efficient DSB repair.

The convoluted history of haem biosynthesis

The capacity of haem to transfer electrons, bind diatomic gases, and catalyse various biochemical reactions makes it one of the essential biomolecules on Earth and one that was likely used by the earliest forms of cellular life. Since the description of haem biosynthesis, our understanding of this multi-step pathway has been almost exclusively derived from a handful of model organisms from narrow taxonomic contexts. Recent advances in genome sequencing and functional studies of diverse and previously neglected groups have led to discoveries of alternative routes of haem biosynthesis that deviate from the 'classical' pathway. In this review, we take an evolutionarily broad approach to illuminate the remarkable diversity and adaptability of haem synthesis, from prokaryotes to eukaryotes, showing the range of strategies that organisms employ to obtain and utilise haem. In particular, the complex evolutionary histories of eukaryotes that involve multiple endosymbioses and horizontal gene transfers are reflected in the mosaic origin of numerous metabolic pathways with haem biosynthesis being a striking case. We show how different evolutionary trajectories and distinct life strategies resulted in pronounced tensions and differences in the spatial organisation of the haem biosynthesis pathway, in some cases leading to a complete loss of a haem-synthesis capacity and, rarely, even loss of a requirement for haem altogether.

The metabolic footprint of Clostridia and Erysipelotrichia reveals their role in depleting sugar alcohols in the cecum

BACKGROUND

The catabolic activity of the microbiota contributes to health by aiding in nutrition, immune education, and niche protection against pathogens. However, the nutrients consumed by common taxa within the gut microbiota remain incompletely understood.

METHODS

Here we combined microbiota profiling with an un-targeted metabolomics approach to determine whether depletion of small metabolites in the cecum of mice correlated with the presence of specific bacterial taxa. Causality was investigated by engrafting germ-free or antibiotic-treated mice with complex or defined microbial communities.

RESULTS

We noted that a depletion of Clostridia and Erysipelotrichia from the gut microbiota triggered by antibiotic treatment was associated with an increase in the cecal concentration of sugar acids and sugar alcohols (polyols). Notably, when we inoculated germ-free mice with a defined microbial community of 14 Clostridia and 3 Erysipelotrichia isolates, we observed the inverse, with a marked decrease in the concentrations of sugar acids and polyols in cecal contents. The carbohydrate footprint produced by the defined microbial community was similar to that observed in gnotobiotic mice receiving a cecal microbiota transplant from conventional mice. Supplementation with sorbitol, a polyol used as artificial sweetener, increased cecal sorbitol concentrations in antibiotic-treated mice, which was abrogated after inoculation with a Clostridia isolate able to grow on sorbitol in vitro.

CONCLUSIONS

We conclude that consumption of sugar alcohols by Clostridia and Erysipelotrichia species depletes these metabolites from the intestinal lumen during homeostasis. Video abstract.

Natural selection shapes maintenance of orthologous sRNAs in divergent host-restricted bacterial genomes

Historically it has been difficult to study the evolution of bacterial small RNAs (sRNAs) across distantly related species. For example, identifying homologs of sRNAs is often difficult in genomes that have undergone multiple structural rearrangements. Also, some types of regulatory sRNAs evolve at rapid rates. The high degree of genomic synteny among divergent host-restricted bacterial lineages, including intracellular symbionts, is conducive to sRNA maintenance and homolog identification. In turn, symbiont genomes can provide us with novel insights into sRNA evolution. Here, we examine the sRNA expression profile of the obligate symbiont of psyllids, Carsonella ruddii, which has one of the smallest cellular genomes described. Using RNA-seq, we identified 36 and 32 antisense sRNAs (asRNAs) expressed by Carsonella from the psyllids Bactericera cockerelli (Carsonella-BC) and Diaphorina citri (Carsonella-DC), respectively. The majority of these asRNAs were associated with genes that are involved in essential amino acid biosynthetic pathways. Eleven of the asRNAs were conserved in both Carsonella lineages and the majority were maintained by selection. Notably, five of the corresponding coding sequences are also the targets of conserved asRNAs in a distantly related insect symbiont, Buchnera. We detected differential expression of two asRNAs for genes involved in arginine and leucine biosynthesis occurring between two distinct Carsonella-BC life stages. Using asRNAs identified in Carsonella, Buchnera, and Profftella which are all endosymbionts, and Escherichia coli, we determined that regions upstream of these asRNAs encode unique conserved patterns of AT/GC richness, GC skew, and sequence motifs which may be involved in asRNA regulation.

Sequence and evolutionary analysis of bacterial ribosomal S1 proteins

The multi-domain bacterial S1 protein is the largest and most functionally important ribosomal protein of the 30S subunit, which interacts with both mRNA and proteins. The family of ribosomal S1 proteins differs in the classical sense from a protein with tandem repeats and has a "bead-on-string" organization, where each repeat is folded into a globular domain. Based on our recent data, the study of evolutionary relationships for the bacterial phyla will provide evidence for one of the proposed theories of the evolutionary development of proteins with structural repeats: from multiple repeats of assembles to single repeats, or vice versa. In this comparative analysis of 1333 S1 sequences that were identified in 24 different phyla, we demonstrate how such phyla can form independently/dependently during evolution. To the best of our knowledge, this work is the first study of the evolutionary history of bacterial ribosomal S1 proteins. The collected and structured data can be useful to computer biologists as a resource for determining percent identity, amino acid composition and logo motifs, as well as dN/dS ratio in bacterial S1 protein. The obtained research data indicate that the evolutionary development of bacterial ribosomal S1 proteins evolved from multiple assemblies to single repeat. The presented data are integrated into the server, which can be accessed at http://oka.protres.ru:4200.

Reclassification of Catabacter hongkongensis as Christensenella hongkongensis comb. nov. based on whole genome analysis

The genera Catabacter (family 'Catabacteraceae') and Christensenella (family Christensenellaceae) are close relatives within the phylum Firmicutes. Members of these genera are strictly anaerobic, non-spore-forming and short straight rods with diverse phenotypes. Phylogenetic analysis of 16S rRNA genes suggest that Catabacter splits Christensenella into a polyphyletic clade. In an effort to ensure that family/genus names represent monophyletic clades, we performed a whole-genome based analysis of the genomes available for the cultured representatives of these genera: four species of Christensenella and two strains of Catabacter hongkongensis. A concatenated alignment of 135 shared protein sequences of single-copy core genes present in the included strains indicates that C. hongkongensis is indeed nested within the Christensenella clade. Based on their evolutionary relationship, we propose the transfer of Catabacter hongkongensis to the genus Christensenella as Christensenella hongkongensis comb. nov.

A Faster Algorithm for Computing the Kernel of Maximum Agreement Subtrees

The maximum agreement subtree method determines the consensus of a collection of phylogenetic trees by identifying maximum cardinality subsets of leaves for which all input trees agree. The trees induced by these maximum cardinality subsets are maximum agreement subtrees (MASTs). A single MAST may be misleading, since there can exist two MASTs which share almost no leaves; nevertheless, it may be impossible to inspect all MASTs, since the number of MASTs can be exponential in the number of leaves. To overcome this drawback, Swenson et al. suggested to further summarize the information common to all MASTs by their intersection, which is called the kernel agreement subtree (KAST). The construction of the KAST is the focus of this paper. Swenson et al. had an O(kn³ + n⁴ + n^{d + 1}) time algorithm for computing the KAST of k trees on n leaves, in which at least one tree has maximum degree d. In this paper, an O(kn³ + n^d)-time algorithm is presented. We demonstrate the efficiency of our algorithm on simulated trees as well as on ribosomal RNA alignments, where trees with 13,000 taxa took only hours to process, whereas the previous algorithm did not terminate after a week of computation.

In Vivo Colonization with Candidate Oral Probiotics Attenuates Colonization and Virulence of Streptococcus mutans

A collection of 113 Streptococcus strains from supragingival dental plaque of caries-free individuals were recently tested in vitro for direct antagonism of the dental caries pathogen Streptococcus mutans, and for their capacity for arginine catabolism via the arginine deiminase system (ADS). To advance their evaluation as potential probiotics, twelve strains of commensal oral streptococci with various antagonistic and ADS potentials were assessed in a mouse model for oral (i.e., oral mucosal pellicles and saliva) and dental colonization under four diets (healthy or high-sucrose, with or without prebiotic arginine). Colonization by autochthonous bacteria was also monitored. One strain failed to colonize, whereas oral colonization by the other eleven strains varied by 3 log units. Dental colonization was high for five strains regardless of diet, six strains increased colonization with at least one high-sucrose diet, and added dietary arginine decreased dental colonization of two strains. Streptococcus sp. A12 (high in vitro ADS activity and antagonism) and two engineered mutants lacking the ADS (ΔarcADS) or pyruvate oxidase-mediated H₂O₂ production (ΔspxB) were tested for competition against S. mutans UA159. A12 wild type and ΔarcADS colonized only transiently, whereas ΔspxB persisted, but without altering oral or dental colonization by S. mutans In testing four additional candidates, S. sanguinis BCC23 markedly attenuated S. mutans' oral and dental colonization, enhanced colonization of autochthonous bacteria, and decreased severity of smooth surface caries under highly cariogenic conditions. Results demonstrate the utility of the mouse model to evaluate potential probiotics, revealing little correlation between in vitro antagonism and competitiveness against S. mutans in vivo IMPORTANCE Our results demonstrate in vivo testing of potential oral probiotics can be accomplished and can yield information to facilitate the ultimate design and optimization of novel anti-caries probiotics. We show human oral commensals associated with dental health are an important source of potential probiotics that may be used to colonize patients under dietary conditions of highly varying cariogenicity. Assessment of competitiveness against dental caries pathogen Streptococcus mutans and impact on caries identified strains or genetic elements for further study. Results also uncovered strains that enhanced oral and dental colonization by autochthonous bacteria when challenged with S. mutans, suggesting cooperative interactions for future elucidation. Distinguishing a rare strain that effectively compete with S. mutans under conditions that promote caries further validates our systematic approach to more critically evaluate probiotics for use in humans.

Pan-genome analysis of Exiguobacterium reveals species delineation and genomic similarity with Exiguobacterium profundum PHM 11

The stint of the bacterial species is convoluting, but the new algorithms to calculate genome-to-genome distance (GGD) and DNA-DNA hybridization (DDH) for comparative genome analysis have rejuvenated the exploration of species and sub-species characterization. The present study reports the first whole genome sequence of Exiguobacterium profundum PHM11. PHM11 genome consist of ~ 2.92 Mb comprising 48 contigs, 47.93% G + C content. Functional annotations revealed a total of 3033 protein coding genes and 33 non-protein coding genes. Out of these, only 2316 could be characterized and others reported as hypothetical proteins. The comparative analysis of predicted proteome of PHM11 with five other Exiguobacterium sp. identified 3806 clusters, out of which the PHM11 shared a total of 2723 clusters having 1664 common clusters, 131 singletons and 928 distributed between five species. The pan-genome analysis of 70 different genomic sequences of Exigubacterium strains devoid of a species taxon was done on the basis of GGD and the DDH which identified eight genomes analogous to the PHM11 at species level and may be characterized as E. profundum. The ANI value and phylogenetic tree analysis also support the same. The results regarding pan-genome analysis provide a convincing insight for delineation of these eight strains to species.

Intrinsic laws of k-mer spectra of genome sequences and evolution mechanism of genomes

Background

K-mer spectra of DNA sequences contain important information about sequence composition and sequence evolution. We want to reveal the evolution rules of genome sequences by studying the k-mer spectra of genome sequences.

Results

The intrinsic laws of k-mer spectra of 920 genome sequences from primate to prokaryote were analyzed. We found that there are two types of evolution selection modes in genome sequences, named as CG Independent Selection and TA Independent Selection. There is a mutual inhibition relationship between CG and TA independent selections. We found that the intensity of CG and TA independent selections correlates closely with genome evolution and G + C content of genome sequences. The living habits of species are related closely to the independent selection modes adopted by species genomes. Consequently, we proposed an evolution mechanism of genomes in which the genome evolution is determined by the intensities of the CG and TA independent selections and the mutual inhibition relationship. Besides, by the evolution mechanism of genomes, we speculated the evolution modes of prokaryotes in mild and extreme environments in the anaerobic age and the evolving process of prokaryotes from anaerobic to aerobic environment on earth as well as the originations of different eukaryotes.

Conclusion

We found that there are two independent selection modes in genome sequences. The evolution of genome sequence is determined by the two independent selection modes and the mutual inhibition relationship between them.

Proposal to reclassify the proteobacterial classes Deltaproteobacteria and Oligoflexia, and the phylum Thermodesulfobacteria into four phyla reflecting major functional capabilities

The class Deltaproteobacteria comprises an ecologically and metabolically diverse group of bacteria best known for dissimilatory sulphate reduction and predatory behaviour. Although this lineage is the fourth described class of the phylum Proteobacteria , it rarely affiliates with other proteobacterial classes and is frequently not recovered as a monophyletic unit in phylogenetic analyses. Indeed, one branch of the class Deltaproteobacteria encompassing Bdellovibrio-like predators was recently reclassified into a separate proteobacterial class, the Oligoflexia . Here we systematically explore the phylogeny of taxa currently assigned to these classes using 120 conserved single-copy marker genes as well as rRNA genes. The overwhelming majority of markers reject the inclusion of the classes Deltaproteobacteria and Oligoflexia in the phylum Proteobacteria . Instead, the great majority of currently recognized members of the class Deltaproteobacteria are better classified into four novel phylum-level lineages. We propose the names Desulfobacterota phyl. nov. and Myxococcota phyl. nov. for two of these phyla, based on the oldest validly published names in each lineage, and retain the placeholder name SAR324 for the third phylum pending formal description of type material. Members of the class Oligoflexia represent a separate phylum for which we propose the name Bdellovibrionota phyl. nov. based on priority in the literature and general recognition of the genus Bdellovibrio. Desulfobacterota phyl. nov. includes the taxa previously classified in the phylum Thermodesulfobacteria , and these reclassifications imply that the ability of sulphate reduction was vertically inherited in the Thermodesulfobacteria rather than laterally acquired as previously inferred. Our analysis also indicates the independent acquisition of predatory behaviour in the phyla Myxococcota and Bdellovibrionota, which is consistent with their distinct modes of action. This work represents a stable reclassification of one of the most taxonomically challenging areas of the bacterial tree and provides a robust framework for future ecological and systematic studies.

Understanding LrgAB Regulation of Streptococcus mutans Metabolism

Lack of LrgAB renders cariogenic Streptococcus mutans more sensitive to oxidative stress, as well as limits the capacity of this organism to re-uptake pyruvate upon starvation. This study was aimed at investigating the ecological and metabolic contribution of LrgAB to competitive fitness, using S. mutans strains, that either lack or overexpress lrgAB. These experiments revealed that impaired aerobic growth of the ΔlrgAB mutant can be effectively restored by supplementation of pyruvate, and that perturbated expression of lrgAB significantly affects pyruvate flux and the conversion of pyruvate to acetyl-CoA by the Pdh pathway, verifying that LrgAB is closely linked to pyruvate catabolism. In vitro competition assays revealed that LrgAB plays an important role in S. mutans competition with H₂O₂-producing S. gordonii, an interaction which can also be modulated by external pyruvate. However, no obvious competitive disadvantage was observed against S. gordonii by either the S. mutans lrgAB mutant or lrgAB overexpression strain in vivo using a mouse caries model. Organic acid analysis of mouse dental biofilms revealed that metabolites produced by the host and/or dental plaque microbiota could complement the deficiency of a lrgAB mutant, and favored S. mutans establishment compared to S. gordonii. Collectively, these results reinforce the importance of the oral microbiota and the metabolic environment in the oral cavity battleground, and highlight that pyruvate uptake through LrgAB may be crucial for interspecies competition that drives niche occupancy.

Choice of species affects phylogenetic stability of deep nodes: an empirical example in Terrabacteria

MOTIVATION

The promise of higher phylogenetic stability through increased dataset sizes within tree of life (TOL) reconstructions has not been fulfilled. Among the many possible causes are changes in species composition (taxon sampling) that could influence phylogenetic accuracy of the methods by altering the relative weight of the evolutionary histories of each individual species. This effect would be stronger in clades that are represented by few lineages, which is common in many prokaryote phyla. Indeed, phyla with fewer taxa showed the most discordance among recent TOL studies. We implemented an approach to systematically test how the identity of taxa among a larger dataset and the number of taxa included affected the accuracy of phylogenetic reconstruction.

RESULTS

Utilizing an empirical dataset within Terrabacteria we found that even within scenarios consisting of the same number of taxa, the species used strongly affected phylogenetic stability. Furthermore, we found that trees with fewer species were more dissimilar to the tree produced from the full dataset. These results hold even when the tree is composed by many phyla and only one of them is being altered. Thus, the effect of taxon sampling in one group does not seem to be buffered by the presence of many other clades, making this issue relevant even to very large datasets. Our results suggest that a systematic evaluation of phylogenetic stability through taxon resampling is advisable even for very large datasets.

AVAILABILITY AND IMPLEMENTATION

https://github.com/BlabOaklandU/PATS.git.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Impact of C-terminal amino acid composition on protein expression in bacteria

The C-terminal sequence of a protein is involved in processes such as efficiency of translation termination and protein degradation. However, the general relationship between features of this C-terminal sequence and levels of protein expression remains unknown. Here, we identified C-terminal amino acid biases that are ubiquitous across the bacterial taxonomy (1,582 genomes). We showed that the frequency is higher for positively charged amino acids (lysine, arginine), while hydrophobic amino acids and threonine are lower. We then studied the impact of C-terminal composition on protein levels in a library of Mycoplasma pneumoniae mutants, covering all possible combinations of the two last codons. We found that charged and polar residues, in particular lysine, led to higher expression, while hydrophobic and aromatic residues led to lower expression, with a difference in protein levels up to fourfold. We further showed that modulation of protein degradation rate could be one of the main mechanisms driving these differences. Our results demonstrate that the identity of the last amino acids has a strong influence on protein expression levels.

Genomic and Metabolomic Analyses of Natural Products in Nodularia spumigena Isolated from a Shrimp Culture Pond

The bloom-forming cyanobacterium Nodularia spumigena CENA596 encodes the biosynthetic gene clusters (BGCs) of the known natural products nodularins, spumigins, anabaenopeptins/namalides, aeruginosins, mycosporin-like amino acids, and scytonemin, along with the terpenoid geosmin. Targeted metabolomics confirmed the production of these metabolic compounds, except for the alkaloid scytonemin. Genome mining of N. spumigena CENA596 and its three closely related Nodularia strains—two planktonic strains from the Baltic Sea and one benthic strain from Japanese marine sediment—revealed that the number of BGCs in planktonic strains was higher than in benthic one. Geosmin—a volatile compound with unpleasant taste and odor—was unique to the Brazilian strain CENA596. Automatic annotation of the genomes using subsystems technology revealed a related number of coding sequences and functional roles. Orthologs from the Nodularia genomes are involved in the primary and secondary metabolisms. Phylogenomic analysis of N. spumigena CENA596 based on 120 conserved protein sequences positioned this strain close to the Baltic Nodularia. Phylogeny of the 16S rRNA genes separated the Brazilian CENA596 strain from those of the Baltic Sea, despite their high sequence identities (99% identity, 100% coverage). The comparative analysis among planktic Nodularia strains showed that their genomes were considerably similar despite their geographically distant origin.

Carbon Use Efficiency and Its Temperature Sensitivity Covary in Soil Bacteria

The strategy that microbial decomposers take with respect to using substrate for growth versus maintenance is one essential biological determinant of the propensity of carbon to remain in soil. To quantify the environmental sensitivity of this key physiological trade-off, we characterized the carbon use efficiency (CUE) of 23 soil bacterial isolates across seven phyla at three temperatures and with up to four substrates. Temperature altered CUE in both an isolate-specific manner and a substrate-specific manner. We searched for genes correlated with the temperature sensitivity of CUE on glucose and deemed those functional genes which were similarly correlated with CUE on other substrates to be validated as markers of CUE. Ultimately, we did not identify any such robust functional gene markers of CUE or its temperature sensitivity. However, we found a positive correlation between rRNA operon copy number and CUE, opposite what was expected. We also found that inefficient taxa increased their CUE with temperature, while those with high CUE showed a decrease in CUE with temperature. Together, our results indicate that CUE is a flexible parameter within bacterial taxa and that the temperature sensitivity of CUE is better explained by observed physiology than by genomic composition across diverse taxa. We conclude that the bacterial CUE response to temperature and substrate is more variable than previously thought.IMPORTANCE Soil microbes respond to environmental change by altering how they allocate carbon to growth versus respiration-or carbon use efficiency (CUE). Ecosystem and Earth System models, used to project how global soil C stocks will continue to respond to the climate crisis, often assume that microbes respond homogeneously to changes in the environment. In this study, we quantified how CUE varies with changes in temperature and substrate quality in soil bacteria and evaluated why CUE characteristics may differ between bacterial isolates and in response to altered growth conditions. We found that bacterial taxa capable of rapid growth were more efficient than those limited to slow growth and that taxa with high CUE were more likely to become less efficient at higher temperatures than those that were less efficient to begin with. Together, our results support the idea that the CUE temperature response is constrained by both growth rate and CUE and that this partly explains how bacteria acclimate to a warming world.

SWeeP: representing large biological sequences datasets in compact vectors

Vectoral and alignment-free approaches to biological sequence representation have been explored in bioinformatics to efficiently handle big data. Even so, most current methods involve sequence comparisons via alignment-based heuristics and fail when applied to the analysis of large data sets. Here, we present “Spaced Words Projection (SWeeP)”, a method for representing biological sequences using relatively small vectors while preserving intersequence comparability. SWeeP uses spaced-words by scanning the sequences and generating indices to create a higher-dimensional vector that is later projected onto a smaller randomly oriented orthonormal base. We constructed phylogenetic trees for all organisms with mitochondrial and bacterial protein data in the NCBI database. SWeeP quickly built complete and accurate trees for these organisms with low computational cost. We compared SWeeP to other alignment-free methods and Sweep was 10 to 100 times quicker than the other techniques. A tool to build SWeeP vectors is available at https://sourceforge.net/projects/spacedwordsprojection/.

Characterization of Bradyrhizobium strains indigenous to Western Australia and South Africa indicates remarkable genetic diversity and reveals putative new species

Bradyrhizobium are N₂-fixing microsymbionts of legumes with relevant applications in agricultural sustainability, and we investigated the phylogenetic relationships of conserved and symbiotic genes of 21 bradyrhizobial strains. The study included strains from Western Australia (WA), isolated from nodules of Glycine spp. the country is one genetic center for the genus and from nodules of other indigenous legumes grown in WA, and strains isolated from forage Glycine sp. grown in South Africa. The 16S rRNA phylogeny divided the strains in two superclades, of B. japonicum and B. elkanii, but with low discrimination among the species. The multilocus sequence analysis (MLSA) with four protein-coding housekeeping genes (dnaK, glnII, gyrB and recA) pointed out seven groups as putative new species, two within the B. japonicum, and five within the B. elkanii superclades. The remaining eleven strains showed higher similarity with six species, B. lupini, B. liaoningense, B. yuanmingense, B. subterraneum, B. brasilense and B. retamae. Phylogenetic analysis of the nodC symbiotic gene clustered 13 strains in three different symbiovars (sv. vignae, sv. genistearum and sv. retamae), while seven others might compose new symbiovars. The genetic profiles of the strains evaluated by BOX-PCR revealed high intra- and interspecific diversity. The results point out the high level of diversity still to be explored within the Bradyrhizobium genus, and further studies might confirm new species and symbiovars.

Sphingomonas: from diversity and genomics to functional role in environmental remediation and plant growth

The species belonging to the Sphingomonas genus possess multifaceted functions ranging from remediation of environmental contaminations to producing highly beneficial phytohormones, such as sphingan and gellan gum. Recent studies have shown an intriguing role of Sphingomonas species in the degradation of organometallic compounds. However, the actual biotechnological potential of this genus requires further assessment. Some of the species from the genus have also been noted to improve plant-growth during stress conditions such as drought, salinity, and heavy metals in agricultural soil. This role has been attributed to their potential to produce plant growth hormones e.g. gibberellins and indole acetic acid. However, the current literature is scattered, and some of the important areas, such as taxonomy, phylogenetics, genome mapping, and cellular transport systems, are still being overlooked in terms of elucidation of the mechanisms behind stress-tolerance and bioremediation. In this review, we elucidated the prospective role and function of this genus for improved utilization during environmental biotechnology.

Comparative Genomic Analysis of Soil Dwelling Bacteria Utilizing a Combinational Codon Usage and Molecular Phylogenetic Approach Accentuating on Key Housekeeping Genes

Soil is a diversified and complex ecological niche, home to a myriad of microorganisms particularly bacteria. The physico-chemical complexities of soil results in a plethora of physiological variations to exist within the different types of soil dwelling bacteria, giving rise to a wide variation in genome structure and complexity. This serves as an attractive proposition to analyze and compare the genome of a large number soil bacteria to comprehend their genome complexity and evolution. In this study a combination of codon usage and molecular phylogenetics of the whole genome and key housekeeping genes like infB (translation initiation factor 2), trpB (tryptophan synthase, beta subunit), atpD (ATP synthase, beta subunit), and rpoB (RNA polymerase, beta subunit) of 92 soil bacterial species spread across the entire eubacterial domain and residing in different soil types was performed. The results indicated the direct relationship of genome size with codon bias and coding frequency in the studied bacteria. The codon usage profile demonstrated by the gene trpB was found to be relatively different from the rest of the housekeeping genes with a large number of bacteria having a greater percentage of genes with Nc values less than the Nc of trpB. The results from the overall codon usage bias profile also depicted that the codon usage bias in the key housekeeping genes of soil bacteria was majorly due to selectional pressure and not mutation. The analysis of hydrophobicity of the gene product encoded by the rpoB coding sequences demonstrated tight clustering across all the soil bacteria suggesting conservation of protein structure for maintenance of form and function. The phylogenetic affinities inferred using 16S rRNA gene and the housekeeping genes demonstrated conflicting signals with trpB gene being the noisiest one. The housekeeping gene atpD was found to depict the least amount of evolutionary change in the soil bacteria considered in this study except in two Clostridium species. The phylogenetic and codon usage analysis of the soil bacteria consistently demonstrated the relatedness of Azotobacter chroococcum with different species of the genus Pseudomonas.

The Environmental Exposures and Inner- and Intercity Traffic Flows of the Metro System May Contribute to the Skin Microbiome and Resistome

The skin functions as the primary interface between the human body and the external environment. To understand how the microbiome varies within urban mass transit and influences the skin microbiota, we profiled the human palm microbiome after contact with handrails within the Hong Kong Mass Transit Railway (MTR) system. Intraday sampling time was identified as the primary determinant of the variation and recurrence of the community composition, whereas human-associated species and clinically important antibiotic resistance genes (ARGs) were captured as p.m. signatures. Line-specific signatures were notably correlated with line-specific environmental exposures and city characteristics. The sole cross-border line appeared as an outlier in most analyses and showed high relative abundance and a significant intraday increment of clinically important ARGs (24.1%), suggesting potential cross-border ARG transmission, especially for tetracycline and vancomycin resistance. Our study provides an important reference for future public health strategies to mitigate intracity and cross-border pathogen and ARG transmission.

Genome-Based Metabolic Reconstruction of a Novel Uncultivated Freshwater Magnetotactic coccus "Ca. Magnetaquicoccus inordinatus" UR-1, and Proposal of a Candidate Family "Ca. Magnetaquicoccaceae"

Magnetotactic bacteria are widely represented microorganisms that have the ability to synthesize magnetosomes. The magnetotactic cocci of the order Magnetococcales are the most frequently identified, but their classification remains unclear due to the low number of cultivated representatives. This paper reports the analysis of an uncultivated magnetotactic coccus UR-1 collected from the Uda River (in eastern Siberia). Genome analyses of this bacterium and comparison to the available Magnetococcales genomes identified a novel species called "Ca. Magnetaquicoccus inordinatus," and a delineated candidate family "Ca. Magnetaquicoccaceae" within the order Magnetococcales is proposed. We used average amino acid identity values <55-56% and <64-65% as thresholds for the separation of families and genera, respectively, within the order Magnetococcales. Analyses of the genome sequence of UR-1 revealed a potential ability for a chemolithoautotrophic lifestyle, with the oxidation of a reduced sulfur compound and carbon assimilation by rTCA. A nearly complete magnetosome genome island, containing a set of mam and mms genes, was also identified. Further comparative analyses of the magnetosome genes showed vertical inheritance as well as horizontal gene transfer as the evolutionary drivers of magnetosome biomineralization genes in strains of the order Magnetococcales.

Fast and accurate genome comparison using genome images: The Extended Natural Vector Method

Using numerical methods for genome comparison has always been of importance in bioinformatics. The Chaos Game Representation (CGR) is an effective genome sequence mapping technology, which converts genome sequences to CGR images. To each CGR image, we associate a vector called an Extended Natural Vector (ENV). The ENV is based on the distribution of intensity values. This mapping produces a one-to-one correspondence between CGR images and their ENVs. We define the distance between two DNA sequences as the distance between their associated ENVs. We cluster and classify several datasets including Influenza A viruses, Bacillus genomes, and Conoidea mitochondrial genomes to build their phylogenetic trees. Results show that our ENV combining CGR method (CGR-ENV) compares favorably in classification accuracy and efficiency against the multiple sequence alignment (MSA) method and other alignment-free methods. The research provides significant insights into the study of phylogeny, evolution, and efficient DNA comparison algorithms for large genomes.

Visualization and prediction of CRISPR incidence in microbial trait-space to identify drivers of antiviral immune strategy

Bacteria and archaea are locked in a near-constant battle with their viral pathogens. Despite previous mechanistic characterization of numerous prokaryotic defense strategies, the underlying ecological drivers of different strategies remain largely unknown and predicting which species will take which strategies remains a challenge. Here, we focus on the CRISPR immune strategy and develop a phylogenetically-corrected machine learning approach to build a predictive model of CRISPR incidence using data on over 100 traits across over 2600 species. We discover a strong but hitherto-unknown negative interaction between CRISPR and aerobicity, which we hypothesize may result from interference between CRISPR-associated proteins and non-homologous end-joining DNA repair due to oxidative stress. Our predictive model also quantitatively confirms previous observations of an association between CRISPR and temperature. Finally, we contrast the environmental associations of different CRISPR system types (I, II, III) and restriction modification systems, all of which act as intracellular immune systems.

Phylogenetic determinants of toxin gene distribution in genomes of Brevibacillus laterosporus

Brevibacillus laterosporus is a globally ubiquitous, spore forming bacterium, strains of which have shown toxic activity against invertebrates and microbes and several have been patented due to their commercial potential. Relatively little is known about this bacterium. Here, we examined the genomes of six published and five newly determined genomes of B. laterosporus, with an emphasis on the relationships between known and putative toxin encoding genes, as well as the phylogenetic relationships between strains. Phylogenetically, strain relationships are similar using average nucleotide identity (ANI) values and multi-gene approaches, although PacBio sequencing revealed multiple copies of the 16S rDNA gene which lessened utility at the strain level. Based on ANI values, the New Zealand isolates were distant from other isolates and may represent a new species. While all of the genomes examined shared some putative toxicity or virulence related proteins, many specific genes were only present in a subset of strains.

•

We examined genomes of 11 Brevibacillus laterosporus, a bacterium which is antagonistic to invertebrates and/or microbes

•

Multiple phylogenetic methods showed New Zealand isolates more distant than all other isolates

•

Each genome could contain 11–13 copies of the 16S rDNA gene, some of which were not identical

•

Many putative toxin encoding genes were present in the genomes, but the toxin complement varied from isolate to isolate

•

Variation in occurrence of toxin-encoding genes indicates the potential to find strains with new combinations of activities

Marker genes as predictors of shared genomic function

Background

Although high-throughput marker gene studies provide valuable insight into the diversity and relative abundance of taxa in microbial communities, they do not provide direct measures of their functional capacity. Recently, scientists have shown a general desire to predict functional profiles of microbial communities based on phylogenetic identification inferred from marker genes, and recent tools have been developed to link the two. However, to date, no large-scale examination has quantified the correlation between the marker gene based taxonomic identity and protein coding gene conservation. Here we utilize 4872 representative prokaryotic genomes from NCBI to investigate the relationship between marker gene identity and shared protein coding gene content.

Results

Even at 99–100% marker gene identity, genomes share on average less than 75% of their protein coding gene content. This occurs regardless of the marker gene(s) used: V4 region of the 16S rRNA, complete 16S rRNA, or single copy orthologs through a multi-locus sequence analysis. An important aspect related to this observation is the intra-organism variation of 16S copies from a single genome. Although the majority of 16S copies were found to have high sequence similarity (> 99%), several genomes contained copies that were highly diverged (< 97% identity).

Conclusions

This is the largest comparison between marker gene similarity and shared protein coding gene content to date. The study highlights the limitations of inferring a microbial community’s functions based on marker gene phylogeny. The data presented expands upon the results of previous studies that examined one or few bacterial species and supports the hypothesis that 16S rRNA and other marker genes cannot be directly used to fully predict the functional potential of a bacterial community.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5641-1) contains supplementary material, which is available to authorized users.

Genome-informed Bradyrhizobium taxonomy: where to from here?

Bradyrhizobium is thought to be the largest and most diverse rhizobial genus, but this is not reflected in the number of described species. Although it was one of the first rhizobial genera recognised, its taxonomy remains complex. Various contemporary studies are showing that genome sequence information may simplify taxonomic decisions. Therefore, the growing availability of genomes for Bradyrhizobium will likely aid in the delineation and characterization of new species. In this study, we addressed two aims: first, we reviewed the availability and quality of available genomic resources for Bradyrhizobium. This was achieved by comparing genome sequences in terms of sequencing technologies used and estimated level of completeness for inclusion in genome-based phylogenetic analyses. Secondly, we utilized these genomes to investigate the taxonomic standing of Bradyrhizobium in light of its diverse lifestyles. Although genome sequences differed in terms of their quality and completeness, our data indicate that the use of these genome sequences is adequate for taxonomic purposes. By using these resources, we inferred a fully resolved, well-supported phylogeny. It separated Bradyrhizobium into seven lineages, three of which corresponded to the so-called supergroups known for the genus. Wide distribution of key lifestyle traits such as nodulation, nitrogen fixation and photosynthesis revealed that these traits have complicated evolutionary histories. We present the first robust Bradyrhizobium species phylogeny based on genome sequence information for investigating the evolution of this important assemblage of bacteria. Furthermore, this study provides the basis for using genome sequence information as a resource to make important taxonomic decisions, particularly at the species and genus levels.

Prot-SpaM: fast alignment-free phylogeny reconstruction based on whole-proteome sequences

Word-based or 'alignment-free' sequence comparison has become an active research area in bioinformatics. While previous word-frequency approaches calculated rough measures of sequence similarity or dissimilarity, some new alignment-free methods are able to accurately estimate phylogenetic distances between genomic sequences. One of these approaches is Filtered Spaced Word Matches. Here, we extend this approach to estimate evolutionary distances between complete or incomplete proteomes; our implementation of this approach is called Prot-SpaM. We compare the performance of Prot-SpaM to other alignment-free methods on simulated sequences and on various groups of eukaryotic and prokaryotic taxa. Prot-SpaM can be used to calculate high-quality phylogenetic trees for dozens of whole-proteome sequences in a matter of seconds or minutes and often outperforms other alignment-free approaches. The source code of our software is available through Github: https://github.com/jschellh/ProtSpaM.

Cyanobacterial phylogenetic analysis based on phylogenomics approaches render evolutionary diversification and adaptation: an overview of representative orders

Phylogenetic studies based on a definite set of marker genes usually reconstruct evolutionary relationships among the prokaryotic species. Based on specific target sequences, such studies represent variations and allow identification of similarities or dissimilarities in organisms. With the advent of completely sequenced genomes and accumulation of information on whole prokaryotic genomes, phylogenetic reconstructions should be considered more reliable if they are ideally based on entire genomes to resolve phylogenetic interest. We applied phylogenomics approaches taking into account completely sequenced cyanobacterial genomes to reconstruct underlying species that represented major taxonomic classes and belonged to distinctly different habitats (freshwater, marine, soils, and rocks). We did not rely on describing phylogeny of all representative class of cyanobacterial species on the basis of only ribosomal gene, 16S rDNA gene. In contrast, we analyzed combined molecular marker and phylogenomics approaches (genome alignment, gene content and gene order, composition vector and protein domain content) for accurately inferring phylogenetic relationship of species. We have shown that this approach reflects the impact of evolution on the organisms and considers connects with the ecological adaptation in cyanobacteria in different habitats. Analysis revealed that the members from marine habitat occupy different profile than those from freshwater. Impact of GC content and genomic repetitiveness over the diversification of cyanobacterial species and their possible role in adaptation was also reflected. Members occupying similar habitats cover more evolutionary distance together and also evolve various strategies for adaptation and survival either through genomic repetitiveness or preferences for genes of particular functions or modified GC content. Genomes undergo different changes for their adaptation in diverse habitats.

A phylogenomic and ecological analysis of the globally abundant Marine Group II archaea (Ca. Poseidoniales ord. nov.)

Marine Group II (MGII) archaea represent the most abundant planktonic archaeal group in ocean surface waters, but our understanding of the group has been limited by a lack of cultured representatives and few sequenced genomes. Here, we conducted a comparative phylogenomic analysis of 270 recently available MGII metagenome-assembled genomes (MAGs) to investigate their evolution and ecology. Based on a rank-normalised genome phylogeny, we propose that MGII is an order-level lineage for which we propose the name Candidatus Poseidoniales (after Gr. n. Poseidon, God of the sea), comprising the families Candidatus Poseidonaceae fam. nov. (formerly subgroup MGIIa) and Candidatus Thalassarchaeaceae fam. nov. (formerly subgroup MGIIb). Within these families, 21 genera could be resolved, many of which had distinct biogeographic ranges and inferred nutrient preferences. Phylogenetic analyses of key metabolic functions suggest that the ancestor of Ca. Poseidoniales was a surface water-dwelling photoheterotroph that evolved to occupy multiple related ecological niches based primarily on spectral tuning of proteorhodopsin genes. Interestingly, this adaptation appears to involve an overwrite mechanism whereby an existing single copy of the proteorhodopsin gene is replaced by a horizontally transferred copy, which in many instances should allow an abrupt change in light absorption capacity. Phototrophy was lost entirely from five Ca. Poseidoniales genera coinciding with their adaptation to deeper aphotic waters. We also report the first instances of nitrate reductase in two genera acquired via horizontal gene transfer (HGT), which was a potential adaptation to oxygen limitation. Additional metabolic traits differentiating families and genera include flagellar-based adhesion, transporters, and sugar, amino acid, and peptide degradation. Our results suggest that HGT has shaped the evolution of Ca. Poseidoniales to occupy a variety of ecological niches and to become the most successful archaeal lineage in ocean surface waters.

The Genome Sequence of "Candidatus Fokinia solitaria": Insights on Reductive Evolution in Rickettsiales

"Candidatus Fokinia solitaria" is an obligate intracellular endosymbiont of a unicellular eukaryote, a ciliate of the genus Paramecium. Here, we present the genome sequence of this bacterium and subsequent analysis. Phylogenomic analysis confirmed the previously reported positioning of the symbiont within the "Candidatus Midichloriaceae" family (order Rickettsiales), as well as its high sequence divergence from other members of the family, indicative of fast sequence evolution. Consistently with this high evolutionary rate, a comparative genomic analysis revealed that the genome of this symbiont is the smallest of the Rickettsiales to date. The reduced genome does not present flagellar genes, nor the pathway for the biosynthesis of lipopolysaccharides (present in all the other so far sequenced members of the family "Candidatus Midichloriaceae") or genes for the Krebs cycle (present, although not always complete, in Rickettsiales). These results indicate an evolutionary trend toward a stronger dependence on the host, in comparison with other members of the family. Two alternative scenarios are compatible with our results; "Candidatus Fokinia solitaria" could be either a recently evolved, vertically transmitted mutualist, or a parasite with a high host-specificity.

A standardized bacterial taxonomy based on genome phylogeny substantially revises the tree of life

Taxonomy is an organizing principle of biology and is ideally based on evolutionary relationships among organisms. Development of a robust bacterial taxonomy has been hindered by an inability to obtain most bacteria in pure culture and, to a lesser extent, by the historical use of phenotypes to guide classification. Culture-independent sequencing technologies have matured sufficiently that a comprehensive genome-based taxonomy is now possible. We used a concatenated protein phylogeny as the basis for a bacterial taxonomy that conservatively removes polyphyletic groups and normalizes taxonomic ranks on the basis of relative evolutionary divergence. Under this approach, 58% of the 94,759 genomes comprising the Genome Taxonomy Database had changes to their existing taxonomy. This result includes the description of 99 phyla, including six major monophyletic units from the subdivision of the Proteobacteria, and amalgamation of the Candidate Phyla Radiation into a single phylum. Our taxonomy should enable improved classification of uncultured bacteria and provide a sound basis for ecological and evolutionary studies.

Evaluation of the infB and rpsB gene fragments as genetic markers intended for identification and phylogenetic analysis of particular representatives of the order Lactobacillales

Detailed differentiation, classification, and phylogenetic analysis of the order Lactobacillales are performed using molecular techniques that involve the comparison of whole genomes, multilocus sequence analysis, DNA-DNA hybridisation, and 16S rRNA sequencing. Despite the wide application of the latter two techniques, issues associated with them are extensively discussed. Although complete genomic analyses are the most appropriate for phylogenetic studies, they are time-consuming and require high levels of expertise. Many phylogenetic/identification markers have been proposed for enterococci, lactobacilli, streptococci, and lactobacilli. However, none have been established for vagococci and some genera within the order Lactobacillales. The objective of the study was to find novel alternative housekeeping genes for classification, typing, and phylogenetic analysis of selected genera within the order Lactobacillales. We designed primers flanking variable regions of the infB (504 nt) and rpsB (333 nt) genes and amplified and sequenced them in 56 strains of different genera within the order Lactobacillales. Statistical analysis and characteristics of the gene regions suggested that they could be used for taxonomic purposes. Phylogenetic analyses, including assessment of (in)congruence between individual phylogenetic trees indicated the possibility of using the concatenation of the two genes as an alternative tool for the evaluation of phylogeny compared with the 16S rRNA gene representing the standard phylogenetic marker of prokaryotes. Moreover, infB, rpsB regions and their concatenate were phylogenetically consistent with two widely applied alternative genetic markers in taxonomy of particular Lactobacillales genera encoding the 60 kDa chaperonin protein (GroEL-hsp60) and phenylalanyl-tRNA synthetase, alpha subunit (pheS).

Phylogenomic analysis on the exceptionally diverse fish clade Gobioidei (Actinopterygii: Gobiiformes) and data-filtering based on molecular clocklikeness

The use of genome-scale data to infer phylogenetic relationships has gained in popularity in recent years due to the progress made in target-gene capture and sequencing techniques. Data filtering, the approach of excluding data inconsistent with the model from analyses, presumably could alleviate problems caused by systematic errors in phylogenetic inference. Different data filtering criteria, such as those based on evolutionary rate and molecular clocklikeness as well as others have been proposed for selecting useful phylogenetic markers, yet few studies have tested these criteria using phylogenomic data. We developed a novel set of single-copy nuclear coding markers to capture thousands of target genes in gobioid fishes, a species-rich lineages of vertebrates, and tested the effects of data-filtering methods based on substitution rate and molecular clocklikeness while attempting to control for the compounding effects of missing data and variation in locus length. We found that molecular clocklikeness was a better predictor than overall substitution rate for phylogenetic usefulness of molecular markers in our study. In addition, when the 100 best ranked loci for our predictors were concatenated and analyzed using maximum likelihood, or combined in a coalescent-based species-tree analysis, the resulting trees showed a well-resolved topology of Gobioidei that mostly agrees with previous studies. However, trees generated from the 100 least clocklike frequently recovered conflicting, and in some cases clearly erroneous topologies with strong support, thus indicating strong systematic biases in those datasets. Collectively these results suggest that data filtering has the potential improve the performance of phylogenetic inference when using both a concatenation approach as well as methods that rely on input from individual gene trees (i.e. coalescent species-tree approaches), which may be preferred in scenarios where incomplete lineage sorting is likely to be an issue.

The importance of designating type material for uncultured taxa

Naming of uncultured Bacteria and Archaea is often inconsistent with the International Code of Nomenclature of Prokaryotes. The recent practice of proposing names for higher taxa without designation of lower ranks and nomenclature types is one of the most important inconsistencies that needs to be addressed to avoid nomenclatural instability. The Code requires names of higher taxa up to the rank of class to be derived from the type genus name, with a proposal pending to formalise this requirement for the rank of phylum. Designation of nomenclature types is crucial for providing priority to names and ensures their uniqueness and stability. However, only legitimate names proposed for axenic cultures can be used for this purpose. Candidatus names reserved for taxa lacking cultured representatives may be granted this right if recent proposals to use genome sequences as type material are endorsed, thereby allowing the Code to be fully applied to lineages represented by metagenome-assembled genomes (MAGs) or single amplified genomes (SAGs). Genome quality standards need to be considered to ensure unambiguous assignment of type material. Here, we illustrate the recommended practice by proposing nomenclature type material for four major uncultured prokaryotic lineages based on high-quality MAGs in accordance with the Code.

Proposal for a new classification of a deep branching bacterial phylogenetic lineage: transfer of Coprothermobacter proteolyticus and Coprothermobacter platensis to Coprothermobacteraceae fam. nov., within Coprothermobacterales ord. nov., Coprothermobacteria classis nov. and Coprothermobacterota phyl. nov. and emended description of the family Thermodesulfobiaceae

The genus Coprothermobacter (initially named Thermobacteroides) is currently placed within the phylum Firmicutes. Early 16S rRNA gene based phylogenetic studies pointed out the great differences between Coprothermobacter and other members of the Firmicutes, revealing that it constitutes a new deep branching lineage. Over the years, several studies based on 16S rRNA gene and whole genome sequences have indicated that Coprothermobacter is very distant phylogenetically to all other bacteria, supporting its placement in a distinct deeply rooted novel phylum. In view of this, we propose its allocation to the new family Coprothermobacteraceae within the novel order Coprothermobacterales, the new class Coprothermobacteria, and the new phylum Coprothermobacterota, and an emended description of the family Thermodesulfobiaceae.

Gene Co-occurrence Networks Reflect Bacteriophage Ecology and Evolution

Bacteriophages are the most abundant and diverse biological entities on the planet, and new phage genomes are being discovered at a rapid pace. As more phage genomes are published, new methods are needed for placing these genomes in an ecological and evolutionary context. Phages are difficult to study by phylogenetic methods, because they exchange genes regularly, and no single gene is conserved across all phages. Here, we demonstrate how gene-level networks can provide a high-resolution view of phage genetic diversity and offer a novel perspective on virus ecology. We focus our analyses on virus host range and show how network topology corresponds to host relatedness, how to find groups of genes with the strongest host-specific signatures, and how this perspective can complement phage host prediction tools. We discuss extensions of gene network analysis to predicting the emergence of phages on new hosts, as well as applications to features of phage biology beyond host range.

Bacteriophages (phages) are viruses that infect bacteria, and they are critical drivers of bacterial evolution and community structure. It is generally difficult to study phages by using tree-based methods, because gene exchange is common, and no single gene is shared among all phages. Instead, networks offer a means to compare phages while placing them in a broader ecological and evolutionary context. In this work, we build a network that summarizes gene sharing across phages and test how a key constraint on phage ecology, host range, corresponds to the structure of the network. We find that the network reflects the relatedness among phage hosts, and phages with genes that are closer in the network are likelier to infect similar hosts. This approach can also be used to identify genes that affect host range, and we discuss possible extensions to analyze other aspects of viral ecology.