Genomic changes following host restriction in bacteria

Pseudogenes in plasmid genomes reveal past transitions in plasmid mobility

Evidence for gene non-functionalization due to mutational processes is found in genomes in the form of pseudogenes. Pseudogenes are known to be rare in prokaryote chromosomes, with the exception of lineages that underwent an extreme genome reduction (e.g. obligatory symbionts). Much less is known about the frequency of pseudogenes in prokaryotic plasmids; those are genetic elements that can transfer between cells and may encode beneficial traits for their host. Non-functionalization of plasmid-encoded genes may alter the plasmid characteristics, e.g. mobility, or their effect on the host. Analyzing 10 832 prokaryotic genomes, we find that plasmid genomes are characterized by threefold-higher pseudogene density compared to chromosomes. The majority of plasmid pseudogenes correspond to deteriorated transposable elements. A detailed analysis of enterobacterial plasmids furthermore reveals frequent gene non-functionalization events associated with the loss of plasmid self-transmissibility. Reconstructing the evolution of closely related plasmids reveals that non-functionalization of the conjugation machinery led to the emergence of non-mobilizable plasmid types. Examples are virulence plasmids in Escherichia and Salmonella. Our study highlights non-functionalization of core plasmid mobility functions as one route for the evolution of domesticated plasmids. Pseudogenes in plasmids supply insights into past transitions in plasmid mobility that are akin to transitions in bacterial lifestyle.

Adaptation to skin mycobiota promotes antibiotic tolerance in Staphylococcus aureus

The microbiota can promote host health by inhibiting pathogen colonization, yet how host-resident fungi, or the mycobiota, contribute to this process remains unclear. The human skin mycobiota is uniquely stable compared to other body sites and dominated by yeasts of the genus Malassezia . We observe that colonization of human skin by Malassezia sympodialis significantly reduces subsequent colonization by the prominent bacterial pathogen Staphylococcus aureus . M. sympodialis secreted products possess potent bactericidal activity against S. aureus and are sufficient to impair S. aureus skin colonization. This bactericidal activity requires an acidic environment and is exacerbated by free fatty acids, demonstrating a unique synergy with host-derived epidermal defenses. Leveraging experimental evolution to pinpoint mechanisms of S. aureus adaptation in response to the skin mycobiota, we identified multiple mutations in the stringent response regulator Rel that promote survival against M. sympodialis . Similar Rel alleles have been reported in S. aureus clinical isolates, and natural Rel variants are sufficient for tolerance to M. sympodialis antagonism. Partial stringent response activation underlies tolerance to clinical antibiotics, with both laboratory-evolved and natural Rel variants conferring multidrug tolerance. These findings demonstrate the ability of the mycobiota to mediate pathogen colonization resistance, identify new mechanisms of bacterial adaptation in response to fungal antagonism, and reveal the potential for microbiota-driven evolution to shape pathogen antibiotic susceptibility.

Highlights

- M. sympodialis reduces colonization of human skin by S. aureus - Bactericidal activity of M. sympodialis is exacerbated by features of the skin niche - S. aureus Rel variants are sufficient for tolerance to Malassezia antagonism - Evolved tolerance to yeast antagonism coincides with S. aureus multidrug tolerance.

Investigating the impact of insertion sequences and transposons in the genomes of the most significant phytopathogenic bacteria

Genetic variability in phytopathogens is one of the main problems encountered for effective plant disease control. This fact may be related to the presence of transposable elements (TEs), but little is known about their role in host genomes. Here, we performed the most comprehensive analysis of insertion sequences (ISs) and transposons (Tns) in the genomes of the most important bacterial plant pathogens. A total of 35 692 ISs and 71 transposons were identified in 270 complete genomes. The level of pathogen-host specialization was found to be a significant determinant of the element distribution among the species. Some Tns were identified as carrying virulence factors, such as genes encoding effector proteins of the type III secretion system and resistance genes for the antimicrobial streptomycin. Evidence for IS-mediated ectopic recombination was identified in Xanthomonas genomes. Moreover, we found that IS elements tend to be inserted in regions near virulence and fitness genes, such ISs disrupting avirulence genes in X. oryzae genomes. In addition, transcriptome analysis under different stress conditions revealed differences in the expression of genes encoding transposases in the Ralstonia solanacearum, X. oryzae, and P. syringae species. Lastly, we also investigated the role of Tns in regulation via small noncoding regulatory RNAs and found these elements may target plant-cell transcriptional activators. Taken together, the results indicate that TEs may have a fundamental role in variability and virulence in plant pathogenic bacteria.

Increasing transposase abundance with ocean depth correlates with a particle-associated lifestyle

Transposases are mobile genetic elements that move within and between genomes, promoting genomic plasticity in microorganisms. In marine microbial communities, the abundance of transposases increases with depth, but the reasons behind this trend remain unclear. Our analysis of metagenomes from the Tara Oceans and Malaspina Expeditions suggests that a particle-associated lifestyle is the main covariate for the high occurrence of transposases in the deep ocean, and this trend holds true for individual genomes as well as in a community-wide sense. We observed a strong and depth-independent correlation between transposase abundance and the presence of biofilm-associated genes, as well as the prevalence of secretory enzymes. This suggests that mobile genetic elements readily propagate among microbial communities within crowded biofilms. Furthermore, we show that particle association positively correlates with larger genome size, which is in turn associated with higher transposase abundance. Cassette sequences associated with transposons are enriched with genes related to defense mechanisms, which are more highly expressed in the deep sea. Thus, while transposons spread at the expense of their microbial hosts, they also introduce novel genes and potentially benefit the hosts in helping to compete for limited resources. Overall, our results suggest a new understanding of deep ocean particles as highways for gene sharing among defensively oriented microbial genomes.IMPORTANCEGenes can move within and between microbial genomes via mobile genetic elements, which include transposases and transposons. In the oceans, there is a puzzling increase in transposase abundance in microbial genomes as depth increases. To gain insight into this trend, we conducted an extensive analysis of marine microbial metagenomes and metatranscriptomes. We found a significant correlation between transposase abundance and a particle-associated lifestyle among marine microbes at both the metagenome and genome-resolved levels. We also observed a link between transposase abundance and genes related to defense mechanisms. These results suggest that as microbes become densely packed into crowded particles, mobile genes are more likely to spread and carry genetic material that provides a competitive advantage in crowded habitats. This may enable deep sea microbes to effectively compete in such environments.

Comparative genomics and host range analysis of four Ralstonia pseudosolanacearum strains isolated from sunflower reveals genomic and phenotypic differences

BACKGROUND

Bacterial wilt caused by Ralstonia solanacearum species complex (RSSC) is one of the devastating diseases in crop production, seriously reducing the yield of crops. R. pseudosolanacearum, is known for its broad infrasubspecific diversity and comprises 36 sequevars that are currently known. Previous studies found that R. pseudosolanacearum contained four sequevars (13, 14, 17 and 54) isolated from sunflowers sown in the same field.

RESULTS

Here, we provided the complete genomes and the results of genome comparison of the four sequevars strains (RS639, RS642, RS647, and RS650). Four strains showed different pathogenicities to the same cultivars and different host ranges. Their genome sizes were about 5.84 ~ 5.94 Mb, encoding 5002 ~ 5079 genes and the average G + C content of 66.85% ~ 67%. Among the coding genes, 146 ~ 159 specific gene families (contained 150 ~ 160 genes) were found in the chromosomes and 34 ~ 77 specific gene families (contained 34 ~ 78 genes) in the megaplasmids from four strains. The average nucleotide identify (ANI) values between any two strains ranged from 99.05% ~ 99.71%, and the proportion of the total base length of collinear blocks accounts for the total gene length of corresponding genome was all more than 93.82%. Then, we performed a search for genomic islands, prophage sequences, the gene clusters macromolecular secretion systems, type III secreted effectors and other virulence factors in these strains, which provided detailed comparison results of their presence and distinctive features compared to the reference strain GMI1000. Among them, the number and types of T2SS gene clusters were different in the four strains, among which RS650 included all five types. T4SS gene cluster of RS639 and RS647 were missed. In the T6SS gene cluster, several genes were inserted in the RS639, RS647, and RS650, and gene deletion was also detected in the RS642. A total of 78 kinds of type III secreted effectors were found, which included 52 core and 9 specific effectors in four strains.

CONCLUSION

This study not only provided the complete genomes of multiple R. pseudosolanacearum strains isolated from a new host, but also revealed the differences in their genomic levels through comparative genomics. Furthermore, these findings expand human knowledge about the range of hosts that Ralstonia can infect, and potentially contribute to exploring rules and factors of the genetic evolution and analyzing its pathogenic mechanism.

Multiple independent losses of sporulation and peptidoglycan in the Mycoplasmatales and related orders of the class Bacilli

