Genomic islands in pathogenic and environmental microorganisms

Transcriptional noise and exaptation as sources for bacterial sRNAs

Understanding how new genes originate and integrate into cellular networks is key to understanding evolution. Bacteria present unique opportunities for both the natural history and experimental study of gene origins, due to their large effective population sizes, rapid generation times, and ease of genetic manipulation. Bacterial small non-coding RNAs (sRNAs), in particular, many of which operate through a simple antisense regulatory logic, may serve as tractable models for exploring processes of gene origin and adaptation. Understanding how and on what timescales these regulatory molecules arise has important implications for understanding the evolution of bacterial regulatory networks, in particular, for the design of comparative studies of sRNA function. Here, we introduce relevant concepts from evolutionary biology and review recent work that has begun to shed light on the timescales and processes through which non-functional transcriptional noise is co-opted to provide regulatory functions. We explore possible scenarios for sRNA origin, focusing on the co-option, or exaptation, of existing genomic structures which may provide protected spaces for sRNA evolution.

Distribution and characteristics of SGI1/PGI2 genomic island from Proteus strains in China

The emergence of multidrug-resistant Salmonella genomic island 1 (SGI1) and Proteus genomic island (PGI) bearing P. mirabilis present a serious threat to public health. In this study, we screened 288 Proteus isolates recovered from seven provinces in China. Fourteen strains (4.9%) all belonged to P. mirabilis were positive for SGI1/PGI2, including twelve from clinical samples (5.3%) and two from food (3.3%). A Blastn search against GenBank and phylogenetic analyses identified eight different SGI1 variants and one PGI2 variant from the fourteen SGI1/PGI2 variants. All SGI1 variants shared a common backbone and harbored different resistance gene(s), except the sul1 gene at its multidrug-resistant (MDR) region. Among the variants, three novel SGI1 variants, designated as SGI1-PmCA11, SGI1-PmCA14 and SGI1-PmCA46, contained different gene cassettes, which were similar to sequences in plasmids or class 1 integrons of Klebsiella pneumoniae, P. mirabilis, Escherichia coli and Salmonella. Moreover, one novel PGI2, designated as PGI2-PmCA72, had an identical gene cassette to the first class 1 integron from PGI2 (GenBank accession no. MG201402.1) in P. mirabilis, but varied due to missing, replaced, inserted and inverted gene clusters. The four novel SGI1/PGI2 variants contained the cmlA5, dfrA14, bla_OXA-10, aadA15, bla_OXA-1, catB3 and dfrA16 resistance genes, which have never been reported in SGI1/PGI2 variants. Phenotypically, all fourteen SGI1/PGI2-containing strains showed multidrug resistance. All except four strains were resistant to the first, or the second and/or-third generation cephalosporins. Considering the increasing number and the emergence of new SGI1/PGI2 variants, further surveillance is needed to prevent the spreading of the MDR genomic islands among Proteus isolates from human and food.

Genomic Insights of Halophilic Planococcus maritimus SAMP MCC 3013 and Detail Investigation of Its Biosurfactant Production

Moderate halophilic bacteria thrive in saline conditions and produce biosurfactant (BS) which facilitates the oil scavenging activity in the oil polluted surroundings. Production of such unusual bioactive molecules plays a vital role for their survival in an extreme and adverse environment. Current research deals with isolation of Planococcus maritimus strain SAMP MCC 3013 from Indian Arabian coastline sea water for BS production. The bacterium tolerated up to 2.7 M NaCl demonstrating osmotic stress bearable physiological systems. We used integrated approach to explore the genomic insight of the strain SAMP and displayed the presence of gene for BS biosynthesis. The genome analysis revealed this potential to be intrinsic to the strain. Preliminary screening techniques viz., surface tension (SFT), drop collapse (DC) and oil displacement (OD) showed SAMP MCC 3013 as a potent BS producer. BS reduced SFT of phosphate buffer saline (PBS) pH: 7.0 from 72 to 30 mN/m with a critical micelle concentration (CMC) value of 1.3 mg/mL. Subsequent investigation on chemical characterization, using thin layer chromatography (TLC), Fourier transform infrared spectroscopy (FT-IR), nuclear magnetic resonance (1H NMR and 13C NMR) and liquid chromatography mass spectrometry (LC-MS) revealed terpene containing BS having sugar, lipid moieties. The genomic sequence analysis of P. maritimus SAMP showed complete genes in the pathway for the synthesis of terpenoid. Probably terpenoid is the accountable backbone molecule for the BS production, but the later stages of terpenoid conversion to the BS could not be found. Moreover, it is important to highlight that till today; no single report documents the in-detailed physico-chemical characterization of BS from Planococcus sp. Based on genomic and functional properties, the term terpene containing BS is denoted for the surfactant produced by P. maritimus.

An efficient genomic signature ranking method for genomic island prediction from a single genome

Genomic islands that are associated with microbial adaptations and carry genomic signatures different from that of the host, and thus many methods have been proposed to select the informative genomic signatures from a range of organisms and discriminate genomic islands from the rest of the genome in terms of these signature biases. However, they are of limited use when closely related genomes are unavailable. In the present work, we proposed a kurtosis-based ranking method to select the informative genomic signatures from a single genome. In simulations with alien fragments from artificial and real genomes, the proposed kurtosis-based ranking method efficiently selected the informative genomic signatures from a single genome, without annotated information of genomes or prior knowledge from other datasets. This understanding can be useful to design more powerful method for genomic island detection.

Virulence factor landscape of a Staphylococcus aureus sequence type 45 strain, MCRF184

Background

We describe the virulence factors of a methicillin-sensitive Staphylococcus aureus sequence type (ST) 45 strain, MCRF184, (spa type t917), that caused severe necrotizing fasciitis in a 72-year-old diabetic male. The genome of MCRF184 possesses three genomic islands: a relatively large type III νSaα with 42 open reading frames (ORFs) that includes superantigen- and lipoprotein-like genes, a truncated νSaβ that consists mostly of the enterotoxin gene cluster (egc), and a νSaγ island with 18 ORFs including α-toxin. Additionally, the genome has two phage-related regions: phage φSa3 with three genes of the immune evasion cluster (IEC), and an incomplete phage that is distinct from other S. aureus phages. Finally, the region between orfX and orfY harbors a putative efflux pump, acetyltransferase, regulators, and mobilization genes instead of genes of SCCmec.

Results

Virulence factors included phenol soluble modulins (PSMs) α1 through α4 and PSMs β1 and β2. Ten ORFs identified in MCRF184 had not been reported in previously sequenced S. aureus strains.

Conclusion

The dire clinical outcome in the patient and the described virulence factors all suggest that MCRF184, a ST45 strain is a highly virulent strain of S. aureus.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5394-2) contains supplementary material, which is available to authorized users.

The Type IV Secretion System of ICEAfe1: Formation of a Conjugative Pilus in Acidithiobacillus ferrooxidans

