Host-bacteria interactions: ecological and evolutionary insights from ancient, professional endosymbionts.

Interactions between eukaryotic hosts and their bacterial symbionts drive key ecological and evolutionary processes, from regulating ecosystems to the evolution of complex molecular machines and processes. Over time, endosymbionts generally evolve reduced genomes, and their relationship with their host tends to stabilize. However, host-bacteria relationships may be heavily influenced by environmental changes. Here, we review these effects on one of the most ancient and diverse endosymbiotic groups, formed by-among others-Legionellales, Francisellaceae, and Piscirickettsiaceae. This group is referred to as Deep-branching Intracellular Gammaproteobacteria (DIG), whose last common ancestor presumably emerged about 2 Ga ago. We show that DIGs are globally distributed, but generally at very low abundance, and are mainly identified in aquatic biomes. Most DIGs harbour a type IVB secretion system, critical for host-adaptation, but its structure and composition vary. Finally, we review the different types of microbial interactions that can occur in diverse environments, with direct or indirect effects on DIG populations. The increased use of omics technologies on environmental samples will allow a better understanding of host-bacterial interactions and help unravel the definition of DIGs as a group from an ecological, molecular, and evolutionary perspective.

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.

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.

The evolution of intramitochondriality in Midichloria bacteria.

Midichloria spp. are intracellular bacterial symbionts of ticks. Representatives of this genus colonise mitochondria in the cells of their hosts. To shed light on this unique interaction we evaluated the presence of an intramitochondrial localization for three Midichloria in the respective tick host species and generated eight high-quality draft genomes and one closed genome, showing that this trait is non-monophyletic, either due to losses or multiple acquisitions. Comparative genomics supports the first hypothesis, as the genomes of non-mitochondrial symbionts are reduced subsets of those capable of colonising the organelles. We detect genomic signatures of mitochondrial tropism, including the differential presence of type IV secretion system and flagellum, which could allow the secretion of unique effectors and/or direct interaction with mitochondria. Other genes, including adhesion molecules, proteins involved in actin polymerisation, cell wall and outer membrane proteins, are only present in mitochondrial symbionts. The bacteria could use these to manipulate host structures, including mitochondrial membranes, to fuse with the organelles or manipulate the mitochondrial network.

Rapid adaptations of Legionella pneumophila to the human host.

Legionella pneumophila are host-adapted bacteria that infect and reproduce primarily in amoeboid protists. Using similar infection mechanisms, they infect human macrophages, and cause Legionnaires' disease, an atypical pneumonia, and the milder Pontiac fever. We hypothesized that, despite the similarities in infection mechanisms, the hosts are different enough that there exist high-selective value mutations that would dramatically increase the fitness of Legionella inside the human host. By comparing a large number of isolates from independent infections, we identified two genes, mutated in three unrelated patients, despite the short duration of the incubation period (2-14 days). One is a gene coding for an outer membrane protein (OMP) belonging to the OmpP1/FadL family. The other is a gene coding for an EAL-domain-containing protein involved in cyclic-di-GMP regulation, which in turn modulates flagellar activity. The clinical strain, carrying the mutated EAL-domain-containing homologue, grows faster in macrophages than the wild-type strain, and thus appears to be better adapted to the human host. As human-to-human transmission is very rare, fixation of these mutations into the population and spread into the environment is unlikely. Therefore, parallel evolution - here mutations in the same genes observed in independent human infections - could point to adaptations to the accidental human host. These results suggest that despite the ability of L. pneumophila to infect, replicate in and exit from macrophages, its human-specific adaptations are unlikely to be fixed in the population.

Early divergence of translation initiation and elongation factors.