Many peptidoglycan-deficient bacteria such as the Mycoplasmatales are known host-associated lineages, lacking the environmental resistance mechanisms and metabolic capabilities necessary for a free-living lifestyle. Several peptidoglycan-deficient and non-sporulating orders of interest are thought to be descended from Gram-positive sporulating Bacilli through reductive evolution. Here we annotate 2650 genomes belonging to the class Bacilli, according to the Genome Taxonomy Database, to predict the peptidoglycan and sporulation phenotypes of three novel orders, RFN20, RF39 and ML615J-28, known only through environmental sequence surveys. These lineages are interspersed between peptidoglycan-deficient non-sporulating orders including the Mycoplasmatales and Acholeplasmatales, and more typical Gram-positive orders such as the Erysipelotrichales and Staphylococcales. We use the extant genotypes to perform ancestral state reconstructions. The novel orders are predicted to have small genomes with minimal metabolic capabilities and to comprise a mix of peptidoglycan-deficient and/or non-sporulating species. In contrast to expectations based on cultured representatives, the order Erysipelotrichales lacks many of the genes involved in peptidoglycan and endospore formation. The reconstructed evolutionary history of these traits suggests multiple independent whole-genome reductions and loss of phenotype via intermediate transition states that continue into the present. We suggest that the evolutionary history of the reduced-genome lineages within the class Bacilli is one driven by multiple independent transitions to host-associated lifestyles, with the degree of reduction in environmental resistance and metabolic capabilities correlated with degree of host association.

Two male-killing Wolbachia from Drosophila birauraia that are closely related but distinct in genome structure

Insects harbour diverse maternally inherited bacteria and viruses, some of which have evolved to kill the male progeny of their hosts (male killing: MK). The fly species Drosophila biauraria carries a maternally transmitted MK-inducing partiti-like virus, but it was unknown if it carries other MK-inducing endosymbionts. Here, we identified two male-killing Wolbachia strains (wBiau1 and wBiau2) from D. biauraria and compared their genomes to elucidate their evolutionary processes. The two strains were genetically closely related but had exceptionally different genome structures with considerable rearrangements compared with combinations of other Wolbachia strains. Despite substantial changes in the genome structure, the two Wolbachia strains did not experience gene losses that would disrupt the male-killing expression or persistence in the host population. The two Wolbachia-infected matrilines carried distinct mitochondrial haplotypes, suggesting that wBiau1 and wBiau2 have invaded D. biauraria independently and undergone considerable genome changes owing to unknown selective pressures in evolutionary history. This study demonstrated the presence of three male-killers from two distinct origins in one fly species and highlighted the diverse and rapid genome evolution of MK Wolbachia in the host.

Ongoing diversification of the global fish pathogen Piscirickettsia salmonis through genetic isolation and transposition bursts

The management of bacterial pathogens remains a key challenge of aquaculture. The marine gammaproteobacterium Piscirickettsia salmonis is the etiological agent of piscirickettsiosis and causes multi-systemic infections in different salmon species, resulting in considerable mortality and substantial commercial losses. Here, we elucidate its global diversity, evolution, and selection during human interventions. Our comprehensive analysis of 73 closed, high quality genome sequences covered strains from major outbreaks and was supplemented by an analysis of all P. salmonis 16S rRNA gene sequences and metagenomic reads available in public databases. Genome comparison showed that Piscirickettsia comprises at least three distinct, genetically isolated species of which two showed evidence for continuing speciation. However, at least twice the number of species exist in marine fish or seawater. A hallmark of Piscirickettsia diversification is the unprecedented amount and diversity of transposases which are particularly active in subgroups undergoing rapid speciation and are key to the acquisition of novel genes and to pseudogenization. Several group-specific genes are involved in surface antigen synthesis and may explain the differences in virulence between strains. However, the frequent failure of antibiotic treatment of piscirickettsiosis outbreaks cannot be explained by horizontal acquisition of resistance genes which so far occurred only very rarely. Besides revealing a dynamic diversification of an important pathogen, our study also provides the data for improving its surveillance, predicting the emergence of novel lineages, and adapting aquaculture management, and thereby contributes towards the sustainability of salmon farming.

Mosaic environment-driven evolution of the deep-sea mussel Gigantidas platifrons bacterial endosymbiont

BACKGROUND

The within-species diversity of symbiotic bacteria represents an important genetic resource for their environmental adaptation, especially for horizontally transmitted endosymbionts. Although strain-level intraspecies variation has recently been detected in many deep-sea endosymbionts, their ecological role in environmental adaptation, their genome evolution pattern under heterogeneous geochemical environments, and the underlying molecular forces remain unclear.

RESULTS

Here, we conducted a fine-scale metagenomic analysis of the deep-sea mussel Gigantidas platifrons bacterial endosymbiont collected from distinct habitats: hydrothermal vent and methane seep. Endosymbiont genomes were assembled using a pipeline that distinguishes within-species variation and revealed highly heterogeneous compositions in mussels from different habitats. Phylogenetic analysis separated the assemblies into three distinct environment-linked clades. Their functional differentiation follows a mosaic evolutionary pattern. Core genes, essential for central metabolic function and symbiosis, were conserved across all clades. Clade-specific genes associated with heavy metal resistance, pH homeostasis, and nitrate utilization exhibited signals of accelerated evolution. Notably, transposable elements and plasmids contributed to the genetic reshuffling of the symbiont genomes and likely accelerated adaptive evolution through pseudogenization and the introduction of new genes.

CONCLUSIONS

The current study uncovers the environment-driven evolution of deep-sea symbionts mediated by mobile genetic elements. Its findings highlight a potentially common and critical role of within-species diversity in animal-microbiome symbioses. Video Abstract.

Genomic distribution of the insertion sequence IS711 reveal a potential role in Brucella genome plasticity and host preference

The Insertion Sequence 711 (IS711) is linked to the Brucella genus. Mapping the genomic distribution of IS711 can help understand this insertion element's biological and evolutionary role. This work aimed to delineate the genomic distribution of the IS711 element and to study its association with Brucella evolution. A total of 124 genomes representing 9 Brucella species were searched using BLASTn sequence alignment tool to identify complete and truncated copies of IS711. Based on the genomic context, each IS711 locus was assigned a code using the initial letters of its neighboring genes. Various tools were used to annotate the neighboring genes and determine the shared synteny around orthologous IS711 loci. The tool Islandviewer 4 was used to scan for genomic islands. The Codon Tree method was used to build phylogenetic trees of B. melitensis, B. abortus, and B. suis genomes. The phylogenetic trees of the three species were analyzed, taking into account the genomic distribution patterns of IS711. The result of IS711 frequency analysis showed a relatively conserved number of copies/genome for the different species and for some biovars. The analysis showed that Brucella species with a relatively low IS711 copy number (4-8 copies/genome) are linked to domestic animals as primary hosts and have potential for zoonotic transmission. However, species with a relatively higher copy number (12-30 copies/genome) are less zoonotic and tend to be linked with wild animals as primary hosts. Analyzing the genomic distribution map of IS711 loci showed several unique patterns of IS711 distribution that are correlated with the evolution of Brucella species and biovars. The results also showed that 46.2% of the conserved IS711 elements are located within genomic islands. Based on our results and previous data, we postulate a model explaining the IS711 role in Brucella evolution. We assume that during the transition from a free-living to an intracellular lifestyle, a descendant of the Brucella genus had acquired a progenitor sequence of the IS711. Subsequently, a burst in IS711 transposition occurred. This parasitic expansion can be deleterious and has to be counteracted by evolutionary forces to prevent lineage extension and to promote adaptation to host. Similar to other plasmid-free pathogenic α-Proteobacteria bacteria, the balance of expansion and reduction of insertion elements could be one of the mechanisms to control genome reduction and streamlining. We hypothesize that the IS711-mediated genomic changes and other small sequence nucleotide changes in specific orthologous genes could significantly contribute to Brucella's evolution and adaptation to different animal hosts.

Cyclitol metabolism is a central feature of Burkholderia leaf symbionts

The symbioses between plants of the Rubiaceae and Primulaceae families with Burkholderia bacteria represent unique and intimate plant-bacterial relationships. Many of these interactions have been identified through PCR-dependent typing methods, but there is little information available about their functional and ecological roles. We assembled 17 new endophyte genomes representing endophytes from 13 plant species, including those of two previously unknown associations. Genomes of leaf endophytes belonging to Burkholderia s.l. show extensive signs of genome reduction, albeit to varying degrees. Except for one endophyte, none of the bacterial symbionts could be isolated on standard microbiological media. Despite their taxonomic diversity, all endophyte genomes contained gene clusters linked to the production of specialized metabolites, including genes linked to cyclitol sugar analog metabolism and in one instance non-ribosomal peptide synthesis. These genes and gene clusters are unique within Burkholderia s.l. and are likely horizontally acquired. We propose that the acquisition of secondary metabolite gene clusters through horizontal gene transfer is a prerequisite for the evolution of a stable association between these endophytes and their hosts.

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

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

Reduced and Nonreduced Genomes in Paraburkholderia Symbionts of Social Amoebas