The dispersal of mobile genetic elements and their gene cargo relies on type IV secretion systems (T4SS). In this work the ICEAfe1 Tra-type T4SS nanomachine, encoded in the publicly available genome of Acidithiobacillus ferrooxidans ATCC 23270^TY, was characterized in terms of its organization, conservation, expression and mating bridge formation. Twenty-one conjugative genes grouped in four genetic clusters encode the ICEAfe1 T4SS, containing all the indispensable functions for the formation and stabilization of the pili and for DNA processing. The clusters' organization resembles that of other mobile genetic elements (such as plasmids and integrative and conjugative elements-ICEs). Sequence conservation, genetic organization and distribution of the tra system in the genomes of other sequenced Acidithiobacillus spp. suggests that the ICEAfe1 T4SS could mediate the lateral gene transfer between related bacteria. All ICEAfe1 T4SS genes are transcriptionally active and expressed from four independent operons. The transcriptional levels of selected marker genes increase in response to Mitomycin C treatment, a DNA damage elicitor that has acknowledged stimulatory effects on excision rates and gene expression of other ICEs, including ICEAfe1. Using a tailor-made pilin-antiserum against ICEAfe1 T4SS TraA pilin and epifluorescence microscopy, the presence of the conjugative pili on the cell surface of A. ferrooxidans could be demonstrated. Additionally, immunodetection assays, by immunogold, allowed the identification of pili-like extracellular structures. Together, the results obtained in this work demonstrate that the ICEAfe1 T4SS is phylogenetically conserved within the taxon, is expressed at mRNA and protein levels in vivo in the A. ferrooxidans type strain, and produces a pili-like structure of extracellular and intercellular localization in this model acidophile, supporting its functionality. Additional efforts will be required to prove conjugation of the ICEAfe1 or parts of this element through the cognate T4SS.

Novel large-scale chromosomal transfer in Bacteroides fragilis contributes to its pan-genome and rapid environmental adaptation

Bacteroides fragilis, an important component of the human gastrointestinal microbiota, can cause lethal extra-intestinal infection upon escape from the gastrointestinal tract. We demonstrated transfer and recombination of large chromosomal segments from B. fragilis HMW615, a multidrug resistant clinical isolate, to B. fragilis 638R. In one example, the transfer of a segment of ~435 Kb/356 genes replaced ~413 Kb/326 genes of the B. fragilis 638R chromosome. In addition to transfer of antibiotic resistance genes, these transfers (1) replaced complete divergent polysaccharide biosynthesis loci; (2) replaced DNA inversion-controlled intergenic shufflons (that control expression of genes encoding starch utilization system outer membrane proteins) with more complex, divergent shufflons; and (3) introduced additional intergenic shufflons encoding divergent Type 1 restriction/modification systems. Conjugative transposon-like genes within a transferred segment and within a putative integrative conjugative element (ICE5) ~45 kb downstream from the transferred segment both encode proteins that may be involved in the observed transfer. These data indicate that chromosomal transfer is a driver of antigenic diversity and nutrient adaptation in Bacteroides that (1) contributes to the dissemination of the extensive B. fragilis pan-genome, (2) allows rapid adaptation to a changing environment and (3) can confer pathogenic characteristics to host symbionts.

Engineering root microbiomes for healthier crops and soils using beneficial, environmentally safe bacteria

The Green Revolution developed new crop varieties, which greatly improved food security worldwide. However, the growth of these plants relied heavily on chemical fertilizers and pesticides, which have led to an overuse of synthetic fertilizers, insecticides, and herbicides with serious environmental consequences and negative effects on human health. Environmentally friendly plant-growth-promoting methods to replace our current reliance on synthetic chemicals and to develop more sustainable agricultural practices to offset the damage caused by many agrochemicals are proposed herein. The increased use of bioinoculants, which consist of microorganisms that establish synergies with target crops and influence production and yield by enhancing plant growth, controlling disease, and providing critical mineral nutrients, is a potential solution. The microorganisms found in bioinoculants are often bacteria or fungi that reside within either external or internal plant microbiomes. However, before they can be used routinely in agriculture, these microbes must be confirmed as nonpathogenic strains that promote plant growth and survival. In this article, besides describing approaches for discovering plant-growth-promoting bacteria in various environments, including phytomicrobiomes and soils, we also discuss methods to evaluate their safety for the environment and for human health.

Impact of Homologous Recombination on the Evolution of Prokaryotic Core Genomes

Microbial populations exchange genetic material through a process called homologous recombination. Although this process has been studied in particular organisms, we lack an understanding of its differential impact over the genome and across microbes with different life-styles. We used a common analytical framework to assess this process in a representative set of microorganisms. Our results uncovered important trends. First, microbes with different lifestyles are differentially impacted, with endosymbionts and obligate pathogens being those less prone to undergo this process. Second, certain genetic elements such as restriction-modification systems seem to be associated with higher rates of recombination. Most importantly, recombined genomes show the footprints of natural selection in which recombined regions preferentially contain genes that can be related to specific ecological adaptations. Taken together, our results clarify the relative contributions of factors modulating homologous recombination and show evidence for a clear a role of this process in shaping microbial genomes and driving ecological adaptations.

Homologous recombination (HR) enables the exchange of genetic material between and within species. Recent studies suggest that this process plays a major role in the microevolution of microbial genomes, contributing to core genome homogenization and to the maintenance of cohesive population structures. However, we still have a very poor understanding of the possible adaptive roles of intraspecific HR and of the factors that determine its differential impact across clades and lifestyles. Here we used a unified methodological framework to assess HR in 338 complete genomes from 54 phylogenetically diverse and representative prokaryotic species, encompassing different lifestyles and a broad phylogenetic distribution. Our results indicate that lifestyle and presence of restriction-modification (RM) machineries are among the main factors shaping HR patterns, with symbionts and intracellular pathogens having the lowest HR levels. Similarly, the size of exchanged genomic fragments correlated with the presence of RM and competence machineries. Finally, genes exchanged by HR showed functional enrichments which could be related to adaptations to different environments and ecological strategies. Taken together, our results clarify the factors underlying HR impact and suggest important adaptive roles of genes exchanged through this mechanism. Our results also revealed that the extent of genetic exchange correlated with lifestyle and some genomic features. Moreover, the genes in exchanged regions were enriched for functions that reflected specific adaptations, supporting identification of HR as one of the main evolutionary mechanisms shaping prokaryotic core genomes.

Comparative analysis reveals the Genomic Islands in Pasteurella multocida population genetics: on Symbiosis and adaptability

BACKGROUND

Pasteurella multocida (P. multocida) is a widespread opportunistic pathogen that infects human and various animals. Genomic Islands (GIs) are one of the most important mobile components that quickly help bacteria acquire large fragments of foreign genes. However, the effects of GIs on P. multocida are unknown in the evolution of bacterial populations.

RESULTS

Ten avian-sourced P. multocida obtained through high-throughput sequencing together with 104 publicly available P. multocida genomes were used to analyse their population genetics, thus constructed a pan-genome containing 3948 protein-coding genes. Through the pan-genome, the open evolutionary pattern of P. multocida was revealed, and the functional components of 944 core genes, 2439 accessory genes and 565 unique genes were analysed. In addition, a total of 280 GIs were predicted in all strains. Combined with the pan-genome of P. multocida, the GIs accounted for 5.8% of the core genes in the pan-genome, mainly related to functional metabolic activities; the accessory genes accounted for 42.3%, mainly for the enrichment of adaptive genes; and the unique genes accounted for 35.4%, containing some defence mechanism-related genes.

CONCLUSIONS

The effects of GIs on the population genetics of P. multocida evolution and adaptation to the environment are reflected by the proportion and function of the pan-genome acquired from GIs, and the large quantities of GI data will aid in additional population genetics studies.

Potentially Mobile Denitrification Genes Identified in Azospirillum sp. Strain TSH58