Protein translation is a foundational attribute of all living cells. The translation function carried out by the ribosome critically depends on an assortment of protein interaction partners, collectively referred to as the translation machinery. Various studies suggest that the diversification of the translation machinery occurred prior to the last universal common ancestor, yet it is unclear whether the predecessors of the extant translation machinery factors were functionally distinct from their modern counterparts. Here we reconstructed the shared ancestral trajectory and subsequent evolution of essential translation factor GTPases, elongation factor EF-Tu (aEF-1A/eEF-1A), and initiation factor IF2 (aIF5B/eIF5B). Based upon their similar functions and structural homologies, it has been proposed that EF-Tu and IF2 emerged from an ancient common ancestor. We generated the phylogenetic tree of IF2 and EF-Tu proteins and reconstructed ancestral sequences corresponding to the deepest nodes in their shared evolutionary history, including the last common IF2 and EF-Tu ancestor. By identifying the residue and domain substitutions, as well as structural changes along the phylogenetic history, we developed an evolutionary scenario for the origins, divergence and functional refinement of EF-Tu and IF2 proteins. Our analyses suggest that the common ancestor of IF2 and EF-Tu was an IF2-like GTPase protein. Given the central importance of the translation machinery to all cellular life, its earliest evolutionary constraints and trajectories are key to characterizing the universal constraints and capabilities of cellular evolution.

Host Adaptation in Legionellales Is 1.9 Ga, Coincident with Eukaryogenesis.

Bacteria adapting to living in a host cell caused the most salient events in the evolution of eukaryotes, namely the seminal fusion with an archaeon, and the emergence of both mitochondrion and chloroplast. A bacterial clade that may hold the key to understanding these events is the deep-branching gammaproteobacterial order Legionellales-containing among others Coxiella and Legionella-of which all known members grow inside eukaryotic cells. Here, by analyzing 35 novel Legionellales genomes mainly acquired through metagenomics, we show that this group is much more diverse than previously thought, and that key host-adaptation events took place very early in its evolution. Crucial virulence factors like the Type IVB secretion (Dot/Icm) system and two shared effector proteins were gained in the last Legionellales common ancestor (LLCA). Many metabolic gene families were lost in LLCA and its immediate descendants, including functions directly and indirectly related to molybdenum metabolism. On the other hand, genome sizes increased in the ancestors of the Legionella genus. We estimate that LLCA lived approximately 1.89 Ga, probably predating the last eukaryotic common ancestor by approximately 0.4-1.0 Gy. These elements strongly indicate that host adaptation arose only once in Legionellales, and that these bacteria were using advanced molecular machinery to exploit and manipulate host cells early in eukaryogenesis.

Heimdallarchaea encodes profilin with eukaryotic-like actin regulation and polyproline binding.

It is now widely accepted that the first eukaryotic cell emerged from a merger of an archaeal host cell and an alphaproteobacterium. However, the exact sequence of events and the nature of the cellular biology of both partner cells is still contentious. Recently the structures of profilins from some members of the newly discovered Asgard superphylum were determined. In addition, it was found that these profilins inhibit eukaryotic rabbit actin polymerization and that this reaction is regulated by phospholipids. However, the interaction with polyproline repeats which are known to be crucial for the regulation of profilin:actin polymerization was found to be absent for these profilins and was thus suggested to have evolved later in the eukaryotic lineage. Here, we show that Heimdallarchaeota LC3, a candidate phylum within the Asgard superphylum, encodes a putative profilin (heimProfilin) that interacts with PIP2 and its binding is regulated by polyproline motifs, suggesting an origin predating the rise of the eukaryotes. More precisely, we determined the 3D-structure of Heimdallarchaeota LC3 profilin and show that this profilin is able to: i) inhibit eukaryotic actin polymerization in vitro; ii) bind to phospholipids; iii) bind to polyproline repeats from enabled/vasodilator-stimulated phosphoprotein; iv) inhibit actin from Heimdallarchaeota from polymerizing into filaments. Our results therefore provide hints of the existence of a complex cytoskeleton already in last eukaryotic common ancestor.

Chromosomal Location Determines the Rate of Intrachromosomal Homologous Recombination in Salmonella.