The social amoeba Dictyostelium discoideum is a predatory soil protist frequently used for studying host-pathogen interactions. A subset of D. discoideum strains isolated from soil persistently carry symbiotic Paraburkholderia, recently formally described as P. agricolaris, P. bonniea, and P. hayleyella. The three facultative symbiont species of D. discoideum present a unique opportunity to study a naturally occurring symbiosis in a laboratory model protist. There is a large difference in genome size between P. agricolaris (8.7 million base pairs [Mbp]) versus P. hayleyella and P. bonniea (4.1 Mbp). We took a comparative genomics approach and compared the three genomes of D. discoideum symbionts to 12 additional Paraburkholderia genomes to test for genome evolution patterns that frequently accompany host adaptation. Overall, P. agricolaris is difficult to distinguish from other Paraburkholderia based on its genome size and content, but the reduced genomes of P. bonniea and P. hayleyella display characteristics indicative of genome streamlining rather than deterioration during adaptation to their protist hosts. In addition, D. discoideum-symbiont genomes have increased secretion system and motility genes that may mediate interactions with their host. Specifically, adjacent BurBor-like type 3 and T6SS-5-like type 6 secretion system operons shared among all three D. discoideum-symbiont genomes may be important for host interaction. Horizontal transfer of these secretion system operons within the amoeba host environment may have contributed to the unique ability of these symbionts to establish and maintain a symbiotic relationship with D. discoideum. IMPORTANCE Protists are a diverse group of typically single cell eukaryotes. Bacteria and archaea that form long-term symbiotic relationships with protists may evolve in additional ways than those in relationships with multicellular eukaryotes such as plants, animals, or fungi. Social amoebas are a predatory soil protist sometimes found with symbiotic bacteria living inside their cells. They present a unique opportunity to explore a naturally occurring symbiosis in a protist frequently used for studying host-pathogen interactions. We show that one amoeba-symbiont species is similar to other related bacteria in genome size and content, while the two reduced-genome-symbiont species show characteristics of genome streamlining rather than deterioration during adaptation to their host. We also identify sets of genes present in all three amoeba-symbiont genomes that are potentially used for host-symbiont interactions. Because the amoeba symbionts are distantly related, the amoeba host environment may be where these genes were shared among symbionts.

Metagenomic methylation patterns resolve bacterial genomes of unusual size and structural complexity

The plasticity of bacterial and archaeal genomes makes examining their ecological and evolutionary dynamics both exciting and challenging. The same mechanisms that enable rapid genomic change and adaptation confound current approaches for recovering complete genomes from metagenomes. Here, we use strain-specific patterns of DNA methylation to resolve complex bacterial genomes from long-read metagenomic data of a marine microbial consortium, the “pink berries” of the Sippewissett Marsh (USA). Unique combinations of restriction-modification (RM) systems encoded by the bacteria produced distinctive methylation profiles that were used to accurately bin and classify metagenomic sequences. Using this approach, we finished the largest and most complex circularized bacterial genome ever recovered from a metagenome (7.9 Mb with >600 transposons), the finished genome of Thiohalocapsa sp. PB-PSB1 the dominant bacteria in the consortia. From genomes binned by methylation patterns, we identified instances of horizontal gene transfer between sulfur-cycling symbionts (Thiohalocapsa sp. PB-PSB1 and Desulfofustis sp. PB-SRB1), phage infection, and strain-level structural variation. We also linked the methylation patterns of each metagenome-assembled genome with encoded DNA methyltransferases and discovered new RM defense systems, including novel associations of RM systems with RNase toxins.

Ecology and evolution of chlamydial symbionts of arthropods

The phylum Chlamydiae consists of obligate intracellular bacteria including major human pathogens and diverse environmental representatives. Here we investigated the Rhabdochlamydiaceae, which is predicted to be the largest and most diverse chlamydial family, with the few described members known to infect arthropod hosts. Using published 16 S rRNA gene sequence data we identified at least 388 genus-level lineages containing about 14 051 putative species within this family. We show that rhabdochlamydiae are mainly found in freshwater and soil environments, suggesting the existence of diverse, yet unknown hosts. Next, we used a comprehensive genome dataset including metagenome assembled genomes classified as members of the family Rhabdochlamydiaceae, and we added novel complete genome sequences of Rhabdochlamydia porcellionis infecting the woodlouse Porcellio scaber, and of 'Candidatus R. oedothoracis' associated with the linyphiid dwarf spider Oedothorax gibbosus. Comparative analysis of basic genome features and gene content with reference genomes of well-studied chlamydial families with known host ranges, namely Parachlamydiaceae (protist hosts) and Chlamydiaceae (human and other vertebrate hosts) suggested distinct niches for members of the Rhabdochlamydiaceae. We propose that members of the family represent intermediate stages of adaptation of chlamydiae from protists to vertebrate hosts. Within the genus Rhabdochlamydia, pronounced genome size reduction could be observed (1.49-1.93 Mb). The abundance and genomic distribution of transposases suggests transposable element expansion and subsequent gene inactivation as a mechanism of genome streamlining during adaptation to new hosts. This type of genome reduction has never been described before for any member of the phylum Chlamydiae. This study provides new insights into the molecular ecology, genomic diversity, and evolution of representatives of one of the most divergent chlamydial families.

High Levels of Multiple Phage WO Infections and Its Evolutionary Dynamics Associated With Wolbachia-Infected Butterflies

Wolbachia is a maternally inherited bacterium that is widely distributed among arthropods, in which it manipulates the reproduction of its hosts. Phage WO is the only bacteriophage known to infect Wolbachia, and may provide benefit to its host or arthropods. We screened for the presence of phage WO in Wolbachia-infected butterfly species for the first time, to investigate their diversity and evolutionary dynamics. All Wolbachia-infected butterfly species, including members of the families Hesperiidae, Lycaenidae, Nymphalidae, Papilionidae, and Pieridae, were found to harbor phage WO. Interestingly, 84% of 19 butterfly species, which were infected with a single Wolbachia strain harbored high levels of multiple phage types (ranging from 3 to 17 types), another three species harbored one or two phage types. For Wolbachia strains (ST-41, ST-19, ST-125 and ST-374) shared among various butterfly species, their host insects all harbored multiple phage types, while two Wolbachia strains (ST-297 and ST-wPcau) were found to infect one butterfly species, whose insect hosts harbored a single phage type, suggesting that horizontal transfer of Wolbachia between insects increased the likelihood of exposure to phages, resulting in increased phage genetic diversity. Twelve horizontal transmission events of phage WO were found, which shared common phage WO types among different Wolbachia strains associated with butterflies. Most horizontal transfer events involved different Wolbachia supergroups (A and B). Horizontal acquisition of phage WO might also occur between eukaryotes without Wolbachia transfer. Furthermore, 22 putative recombination events were identified in 13 of 16 butterfly species which harbored multiple phage types. These results showed that horizontal transfer of Wolbachia caused it to be exposed to the phage gene pool, and that horizontal transmission of phage WO, as well as intragenic recombination were important dynamics for phage WO genome evolution, which effectively promoted the high level of phage WO diversity associated with butterflies.

Horizontal gene transfer and adaptive evolution in bacteria

Horizontal gene transfer (HGT) is arguably the most conspicuous feature of bacterial evolution. Evidence for HGT is found in most bacterial genomes. Although HGT can considerably alter bacterial genomes, not all transfer events may be biologically significant and may instead represent the outcome of an incessant evolutionary process that only occasionally has a beneficial purpose. When adaptive transfers occur, HGT and positive selection may result in specific, detectable signatures in genomes, such as gene-specific sweeps or increased transfer rates for genes that are ecologically relevant. In this Review, we first discuss the various mechanisms whereby HGT occurs, how the genetic signatures shape patterns of genomic variation and the distinct bioinformatic algorithms developed to detect these patterns. We then discuss the evolutionary theory behind HGT and positive selection in bacteria, and discuss the approaches developed over the past decade to detect transferred DNA that may be involved in adaptation to new environments.

A high-quality reference genome for the fish pathogen Streptococcus iniae

Fish mortality caused by Streptococcus iniae is a major economic problem in aquaculture in warm and temperate regions globally. There is also risk of zoonotic infection by S. iniae through handling of contaminated fish. In this study, we present the complete genome sequence of S. iniae strain QMA0248, isolated from farmed barramundi in South Australia. The 2.12 Mb genome of S. iniae QMA0248 carries a 32 kb prophage, a 12 kb genomic island and 92 discrete insertion sequence (IS) elements. These include nine novel IS types that belong mostly to the IS3 family. Comparative and phylogenetic analysis between S. iniae QMA0248 and publicly available complete S. iniae genomes revealed discrepancies that are probably due to misassembly in the genomes of isolates ISET0901 and ISNO. Long-range PCR confirmed five rRNA loci in the PacBio assembly of QMA0248, and, unlike S. iniae 89353, no tandemly repeated rRNA loci in the consensus genome. However, we found sequence read evidence that the tandem rRNA repeat existed within a subpopulation of the original QMA0248 culture. Subsequent nanopore sequencing revealed that the tandem rRNA repeat was the most prevalent genotype, suggesting that there is selective pressure to maintain fewer rRNA copies under uncertain laboratory conditions. Our study not only highlights assembly problems in existing genomes, but provides a high-quality reference genome for S. iniae QMA0248, including manually curated mobile genetic elements, that will assist future S. iniae comparative genomic and evolutionary studies.

Evolutionarily recent dual obligatory symbiosis among adelgids indicates a transition between fungus- and insect-associated lifestyles

