Phylogeny and Expansion of Serine/Threonine Kinases in Phagocytotic Bacteria in the Phylum Planctomycetota.

The recently isolated bacterium "Candidatus Uabimicrobium amorphum" is the only known prokaryote that can engulf other bacterial cells. Its proteome contains a high fraction of proteins involved in signal transduction systems, which is a feature normally associated with multicellularity in eukaryotes. Here, we present a protein-based phylogeny which shows that "Ca. Uabimicrobium amorphum" represents an early diverging lineage that clusters with the Saltatorellus clade within the phylum Planctomycetota. A gene flux analysis indicated a gain of 126 protein families for signal transduction functions in "Ca. Uabimicrobium amorphum", of which 66 families contained eukaryotic-like Serine/Threonine kinases with Pkinase domains. In total, we predicted 525 functional Serine/Threonine kinases in "Ca. Uabimicrobium amorphum", which represent 8% of the proteome and is the highest fraction of Serine/Threonine kinases in a bacterial proteome. The majority of Serine/Threonine kinases in this species are membrane proteins and 30% contain long, tandem arrays of WD40 or TPR domains. The pKinase domain was predicted to be located in the cytoplasm, while the WD40 and TPR domains were predicted to be located in the periplasm. Such domain combinations were also identified in the Serine/Threonine kinases of other species in the Planctomycetota, although in much lower abundances. A phylogenetic analysis of the Serine/Threonine kinases in the Planctomycetota inferred from the Pkinase domain alone provided support for lineage-specific expansions of the Serine/Threonine kinases in "Ca. Uabimicrobium amorphum". The results imply that expansions of eukaryotic-like signal transduction systems are not restricted to multicellular organisms, but have occurred in parallel in prokaryotes with predatory lifestyles and phagocytotic-like behaviors.

Folded Alpha Helical Putative New Proteins from Apilactobacillus kunkeei.

The emergence of new proteins is a central question in biology. Most tertiary protein folds known to date appear to have an ancient origin, but it is clear from bioinformatic analyses that new proteins continuously emerge in all organismal groups. However, there is a paucity of experimental data on new proteins regarding their structure and biophysical properties. We performed a detailed phylogenetic analysis and identified 48 putative open reading frames in the honeybee-associated bacterium Apilactobacillus kunkeei for which no or few homologs could be identified in closely-related species, suggesting that they could be relatively new on an evolutionary time scale and represent recently evolved proteins. Using circular dichroism-, fluorescence- and nuclear magnetic resonance (NMR) spectroscopy we investigated six of these proteins and show that they are not intrinsically disordered, but populate alpha-helical dominated folded states with relatively low thermodynamic stability (0-3 kcal/mol). The NMR and biophysical data demonstrate that small new proteins readily adopt simple folded conformations suggesting that more complex tertiary structures can be continuously re-invented during evolution by fusion of such simple secondary structure elements. These findings have implications for the general view on protein evolution, where de novo emergence of folded proteins may be a common event.

TADA: taxonomy-aware dataset aggregator.

SUMMARY: The profusion of sequenced genomes across the bacterial and archeal domains offers unprecedented possibilities for phylogenetic and comparative genomic analyses. In general, phylogenetic reconstruction is improved by the use of more data. However, including all available data is (i) not computationally tractable, and (ii) prone to biases, as the abundance of genomes is very unequally distributed over the biological diversity. Thus, in most cases, subsampling taxa to build a phylogeny is necessary. Currently, though, there is no available software to perform that handily. Here we present TADA, a taxonomic-aware dataset selection workflow that allows sampling across user-defined portions of the prokaryotic diversity with variable granularity, while setting constraints on genome quality and balance between branches. AVAILABILITY AND IMPLEMENTATION: TADA is implemented as a snakemake workflow and is freely available at https://github.com/emilhaegglund/TADA.

Apilactobacillus kunkeei releases RNA-associated membrane vesicles and proteinaceous nanoparticles.