Denitrification ability is sporadically distributed among diverse bacteria, archaea, and fungi. In addition, disagreement has been found between denitrification gene phylogenies and the 16S rRNA gene phylogeny. These facts have suggested potential occurrences of horizontal gene transfer (HGT) for the denitrification genes. However, evidence of HGT has not been clearly presented thus far. In this study, we identified the sequences and the localization of the nitrite reductase genes in the genomes of 41 denitrifying Azospirillum sp. strains and searched for mobile genetic elements that contain denitrification genes. All Azospirillum sp. strains examined in this study possessed multiple replicons (4 to 11 replicons), with their sizes ranging from 7 to 1,031 kbp. Among those, the nitrite reductase gene nirK was located on large replicons (549 to 941 kbp). Genome sequencing showed that Azospirillum strains that had similar nirK sequences also shared similar nir-nor gene arrangements, especially between the TSH58, Sp7T, and Sp245 strains. In addition to the high similarity between nir-nor gene clusters among the three Azospirillum strains, a composite transposon structure was identified in the genome of strain TSH58, which contains the nir-nor gene cluster and the novel IS6 family insertion sequences (ISAz581 and ISAz582). The nirK gene within the composite transposon system was actively transcribed under denitrification-inducing conditions. Although not experimentally verified in this study, the composite transposon system containing the nir-nor gene cluster could be transferred to other cells if it is moved to a prophage region and the phage becomes activated and released outside the cells. Taken together, strain TSH58 most likely acquired its denitrification ability by HGT from closely related Azospirillum sp. denitrifiers.IMPORTANCE The evolutionary history of denitrification is complex. While the occurrence of horizontal gene transfer has been suggested for denitrification genes, most studies report circumstantial evidences, such as disagreement between denitrification gene phylogenies and the 16S rRNA gene phylogeny. Based on the comparative genome analyses of Azospirillum sp. denitrifiers, we identified denitrification genes, including nirK and norCBQD, located on a mobile genetic element in the genome of Azospirillum sp. strain TSH58. The nirK was actively transcribed under denitrification-inducing conditions. Since this gene was the sole nitrite reductase gene in strain TSH58, this strain most likely benefitted by acquiring denitrification genes via horizontal gene transfer. This finding will significantly advance our scientific knowledge regarding the ecology and evolution of denitrification.

Evolutionary dynamics of origin and loss in the deep history of phospholipase D toxin genes

Background

Venom-expressed sphingomyelinase D/phospholipase D (SMase D/PLD) enzymes evolved from the ubiquitous glycerophosphoryl diester phosphodiesterases (GDPD). Expression of GDPD-like SMaseD/PLD toxins in both arachnids and bacteria has inspired consideration of the relative contributions of lateral gene transfer and convergent recruitment in the evolutionary history of this lineage. Previous work recognized two distinct lineages, a SicTox-like (ST-like) clade including the arachnid toxins, and an Actinobacterial-toxin like (AT-like) clade including the bacterial toxins and numerous fungal homologs.

Results

Here we expand taxon sampling by homology detection to discover new GDPD-like SMase D/PLD homologs. The ST-like clade now includes homologs in a wider variety of arthropods along with a sister group in Cnidaria; the AT-like clade now includes additional fungal phyla and proteobacterial homologs; and we report a third clade expressed in diverse aquatic metazoan taxa, a few single-celled eukaryotes, and a few aquatic proteobacteria. GDPD-like SMaseD/PLDs have an ancient presence in chelicerates within the ST-like family and ctenophores within the Aquatic family. A rooted phylogenetic tree shows that the three clades derived from a basal paraphyletic group of proteobacterial GDPD-like SMase D/PLDs, some of which are on mobile genetic elements. GDPD-like SMase D/PLDs share a signature C-terminal motif and a shortened βα1 loop, features that distinguish them from GDPDs. The three major clades also have active site loop signatures that distinguish them from GDPDs and from each other. Analysis of molecular phylogenies with respect to organismal relationships reveals a dynamic evolutionary history including both lateral gene transfer and gene duplication/loss.

Conclusions

The GDPD-like SMaseD/PLD enzymes derive from a single ancient ancestor, likely proteobacterial, and radiated into diverse organismal lineages at least in part through lateral gene transfer.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1302-2) contains supplementary material, which is available to authorized users.

Occurrence and significance of pathogenicity and fitness islands in environmental vibrios

Pathogenicity islands (PAIs) are large genomic regions that contain virulence genes, which aid pathogens in establishing infections. While PAIs in clinical strains (strains isolated from a human infection) are well-studied, less is known about the occurrence of PAIs in strains isolated from the environment. In this study we describe three PAIs found in environmental Vibrio vulnificus and Vibrio parahaemolyticus strains, as well as a genomic fitness island found in a Vibrio diabolicus strain. All four islands had markedly different GC profiles than the rest of the genome, indicating that all of these islands were acquired via lateral gene transfer. Genes on the PAIs and fitness island were characterized. The PAI found in V. parahaemolyticus contained the tdh gene, a collagenase gene, and genes involved in the type 3 secretion system II (T3SS2). A V. vulnificus environmental strain contained two PAIs, a small 25 kbp PAI and a larger 143 kbp PAI. Both PAIs contained virulence genes. Toxin-antitoxin (TA) genes were found in all three species: on the V. diabolicus fitness island, and on the V. parahaemolyticus and V. vulnificus PAIs.

How Can We Define "Optimal Microbiota?": A Comparative Review of Structure and Functions of Microbiota of Animals, Fish, and Plants in Agriculture

All multicellular organisms benefit from their own microbiota, which play important roles in maintaining the host nutritional health and immunity. Recently, the number of studies on the microbiota of animals, fish, and plants of economic importance is rapidly expanding and there are increasing expectations that productivity and sustainability in agricultural management can be improved by microbiota manipulation. However, optimizing microbiota is still a challenging task because of the lack of knowledge on the dominant microorganisms or significant variations between microbiota, reflecting sampling biases, different agricultural management as well as breeding backgrounds. To offer a more generalized view on microbiota in agriculture, which can be used for defining criteria of “optimal microbiota” as the goal of manipulation, we summarize here current knowledge on microbiota on animals, fish, and plants with emphasis on bacterial community structure and metabolic functions, and how microbiota can be affected by domestication, conventional agricultural practices, and use of antimicrobial agents. Finally, we discuss future tasks for defining “optimal microbiota,” which can improve host growth, nutrition, and immunity and reduce the use of antimicrobial agents in agriculture.

Comparative genome analysis reveals key genetic factors associated with probiotic property in Enterococcus faecium strains

BACKGROUND

Enterococcus faecium though commensal in the human gut, few strains provide a beneficial effect to humans as probiotics while few are responsible for the nosocomial infection. Comparative genomics of E. faecium can decipher the genomic differences responsible for probiotic, pathogenic and non-pathogenic properties. In this study, we compared E. faecium strain 17OM39 with a marketed probiotic, non-pathogenic non-probiotic (NPNP) and pathogenic strains.

RESULTS

E. faecium 17OM39 was found to be closely related with marketed probiotic strain T110 based on core genome analysis. Strain 17OM39 was devoid of known vancomycin, tetracycline resistance and functional virulence genes. Moreover, E. faecium 17OM39 genome was found to be more stable due to the absence of frequently found transposable elements. Genes imparting beneficial functional properties were observed to be present in marketed probiotic T110 and 17OM39 strains. Genes associated with colonization and survival within gastrointestinal tract was also detected across all the strains.

CONCLUSIONS

Beyond shared genetic features; this study particularly identified genes that are responsible for imparting probiotic, non-pathogenic and pathogenic features to the strains of E. faecium. Higher genomic stability, absence of known virulence factors and antibiotic resistance genes and close genomic relatedness with marketed probiotics makes E. faecium 17OM39 a potential probiotic candidate. The work presented here demonstrates that comparative genome analyses can be applied to large numbers of genomes, to find potential probiotic candidates.

Photocatalytic Reactive Ultrafiltration Membrane for Removal of Antibiotic Resistant Bacteria and Antibiotic Resistance Genes from Wastewater Effluent