Adelgids (Insecta: Hemiptera: Adelgidae) form a small group of insects but harbor a surprisingly diverse set of bacteriocyte-associated endosymbionts, which suggest multiple replacement and acquisition of symbionts over evolutionary time. Specific pairs of symbionts have been associated with adelgid lineages specialized on different secondary host conifers. Using a metagenomic approach, we investigated the symbiosis of the Adelges laricis/Adelges tardus species complex containing betaproteobacterial ("Candidatus Vallotia tarda") and gammaproteobacterial ("Candidatus Profftia tarda") symbionts. Genomic characteristics and metabolic pathway reconstructions revealed that Vallotia and Profftia are evolutionary young endosymbionts, which complement each other's role in essential amino acid production. Phylogenomic analyses and a high level of genomic synteny indicate an origin of the betaproteobacterial symbiont from endosymbionts of Rhizopus fungi. This evolutionary transition was accompanied with substantial loss of functions related to transcription regulation, secondary metabolite production, bacterial defense mechanisms, host infection, and manipulation. The transition from fungus to insect endosymbionts extends our current framework about evolutionary trajectories of host-associated microbes.

The physiology and metabolic properties of a novel, low-abundance Psychrilyobacter species isolated from the anoxic Black Sea shed light on its ecological role

Members of the Psychrilyobacter spp. of the phylum Fusobacteria have been recently suggested to be amongst the most significant primary degraders of the detrital organic matter in sulfidic marine habitats, despite representing only a small proportion (<0.1%) of the microbial community. In this study, we have isolated a previously uncultured Psychrilyobacter species (strains SD5^T and BL5; Psychrilyobacter piezotolerans sp. nov.) from the sulfidic waters (i.e., 2000 m depth) of the Black Sea and investigated its physiology and genomic capability in order to better understand potential ecological adaptation strategies. P. piezotolerans utilized a broad range of organic substituents (carbohydrates and proteins) and, remarkably, grew at sulfide concentrations up to 32 mM. These flexible physiological properties were supported by the presence of the respective metabolic pathways in the genomes of both strains. Growth at varying hydrostatic pressure (0.1-50 MPa) was sustained by modifying its membrane lipid composition. Thus, we have isolated a novel member of the 'rare biosphere', which endures the extreme conditions and may play a significant role in the degradation of detrital organic matter sinking into the sulfidic waters of the Black Sea.

Genic Selection Within Prokaryotic Pangenomes

Understanding the evolutionary forces shaping prokaryotic pangenome structure is a major goal of microbial evolution research. Recent work has highlighted that a substantial proportion of accessory genes appear to confer niche-specific adaptations. This work has primarily focused on selection acting at the level of individual cells. Herein, we discuss a lower level of selection that also contributes to pangenome variation: genic selection. This refers to cases where genetic elements, rather than individual cells, are the entities under selection. The clearest examples of this form of selection are selfish mobile genetic elements, which are those that have either a neutral or a deleterious effect on host fitness. We review the major classes of these and other mobile elements and discuss the characteristic features of such elements that could be under genic selection. We also discuss how genetic elements that are beneficial to hosts can also be under genic selection, a scenario that may be more prevalent but not widely appreciated, because disentangling the effects of selection at different levels (i.e., organisms vs. genes) is challenging. Nonetheless, an appreciation for the potential action and implications of genic selection is important to better understand the evolution of prokaryotic pangenomes.

A Unique Reverse Adaptation Mechanism Assists Bordetella pertussis in Resistance to Both Scarcity and Toxicity of Manganese

The ability of bacterial pathogens to acquire essential micronutrients is critical for their survival in the host environment. Manganese plays a complex role in the virulence of a variety of pathogens due to its function as an antioxidant and enzymatic cofactor. Therefore, host cells deprive pathogens of manganese to prevent or attenuate infection. Here, we show that evolution of the human-restricted pathogen Bordetella pertussis has selected for an inhibitory duplication within a manganese exporter of the calcium:cation antiporter superfamily. Intriguingly, upon exposure to toxic levels of manganese, the nonfunctional exporter becomes operative in resister cells due to a unique reverse adaptation mechanism. However, compared with wild-type (wt) cells, the resisters carrying a functional copy of the exporter displayed strongly reduced intracellular levels of manganese and impaired growth under oxidative stress. Apparently, inactivation of the manganese exporter and the resulting accumulation of manganese in the cytosol benefited the pathogen by improving its survival under stress conditions. The inhibitory duplication within the exporter gene is highly conserved among B. pertussis strains, absent from all other Bordetella species and from a vast majority of organisms across all kingdoms of life. Therefore, we conclude that inactivation of the exporter gene represents an exceptional example of a flexible genome decay strategy employed by a human pathogen to adapt to its exclusive host.

Adaptation and genomic erosion in fragmented Pseudomonas aeruginosa populations in the sinuses of people with cystic fibrosis

Pseudomonas aeruginosa notoriously adapts to the airways of people with cystic fibrosis (CF), yet how infection-site biogeography and associated evolutionary processes vary as lifelong infections progress remains unclear. Here we test the hypothesis that early adaptations promoting aggregation influence evolutionary-genetic trajectories by examining longitudinal P. aeruginosa from the sinuses of six adults with CF. Highly host-adapted lineages harbored mutator genotypes displaying signatures of early genome degradation associated with recent host restriction. Using an advanced imaging technique (MiPACT-HCR [microbial identification after passive clarity technique]), we find population structure tracks with genome degradation, with the most host-adapted, genome-degraded P. aeruginosa (the mutators) residing in small, sparse aggregates. We propose that following initial adaptive evolution in larger populations under strong selection for aggregation, P. aeruginosa persists in small, fragmented populations that experience stronger effects of genetic drift. These conditions enrich for mutators and promote degenerative genome evolution. Our findings underscore the importance of infection-site biogeography to pathogen evolution.

Development, Phenotypic Characterization and Genomic Analysis of a Francisella tularensis Panel for Tularemia Vaccine Testing

Francisella tularensis is one of several biothreat agents for which a licensed vaccine is needed to protect against this pathogen. To aid in the development of a vaccine protective against pneumonic tularemia, we generated and characterized a panel of F. tularensis isolates that can be used as challenge strains to assess vaccine efficacy. Our panel consists of both historical and contemporary isolates derived from clinical and environmental sources, including human, tick, and rabbit isolates. Whole genome sequencing was performed to assess the genetic diversity in comparison to the reference genome F. tularensis Schu S4. Average nucleotide identity analysis showed >99% genomic similarity across the strains in our panel, and pan-genome analysis revealed a core genome of 1,707 genes, and an accessory genome of 233 genes. Three of the strains in our panel, FRAN254 (tick-derived), FRAN255 (a type B strain), and FRAN256 (a human isolate) exhibited variation from the other strains. Moreover, we identified several unique mutations within the Francisella Pathogenicity Island across multiple strains in our panel, revealing unexpected diversity in this region. Notably, FRAN031 (Scherm) completely lacked the second pathogenicity island but retained virulence in mice. In contrast, FRAN037 (Coll) was attenuated in a murine pneumonic tularemia model and had mutations in pdpB and iglA which likely led to attenuation. All of the strains, except FRAN037, retained full virulence, indicating their effectiveness as challenge strains for future vaccine testing. Overall, we provide a well-characterized panel of virulent F. tularensis strains that can be utilized in ongoing efforts to develop an effective vaccine against pneumonic tularemia to ensure protection is achieved across a range F. tularensis strains.

Enhanced Mutation Rate, Relaxed Selection, and the "Domino Effect" are associated with Gene Loss in Blattabacterium, A Cockroach Endosymbiont

Intracellular endosymbionts have reduced genomes that progressively lose genes at a timescale of tens of million years. We previously reported that gene loss rate is linked to mutation rate in Blattabacterium, however, the mechanisms causing gene loss are not yet fully understood. Here, we carried out comparative genomic analyses on the complete genome sequences of a representative set of 67 Blattabacterium strains, with sizes ranging between 511 and 645 kb. We found that 200 of the 566 analyzed protein-coding genes were lost in at least one lineage of Blattabacterium, with the most extreme case being one gene that was lost independently in 24 lineages. We found evidence for three mechanisms influencing gene loss in Blattabacterium. First, gene loss rates were found to increase exponentially with the accumulation of substitutions. Second, genes involved in vitamin and amino acid metabolism experienced relaxed selection in Cryptocercus and Mastotermes, possibly triggered by their vertically inherited gut symbionts. Third, we found evidence of epistatic interactions among genes leading to a "domino effect" of gene loss within pathways. Our results highlight the complexity of the process of genome erosion in an endosymbiont.

Design and Testing of Effective Primers for Amplification of the orf7 Gene of Phage WO Associated with Andricus hakonensis

Simple Summary

Andricus hakonensis is thought to contain the most complex and diverse phage types known and should be an ideal model material for studying interactions among bacteriophages, bacteria and eukaryotes. As shown in previous studies, existing primers are not effective enough to amplify all virus groups in A. hakonensis. Based on a comprehensive analysis of all virus groups reported to date, we designed a relatively conservative primer for virus detection. This primer can accurately and efficiently detect the presence of phage WO in arthropod hosts. Using gene alignment, clear evidence was provided for the existence of hitherto unreported base deletions, which are an important cause of diversity in phage WO associated with A. hakonensis.