Extracellularly released particles, including membrane vesicles, have increasingly been recognized as important for bacterial community functions and host-interaction processes, but their compositions and functional roles differ between species and also between strains of the same species. In this study, we have determined the composition of membrane vesicles and protein particles identified in the cell-free pellets of two strains of Apilactobacillus kunkeei, a defensive symbiont of honeybees. The membrane vesicles were separated from the extracellular particles using density gradient ultracentrifugation. The peaks of the RNA and protein distributions were separated from each other and the highest concentration of RNA was observed in the fractions that contained the membrane vesicles while the highest protein concentration coincided with the fractions that contained extracellular particles. A comparative proteomics analysis by LC-MS/MS showed that 37 proteins with type-I signal peptides were consistently identified across the fractionated samples obtained from the cell-free pellets, of which 29 were orthologs detected in both strains. Functional predictions of the extracellular proteins revealed the presence of glycoside hydrolases, glycosyltransferases, giant proteins and peptidases. The extracellular transcriptomes mapped to a broad set of genes with a similar functional profile as the whole cell transcriptome. This study provides insights into the composition of membrane vesicles and extracellular proteins of a bee-associated symbiont.

Genome Evolution of a Symbiont Population for Pathogen Defense in Honeybees.

The honeybee gut microbiome is thought to be important for bee health, but the role of the individual members is poorly understood. Here, we present closed genomes and associated mobilomes of 102 Apilactobacillus kunkeei isolates obtained from the honey crop (foregut) of honeybees sampled from beehives in Helsingborg in the south of Sweden and from the islands Gotland and Åland in the Baltic Sea. Each beehive contained a unique composition of isolates and repeated sampling of similar isolates from two beehives in Helsingborg suggests that the bacterial community is stably maintained across bee generations during the summer months. The sampled bacterial population contained an open pan-genome structure with a high genomic density of transposons. A subset of strains affiliated with phylogroup A inhibited growth of the bee pathogen Melissococcus plutonius, all of which contained a 19.5 kb plasmid for the synthesis of the antimicrobial compound kunkecin A, while a subset of phylogroups B and C strains contained a 32.9 kb plasmid for the synthesis of a putative polyketide antibiotic. This study suggests that the mobile gene pool of A. kunkeei plays a key role in pathogen defense in honeybees, providing new insights into the evolutionary dynamics of defensive symbiont populations.

The Subcellular Proteome of a Planctomycetes Bacterium Shows That Newly Evolved Proteins Have Distinct Fractionation Patterns.

The Planctomycetes bacteria have unique cell architectures with heavily invaginated membranes as confirmed by three-dimensional models reconstructed from FIB-SEM images of Tuwongella immobilis and Gemmata obscuriglobus. The subcellular proteome of T. immobilis was examined by differential solubilization followed by LC-MS/MS analysis, which identified 1569 proteins in total. The Tris-soluble fraction contained mostly cytoplasmic proteins, while inner and outer membrane proteins were found in the Triton X-100 and SDS-soluble fractions, respectively. For comparisons, the subcellular proteome of Escherichia coli was also examined using the same methodology. A notable difference in the overall fractionation pattern of the two species was a fivefold higher number of predicted cytoplasmic proteins in the SDS-soluble fraction in T. immobilis. One category of such proteins is represented by innovations in the Planctomycetes lineage, including unique sets of serine/threonine kinases and extracytoplasmic sigma factors with WD40 repeat domains for which no homologs are present in E. coli. Other such proteins are members of recently expanded protein families in which the newly evolved paralog with a new domain structure is recovered from the SDS-soluble fraction, while other paralogs may have similar domain structures and fractionation patterns as the single homolog in E. coli. The expanded protein families in T. immobilis include enzymes involved in replication-repair processes as well as in rRNA and tRNA modification and degradation. These results show that paralogization and domain shuffling have yielded new proteins with distinct fractionation characteristics. Understanding the molecular intricacies of these adaptive changes might aid in the development of a model for the evolution of cellular complexity.

Evolutionary Remodeling of the Cell Envelope in Bacteria of the Planctomycetes Phylum.

Bacteria of the Planctomycetes phylum have many unique cellular features, such as extensive membrane invaginations and the ability to import macromolecules. These features raise intriguing questions about the composition of their cell envelopes. In this study, we have used microscopy, phylogenomics, and proteomics to examine the composition and evolution of cell envelope proteins in Tuwongella immobilis and other members of the Planctomycetes. Cryo-electron tomography data indicated a distance of 45 nm between the inner and outer membranes in T. immobilis. Consistent with the wide periplasmic space, our bioinformatics studies showed that the periplasmic segments of outer-membrane proteins in type II secretion systems are extended in bacteria of the order Planctomycetales. Homologs of two highly abundant cysteine-rich cell wall proteins in T. immobilis were identified in all members of the Planctomycetales, whereas genes for peptidoglycan biosynthesis and cell elongation have been lost in many members of this bacterial group. The cell wall proteins contain multiple copies of the YTV motif, which is the only domain that is conserved and unique to the Planctomycetales. Earlier diverging taxa in the Planctomycetes phylum contain genes for peptidoglycan biosynthesis but no homologs to the YTV cell wall proteins. The major remodeling of the cell envelope in the ancestor of the Planctomycetales coincided with the emergence of budding and other unique cellular phenotypes. The results have implications for hypotheses about the process whereby complex cellular features evolve in bacteria.