Biological wastewater treatment is not effective in removal of antibiotic resistant bacteria (ARB) and antibiotic resistance genes (ARGs). In this study, we fabricated a photocatalytic reactive membrane by functionalizing polyvinylidene fluoride (PVDF) ultrafiltration (UF) membrane with titanium oxide (TiO₂) nanoparticles for the removal of ARB and ARGs from a secondary wastewater effluent. The TiO₂-modified PVDF membrane provided complete retention of ARB and effective photocatalytic degradation of ARGs and integrons. Specifically, the total removal efficiency of ARGs (i.e., plasmid-mediated floR, sul1, and sul2) with TiO₂-modified PVDF membrane reached ∼98% after exposure to UV irradiation. Photocatalytic degradation of ARGs located in the genome was found to be more efficient than those located in plasmid. Excellent removal of integrons (i.e., intI1, intI2, and intI3) after UV treatment indicated that the horizontal transfer potential of ARGs was effectively controlled by the TiO₂ photocatalytic reaction. We also evaluated the antifouling properties of the TiO₂-UF membrane to demonstrate its potential application in wastewater treatment.

Mutation Landscape of Base Substitutions, Duplications, and Deletions in the Representative Current Cholera Pandemic Strain

Pandemic cholera is a major concern for public health because of its high mortality and morbidity. Mutation accumulation (MA) experiments were performed on a representative strain of the current cholera pandemic. Although the base-pair substitution mutation rates in Vibrio cholerae (1.24 × 10-10 per site per generation for wild-type lines and 3.29 × 10-8 for mismatch repair deficient lines) are lower than that previously reported in other bacteria using MA analysis, we discovered specific high rates (8.31 × 10-8 site/generation for wild-type lines and 1.82 × 10-6 for mismatch repair deficient lines) of base duplication or deletion driven by large-scale copy number variations (CNVs). These duplication-deletions are located in two pathogenic islands, IMEX and the large integron island. Each element of these islands has discrepant rate in rapid integration and excision, which provides clues to the pandemicity evolution of V. cholerae. These results also suggest that large-scale structural variants such as CNVs can accumulate rapidly during short-term evolution. Mismatch repair deficient lines exhibit a significantly increased mutation rate in the larger chromosome (Chr1) at specific regions, and this pattern is not observed in wild-type lines. We propose that the high frequency of GATC sites in Chr1 improves the efficiency of MMR, resulting in similar rates of mutation in the wild-type condition. In addition, different mutation rates and spectra were observed in the MA lines under distinct growth conditions, including minimal media, rich media and antibiotic treatments.

Genome Variation in the Model Halophilic Bacterium Salinibacter ruber

The halophilic bacterium Salinibacter ruber is an abundant and ecologically important member of halophilic communities worldwide. Given its broad distribution and high intraspecific genetic diversity, S. ruber is considered one of the main models for ecological and evolutionary studies of bacterial adaptation to hypersaline environments. However, current insights on the genomic diversity of this species is limited to the comparison of the genomes of two co-isolated strains. Here, we present a comparative genomic analysis of eight S. ruber strains isolated at two different time points in each of two different Mediterranean solar salterns. Our results show an open pangenome with contrasting evolutionary patterns in the core and accessory genomes. We found that the core genome is shaped by extensive homologous recombination (HR), which results in limited sequence variation within population clusters. In contrast, the accessory genome is modulated by horizontal gene transfer (HGT), with genomic islands and plasmids acting as gateways to the rest of the genome. In addition, both types of genetic exchange are modulated by restriction and modification (RM) or CRISPR-Cas systems. Finally, genes differentially impacted by such processes reveal functional processes potentially relevant for environmental interactions and adaptation to extremophilic conditions. Altogether, our results support scenarios that conciliate “Neutral” and “Constant Diversity” models of bacterial evolution.

Comparative and phylogenetic analysis of a novel family of Enterobacteriaceae-associated genomic islands that share a conserved excision/integration module

Genomic Islands (GIs) are DNA regions acquired through horizontal gene transfer that encode advantageous traits for bacteria. Many GIs harbor genes that encode the molecular machinery required for their excision from the bacterial chromosome. Notably, the excision/integration dynamics of GIs may modulate the virulence of some pathogens. Here, we report a novel family of GIs found in plant and animal Enterobacteriaceae pathogens that share genes with those found in ROD21, a pathogenicity island whose excision is involved in the virulence of Salmonella enterica serovar Enteritidis. In these GIs we identified a conserved set of genes that includes an excision/integration module, suggesting that they are excisable. Indeed, we found that GIs within carbapenem-resistant Klebsiella pneumoniae ST258 KP35 and enteropathogenic Escherichia coli O127:H6 E2348/69 are excised from the bacterial genome. In addition to putative virulence factors, these GIs encode conjugative transfer-related proteins and short and full-length homologues of the global transcriptional regulator H-NS. Phylogenetic analyses suggest that the identified GIs likely originated in phytopathogenic bacteria. Taken together, our findings indicate that these GIs are excisable and may play a role in bacterial interactions with their hosts.

GI-Cluster: Detecting genomic islands via consensus clustering on multiple features

The accurate detection of genomic islands (GIs) in microbial genomes is important for both evolutionary study and medical research, because GIs may promote genome evolution and contain genes involved in pathogenesis. Various computational methods have been developed to predict GIs over the years. However, most of them cannot make full use of GI-associated features to achieve desirable performance. Additionally, many methods cannot be directly applied to newly sequenced genomes. We develop a new method called GI-Cluster, which provides an effective way to integrate multiple GI-related features via consensus clustering. GI-Cluster does not require training datasets or existing genome annotations, but it can still achieve comparable or better performance than supervised learning methods in comprehensive evaluations. Moreover, GI-Cluster is widely applicable, either to complete and incomplete genomes or to initial GI predictions from other programs. GI-Cluster also provides plots to visualize the distribution of predicted GIs and related features. GI-Cluster is available at https://github.com/icelu/GI_Cluster.

An antifungal polyketide associated with horizontally acquired genes supports symbiont-mediated defense in Lagria villosa beetles

Microbial symbionts are often a source of chemical novelty and can contribute to host defense against antagonists. However, the ecological relevance of chemical mediators remains unclear for most systems. Lagria beetles live in symbiosis with multiple strains of Burkholderia bacteria that protect their offspring against pathogens. Here, we describe the antifungal polyketide lagriamide, and provide evidence supporting that it is produced by an uncultured symbiont, Burkholderia gladioli Lv-StB, which is dominant in field-collected Lagria villosa. Interestingly, lagriamide is structurally similar to bistramides, defensive compounds found in marine tunicates. We identify a gene cluster that is probably involved in lagriamide biosynthesis, provide evidence for horizontal acquisition of these genes, and show that the naturally occurring symbiont strains on the egg are protective in the soil environment. Our findings highlight the potential of microbial symbionts and horizontal gene transfer as influential sources of ecological innovation.

Streptococcus pneumoniae's Virulence and Host Immunity: Aging, Diagnostics, and Prevention

Streptococcus pneumoniae is an infectious pathogen responsible for millions of deaths worldwide. Diseases caused by this bacterium are classified as pneumococcal diseases. This pathogen colonizes the nasopharynx of its host asymptomatically, but overtime can migrate to sterile tissues and organs and cause infections. Pneumonia is currently the most common pneumococcal disease. Pneumococcal pneumonia is a global health concern and vastly affects children under the age of five as well as the elderly and individuals with pre-existing health conditions. S. pneumoniae has a large selection of virulence factors that promote adherence, invasion of host tissues, and allows it to escape host immune defenses. A clear understanding of S. pneumoniae's virulence factors, host immune responses, and examining the current techniques available for diagnosis, treatment, and disease prevention will allow for better regulation of the pathogen and its diseases. In terms of disease prevention, other considerations must include the effects of age on responses to vaccines and vaccine efficacy. Ongoing work aims to improve on current vaccination paradigms by including the use of serotype-independent vaccines, such as protein and whole cell vaccines. Extending our knowledge of the biology of, and associated host immune response to S. pneumoniae is paramount for our improvement of pneumococcal disease diagnosis, treatment, and improvement of patient outlook.