Homologous recombination is an important mechanism directly involved in the repair, organization, and evolution of prokaryotic and eukaryotic chromosomes. We developed a system, based on two genetic cassettes, that allows the measurement of recombinational repair rates between different locations on the chromosome. Using this system, we analyzed 81 different positional combinations throughout the chromosome to answer the question of how the position and orientation of sequences affect intrachromosomal homologous recombination. Our results show that recombination was possible between any two locations tested in this study and that recombinational repair rates varied by just above an order of magnitude. The observed differences in rate do not correlate with distance between the recombination cassettes or with distance from the origin of replication but could be explained if each location contributes individually to the recombination event. The relative levels of accessibility for recombination vary 5-fold between the various cassette locations, and we found that the nucleoid structure of the chromosome may be the major factor influencing the recombinational accessibility of each chromosomal site. Furthermore, we found that the orientation of the recombination cassettes had a significant impact on recombination. Recombinational repair rates for the cassettes inserted as direct repeats are, on average, 2.2-fold higher than those for the same sets inserted as inverted repeats. These results suggest that the bacterial chromosome is not homogenous with regard to homologous recombination, with regions that are more or less accessible, and that the orientation of genes affects recombination rates. IMPORTANCE Bacterial chromosomes frequently carry multiple copies of genes at separate chromosomal locations. In Salmonella, these include the 7 rrn operons and the duplicate tuf genes. Genes within these families coevolve by homologous recombination, but it is not obvious whether their rates of recombination reflect general rates of intrachromosomal recombination or are an evolved property particularly associated with these conserved genes and locations. Using a novel experimental system, we show that recombination is possible between all tested pairs of locations at rates that vary by just above 1 order of magnitude. Differences in rate do not correlate with distance between the sites or distance to the origin of replication but may be explained if each location contributes individually to the recombination event. Our results suggest the existence of bacterial chromosomal domains that are differentially available for recombination and that gene orientation affects recombination rates.

Molecular Characterization of Multidrug-Resistant Yersinia enterocolitica From Foodborne Outbreaks in Sweden.

The foodborne pathogen Yersinia enterocolitica causes gastrointestinal infections worldwide. In the spring of 2019, the Swedish Public Health Agency and Statens Serum Institut in Denmark independently identified an outbreak caused by Yersinia enterocolitica 4/O:3 that after sequence comparison turned out to be a cross-border outbreak. A trace-back investigation suggested shipments of fresh prewashed spinach from Italy as a common source for the outbreak. Here, we determined the genome sequences of five Y. enterocolitica clinical isolates during the Swedish outbreak using a combination of Illumina HiSeq short-read and Nanopore Technologies' MinION long-read whole-genome sequencing. WGS results showed that all clinical strains have a fully assembled chromosome of approximately 4.6 Mbp in size and a 72-kbp virulence plasmid; one of the strains was carrying an additional 5.7-kbp plasmid, pYE-tet. All strains showed a high pathogen probability score (87.5%) with associated genes for virulence, all of which are closely related to an earlier clinical strain Y11 from Germany. In addition, we identified a chromosomally encoded multidrug-resistance cassette carrying resistance genes against chloramphenicol (catA1), streptomycin (aadA1), sulfonamides (sul1), and a mercury resistance module. This chromosomally encoded Tn2670 transposon has previously been reported associated with IncFII plasmids in Enterobacteriaceae: a Shigella flexneri clinical isolate from Japan in 1950s, a Klebsiella pneumoniae outbreak from Australia in 1997, and Salmonella enterica serovar Typhimurium. Interestingly, we identified an additional 5.7-kbp plasmid with tetB (encoding an ABC transporter), Rep, and its own ORI and ORIt sites, sharing high homology with small tetB-Rep plasmids from Pasteurellaceae. This is the first time that Tn2670 and Pasteurellaceae plasmids have been reported in Y. enterocolitica. Taken together, our study showed that the Swedish Y. enterocolitica outbreak strains acquired multi-antibiotic and metal-resistance genes through horizontal gene transfer, suggesting a potential reservoir of intraspecies dissemination of multidrug-resistance genes among foodborne pathogens. This study also highlights the concern of food-chain contamination of prewashed vegetables as a perpetual hazard against public health.

Expression of the qepA1 gene is induced under antibiotic exposure.