Correction to: Bacteria: back pain, leg pain and Modic sign-a surgical multicenter comparative study.

Unfortunately, the 5th author name was incorrectly published in the original paper. The complete correct name is given below.

Paralogization and New Protein Architectures in Planctomycetes Bacteria with Complex Cell Structures.

Bacteria of the phylum Planctomycetes have a unique cell plan with an elaborate intracellular membrane system, thereby resembling eukaryotic cells. The origin and evolution of these remarkable features is debated. To study the evolutionary genomics of bacteria with complex cell architectures, we have resequenced the 9.2-Mb genome of the model organism Gemmata obscuriglobus and sequenced the 10-Mb genome of G. massiliana Soil9, the 7.9-Mb genome of CJuql4, and the 6.7-Mb genome of Tuwongella immobilis, all of which belong to the family Gemmataceae. A gene flux analysis of the Planctomycetes revealed a massive emergence of novel protein families at multiple nodes within the Gemmataceae. The expanded protein families have unique multidomain architectures composed of domains that are characteristic of prokaryotes, such as the sigma factor domain of extracytoplasmic sigma factors, and domains that have proliferated in eukaryotes, such as the WD40, leucine-rich repeat, tetratricopeptide repeat and Ser/Thr kinase domains. Proteins with identifiable domains in the Gemmataceae have longer lengths and linkers than proteins in most other bacteria, and the analyses suggest that these traits were ancestrally present in the Planctomycetales. A broad comparison of protein length distribution profiles revealed an overlap between the longest proteins in prokaryotes and the shortest proteins in eukaryotes. We conclude that the many similarities between proteins in the Planctomycetales and the eukaryotes are due to convergent evolution and that there is no strict boundary between prokaryotes and eukaryotes with regard to features such as gene paralogy, protein length, and protein domain composition patterns.

Bacteria: back pain, leg pain and Modic sign-a surgical multicentre comparative study.

PURPOSE: To compare bacterial findings in pain-generating degenerated discs in adults operated on for lumbar disc herniation (LDH), and mostly also suffering from low back pain (LBP), with findings in adolescent patients with non-degenerated non-pain-generating discs operated on for scoliosis, and to evaluate associations with Modic signs on magnetic resonance imaging (MRI). Cutibacterium acnes (Propionibacterium acnes) has been found in painful degenerated discs, why it has been suggested treating patients with LDH/LBP with antibiotics. As multidrug-resistant bacteria are a worldwide concern, new indications for using antibiotics should be based on solid scientific evidence. METHODS: Between 2015 and 2017, 40 adults with LDH/LBP (median age 43, IQR 33-49) and 20 control patients with scoliosis (median age 17, IQR 15-20) underwent surgery at seven Swedish hospitals. Samples were cultured from skin, surgical wound, discs and vertebrae. Genetic relatedness of C. acnes isolates was investigated using single-nucleotide polymorphism analysis. DNA samples collected from discs/vertebrae were analysed using 16S rRNA-based PCR sequencing. MRI findings were assessed for Modic changes. RESULTS: No bacterial growth was found in 6/40 (15%) LDH patients, compared with 3/20 (15%) scoliosis patients. Most positive samples in both groups were isolated from the skin and then from subcutis or deep within the wound. Of the four disc and vertebral samples from each of the 60 patients, 235/240 (98%) were DNA negative by bacterial PCR. A single species, C. acnes, was found exclusively in the disc/vertebra from one patient in each group. In the LDH group, 29/40 (72%) patients had at least one sample with growth of C. acnes, compared to 14/20 (70%) in the scoliosis group. Bacterial findings and Modic changes were not associated. CONCLUSIONS: Cutibacterium acnes found in discs and vertebrae during surgery for disc herniation in adults with degenerated discs may be caused by contamination, as findings in this group were similar to findings in a control group of young patients with scoliosis and non-degenerated discs. Furthermore, such findings were almost always combined with bacterial findings on the skin and/or in the wound. There was no association between preoperative Modic changes and bacterial findings. Antibiotic treatment of lumbar disc herniation with sciatica and/or low back pain, without signs of clinical discitis/spondylitis, should be seriously questioned. These slides can be retrieved under Electronic Supplementary Material.