Chalkophores

Copper-binding metallophores, or chalkophores, play a role in microbial copper homeostasis that is analogous to that of siderophores in iron homeostasis. The best-studied chalkophores are members of the methanobactin (Mbn) family-ribosomally produced, posttranslationally modified natural products first identified as copper chelators responsible for copper uptake in methane-oxidizing bacteria. To date, Mbns have been characterized exclusively in those species, but there is genomic evidence for their production in a much wider range of bacteria. This review addresses the current state of knowledge regarding the function, biosynthesis, transport, and regulation of Mbns. While the roles of several proteins in these processes are supported by substantial genetic and biochemical evidence, key aspects of Mbn manufacture, handling, and regulation remain unclear. In addition, other natural products that have been proposed to mediate copper uptake as well as metallophores that have biologically relevant roles involving copper binding, but not copper uptake, are discussed.

Genus-wide comparison of Pseudovibrio bacterial genomes reveal diverse adaptations to different marine invertebrate hosts

Bacteria belonging to the genus Pseudovibrio have been frequently found in association with a wide variety of marine eukaryotic invertebrate hosts, indicative of their versatile and symbiotic lifestyle. A recent comparison of the sponge-associated Pseudovibrio genomes has shed light on the mechanisms influencing a successful symbiotic association with sponges. In contrast, the genomic architecture of Pseudovibrio bacteria associated with other marine hosts has received less attention. Here, we performed genus-wide comparative analyses of 18 Pseudovibrio isolated from sponges, coral, tunicates, flatworm, and seawater. The analyses revealed a certain degree of commonality among the majority of sponge- and coral-associated bacteria. Isolates from other marine invertebrate host, tunicates, exhibited a genetic repertoire for cold adaptation and specific metabolic abilities including mucin degradation in the Antarctic tunicate-associated bacterium Pseudovibrio sp. Tun.PHSC04_5.I4. Reductive genome evolution was simultaneously detected in the flatworm-associated bacteria and the sponge-associated bacterium P. axinellae AD2, through the loss of major secretion systems (type III/VI) and virulence/symbioses factors such as proteins involved in adhesion and attachment to the host. Our study also unraveled the presence of a CRISPR-Cas system in P. stylochi UST20140214-052 a flatworm-associated bacterium possibly suggesting the role of CRISPR-based adaptive immune system against the invading virus particles. Detection of mobile elements and genomic islands (GIs) in all bacterial members highlighted the role of horizontal gene transfer for the acquisition of novel genetic features, likely enhancing the bacterial ecological fitness. These findings are insightful to understand the role of genome diversity in Pseudovibrio as an evolutionary strategy to increase their colonizing success across a wide range of marine eukaryotic hosts.

Mobile genetic elements in Neisseria gonorrhoeae: movement for change

Neisseria gonorrhoeae, the causative agent of the sexually transmitted disease gonorrhoeae, possesses several mobile genetic elements (MGEs). The MGEs such as transposable elements mediate intrachromosomal rearrangements, while plasmids and the gonococcal genetic island are involved in interchromosomal gene transfer. Additionally, gonococcal MGEs serve as hotspots for recombination and integration of other genetic elements such as bacteriophages, contribute to gene regulation or spread genes through gonococcal populations by horizontal gene transfer. In this review, we summarise the literature on the structure and biology of MGEs and discuss how these genetic elements may play a role in the pathogenesis and spread of antimicrobial resistance in N. gonorrhoeae. Although an abundance of information about gonococcal MGEs exists (mainly from whole genome sequencing and bioinformatic analysis), there are still many open questions on how MGEs influence the biology of N. gonorrhoeae.

Male-killing toxin in a bacterial symbiont of Drosophila

Several lineages of symbiotic bacteria in insects selfishly manipulate host reproduction to spread in a population ¹ , often by distorting host sex ratios. Spiroplasma poulsonii^2,3 is a helical and motile, Gram-positive symbiotic bacterium that resides in a wide range of Drosophila species ⁴ . A notable feature of S. poulsonii is male killing, whereby the sons of infected female hosts are selectively killed during development^1,2. Although male killing caused by S. poulsonii has been studied since the 1950s, its underlying mechanism is unknown. Here we identify an S. poulsonii protein, designated Spaid, whose expression induces male killing. Overexpression of Spaid in D. melanogaster kills males but not females, and induces massive apoptosis and neural defects, recapitulating the pathology observed in S. poulsonii-infected male embryos^5-11. Our data suggest that Spaid targets the dosage compensation machinery on the male X chromosome to mediate its effects. Spaid contains ankyrin repeats and a deubiquitinase domain, which are required for its subcellular localization and activity. Moreover, we found a laboratory mutant strain of S. poulsonii with reduced male-killing ability and a large deletion in the spaid locus. Our study has uncovered a bacterial protein that affects host cellular machinery in a sex-specific way, which is likely to be the long-searched-for factor responsible for S. poulsonii-induced male killing.

Coalescent framework for prokaryotes undergoing interspecific homologous recombination

Coalescent process for prokaryote species is theoretically considered. Prokaryotes undergo homologous recombination with individuals of the same species (intraspecific recombination) and with individuals of other species (interspecific recombination). This work particularly focuses on interspecific recombination because intraspecific recombination has been well incorporated in coalescent framework. We present a simulation framework for generating SNP (single-nucleotide polymorphism) patterns that allows external DNA integration into host genome from other species. Using this simulation tool, msPro, we observed that the joint processes of intra- and interspecific recombination generate complex SNP patterns. The direct effect of interspecific recombination includes increased polymorphism. Because interspecific recombination is very rare in nature, it generates regions with exceptionally high polymorphism. Following interspecific recombination, intraspecific recombination cuts the integrated external DNA into small fragments, generating a complex SNP pattern that appears as if external DNA was integrated multiple times. The insight gained from our work using the msPro simulator will be useful for understanding and evaluating the relative contributions of intra- and interspecific recombination events in generating complex SNP patters in prokaryotes.

PGI2 Is a Novel SGI1-Relative Multidrug-Resistant Genomic Island Characterized in Proteus mirabilis

A novel 61,578-bp genomic island named Proteus genomic island 2 (PGI2) was characterized in Proteus mirabilis of swine origin in China. The 23.85-kb backbone of PGI2 is related to those of Salmonella genomic island 1 and Acinetobacter genomic island 1. The multidrug resistance (MDR) region of PGI2 is a complex class 1 integron containing 14 different resistance genes. PGI2 was conjugally mobilized in trans to Escherichia coli in the presence of a conjugative IncC helper plasmid.

Physiological and transcriptome changes induced by Pseudomonas putida acquisition of an integrative and conjugative element