Abstract

Phage WO was first characterized in Wolbachia, an obligate intracellular Rickettsiales known for its ability to regulate the reproduction of arthropod hosts. In this paper, we focus on the study of virus diversity in Andricus hakonensis and the development of highly effective primers. Based on the existing Wolbachia genome sequence, we designed primers (WO-TF and WO-TR) to amplify the full-length orf7 gene of phage WO. Surprisingly, sequencing results showed a high abundance of other phage WO groups in A. hakonensis, in addition to the four groups previously identified. The results also showed that A. hakonensis contained most of the known types of orf7 genes (I, III, IV, V and VI) and the level of diversity of harbored phage WO was very high. Therefore, we speculated that existing primers were not specific enough and that new primers for the detection of phage WO were needed. Based on the existing orf7 gene sequence, we designed specific detection primers (WO-SUF and WO-SUR). Sequencing results showed that the primers effectively amplified all known types of phage WO. In addition to amplifying most of the known sequences, we also detected some new genotypes in A. hakonensis using the new primers. Importantly, all phage WO groups could be efficiently detected. Combined with the results of previous studies, our results suggest that A. hakonensis contains the largest number of phage types (up to 36 types). This study is novel in that it provides practical molecular evidence supporting base deletions, in addition to gene mutations and genetic recombination, as an important cause of phage WO diversity.

Gulf of Mexico blue hole harbors high levels of novel microbial lineages

Exploration of oxygen-depleted marine environments has consistently revealed novel microbial taxa and metabolic capabilities that expand our understanding of microbial evolution and ecology. Marine blue holes are shallow karst formations characterized by low oxygen and high organic matter content. They are logistically challenging to sample, and thus our understanding of their biogeochemistry and microbial ecology is limited. We present a metagenomic and geochemical characterization of Amberjack Hole on the Florida continental shelf (Gulf of Mexico). Dissolved oxygen became depleted at the hole's rim (32 m water depth), remained low but detectable in an intermediate hypoxic zone (40-75 m), and then increased to a secondary peak before falling below detection in the bottom layer (80-110 m), concomitant with increases in nutrients, dissolved iron, and a series of sequentially more reduced sulfur species. Microbial communities in the bottom layer contained heretofore undocumented levels of the recently discovered phylum Woesearchaeota (up to 58% of the community), along with lineages in the bacterial Candidate Phyla Radiation (CPR). Thirty-one high-quality metagenome-assembled genomes (MAGs) showed extensive biochemical capabilities for sulfur and nitrogen cycling, as well as for resisting and respiring arsenic. One uncharacterized gene associated with a CPR lineage differentiated hypoxic from anoxic zone communities. Overall, microbial communities and geochemical profiles were stable across two sampling dates in the spring and fall of 2019. The blue hole habitat is a natural marine laboratory that provides opportunities for sampling taxa with under-characterized but potentially important roles in redox-stratified microbial processes.

At the Gate of Mutualism: Identification of Genomic Traits Predisposing to Insect-Bacterial Symbiosis in Pathogenic Strains of the Aphid Symbiont Serratia symbiotica

Mutualistic associations between insects and heritable bacterial symbionts are ubiquitous in nature. The aphid symbiont Serratia symbiotica is a valuable candidate for studying the evolution of bacterial symbiosis in insects because it includes a wide diversity of strains that reflect the diverse relationships in which bacteria can be engaged with insects, from pathogenic interactions to obligate intracellular mutualism. The recent discovery of culturable strains, which are hypothesized to resemble the ancestors of intracellular strains, provide an opportunity to study the mechanisms underlying bacterial symbiosis in its early stages. In this study, we analyzed the genomes of three of these culturable strains that are pathogenic to aphid hosts, and performed comparative genomic analyses including mutualistic host-dependent strains. All three genomes are larger than those of the host-restricted S. symbiotica strains described so far, and show significant enrichment in pseudogenes and mobile elements, suggesting that these three pathogenic strains are in the early stages of the adaptation to their host. Compared to their intracellular mutualistic relatives, the three strains harbor a greater diversity of genes coding for virulence factors and metabolic pathways, suggesting that they are likely adapted to infect new hosts and are a potential source of metabolic innovation for insects. The presence in their genomes of secondary metabolism gene clusters associated with the production of antimicrobial compounds and phytotoxins supports the hypothesis that S. symbiotia symbionts evolved from plant-associated strains and that plants may serve as intermediate hosts. Mutualistic associations between insects and bacteria are the result of independent transitions to endosymbiosis initiated by the acquisition of environmental progenitors. In this context, the genomes of free-living S. symbiotica strains provide a rare opportunity to study the inventory of genes held by bacterial associates of insects that are at the gateway to a host-dependent lifestyle.

A Case of Intragenic Recombination Dramatically Impacting the Phage WO Genetic Diversity in Gall Wasps

The phage WO was characterized in Wolbachia, a strictly intracellular bacterium causing several reproductive alterations in its arthropod hosts. This study aimed to screen the presence of Wolbachia and phage WO in 15 gall wasp species from six provinces of southern China to investigate their diversity and prevalence patterns. A high incidence of Wolbachia infection was determined in the gall wasp species, with an infection rate of 86.7% (13/15). Moreover, seven species had double or multiple infections. All Wolbachia-infected gall wasp species were found to harbor phage WO. The gall wasp species infected with a single Wolbachia strain were found to harbor a single phage WO type. On the contrary, almost all species with double or multiple Wolbachia infections harbored a high level of phage WO diversity (ranging from three to 27 types). Six horizontal transfer events of phage WO in Wolbachia were found to be associated with gall wasps, which shared identical orf7 sequences among their respective accomplices. The transfer potentially took place through gall inducers and associated inquilines infected with or without Wolbachia. Furthermore, 10 putative recombination events were identified from Andricus hakonensis and Andricus sp2, which harbored multiple phage WO types, suggesting that intragenic recombination was the important evolutionary force, which effectively promoted the high level of phage WO diversity associated with gall wasps.

Radical genome remodelling accompanied the emergence of a novel host-restricted bacterial pathogen

The emergence of new pathogens is a major threat to public and veterinary health. Changes in bacterial habitat such as a switch in host or disease tropism are typically accompanied by genetic diversification. Staphylococcus aureus is a multi-host bacterial species associated with human and livestock infections. A microaerophilic subspecies, Staphylococcus aureus subsp. anaerobius, is responsible for Morel’s disease, a lymphadenitis restricted to sheep and goats. However, the evolutionary history of S. aureus subsp. anaerobius and its relatedness to S. aureus are unknown. Population genomic analyses of clinical S. aureus subsp. anaerobius isolates revealed a highly conserved clone that descended from a S. aureus progenitor about 1000 years ago before differentiating into distinct lineages that contain African and European isolates. S. aureus subsp. anaerobius has undergone limited clonal expansion, with a restricted population size, and an evolutionary rate 10-fold slower than S. aureus. The transition to its current restricted ecological niche involved acquisition of a pathogenicity island encoding a ruminant host-specific effector of abscess formation, large chromosomal re-arrangements, and the accumulation of at least 205 pseudogenes, resulting in a highly fastidious metabolism. Importantly, expansion of ~87 insertion sequences (IS) located largely in intergenic regions provided distinct mechanisms for the control of expression of flanking genes, including a novel mechanism associated with IS-mediated anti-anti-sense decoupling of ancestral gene repression. Our findings reveal the remarkable evolutionary trajectory of a host-restricted bacterial pathogen that resulted from extensive remodelling of the S. aureus genome through an array of diverse mechanisms in parallel.

The emergence of new pathogens is a major threat to public and veterinary health. Some bacteria such as Staphylococcus aureus, have the capacity to infect many different host species including humans and livestock while others such as the closely-related S. aureus subsp. anaerobius, associated with a single type of pathology called Morel’s disease in small ruminants, are highly niche-restricted. However, our understanding of the genetic basis for such differences in bacterial host-tropism is very limited. Here, we discovered that S. aureus subsp. anaerobius evolved from an S. aureus ancestor and underwent an array of extensive changes to its genome that accompanied the transition to its current restricted lifestyle. We observed genome decay involving loss of function of hundreds of genes, large intra-chromosomal rearrangements affecting most of the genome, acquisition of a pathogenicity island, and expansion of large numbers of insertion sequences that are inserted at intergenic sites around the genome. Importantly, we found that IS elements affect the expression of neighbouring genes in different ways including a novel mechanism of IS-enabled disruption of ancestral gene repression. Taken together, we provide a remarkable example of radical genomic changes associated with evolutionary transition from a multi-host to highly restricted host ecology.

Whole-Genome Analysis of Mycobacterium avium subsp. paratuberculosis IS900 -