Genome Evolution of Bartonellaceae Symbionts of Ants at the Opposite Ends of the Trophic Scale.

Many insects rely on bacterial symbionts to supply essential amino acids and vitamins that are deficient in their diets, but metabolic comparisons of closely related gut bacteria in insects with different dietary preferences have not been performed. Here, we demonstrate that herbivorous ants of the genus Dolichoderus from the Peruvian Amazon host bacteria of the family Bartonellaceae, known for establishing chronic or pathogenic infections in mammals. We detected these bacteria in all studied Dolichoderus species, and found that they reside in the midgut wall, that is, the same location as many previously described nutritional endosymbionts of insects. The genomic analysis of four divergent strains infecting different Dolichoderus species revealed genes encoding pathways for nitrogen recycling and biosynthesis of several vitamins and all essential amino acids. In contrast, several biosynthetic pathways have been lost, whereas genes for the import and conversion of histidine and arginine to glutamine have been retained in the genome of a closely related gut bacterium of the carnivorous ant Harpegnathos saltator. The broad biosynthetic repertoire in Bartonellaceae of herbivorous ants resembled that of gut bacteria of honeybees that likewise feed on carbohydrate-rich diets. Taken together, the broad distribution of Bartonellaceae across Dolichoderus ants, their small genome sizes, the specific location within hosts, and the broad biosynthetic capability suggest that these bacteria are nutritional symbionts in herbivorous ants. The results highlight the important role of the host nutritional biology for the genomic evolution of the gut microbiota-and conversely, the importance of the microbiota for the nutrition of hosts.

Specificity in Arabidopsis thaliana recruitment of root fungal communities from soil and rhizosphere.

Biotic and abiotic conditions in soil pose major constraints on growth and reproductive success of plants. Fungi are important agents in plant soil interactions but the belowground mycobiota associated with plants remains poorly understood. We grew one genotype each from Sweden and Italy of the widely-studied plant model Arabidopsis thaliana. Plants were grown under controlled conditions in organic topsoil local to the Swedish genotype, and harvested after ten weeks. Total DNA was extracted from three belowground compartments: endosphere (sonicated roots), rhizosphere and bulk soil, and fungal communities were characterized from each by amplification and sequencing of the fungal barcode region ITS2. Fungal species diversity was found to decrease from bulk soil to rhizosphere to endosphere. A significant effect of plant genotype on fungal community composition was detected only in the endosphere compartment. Despite A. thaliana being a non-mycorrhizal plant, it hosts a number of known mycorrhiza fungi in its endosphere compartment, which is also colonized by endophytic, pathogenic and saprotrophic fungi. Species in the Archaeorhizomycetes were most abundant in rhizosphere samples suggesting an adaptation to environments with high nutrient turnover for some of these species. We conclude that A. thaliana endosphere fungal communities represent a selected subset of fungi recruited from soil and that plant genotype has small but significant quantitative and qualitative effects on these communities.

Origin and Evolution of the Bartonella Gene Transfer Agent.

Gene transfer agents (GTAs) are domesticated bacteriophages that have evolved into molecular machines for the transfer of bacterial DNA. Despite their widespread nature and their biological implications, the mechanisms and selective forces that drive the emergence of GTAs are still poorly understood. Two GTAs have been identified in the Alphaproteobacteria: the RcGTA, which is widely distributed in a broad range of species; and the BaGTA, which has a restricted host range that includes vector-borne intracellular bacteria of the genus Bartonella. The RcGTA packages chromosomal DNA randomly, whereas the BaGTA particles contain a relatively higher fraction of genes for host interaction factors that are amplified from a nearby phage-derived origin of replication. In this study, we compare the BaGTA genes with homologous bacteriophage genes identified in the genomes of Bartonella species and close relatives. Unlike the BaGTA, the prophage genes are neither present in all species, nor inserted into homologous genomic sites. Phylogenetic inferences and substitution frequency analyses confirm codivergence of the BaGTA with the host genome, as opposed to multiple integration and recombination events in the prophages. Furthermore, the organization of segments flanking the BaGTA differs from that of the prophages by a few rearrangement events, which have abolished the normal coordination between phage genome replication and phage gene expression. Based on the results of our comparative analysis, we propose a model for how a prophage may be transformed into a GTA that transfers amplified bacterial DNA segments.