Integrative and conjugative elements (ICEs) comprise ubiquitous large mobile regions in prokaryotic chromosomes that transmit vertically to daughter cells and transfer horizontally to distantly related lineages. Their evolutionary success originates in maximized combined ICE-host fitness trade-offs, but how the ICE impacts on the host metabolism and physiology is poorly understood. Here we investigate global changes in the host genetic network and physiology of Pseudomonas putida with or without an integrated ICEclc, a model ICE widely distributed in proteobacterial genomes. Genome-wide gene expression differences were analyzed by RNA-seq using exponentially growing or stationary phase-restimulated cultures on 3-chlorobenzoate, an aromatic compound metabolizable thanks to specific ICEclc-located genes. We found that the presence of ICEclc imposes a variety of changes in global pathways such as cell cycle and amino acid metabolism, which were more numerous in stationary-restimulated than exponential phase cells. Unexpectedly, ICEclc stimulates cellular motility and leads to more rapid growth on 3-chlorobenzoate than cells carrying only the integrated clc genes. ICEclc also concomitantly activates the P. putida Pspu28-prophage, but this in itself did not provoke measurable fitness effects. ICEclc thus interferes in a number of cellular pathways, inducing both direct benefits as well as indirect costs in P. putida.

Sequential induction of three recombination directionality factors directs assembly of tripartite integrative and conjugative elements

Tripartite integrative and conjugative elements (ICE³) are a novel form of ICE that exist as three separate DNA regions integrated within the genomes of Mesorhizobium spp. Prior to conjugative transfer the three ICE³ regions of M. ciceri WSM1271 ICEMcSym¹²⁷¹ combine and excise to form a single circular element. This assembly requires three coordinated recombination events involving three site-specific recombinases IntS, IntG and IntM. Here, we demonstrate that three excisionases–or recombination directionality factors—RdfS, RdfG and RdfM are required for ICE³ excision. Transcriptome sequencing revealed that expression of ICE³ transfer and conjugation genes was induced by quorum sensing. Quorum sensing activated expression of rdfS, and in turn RdfS stimulated transcription of both rdfG and rdfM. Therefore, RdfS acts as a “master controller” of ICE³ assembly and excision. The dependence of all three excisive reactions on RdfS ensures that ICE³ excision occurs via a stepwise sequence of recombination events that avoids splitting the chromosome into a non-viable configuration. These discoveries expose a surprisingly simple control system guiding molecular assembly of these novel and complex mobile genetic elements and highlight the diverse and critical functions of excisionase proteins in control of horizontal gene transfer.

Bacteria evolve and adapt quickly through the horizontal transfer of DNA. A major mechanism facilitating this transfer is conjugation. Conjugative DNA elements that integrate into the chromosome are termed ‘Integrative and Conjugative Elements’ (ICE). We recently discovered a unique form of ICE that undergoes a complex series of recombination events with the host chromosome to split itself into three separate parts. This tripartite ICE must also precisely order its recombination when leaving the current host to avoid splitting the host chromosome and the ICE into non-viable parts. In this work, we show that the tripartite ICEs use chemical cell-cell communication to stimulate recombination and that recombination events are specifically ordered through cascaded transcriptional activation of small DNA-binding proteins called recombination directionality factors. Despite the inherent complexity of tripartite ICEs this work exposes a surprisingly simple system to stimulate their precise and ordered molecular assembly prior to horizontal transfer.

xenoGI: reconstructing the history of genomic island insertions in clades of closely related bacteria

Background

Genomic islands play an important role in microbial genome evolution, providing a mechanism for strains to adapt to new ecological conditions. A variety of computational methods, both genome-composition based and comparative, have been developed to identify them. Some of these methods are explicitly designed to work in single strains, while others make use of multiple strains. In general, existing methods do not identify islands in the context of the phylogeny in which they evolved. Even multiple strain approaches are best suited to identifying genomic islands that are present in one strain but absent in others. They do not automatically recognize islands which are shared between some strains in the clade or determine the branch on which these islands inserted within the phylogenetic tree.

Results

We have developed a software package, xenoGI, that identifies genomic islands and maps their origin within a clade of closely related bacteria, determining which branch they inserted on. It takes as input a set of sequenced genomes and a tree specifying their phylogenetic relationships. Making heavy use of synteny information, the package builds gene families in a species-tree-aware way, and then attempts to combine into islands those families whose members are adjacent and whose most recent common ancestor is shared. The package provides a variety of text-based analysis functions, as well as the ability to export genomic islands into formats suitable for viewing in a genome browser. We demonstrate the capabilities of the package with several examples from enteric bacteria, including an examination of the evolution of the acid fitness island in the genus Escherichia. In addition we use output from simulations and a set of known genomic islands from the literature to show that xenoGI can accurately identify genomic islands and place them on a phylogenetic tree.

Conclusions

xenoGI is an effective tool for studying the history of genomic island insertions in a clade of microbes. It identifies genomic islands, and determines which branch they inserted on within the phylogenetic tree for the clade. Such information is valuable because it helps us understand the adaptive path that has produced living species.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2038-0) contains supplementary material, which is available to authorized users.

Escherichia coli: an old friend with new tidings

Escherichia coli is one of the most-studied microorganisms worldwide but its characteristics are continually changing. Extraintestinal E. coli infections, such as urinary tract infections and neonatal sepsis, represent a huge public health problem. They are caused mainly by specialized extraintestinal pathogenic E. coli (ExPEC) strains that can innocuously colonize human hosts but can also cause disease upon entering a normally sterile body site. The virulence capability of such strains is determined by a combination of distinctive accessory traits, called virulence factors, in conjunction with their distinctive phylogenetic background. It is conceivable that by developing interventions against the most successful ExPEC lineages or their key virulence/colonization factors the associated burden of disease and health care costs could foreseeably be reduced in the future. On the other hand, one important problem worldwide is the increase of antimicrobial resistance shown by bacteria. As underscored in the last WHO global report, within a wide range of infectious agents including E. coli, antimicrobial resistance has reached an extremely worrisome situation that ‘threatens the achievements of modern medicine’. In the present review, an update of the knowledge about the pathogenicity, antimicrobial resistance and clinical aspects of this ‘old friend’ was presented.

Gene Flow and Molecular Innovation in Bacteria

The small molecules produced by environmental bacteria have been mainstays of both chemical and biological research for decades, and some have led to important therapeutic interventions. These small molecules have been shaped by natural selection as they evolved to fulfill changing functional roles in their native environments. This minireview describes some recent systematic studies providing illustrative examples that involve the acquisition and alteration of genetic information for molecular innovation by bacteria in well-defined environments. Two different bacterial genera are featured, Pseudonocardia and Salinispora, and, although the small-molecule repertoires of both have benefited from horizontal gene transfer, Pseudonocardia spp. have relied on plasmid-based tactics while Salinispora spp. have relied on chromosomally integrated genomic islands.

Metagenomic analysis of soil and freshwater from zoo agricultural area with organic fertilization

Microbial communities drive biogeochemical cycles in agricultural areas by decomposing organic materials and converting essential nutrients. Organic amendments improve soil quality by increasing the load of essential nutrients and enhancing the productivity. Additionally, fresh water used for irrigation can affect soil quality of agricultural soils, mainly due to the presence of microbial contaminants and pathogens. In this study, we investigated how microbial communities in irrigation water might contribute to the microbial diversity and function of soil. Whole-metagenomic sequencing approaches were used to investigate the taxonomic and the functional profiles of microbial communities present in fresh water used for irrigation, and in soil from a vegetable crop, which received fertilization with organic compost made from animal carcasses. The taxonomic analysis revealed that the most abundant genera were Polynucleobacter (~8% relative abundance) and Bacillus (~10%) in fresh water and soil from the vegetable crop, respectively. Low abundance (0.38%) of cyanobacterial groups were identified. Based on functional gene prediction, denitrification appears to be an important process in the soil community analysed here. Conversely, genes for nitrogen fixation were abundant in freshwater, indicating that the N-fixation plays a crucial role in this particular ecosystem. Moreover, pathogenicity islands, antibiotic resistance and potential virulence related genes were identified in both samples, but no toxigenic genes were detected. This study provides a better understanding of the community structure of an area under strong agricultural activity with regular irrigation and fertilization with an organic compost made from animal carcasses. Additionally, the use of a metagenomic approach to investigate fresh water quality proved to be a relevant method to evaluate its use in an agricultural ecosystem.