Mycobacterium avium subsp. paratuberculosis (Map) is the etiological agent of Johne’s disease in ruminants. The IS900 insertion sequence (IS) has been used widely as an epidemiological marker and target for PCR diagnosis. Updated DNA sequencing technologies have led to a rapid increase in available Map genomes, which makes it possible to analyze the distribution of IS900 in this slow-growing bacterium. The objective of this study is to characterize the distribution of the IS900 element and how it affects genomic evolution and gene function of Map. A secondary goal is to develop automated in silico restriction fragment length polymorphism (RFLP) analysis using IS900. Complete genomes from the major phylogenetic lineages known as C-type and S-type (including subtypes I and III), were chosen to represent the genetic diversity of Map. IS900 elements were located in these genomes using BLAST software and the relevant fragments extracted. An in silico RFLP analysis using the BstEII restriction site was performed to obtain exact sizes of the DNA fragments carrying a copy of IS900 and the resulting RFLP profiles were analyzed and compared by digital visualization of the separated restriction fragments. The program developed for this study allowed automated localization of IS900 sequences to identify their position within each genome along with the exact number of copies per genome. The number of IS900 copies ranged from 16 in the C-type isolate to 22 in the S-type subtype I isolate. A loci-by-loci sequence alignment of all IS900 copies within the three genomes revealed new sequence polymorphisms that define three sequevars distinguishing the subtypes. Nine IS900 insertion site locations were conserved across all genomes studied while smaller subsets were unique to a particular lineage. Preferential insertion motif sequences were identified for IS900 along with genes bordering all IS900 insertions. Rarely did IS900 insert within coding sequences as only three genes were disrupted in this way. This study makes it possible to automate IS900 distribution in Map genomes to enrich knowledge on the distribution dynamics of this IS for epidemiological purposes, for understanding Map evolution and for studying the biological implications of IS900 insertions.

Comparison of 15 dinoflagellate genomes reveals extensive sequence and structural divergence in family Symbiodiniaceae and genus Symbiodinium

Background

Dinoflagellates in the family Symbiodiniaceae are important photosynthetic symbionts in cnidarians (such as corals) and other coral reef organisms. Breakdown of the coral-dinoflagellate symbiosis due to environmental stress (i.e. coral bleaching) can lead to coral death and the potential collapse of reef ecosystems. However, evolution of Symbiodiniaceae genomes, and its implications for the coral, is little understood. Genome sequences of Symbiodiniaceae remain scarce due in part to their large genome sizes (1–5 Gbp) and idiosyncratic genome features.

Results

Here, we present de novo genome assemblies of seven members of the genus Symbiodinium, of which two are free-living, one is an opportunistic symbiont, and the remainder are mutualistic symbionts. Integrating other available data, we compare 15 dinoflagellate genomes revealing high sequence and structural divergence. Divergence among some Symbiodinium isolates is comparable to that among distinct genera of Symbiodiniaceae. We also recovered hundreds of gene families specific to each lineage, many of which encode unknown functions. An in-depth comparison between the genomes of the symbiotic Symbiodinium tridacnidorum (isolated from a coral) and the free-living Symbiodinium natans reveals a greater prevalence of transposable elements, genetic duplication, structural rearrangements, and pseudogenisation in the symbiotic species.

Conclusions

Our results underscore the potential impact of lifestyle on lineage-specific gene-function innovation, genome divergence, and the diversification of Symbiodinium and Symbiodiniaceae. The divergent features we report, and their putative causes, may also apply to other microbial eukaryotes that have undergone symbiotic phases in their evolutionary history.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-00994-6.

A horizontally acquired expansin gene increases virulence of the emerging plant pathogen Erwinia tracheiphila

Erwinia tracheiphila is a bacterial plant pathogen that causes a fatal wilt infection in some cucurbit crop plants. Wilt symptoms are thought to be caused by systemic bacterial colonization through xylem that impedes sap flow. However, the genetic determinants of within-plant movement are unknown for this pathogen species. Here, we find that E. tracheiphila has horizontally acquired an operon with a microbial expansin (exlx) gene adjacent to a glycoside hydrolase family 5 (gh5) gene. Plant inoculation experiments with deletion mutants in the individual genes (Δexlx and Δgh5) and the full operon (Δexlx-gh5) resulted in decreased severity of wilt symptoms, decreased mortality rate, and impaired systemic colonization compared to the Wt strain. Co-inoculation experiments with Wt and Δexlx-gh5 rescued the movement defect of the mutant strain, suggesting that expansin and GH5 function extracellularly. Together, these results show that expansin-GH5 contributes to systemic movement through xylem, leading to rapid wilt symptom development and higher rates of plant death. The presence of expansin genes in diverse species of bacterial and fungal wilt-inducing pathogens suggests that microbial expansin proteins may be an under-appreciated virulence factor for many pathogen species.

The effects of taxonomy, diet, and ecology on the microbiota of riverine macroinvertebrates

Freshwater macroinvertebrates play key ecological roles in riverine food webs, such as the transfer of nutrients to consumers and decomposition of organic matter. Although local habitat quality drives macroinvertebrate diversity and abundance, little is known about their microbiota. In most animals, the microbiota provides benefits, such as increasing the rate at which nutrients are metabolized, facilitating immune system development, and defending against pathogenic attack. Our objectives were to identify the bacteria within aquatic invertebrates and determine whether their composition varied with taxonomy, habitat, diet, and time of sample collection. In 2016 and 2017, we collected 264 aquatic invertebrates from the mainstem Saint John (Wolastoq) River in New Brunswick, Canada, representing 15 orders. We then amplified the V3-V4 hypervariable region of the 16S rRNA gene within each individual, which revealed nearly 20,000 bacterial operational taxonomic units (OTUs). The microbiota across all aquatic invertebrates were dominated by Proteobacteria (69.25% of the total sequence reads), but they differed significantly in beta diversity, both among host invertebrate taxa (genus-, family-, and order-levels) and temporally. In contrast to previous work, we observed no microbiota differences among functional feeding groups or traditional feeding habits, and neither water velocity nor microhabitat type structured microbiota variability. Our findings suggest that host invertebrate taxonomy was the most important factor in modulating the composition of the microbiota, likely through a combination of vertical and horizontal bacterial transmission, and evolutionary processes. This is one of the most comprehensive studies of freshwater invertebrate microbiota to date, and it underscores the need for future studies of invertebrate microbiota evolution and linkages to environmental bacteria and physico-chemical conditions.

Pathogenic signature of invasive non-typhoidal Salmonella in Africa: implications for vaccine development

Invasive non-typhoidal Salmonella (iNTS) infections are a leading cause of bacteremia in Sub-Saharan Africa (sSA), thereby representing a major public health threat. Salmonella Typhimurium clade ST313 and Salmonella Enteriditis lineages associated with Western and Central/Eastern Africa are among the iNTS serovars which are of the greatest concern due to their case-fatality rate, especially in children and in the immunocompromised population. Identification of pathogen-associated features and host susceptibility factors that increase the risk for invasive non-typhoidal salmonellosis would be instrumental for the design of targeted prevention strategies, which are urgently needed given the increasing spread of multidrug-resistant iNTS in Africa. This review summarizes current knowledge of bacterial traits and host immune responses associated with iNTS infections in sSA, then discusses how this knowledge can guide vaccine development while providing a summary of vaccine candidates in preclinical and early clinical development.

Impact of insertion sequences on convergent evolution of Shigella species

Shigella species are specialised lineages of Escherichia coli that have converged to become human-adapted and cause dysentery by invading human gut epithelial cells. Most studies of Shigella evolution have been restricted to comparisons of single representatives of each species; and population genomic studies of individual Shigella species have focused on genomic variation caused by single nucleotide variants and ignored the contribution of insertion sequences (IS) which are highly prevalent in Shigella genomes. Here, we investigate the distribution and evolutionary dynamics of IS within populations of Shigella dysenteriae Sd1, Shigella sonnei and Shigella flexneri. We find that five IS (IS1, IS2, IS4, IS600 and IS911) have undergone expansion in all Shigella species, creating substantial strain-to-strain variation within each population and contributing to convergent patterns of functional gene loss within and between species. We find that IS expansion and genome degradation are most advanced in S. dysenteriae and least advanced in S. sonnei; and using genome-scale models of metabolism we show that Shigella species display convergent loss of core E. coli metabolic capabilities, with S. sonnei and S. flexneri following a similar trajectory of metabolic streamlining to that of S. dysenteriae. This study highlights the importance of IS to the evolution of Shigella and provides a framework for the investigation of IS dynamics and metabolic reduction in other bacterial species.

Dynamics of Insect-Microbiome Interaction Influence Host and Microbial Symbiont

Insects share an intimate relationship with their gut microflora and this symbiotic association has developed into an essential evolutionary outcome intended for their survival through extreme environmental conditions. While it has been clearly established that insects, with very few exceptions, associate with several microbes during their life cycle, information regarding several aspects of these associations is yet to be fully unraveled. Acquisition of bacteria by insects marks the onset of microbial symbiosis, which is followed by the adaptation of these bacterial species to the gut environment for prolonged sustenance and successful transmission across generations. Although several insect-microbiome associations have been reported and each with their distinctive features, diversifications and specializations, it is still unclear as to what led to these diversifications. Recent studies have indicated the involvement of various evolutionary processes operating within an insect body that govern the transition of a free-living microbe to an obligate or facultative symbiont and eventually leading to the establishment and diversification of these symbiotic relationships. Data from various studies, summarized in this review, indicate that the symbiotic partners, i.e., the bacteria and the insect undergo several genetic, biochemical and physiological changes that have profound influence on their life cycle and biology. An interesting outcome of the insect-microbe interaction is the compliance of the microbial partner to its eventual genome reduction. Endosymbionts possess a smaller genome as compared to their free-living forms, and thus raising the question what is leading to reductive evolution in the microbial partner. This review attempts to highlight the fate of microbes within an insect body and its implications for both the bacteria and its insect host. While discussion on each specific association would be too voluminous and outside the scope of this review, we present an overview of some recent studies that contribute to a better understanding of the evolutionary trajectory and dynamics of the insect-microbe association and speculate that, in the future, a better understanding of the nature of this interaction could pave the path to a sustainable and environmentally safe way for controlling economically important pests of crop plants.