BACKGROUND: The qepA1 gene encodes an efflux pump that reduces susceptibility to ciprofloxacin. Little is known about the regulation of qepA1 expression. OBJECTIVES: To assess the potential role of ciprofloxacin and other antibiotics in the regulation of qepA1 gene expression. To identify the promoter that drives qepA1 expression and other factors involved in expression regulation. To assess whether the identified features are universal among qepA alleles. METHODS: A translational qepA1-yfp fusion under the control of the qepA1 upstream region was cloned into the Escherichia coli chromosome. Expression of the fusion protein was measured in the presence of various antibiotics. Deletions within the upstream region were introduced to identify regions involved in gene expression and regulation. The qepA1 coding sequence and upstream region were compared with all available qepA sequences. RESULTS: Cellular stress caused by the presence of various antibiotics can induce qepA1 expression. The qepA1 gene is fused to a class I integron and gene expression is driven by the Pc promoter within the integrase gene. A segment within the integron belonging to a truncated dfrB4 gene is essential for the regulation of qepA1 expression. This genetic context is universal among all sequenced qepA alleles. CONCLUSIONS: The fusion of the qepA1 gene to a class I integron has created a novel regulatory unit that enables qepA1 expression to be under the control of antibiotic exposure. This setup mitigates potential negative effects of QepA1 production on bacterial fitness by restricting high-level expression to environmental conditions in which QepA1 is beneficial.

Chlamydia psittaci in fulmars on the Faroe Islands: a causative link to South American psittacines eight decades after a severe epidemic.

A psittacosis epidemic linked to fulmar hunting occurred on the Faroe Islands in the 1930s. This study investigates a plausible explanation to the 20% human mortality in this outbreak. Phylogenetic analysis showed that Chlamydia psittaci isolated from fulmars were closely related to the highly virulent 6BC strains from psittacines and are compatible with an acquisition by fulmars of an ancestor of the 6BC clade in the 1930s. This supports the hypothesis that the outbreak on the Faroe Islands started after naïve fulmars acquired C. psittaci from infected dead parrots thrown overboard when shipped to Europe in the 1930s.

miComplete: weighted quality evaluation of assembled microbial genomes.

SUMMARY: Metagenomics and single-cell genomics have revolutionized the study of microorganisms, increasing our knowledge of microbial genomic diversity by orders of magnitude. A major issue pertaining to metagenome-assembled genomes (MAGs) and single-cell amplified genomes (SAGs) is to estimate their completeness and redundancy. Most approaches rely on counting conserved gene markers. In miComplete, we introduce a weighting strategy, where we normalize the presence/absence of markers by their median distance to the next marker in a set of complete reference genomes. This approach alleviates biases introduced by the presence/absence of shorter DNA pieces containing many markers, e.g. ribosomal protein operons. AVAILABILITY AND IMPLEMENTATION: miComplete is written in Python 3 and released under GPLv3. Source code and documentation are available at https://bitbucket.org/evolegiolab/micomplete. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

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.

Withdrawn as Duplicate: The all-intracellular order Legionellales is unexpectedly diverse, globally distributed and lowly abundant.

Legionellales is an order of the Gammaproteobacteria, only composed of host-adapted, intracellular bacteria, including the accidental human pathogens Legionella pneumophila and Coxiella burnetii. Although the diversity in terms of lifestyle is large across the order, only a few genera have been sequenced, owing to the difficulty to grow intracellular bacteria in pure culture. In particular, we know little about their global distribution and abundance.Here, we analyze 16/18S rDNA amplicons both from tens of thousands of published studies and from two separate sampling campaigns in and around ponds and in a silver mine. We demonstrate that the diversity of the order is much larger than previously thought, with over 450 uncultured genera. We show that Legionellales are found in about half of the samples from freshwater, soil and marine environments, and quasi-ubiquitous in man-made environments. Their abundance is low, typically 0.1%, with few samples up to 1%. Most Legionellales OTUs are globally distributed, while many do not belong to a previously identified species.This study sheds a new light on the ubiquity and diversity of one major group of host-adapted bacteria. It also emphasizes the need to use metagenomics to better understand the role of host-adapted bacteria in all environments.

Genomes of two archaeal endosymbionts show convergent adaptations to an intracellular lifestyle.