Stand-alone ClpG disaggregase confers superior heat tolerance to bacteria

AAA+ disaggregases solubilize aggregated proteins and confer heat tolerance to cells. Their disaggregation activities crucially depend on partner proteins, which target the AAA+ disaggregases to protein aggregates while concurrently stimulating their ATPase activities. Here, we report on two potent ClpG disaggregase homologs acquired through horizontal gene transfer by the species Pseudomonas aeruginosa and subsequently abundant P. aeruginosa clone C. ClpG exhibits high, stand-alone disaggregation potential without involving any partner cooperation. Specific molecular features, including high basal ATPase activity, a unique aggregate binding domain, and almost exclusive expression in stationary phase distinguish ClpG from other AAA+ disaggregases. Consequently, ClpG largely contributes to heat tolerance of P. aeruginosa primarily in stationary phase and boosts heat resistance 100-fold when expressed in Escherichia coli This qualifies ClpG as a potential persistence and virulence factor in P. aeruginosa.

Zisland Explorer: detect genomic islands by combining homogeneity and heterogeneity properties

Genomic islands are genomic fragments of alien origin in bacterial and archaeal genomes, usually involved in symbiosis or pathogenesis. In this work, we described Zisland Explorer, a novel tool to predict genomic islands based on the segmental cumulative GC profile. Zisland Explorer was designed with a novel strategy, as well as a combination of the homogeneity and heterogeneity of genomic sequences. While the sequence homogeneity reflects the composition consistence within each island, the heterogeneity measures the composition bias between an island and the core genome. The performance of Zisland Explorer was evaluated on the data sets of 11 different organisms. Our results suggested that the true-positive rate (TPR) of Zisland Explorer was at least 10.3% higher than that of four other widely used tools. On the other hand, the new tool did not lose overall accuracy with the improvement in the TPR and showed better equilibrium among various evaluation indexes. Also, Zisland Explorer showed better accuracy in the prediction of experimental island data. Overall, the tool provides an alternative solution over other tools, which expands the field of island prediction and offers a supplement to increase the performance of the distinct predicting strategy. We have provided a web service as well as a graphical user interface and open-source code across multiple platforms for Zisland Explorer, which is available at http://cefg.uestc.edu.cn/Zisland_Explorer/ or http://tubic.tju.edu.cn/Zisland_Explorer/.

The nitrate-reduction gene cluster components exert lineage-dependent contributions to optimization of Sinorhizobium symbiosis with soybeans

Receiving nodulation and nitrogen fixation genes does not guarantee rhizobia an effective symbiosis with legumes. Here, variations in gene content were determined for three Sinorhizobium species showing contrasting symbiotic efficiency on soybeans. A nitrate-reduction gene cluster absent in S. sojae was found to be essential for symbiotic adaptations of S. fredii and S. sp. III. In S. fredii, the deletion mutation of the nap (nitrate reductase), instead of nir (nitrite reductase) and nor (nitric oxide reductase), led to defects in nitrogen-fixation (Fix- ). By contrast, none of these core nitrate-reduction genes were required for the symbiosis of S. sp. III. However, within the same gene cluster, the deletion of hemN1 (encoding oxygen-independent coproporphyrinogen III oxidase) in both S. fredii and S. sp. III led to the formation of nitrogen-fixing (Fix+ ) but ineffective (Eff- ) nodules. These Fix+ /Eff- nodules were characterized by significantly lower enzyme activity of glutamine synthetase indicating rhizobial modulation of nitrogen-assimilation by plants. A distant homologue of HemN1 from S. sojae can complement this defect in S. fredii and S. sp. III, but exhibited a more pleotropic role in symbiosis establishment. These findings highlighted the lineage-dependent optimization of symbiotic functions in different rhizobial species associated with the same host.

Transformation Asymmetry and the Evolution of the Bacterial Accessory Genome

Bacterial transformation can insert or delete genomic islands (GIs), depending on the donor and recipient genotypes, if an homologous recombination spans the GI’s integration site and includes sufficiently long flanking homologous arms. Combining mathematical models of recombination with experiments using pneumococci found GI insertion rates declined geometrically with the GI’s size. The decrease in acquisition frequency with length (1.08×10⁻³bp⁻¹) was higher than a previous estimate of the analogous rate at which core genome recombinations terminated. Although most efficient for shorter GIs, transformation-mediated deletion frequencies did not vary consistently with GI length, with removal of 10-kb GIs ∼50% as efficient as acquisition of base substitutions. Fragments of 2 kb, typical of transformation event sizes, could drive all these deletions independent of island length. The strong asymmetry of transformation, and its capacity to efficiently remove GIs, suggests nonmobile accessory loci will decline in frequency without preservation by selection.

Recruitment of a Lineage-Specific Virulence Regulatory Pathway Promotes Intracellular Infection by a Plant Pathogen Experimentally Evolved into a Legume Symbiont

Ecological transitions between different lifestyles, such as pathogenicity, mutualism and saprophytism, have been very frequent in the course of microbial evolution, and often driven by horizontal gene transfer. Yet, how genomes achieve the ecological transition initiated by the transfer of complex biological traits remains poorly known. Here, we used experimental evolution, genomics, transcriptomics and high-resolution phenotyping to analyze the evolution of the plant pathogen Ralstonia solanacearum into legume symbionts, following the transfer of a natural plasmid encoding the essential mutualistic genes. We show that a regulatory pathway of the recipient R. solanacearum genome involved in extracellular infection of natural hosts was reused to improve intracellular symbiosis with the Mimosa pudica legume. Optimization of intracellular infection capacity was gained through mutations affecting two components of a new regulatory pathway, the transcriptional regulator efpR and a region upstream from the RSc0965-0967 genes of unknown functions. Adaptive mutations caused the downregulation of efpR and the over-expression of a downstream regulatory module, the three unknown genes RSc3146-3148, two of which encoding proteins likely associated to the membrane. This over-expression led to important metabolic and transcriptomic changes and a drastic qualitative and quantitative improvement of nodule intracellular infection. In addition, these adaptive mutations decreased the virulence of the original pathogen. The complete efpR/RSc3146-3148 pathway could only be identified in the genomes of the pathogenic R. solanacearum species complex. Our findings illustrate how the rewiring of a genetic network regulating virulence allows a radically different type of symbiotic interaction and contributes to ecological transitions and trade-offs.

Horizontal acquisition of a hypoxia-responsive molybdenum cofactor biosynthesis pathway contributed to Mycobacterium tuberculosis pathoadaptation

The unique ability of the tuberculosis (TB) bacillus, Mycobacterium tuberculosis, to persist for long periods of time in lung hypoxic lesions chiefly contributes to the global burden of latent TB. We and others previously reported that the M. tuberculosis ancestor underwent massive episodes of horizontal gene transfer (HGT), mostly from environmental species. Here, we sought to explore whether such ancient HGT played a part in M. tuberculosis evolution towards pathogenicity. We were interested by a HGT-acquired M. tuberculosis-specific gene set, namely moaA1-D1, which is involved in the biosynthesis of the molybdenum cofactor. Horizontal acquisition of this gene set was striking because homologues of these moa genes are present all across the Mycobacterium genus, including in M. tuberculosis. Here, we discovered that, unlike their paralogues, the moaA1-D1 genes are strongly induced under hypoxia. In vitro, a M. tuberculosis moaA1-D1-null mutant has an impaired ability to respire nitrate, to enter dormancy and to survive in oxygen-limiting conditions. Conversely, heterologous expression of moaA1-D1 in the phylogenetically closest non-TB mycobacterium, Mycobacterium kansasii, which lacks these genes, improves its capacity to respire nitrate and grants it with a marked ability to survive oxygen depletion. In vivo, the M. tuberculosis moaA1-D1-null mutant shows impaired survival in hypoxic granulomas in C3HeB/FeJ mice, but not in normoxic lesions in C57BL/6 animals. Collectively, our results identify a novel pathway required for M. tuberculosis resistance to host-imposed stress, namely hypoxia, and provide evidence that ancient HGT bolstered M. tuberculosis evolution from an environmental species towards a pervasive human-adapted pathogen.