Contrasting patterns of genome-level diversity across distinct co-occurring bacterial populations.

To understand the forces driving differentiation and diversification in wild bacterial populations, we must be able to delineate and track ecologically relevant units through space and time. Mapping metagenomic sequences to reference genomes derived from the same environment can reveal genetic heterogeneity within populations, and in some cases, be used to identify boundaries between genetically similar, but ecologically distinct, populations. Here we examine population-level heterogeneity within abundant and ubiquitous freshwater bacterial groups such as the acI Actinobacteria and LD12 Alphaproteobacteria (the freshwater sister clade to the marine SAR11) using 33 single-cell genomes and a 5-year metagenomic time series. The single-cell genomes grouped into 15 monophyletic clusters (termed "tribes") that share at least 97.9% 16S rRNA identity. Distinct populations were identified within most tribes based on the patterns of metagenomic read recruitments to single-cell genomes representing these tribes. Genetically distinct populations within tribes of the acI Actinobacterial lineage living in the same lake had different seasonal abundance patterns, suggesting these populations were also ecologically distinct. In contrast, sympatric LD12 populations were less genetically differentiated. This suggests that within one lake, some freshwater lineages harbor genetically discrete (but still closely related) and ecologically distinct populations, while other lineages are composed of less differentiated populations with overlapping niches. Our results point at an interplay of evolutionary and ecological forces acting on these communities that can be observed in real time.

Tuwongella immobilis gen. nov., sp. nov., a novel non-motile bacterium within the phylum Planctomycetes.

A gram-negative, budding, catalase negative, oxidase positive and non-motile bacterium (MBLW1T) with a complex endomembrane system has been isolated from a freshwater lake in southeast Queensland, Australia. Phylogeny based on 16S rRNA gene sequence analysis places the strain within the family Planctomycetaceae, related to Zavarzinella formosa (93.3 %), Telmatocola sphagniphila (93.3 %) and Gemmata obscuriglobus (91.9 %). Phenotypic and chemotaxonomic analysis demonstrates considerable differences to the type strains of the related genera. MBLW1T displays modest salt tolerance and grows optimally at pH values of 7.5-8.0 and at temperatures of 32-36 °C. Transmission electron microscopy analysis demonstrates the presence of a complex endomembrane system, however, without the typically condensed nucleoid structure found in related genera. The major fatty acids are 16 : 1 ω5c, 16 : 0 and 18 : 0. Based on discriminatory results from 16S rRNA gene sequence analysis, phenotypic, biochemical and chemotaxonomic analysis, MBLW1T should be considered as a new genus and species, for which the name Tuwongella immobilis gen. nov., sp. nov. is proposed. The type strain is MBLW1T (=CCUG 69661T=DSM 105045T).

Switches in Genomic GC Content Drive Shifts of Optimal Codons under Sustained Selection on Synonymous Sites.

The major codon preference model suggests that codons read by tRNAs in high concentrations are preferentially utilized in highly expressed genes. However, the identity of the optimal codons differs between species although the forces driving such changes are poorly understood. We suggest that these questions can be tackled by placing codon usage studies in a phylogenetic framework and that bacterial genomes with extreme nucleotide composition biases provide informative model systems. Switches in the background substitution biases from GC to AT have occurred in Gardnerella vaginalis (GC = 32%), and from AT to GC in Lactobacillus delbrueckii (GC = 62%) and Lactobacillus fermentum (GC = 63%). We show that despite the large effects on codon usage patterns by these switches, all three species evolve under selection on synonymous sites. In G. vaginalis, the dramatic codon frequency changes coincide with shifts of optimal codons. In contrast, the optimal codons have not shifted in the two Lactobacillus genomes despite an increased fraction of GC-ending codons. We suggest that all three species are in different phases of an on-going shift of optimal codons, and attribute the difference to a stronger background substitution bias and/or longer time since the switch in G. vaginalis. We show that comparative and correlative methods for optimal codon identification yield conflicting results for genomes in flux and discuss possible reasons for the mispredictions. We conclude that switches in the direction of the background substitution biases can drive major shifts in codon preference patterns even under sustained selection on synonymous codon sites.