Endosymbiosis is a widespread phenomenon in the microbial world and can be based on diverse interactions between endosymbiont and host cell. The vast majority of the known endosymbiotic interactions involve bacteria that have invaded eukaryotic host cells. However, methanogenic archaea have been found to thrive in anaerobic, hydrogenosome-containing protists and it was suggested that this symbiosis is based on the transfer of hydrogen. Here, we used culture-independent genomics approaches to sequence the genomes of two distantly related methanogenic endosymbionts that have been acquired in two independent events by closely related anaerobic ciliate hosts Nyctotherus ovalis and Metopus contortus, respectively. The sequences obtained were then validated as originating from the ciliate endosymbionts by in situ probing experiments. Comparative analyses of these genomes and their closest free-living counterparts reveal that the genomes of both endosymbionts are in an early stage of adaptation towards endosymbiosis as evidenced by the large number of genes undergoing pseudogenization. For instance, the observed loss of genes involved in amino acid biosynthesis in both endosymbiont genomes indicates that the endosymbionts rely on their hosts for obtaining several essential nutrients. Furthermore, the endosymbionts appear to have gained significant amounts of genes of potentially secreted proteins, providing targets for future studies aiming to elucidate possible mechanisms underpinning host-interactions. Altogether, our results provide the first genomic insights into prokaryotic endosymbioses from the archaeal domain of life.

Legionella pneumophila recurrently isolated in a Spanish hospital: Two years of antimicrobial resistance surveillance.

OBJECTIVES: The aim of this study was to monitor the spread, persistence and antibiotic resistance patterns of Legionella spp. strains found in a hospital water distribution system. These environmental studies are intended to help detect the presence of antibiotic resistant strains before they infect patients. METHODS: Antimicrobial surveillance tests were performed at 27 different sampling points of the water network of a large Spanish hospital over two years. Water samples were screened for Legionella according to ISO 11731:2007. Legionella spp. isolates were identified by serotyping and by mass spectrometry (MALDI-ToF). Epidemiological molecular typing was done by Pulse-Field Gel Electrophoresis (PFGE) and by Sequence-Based Typing (SBT). Antibiotic susceptibility tests were performed using disk diffusion and ETEST®. RESULTS: Legionella spp. were recurrently isolated for 2 years. All isolates belonged the same group, L. pneumophila serogroups 2-14. Isolates were all attributed by SBT to sequence type (ST) ST328, although PFGE revealed 5 different patterns. No significant change in antibiotic susceptibility could be observed for this study period, irrespectively of the method used. CONCLUSION: Colonization of water systems by Legionella spp. is still occurring, although all the prevention rules were strictly followed. Antibiotic resistance monitoring may help us to find resistance in bacteria with environmental reservoirs but difficult to isolate from patients. The knowledge of the antibiotic susceptibility in environmental strains may help us to predict changes in clinical strains. This study might also help reconsidering Legionnaires' disease (LD) diagnostic methods. L. pneumophila serogroups 2-14 present all along the time of the investigation in the water distribution system can cause LD. However, they may not be detected by routine urine tests run on patients, thereby missing an ongoing LD infection.

Deep mitochondrial origin outside the sampled alphaproteobacteria.

Mitochondria are ATP-generating organelles, the endosymbiotic origin of which was a key event in the evolution of eukaryotic cells 1 . Despite strong phylogenetic evidence that mitochondria had an alphaproteobacterial ancestry 2 , efforts to pinpoint their closest relatives among sampled alphaproteobacteria have generated conflicting results, complicating detailed inferences about the identity and nature of the mitochondrial ancestor. While most studies support the idea that mitochondria evolved from an ancestor related to Rickettsiales3-9, an order that includes several host-associated pathogenic and endosymbiotic lineages10,11, others have suggested that mitochondria evolved from a free-living group12-14. Here we re-evaluate the phylogenetic placement of mitochondria. We used genome-resolved binning of oceanic metagenome datasets and increased the genomic sampling of Alphaproteobacteria with twelve divergent clades, and one clade representing a sister group to all Alphaproteobacteria. Subsequent phylogenomic analyses that specifically address long branch attraction and compositional bias artefacts suggest that mitochondria did not evolve from Rickettsiales or any other currently recognized alphaproteobacterial lineage. Rather, our analyses indicate that mitochondria evolved from a proteobacterial lineage that branched off before the divergence of all sampled alphaproteobacteria. In light of this new result, previous hypotheses on the nature of the mitochondrial ancestor6,15,16 should be re-evaluated.