Mycobacterium tuberculosis, the etiological agent of tuberculosis (TB), can persist for years and even decades in the lungs of its human host. Here we report that a unique M. tuberculosis gene cluster involved in the synthesis of the molybdenum cofactor, a cofactor for several oxidoreductases including the nitrate reductase, allows this major pathogen to respire nitrate and to persist in a dormant state under hypoxia, a stress condition encountered in lung TB lesions. Strikingly the M. tuberculosis ancestor, which most likely was an environmental harmless bacterium, acquired this gene cluster, together with its hypoxia-responsive transcriptional regulator, horizontally from neighboring bacteria. Our results uncover a key step in M. tuberculosis evolution towards pathogenicity.

Dissemination and loss of a biofilm-related genomic island in marine Pseudoalteromonas mediated by integrative and conjugative elements

Acquisition of genomic islands (GIs) plays a central role in the diversification and adaptation of bacteria. Some GIs can be mobilized in trans by integrative and conjugative elements (ICEs) or conjugative plasmids if the GIs carry specific transfer-related sequences. However, the transfer mechanism of GIs lacking such elements remains largely unexplored. Here, we investigated the transmissibility of a GI found in a coral-associated marine bacterium. This GI does not carry genes with transfer functions, but it carries four genes required for robust biofilm formation. Notably, this GI is inserted in the integration site for SXT/R391 ICEs. We demonstrated that acquisition of an SXT/R391 ICE results in either a tandem GI/ICE arrangement or the complete displacement of the GI. The GI displacement by the ICE greatly reduces biofilm formation. In contrast, the tandem integration of the ICE with the GI in cis allows the GI to hijack the transfer machinery of the ICE to excise, transfer and re-integrate into a new host. Collectively, our findings reveal that the integration of an ICE into a GI integration site enables rapid genome dynamics and a new mechanism by which SXT/R391 ICEs can augment genome plasticity.

Phenetic Comparison of Prokaryotic Genomes Using k-mers

Bacterial genomics studies are getting more extensive and complex, requiring new ways to envision analyses. Using the Ray Surveyor software, we demonstrate that comparison of genomes based on their k-mer content allows reconstruction of phenetic trees without the need of prior data curation, such as core genome alignment of a species. We validated the methodology using simulated genomes and previously published phylogenomic studies of Streptococcus pneumoniae and Pseudomonas aeruginosa. We also investigated the relationship of specific genetic determinants with bacterial population structures. By comparing clusters from the complete genomic content of a genome population with clusters from specific functional categories of genes, we can determine how the population structures are correlated. Indeed, the strain clustering based on a subset of k-mers allows determination of its similarity with the whole genome clusters. We also applied this methodology on 42 species of bacteria to determine the correlational significance of five important bacterial genomic characteristics. For example, intrinsic resistance is more important in P. aeruginosa than in S. pneumoniae, and the former has increased correlation of its population structure with antibiotic resistance genes. The global view of the pangenome of bacteria also demonstrated the taxa-dependent interaction of population structure with antibiotic resistance, bacteriophage, plasmid, and mobile element k-mer data sets.

Identification of variable genomic regions related to stress response in Oenococcus oeni

The lactic acid bacterium Oenococcus oeni is the most important species involved in malolactic fermentation due to its capability to survive in presence of ethanol and in the acidic environment of wine. In order to identify novel genes involved in adaptation to wine, a new approach using genome-wide analysis based on stress-related genes was performed in strain O. oeni PSU-1, and 106 annotated stress genes were identified. The in silico analysis revealed the high similarity of all those genes through 57 O. oeni genomes; however, seven variable regions of genomic plasticity could be determined for their different presence observed among these strains. Regions 3 and 5 had the typical hallmarks of horizontal transfer, suggesting that the strategy of acquiring genes from other bacteria enhanced the fitness of O. oeni strains. Certain genes related to stress resistance were described in these regions, and similarities of putative acquired regions with other lactic acid bacteria species were found. Some genomic fragments present in all the strains were described and another new genomic island harbouring a threonine dehydrogenase was found. The association of selected sequences with adaptation to wine was assessed by screening 31 O. oeni strains using PCR of single genes, but no sequences were found to be exclusive to highly performing malolactic fermentation strains. This study provides new information about the genomic variability of O. oeni strains contributing to a further understanding of this species and the relationship of its genomic traits with the ability to adapt to stress conditions.

Genomic insights into specialized metabolism in the marine actinomycete Salinispora

Comparative genomics is providing new opportunities to address the diversity and distributions of genes encoding the biosynthesis of specialized metabolites. An analysis of 119 genome sequences representing three closely related species of the marine actinomycete genus Salinispora reveals extraordinary biosynthetic diversity in the form of 176 distinct biosynthetic gene clusters (BGCs) of which only 24 have been linked to their products. Remarkably, more than half of the BGCs were observed in only one or two strains, suggesting they were acquired relatively recently in the evolutionary history of the genus. These acquired gene clusters are concentrated in specific genomic islands, which represent hot spots for BGC acquisition. While most BGCs are stable in terms of their chromosomal position, others migrated to different locations or were exchanged with unrelated gene clusters suggesting a plug and play type model of evolution that provides a mechanism to test the relative fitness effects of specialized metabolites. Transcriptome analyses were used to address the relationships between BGC abundance, chromosomal position and product discovery. The results indicate that recently acquired BGCs can be functional and that complex evolutionary processes shape the micro-diversity of specialized metabolism observed in closely related environmental bacteria.

Variation of genomic islands and flanking fragments in Vibrio parahaemolyticus isolates from environmental and clinical sources in Taiwan

Vibrio parahaemolyticus is a halophilic foodborne pathogenic bacterium that causes gastroenteritis; it has become an issue of global concern since the emergence and spread of pandemic O3:K6 strains. This study evaluated the role of Vibrio pathogenicity island (VPaI)-associated fragments in the genetic variation and grouping of this pathogen. Distribution of some VPaI fragments and flanking fragments (VPaI-1, VPaI-4, VPaI-5, VPaI-6 and VPaI-7) was determined in a total of 53 V. parahaemolyticus isolates from environmental and clinical sources in Taiwan, and supported by the sequences of seven fragments of VPaI-4 and its flanking fragment VP2145. As determined from the distribution of these VPaI-associated fragments, the clinical pandemic isolates were closely related in a single cluster; the clinical nonpandemic isolates were grouped into several clusters, while the environmental isolates were comparatively highly diversified. The profiles of virulence-associated genes of environmental pathogenic isolates varied, and were closer to those of clinical nonpandemic isolates than those of pandemic isolates. Isolates with atypical profiles of the VPaI-associated fragments and virulence-associated genes were identified. Sequences of VP2145 exhibited a close phylogenetic relationship among these local isolates, which were distinct from most V. parahaemolyticus strains from other geographic regions. This investigation demonstrated the application of VPaI-associated fragments in studying the genetic variation and clustering of V. parahaemolyticus isolates from different sources.