Comparative genomics of whole-cell pertussis vaccine strains from India

BACKGROUND

Despite high vaccination coverage using acellular (ACV) and whole-cell pertussis (WCV) vaccines, the resurgence of pertussis is observed globally. Genetic divergence in circulating strains of Bordetella pertussis has been reported as one of the contributing factors for the resurgence of the disease. Our current knowledge of B. pertussis genetic evolution in circulating strains is mostly based on studies conducted in countries using ACVs targeting only a few antigens used in the production of ACVs. To better understand the adaptation to vaccine-induced selection pressure, it will be essential to study B. pertussis populations in developing countries which are using WCVs. India is a significant user and global supplier of WCVs. We report here comparative genome analyses of vaccine and clinical isolates reported from India. Whole-genome sequences obtained from vaccine strains: WCV (J445, J446, J447 and J448), ACV (BP165) were compared with Tohama-I reference strain and recently reported clinical isolates from India (BPD1, BPD2). Core genome-based phylogenetic analysis was also performed using 166 isolates reported from countries using ACV.

RESULTS

Whole-genome analysis of vaccine and clinical isolates reported from India revealed high genetic similarity and conserved genome among strains. Phylogenetic analysis showed that clinical and vaccine strains share genetic closeness with reference strain Tohama-I. The allelic profile of vaccine strains (J445:ptxP1/ptxA2/prn1/fim2-1/fim3-1; J446: ptxP2/ptxA4/prn7/fim2-2/fim3-1; J447 and J448: ptxP1/ptxA1/ prn1/fim2-1/fim3-1), which matched entirely with clinical isolates (BPD1:ptxP1/ptxA1/prn1/fim2-1 and BPD2: ptxP1/ptxA1/prn1/fim2-1) reported from India. Multi-locus sequence typing (MLST) demonstrated the presence of dominant sequence types ST2 and primitive ST1 in vaccine strains which will allow better coverage against circulating strains of B. pertussis.

CONCLUSIONS

The study provides a detailed characterization of vaccine and clinical strains reported from India, which will further facilitate epidemiological studies on genetic shifts in countries which are using WCVs in their immunization programs.

Comparative genomics reveals different population structures associated with host and geographic origin in antimicrobial-resistant Salmonella enterica

Genetic variation in a pathogen, including the causative agent of salmonellosis, Salmonella enterica, can occur as a result of eco-evolutionary forces triggered by dissimilarities of ecological niches. Here, we applied comparative genomics to study 90 antimicrobial resistant (AMR) S. enterica isolates from bovine and human hosts in New York and Washington states to understand host- and geographic-associated population structure. Results revealed distinct presence/absence profiles of functional genes and pseudogenes (e.g., virulence genes) associated with bovine and human isolates. Notably, bovine isolates contained significantly more transposase genes but fewer transposase pseudogenes than human isolates, suggesting the occurrence of large-scale transposition in genomes of bovine and human isolates at different times. The high correlation between transposase genes and AMR genes, as well as plasmid replicons, highlights the potential role of horizontally transferred transposons in promoting adaptation to antibiotics. By contrast, a number of potentially geographic-associated single-nucleotide polymorphisms (SNPs), rather than geographic-associated genes, were identified. Interestingly, 38% of these SNPs were in genes annotated as cell surface protein-encoding genes, including some essential for antibiotic resistance and host colonization. Overall, different evolutionary forces and limited recent inter-population transmission appear to shape AMR S. enterica population structure in different hosts and geographic origins.

Cultivation-Assisted Genome of Candidatus Fukatsuia symbiotica; the Enigmatic "X-Type" Symbiont of Aphids

Heritable symbionts are common in terrestrial arthropods and often provide beneficial services to hosts. Unlike obligate, nutritional symbionts that largely persist under strict host control within specialized host cells, heritable facultative symbionts exhibit large variation in within-host lifestyles and services rendered with many retaining the capacity to transition among roles. One enigmatic symbiont, Candidatus Fukatsuia symbiotica, frequently infects aphids with reported roles ranging from pathogen, defensive symbiont, mutualism exploiter, and nutritional co-obligate symbiont. Here, we used an in vitro culture-assisted protocol to sequence the genome of a facultative strain of Fukatsuia from pea aphids (Acyrthosiphon pisum). Phylogenetic and genomic comparisons indicate that Fukatsuia is an aerobic heterotroph, which together with Regiella insecticola and Hamiltonella defensa form a clade of heritable facultative symbionts within the Yersiniaceae (Enterobacteriales). These three heritable facultative symbionts largely share overlapping inventories of genes associated with housekeeping functions, metabolism, and nutrient acquisition, while varying in complements of mobile DNA. One unusual feature of Fukatsuia is its strong tendency to occur as a coinfection with H. defensa. However, the overall similarity of gene inventories among aphid heritable facultative symbionts suggests that metabolic complementarity is not the basis for coinfection, unless playing out on a H. defensa strain-specific basis. We also compared the pea aphid Fukatsuia with a strain from the aphid Cinara confinis (Lachninae) where it is reported to have transitioned to co-obligate status to support decaying Buchnera function. Overall, the two genomes are very similar with no clear genomic signatures consistent with such a transition, which suggests co-obligate status in C. confinis was a recent event.

Aphids in focus: unravelling their complex ecology and evolution using genetic and molecular approaches

Aphids are renowned plant parasites of agriculture, horticulture and forestry, causing direct physical damage by sucking phloem and especially by transmission of plant pathogenic viruses. The huge yield loss they cause amounts to hundreds of millions of dollars globally, and because of this damage and the intense efforts expended on control, some 20 species are now resistant to pesticides worldwide. Aphids represent an ancient, mainly northern temperate group, although some species occur in the tropics, often as obligate asexual lineages or even asexual ‘species’. However, besides their notoriety as enemies of plant growers, aphids are also extremely interesting scientifically, especially at the molecular and genetic levels. They reproduce mainly asexually, one female producing 10–90 offspring in 7–10 days and therefore, theoretically, could produce billions of offspring in one growing season in the absence of mortality factors (i.e. climate/weather and antagonists). In this overview, we provide examples of what molecular and genetic studies of aphids have revealed concerning a range of topics, especially fine-grained ecological processes. Aphids, despite their apparently limited behavioural repertoire, are in fact masters (or, perhaps more accurately, mistresses) of adaptation and evolutionary flexibility and continue to flourish in a variety of ecosystems, including the agro-ecosystem, regardless of our best efforts to combat them.

Microsporidia with Vertical Transmission Were Likely Shaped by Nonadaptive Processes

Microsporidia have the leanest genomes among eukaryotes, and their physiological and genomic simplicity has been attributed to their intracellular, obligate parasitic life-style. However, not all microsporidia genomes are small or lean, with the largest dwarfing the smallest ones by at least an order of magnitude. To better understand the evolutionary mechanisms behind this genomic diversification, we explore here two clades of microsporidia with distinct life histories, Ordospora and Hamiltosporidium, parasitizing the same host species, Daphnia magna. Based on seven newly assembled genomes, we show that mixed-mode transmission (the combination of horizontal and vertical transmission), which occurs in Hamiltosporidium, is found to be associated with larger and AT-biased genomes, more genes, and longer intergenic regions, as compared with the exclusively horizontally transmitted Ordospora. Furthermore, the Hamiltosporidium genome assemblies contain a variety of repetitive elements and long segmental duplications. We show that there is an excess of nonsynonymous substitutions in the microsporidia with mixed-mode transmission, which cannot be solely attributed to the lack of recombination, suggesting that bursts of genome size in these microsporidia result primarily from genetic drift. Overall, these findings suggest that the switch from a horizontal-only to a mixed mode of transmission likely produces population bottlenecks in Hamiltosporidium species, therefore reducing the effectiveness of natural selection, and allowing their genomic features to be largely shaped by nonadaptive processes.

Mycobacterium leprae's evolution and environmental adaptation

Leprosy is an ancient disease caused by the acid-fast bacillus Mycobacterium leprae, also known as Hansen's bacillus. M. leprae is an obligate intracellular microorganism with a marked Schwann cell tropism and is the only human pathogen capable of invading the superficial peripheral nerves. The transmission mechanism of M. leprae is not fully understood; however, the nasal mucosa is accepted as main route of M. leprae entry to the human host. The complete sequencing and the comparative genome analysis show that M. leprae underwent a genome reductive evolution process, as result of lifestyle change and adaptation to different environments; some of lost genes are homologous to those of host cells. Thus, M. leprae reduced its genome size to 3.3 Mbp, contributing to obtain the lowest GC content (approximately 58%) among mycobacteria. The M. leprae genome contains 1614 open reading frames coding for functional proteins, and 1310 pseudogenes corresponding to 41% of the genome, approximately. Comparative analyses to different microorganisms showed that M. leprae possesses the highest content of pseudogenes among pathogenic and non-pathogenic bacteria and archaea. The pathogen adaptation into host cells, as the Schwann cells, brought about the reduction of the genome and induced multiple gene inactivation. The present review highlights the characteristics of genome's reductive evolution that M. leprae experiences in the genetic aspects compared with other pathogens. The possible mechanisms of pseudogenes formation are discussed.

The population structure of Ochrobactrum isolated from entomopathogenic nematodes indicates interactions with the symbiotic system

Entomopathogenic nematodes (EPNs) form specific mutualistic associations with bioluminescent enterobacteria. In Heterorhabditidis indica, Ochrobactrum spp. was identified beside the symbiont Photorhabdus luminescens but its involvement in the symbiotic association in the EPNs remains unclear. This study describe the population structure and the diversity in Ochrobactrum natural populations isolated from EPNs in the Caribbean basin in order to question the existence of EPN-specialized clones and to gain a better insight into Ochrobactrum-EPNs relationships. EPN-associated Ochrobactrum and Photorhabdus strains were characterized by multi-locus sequence typing, Pulsed-Field Gel Electrophoresis fingerprinting and phenotypic traits. Population study showed the absence of EPN-specialized clones in O. intermedium and O. anthropi but suggested the success of some particular lineages. A low level of genetic and genomic diversification of Ochrobactrum isolated from the natural population of Caribbean nematodes was observed comparatively to the diversity of human-associated Ochrobactrum strains. Correspondences between Ochrobactrum and P. luminescens PFGE clusters have been observed, particularly in the case of nematodes from Dominican Republic and Puerto Rico. O. intermedium and O. anthropi associated to EPNs formed less biofilm than human-associated strains. These results evoke interactions between Ochrobactrum and the EPN symbiotic system rather than transient contamination. The main hypothesis to investigate is a toxic/antitoxic relationship because of the ability of Ochrobactrum to resist to antimicrobial and toxic compounds produced by Photorhabdus.

Fueled by methane: deep-sea sponges from asphalt seeps gain their nutrition from methane-oxidizing symbionts

Sponges host a remarkable diversity of microbial symbionts, however, the benefit their microbes provide is rarely understood. Here, we describe two new sponge species from deep-sea asphalt seeps and show that they live in a nutritional symbiosis with methane-oxidizing (MOX) bacteria. Metagenomics and imaging analyses revealed unusually high amounts of MOX symbionts in hosts from a group previously assumed to have low microbial abundances. These symbionts belonged to the Marine Methylotrophic Group 2 clade. They are host-specific and likely vertically transmitted, based on their presence in sponge embryos and streamlined genomes, which lacked genes typical of related free-living MOX. Moreover, genes known to play a role in host-symbiont interactions, such as those that encode eukaryote-like proteins, were abundant and expressed. Methane assimilation by the symbionts was one of the most highly expressed metabolic pathways in the sponges. Molecular and stable carbon isotope patterns of lipids confirmed that methane-derived carbon was incorporated into the hosts. Our results revealed that two species of sponges, although distantly related, independently established highly specific, nutritional symbioses with two closely related methanotrophs. This convergence in symbiont acquisition underscores the strong selective advantage for these sponges in harboring MOX bacteria in the food-limited deep sea.

A New Species of the γ-Proteobacterium Francisella, F. adeliensis Sp. Nov., Endocytobiont in an Antarctic Marine Ciliate and Potential Evolutionary Forerunner of Pathogenic Species

The study of the draft genome of an Antarctic marine ciliate, Euplotes petzi, revealed foreign sequences of bacterial origin belonging to the γ-proteobacterium Francisella that includes pathogenic and environmental species. TEM and FISH analyses confirmed the presence of a Francisella endocytobiont in E. petzi. This endocytobiont was isolated and found to be a new species, named F. adeliensis sp. nov.. F. adeliensis grows well at wide ranges of temperature, salinity, and carbon dioxide concentrations implying that it may colonize new organisms living in deeply diversified habitats. The F. adeliensis genome includes the igl and pdp gene sets (pdpC and pdpE excepted) of the Francisella pathogenicity island needed for intracellular growth. Consistently with an F. adeliensis ancient symbiotic lifestyle, it also contains a single insertion-sequence element. Instead, it lacks genes for the biosynthesis of essential amino acids such as cysteine, lysine, methionine, and tyrosine. In a genome-based phylogenetic tree, F. adeliensis forms a new early branching clade, basal to the evolution of pathogenic species. The correlations of this clade with the other clades raise doubts about a genuine free-living nature of the environmental Francisella species isolated from natural and man-made environments, and suggest to look at F. adeliensis as a pioneer in the Francisella colonization of eukaryotic organisms.

Comparative Genome Analysis of an Extensively Drug-Resistant Isolate of Avian Sequence Type 167 Escherichia coli Strain Sanji with Novel In Silico Serotype O89b:H9

E. coli strain Sanji is the first sequenced and analyzed genome of the recently emerged pathogenic XDR strains with sequence type ST167 and novel in silico serotype O89b:H9. Comparison of the genomes of Sanji with other ST167 strains revealed distinct sets of different plasmids, mobile IS elements, and antibiotic resistance genes in each genome, indicating that there exist multiple paths toward achieving XDR. The emergence of these pathogenic ST167 E. coli strains with diverse XDR capabilities highlights the difficulty of preventing or mitigating the development of XDR properties in bacteria and points to the importance of better understanding of the shared underlying virulence mechanisms and physiology of pathogenic bacteria.

Extensive drug resistance (XDR) is an escalating global problem. Escherichia coli strain Sanji was isolated from an outbreak of pheasant colibacillosis in Fujian province, China, in 2011. This strain has XDR properties, exhibiting sensitivity to carbapenems but no other classes of known antibiotics. Whole-genome sequencing revealed a total of 32 known antibiotic resistance genes, many associated with insertion sequence 26 (IS26) elements. These were found on the Sanji chromosome and 2 of its 6 plasmids, pSJ_255 and pSJ_82. The Sanji chromosome also harbors a type 2 secretion system (T2SS), a type 3 secretion system (T3SS), a type 6 secretion system (T6SS), and several putative prophages. Sanji and other ST167 strains have a previously uncharacterized O-antigen (O89b) that is most closely related to serotype O89 as determined on the basis of analysis of the wzm-wzt genes and in silico serotyping. This O89b-antigen gene cluster was also found in the genomes of a few other pathogenic sequence type 617 (ST617) and ST10 complex strains. A time-scaled phylogeny inferred from comparative single nucleotide variant analysis indicated that development of these O89b-containing lineages emerged about 30 years ago. Comparative sequence analysis revealed that the core genome of Sanji is nearly identical to that of several recently sequenced strains of pathogenic XDR E. coli belonging to the ST167 group. Comparison of the mobile elements among the different ST167 genomes revealed that each genome carries a distinct set of multidrug resistance genes on different types of plasmids, indicating that there are multiple paths toward the emergence of XDR in E. coli.

IMPORTANCEE. coli strain Sanji is the first sequenced and analyzed genome of the recently emerged pathogenic XDR strains with sequence type ST167 and novel in silico serotype O89b:H9. Comparison of the genomes of Sanji with other ST167 strains revealed distinct sets of different plasmids, mobile IS elements, and antibiotic resistance genes in each genome, indicating that there exist multiple paths toward achieving XDR. The emergence of these pathogenic ST167 E. coli strains with diverse XDR capabilities highlights the difficulty of preventing or mitigating the development of XDR properties in bacteria and points to the importance of better understanding of the shared underlying virulence mechanisms and physiology of pathogenic bacteria.

Author Video: An author video summary of this article is available.

The Domestication of a Large DNA Virus by the Wasp Venturia canescens Involves Targeted Genome Reduction through Pseudogenization

Polydnaviruses (PDVs) are compelling examples of viral domestication, in which wasps express a large set of genes originating from a chromosomally integrated virus to produce particles necessary for their reproductive success. Parasitoid wasps generally use PDVs as a virulence gene delivery system allowing the protection of their progeny in the body of parasitized host. However, in the wasp Venturia canescens an independent viral domestication process led to an alternative strategy as the wasp incorporates virulence proteins in viral liposomes named virus-like particles (VLPs), instead of DNA molecules. Proteomic analysis of purified VLPs and transcriptome sequencing revealed the loss of some viral functions. In particular, the genes coding for capsid components are no longer expressed, which explains why VLPs do not incorporate DNA. Here a thorough examination of V. canescens genome revealed the presence of the pseudogenes corresponding to most of the genes involved in lost functions. This strongly suggests that an accumulation of mutations that leads to gene specific pseudogenization precedes the loss of viral genes observed during virus domestication. No evidence was found for block loss of collinear genes, although extensive gene order reshuffling of the viral genome was identified from comparisons between endogenous and exogenous viruses. These results provide the first insights on the early stages of large DNA virus domestication implicating massive genome reduction through gene-specific pseudogenization, a process which differs from the large deletions described for bacterial endosymbionts.