Molecular archaeology of the Escherichia coli genome

Discovery of mono-ADP ribosylating toxins with high structural homology to Pseudomonas exotoxin A

Mono-ADP-ribosyl transferase (mART) proteins are secreted virulence factors produced by several human pathogens, the founding member being diphtheria toxin (DT). Pseudomonas aeruginosa can also secrete a mART toxin, known as exotoxin A (PE), but with an organization of its three functional domains (receptor, translocation, and enzymatic elements) that is opposite to DT. Two additional PE-like toxins (PLTs) have been identified from Vibrio cholerae and Aeromonas hydrophila, suggesting more PLT family members may exist. Database mining discovered six additional putative homologues, considerably extending this group of PLTs across a wide range of bacterial species. Here, we examine sequence and structural information for these new family members with respect to previously identified PLTs. The X-ray crystal structures of four new homologues show the conservation of critical features responsible for structure and function. This study shows the potential of these newly described toxins for the development of novel drug delivery platforms. Additionally, genomic analysis suggests horizontal gene transfer to account for the wide distribution of PLTs across a range of eubacteria species, highlighting the need to monitor emerging pathogens and their virulence factors.

A transcription factor from the cryptic Escherichia coli Rac prophage controls both phage and host operons

Bacterial genomes are shaped by cryptic prophages, which are viral genomes integrated into the bacterial chromosome. Escherichia coli genomes have 10 prophages on average. Though usually inactive, prophage genes can profoundly impact host cell physiology. Among the phage genes in the E. coli chromosome, there are several putative transcription factors (TFs). These prophage TFs are predicted to control only phage promoters; however, their regulatory functions are not well characterized. The cohabitation of prophages and bacteria has led to conditions under which the majority of prophage genes are unexpressed, at least under normal growth conditions. We characterized a Rac prophage TF, YdaT, expression of which is normally inhibited by Rac TFs and, surprisingly, by the host global regulator OxyR. YdaT, when expressed, leads to a toxic phenotype manifested by drastic cell filamentation and cell death. We determined the binding sites and regulatory action for YdaT, finding two sites within the Rac locus, and one upstream of the host rcsA gene, which codes for the global regulator RcsA. The resulting increase in RcsA strongly impacts the bacterial RcsA/B regulon, which includes operons related to motility, capsule biosynthesis, colanic acid production, biofilm formation, and cell division. Our results provide novel insights into the host's genetic network, which appears to integrate YdaT in a complex manner, favoring its maintenance in the silenced state. The fact that the potentially toxic YdaT locus remains unmutated suggests its importance and potential benefits for the host, which may appear under stress conditions that are not yet known.

Role of the LuxR solo, SdiA, in eavesdropping on foreign bacteria

Bacteria can cooperate by coordinating their gene expression through the production, release, and detection of small molecules, a phenomenon known as quorum sensing (QS). One type of QS commonly found in Gram-negative bacteria utilizes a LuxI-type enzyme to produce a signaling molecule of the N-acyl-homoserine lactone (AHL) family, and a transcription factor of the LuxR family to detect and respond to the AHL. In a subset of Enterobacteriaceae, including Escherichia coli and Salmonella, no LuxI family member is present and no AHLs are synthesized. However, they encode a LuxR family member, SdiA, that is used to detect the QS molecules of other bacterial species, a behavior known as eavesdropping. Despite significant research on the topic, the overall role of SdiA-mediated eavesdropping in these bacteria remains unclear. In this review, we discuss the phenotypes and regulons of SdiA in the Enterobacteriaceae.

A new selective force driving metabolic gene clustering

The evolution of operons has puzzled evolutionary biologists since their discovery, and many theories exist to explain their emergence, spreading, and evolutionary conservation. In this work, we suggest that DNA replication introduces a selective force for the clustering of functionally related genes on chromosomes, which we interpret as a preliminary and necessary step in operon formation. Our reasoning starts from the observation that DNA replication produces copy number variations of genomic regions, and we propose that such changes perturb metabolism. The formalization of this effect by exploiting concepts from metabolic control analysis suggests that the minimization of such perturbations during evolution could be achieved through the formation of gene clusters and operons. We support our theoretical derivations with simulations based on a realistic metabolic network, and we confirm that present-day genomes have a degree of compaction of functionally related genes, which is significantly correlated to the proposed perturbations introduced by replication. The formation of clusters of functionally related genes in microbial genomes has puzzled microbiologists since their first discovery. Here, we suggest that replication, and the copy number variations due to the replisome passage, might play a role in the process through a perturbation in metabolite homeostasis. We provide theoretical support to this hypothesis, and we found that both simulations and genomic analysis support our hypothesis.

IMPORTANCE

The formation of clusters of functionally related genes in microbial genomes has puzzled microbiologists since their discovery. Here, we suggest that replication, and the copy number variations due to the replisome passage, might play a role in the process through a perturbation in metabolite homeostasis. We provide theoretical support to this hypothesis, and we found that both simulations and genomic analysis support our hypothesis.

Genomic diversity, antibiotic resistance, and virulence in South African Enterococcus faecalis and Enterococcus lactis isolates

This study presents the empirical findings of an in-depth genomic analysis of Enterococcus faecalis and Enterococcus lactis isolates from South Africa. It offers valuable insights into their genetic characteristics and their significant implications for public health. The study uncovers nuanced variations in the gene content of these isolates, despite their similar GC contents, providing a comprehensive view of the evolutionary diversity within the species. Genomic islands are identified, particularly in E. faecalis, emphasizing its propensity for horizontal gene transfer and genetic diversity, especially in terms of antibiotic resistance genes. Pangenome analysis reveals the existence of a core genome, accounting for a modest proportion of the total genes, with 2157 core genes, 1164 shell genes, and 4638 cloud genes out of 7959 genes in 52 South African E. faecalis genomes (2 from this study, 49 south Africa genomes downloaded from NCBI, and E. faecalis reference genome). Detecting large-scale genomic rearrangements, including chromosomal inversions, underscores the dynamic nature of bacterial genomes and their role in generating genetic diversity. The study uncovers an array of antibiotic resistance genes, with trimethoprim, tetracycline, glycopeptide, and multidrug resistance genes prevalent, raising concerns about the effectiveness of antibiotic treatment. Virulence gene profiling unveils a diverse repertoire of factors contributing to pathogenicity, encompassing adhesion, biofilm formation, stress resistance, and tissue damage. These empirical findings provide indispensable insights into these bacteria's genomic dynamics, antibiotic resistance mechanisms, and virulence potential, underlining the pressing need to address antibiotic resistance and implement robust control measures.

Soil Giant Phage: Genome and Biological Characteristics of Sinorhizobium Jumbo Phage

This paper presents the first in-depth research on the biological and genomic properties of lytic rhizobiophage AP-J-162 isolated from the soils of the mountainous region of Dagestan (North Caucasus), which belongs to the centers of origin of cultivated plants, according to Vavilov N.I. The rhizobiophage host strains are nitrogen-fixing bacteria of the genus Sinorhizobium spp., symbionts of leguminous forage grasses. The phage particles have a myovirus virion structure. The genome of rhizobiophage AP-J-162 is double-stranded DNA of 471.5 kb in length; 711 ORFs are annotated and 41 types of tRNAs are detected. The closest phylogenetic relative of phage AP-J-162 is Agrobacterium phage Atu-ph07, but no rhizobiophages are known. The replicative machinery, capsid, and baseplate proteins of phage AP-J-162 are structurally similar to those of Escherichia phage T4, but there is no similarity between their tail protein subunits. Amino acid sequence analysis shows that 339 of the ORFs encode hypothetical or functionally relevant products, while the remaining 304 ORFs are unique. Additionally, 153 ORFs are similar to those of Atu_ph07, with one-third of the ORFs encoding different enzymes. The biological properties and genomic characteristics of phage AP-J-162 distinguish it as a unique model for exploring phage-microbe interactions with nitrogen-fixing symbiotic microorganisms.

Specificities and commonalities of the Planctomycetes plasmidome

Plasmids, despite their critical role in antibiotic resistance and modern biotechnology, are understood in only a few bacterial groups in terms of their natural ecological dynamics. The bacterial phylum Planctomycetes, known for its unique molecular and cellular biology, has a largely unexplored plasmidome. This study offers a thorough exploration of the diversity of natural plasmids within Planctomycetes, which could serve as a foundation for developing various genetic research tools for this phylum. Planctomycetes plasmids encode a broad range of biological functions and appear to have coevolved significantly with their host chromosomes, sharing many homologues. Recent transfer events of insertion sequences between cohabiting chromosomes and plasmids were also observed. Interestingly, 64% of plasmid genes are distantly related to either chromosomally encoded genes or have homologues in plasmids from other bacterial groups. The planctomycetal plasmidome is composed of 36% exclusive proteins. Most planctomycetal plasmids encode a replication initiation protein from the Replication Protein A family near a putative iteron-containing replication origin, as well as active type I partition systems. The identification of one conjugative and three mobilizable plasmids suggests the occurrence of horizontal gene transfer via conjugation within this phylum. This comprehensive description enhances our understanding of the plasmidome of Planctomycetes and its potential implications in antibiotic resistance and biotechnology.

Origin, evolution, and maintenance of gene-strand bias in bacteria

Gene-strand bias is a characteristic feature of bacterial genome organization wherein genes are preferentially encoded on the leading strand of replication, promoting co-orientation of replication and transcription. This co-orientation bias has evolved to protect gene essentiality, expression, and genomic stability from the harmful effects of head-on replication-transcription collisions. However, the origin, variation, and maintenance of gene-strand bias remain elusive. Here, we reveal that the frequency of inversions that alter gene orientation exhibits large variation across bacterial populations and negatively correlates with gene-strand bias. The density, distance, and distribution of inverted repeats show a similar negative relationship with gene-strand bias explaining the heterogeneity in inversions. Importantly, these observations are broadly evident across the entire bacterial kingdom uncovering inversions and inverted repeats as primary factors underlying the variation in gene-strand bias and its maintenance. The distinct catalytic subunits of replicative DNA polymerase have co-evolved with gene-strand bias, suggesting a close link between replication and the origin of gene-strand bias. Congruently, inversion frequencies and inverted repeats vary among bacteria with different DNA polymerases. In summary, we propose that the nature of replication determines the fitness cost of replication-transcription collisions, establishing a selection gradient on gene-strand bias by fine-tuning DNA sequence repeats and, thereby, gene inversions.

Niche availability and competitive loss by facilitation control proliferation of bacterial strains intended for soil microbiome interventions

Microbiome engineering - the targeted manipulation of microbial communities - is considered a promising strategy to restore ecosystems, but experimental support and mechanistic understanding are required. Here, we show that bacterial inoculants for soil microbiome engineering may fail to establish because they inadvertently facilitate growth of native resident microbiomes. By generating soil microcosms in presence or absence of standardized soil resident communities, we show how different nutrient availabilities limit outgrowth of focal bacterial inoculants (three Pseudomonads), and how this might be improved by adding an artificial, inoculant-selective nutrient niche. Through random paired interaction assays in agarose micro-beads, we demonstrate that, in addition to direct competition, inoculants lose competitiveness by facilitating growth of resident soil bacteria. Metatranscriptomics experiments with toluene as selective nutrient niche for the inoculant Pseudomonas veronii indicate that this facilitation is due to loss and uptake of excreted metabolites by resident taxa. Generation of selective nutrient niches for inoculants may help to favor their proliferation for the duration of their intended action while limiting their competitive loss.

Pathogenomes and virulence profiles of representative big six non-O157 serogroup Shiga toxin-producing Escherichia coli

Shiga toxin (Stx)-producing Escherichia coli (STEC) of non-O157:H7 serotypes are responsible for global and widespread human food-borne disease. Among these serogroups, O26, O45, O103, O111, O121, and O145 account for the majority of clinical infections and are colloquially referred to as the "Big Six." The "Big Six" strain panel we sequenced and analyzed in this study are reference type cultures comprised of six strains representing each of the non-O157 STEC serogroups curated and distributed by the American Type Culture Collection (ATCC) as a resource to the research community under panel number ATCC MP-9. The application of long- and short-read hybrid sequencing yielded closed chromosomes and a total of 14 plasmids of diverse functions. Through high-resolution comparative phylogenomics, we cataloged the shared and strain-specific virulence and resistance gene content and established the close relationship of serogroup O26 and O103 strains featuring flagellar H-type 11. Virulence phenotyping revealed statistically significant differences in the Stx-production capabilities that we found to be correlated to the strain's individual stx-status. Among the carried Stx1a, Stx2a, and Stx2d phages, the Stx2a phage is by far the most responsive upon RecA-mediated phage mobilization, and in consequence, stx2a + isolates produced the highest-level of toxin in this panel. The availability of high-quality closed genomes for this "Big Six" reference set, including carried plasmids, along with the recorded genomic virulence profiles and Stx-production phenotypes will provide a valuable foundation to further explore the plasticity in evolutionary trajectories in these emerging non-O157 STEC lineages, which are major culprits of human food-borne disease.

A universal and constant rate of gene content change traces pangenome flux to LUCA

Prokaryotic genomes constantly undergo gene flux via lateral gene transfer, generating a pangenome structure consisting of a conserved core genome surrounded by a more variable accessory genome shell. Over time, flux generates change in genome content. Here, we measure and compare the rate of genome flux for 5655 prokaryotic genomes as a function of amino acid sequence divergence in 36 universally distributed proteins of the informational core (IC). We find a clock of gene content change. The long-term average rate of gene content flux is remarkably constant across all higher prokaryotic taxa sampled, whereby the size of the accessory genome-the proportion of the genome harboring gene content difference for genome pairs-varies across taxa. The proportion of species-level accessory genes per genome, varies from 0% (Chlamydia) to 30%-33% (Alphaproteobacteria, Gammaproteobacteria, and Clostridia). A clock-like rate of gene content change across all prokaryotic taxa sampled suggest that pangenome structure is a general feature of prokaryotic genomes and that it has been in existence since the divergence of bacteria and archaea.

Comparative genomics analysis and characterization of Shiga toxin-producing Escherichia coli O157:H7 strains reveal virulence genes, resistance genes, prophages and plasmids

Escherichia coli O157:H7 is a foodborne pathogen that has been linked to global disease outbreaks. These diseases include hemorrhagic colitis and hemolytic uremic syndrome. It is vital to know the features that make this strain pathogenic to understand the development of disease outbreaks. In the current study, a comparative genomic analysis was carried out to determine the presence of structural and functional features of O157:H7 strains obtained from 115 National Center for Biotechnology Information database. These strains of interest were analysed in the following programs: BLAST Ring Image Generator, PlasmidFinder, ResFinder, VirulenceFinder, IslandViewer 4 and PHASTER. Five strains (ECP19-198, ECP19-798, F7508, F8952, H2495) demonstrated a great homology with Sakai because of a few regions missing. Five resistant genes were identified, however, Macrolide-associated resistance gene mdf(A) was commonly found in all genomes. Majority of the strains (97%) were positive for 15 of the virulent genes (espA, espB, espF, espJ, gad, chuA, eae, iss, nleA, nleB, nleC, ompT, tccP, terC and tir). The plasmid analysis demonstrated that the IncF group was the most prevalent in the strains analysed. The prophage and genomic island analysis showed a distribution of bacteriophages and genomic islands respectively. The results indicated that structural and functional features of the many O157:H7 strains differ and may be a result of obtaining mobile genetic elements via horizontal gene transfer. Understanding the evolution of O157:H7 strains pathogenicity in terms of their structural and functional features will enable the development of detection and control of transmission strategies.

Sequence-discrete species for prokaryotes and other microbes: A historical perspective and pending issues

Whether prokaryotes, and other microorganisms, form distinct clusters that can be recognized as species remains an issue of paramount theoretical as well as practical consequence in identifying, regulating, and communicating about these organisms. In the past decade, comparisons of thousands of genomes of isolates and hundreds of metagenomes have shown that prokaryotic diversity may be predominantly organized in such sequence-discrete clusters, albeit organisms of intermediate relatedness between the identified clusters are also frequently found. Accumulating evidence suggests, however, that the latter "intermediate" organisms show enough ecological and/or functional distinctiveness to be considered different species. Notably, the area of discontinuity between clusters often-but not always-appears to be around 85%-95% genome-average nucleotide identity, consistently among different taxa. More recent studies have revealed remarkably similar diversity patterns for viruses and microbial eukaryotes as well. This high consistency across taxa implies a specific mechanistic process that underlies the maintenance of the clusters. The underlying mechanism may be a substantial reduction in the efficiency of homologous recombination, which mediates (successful) horizontal gene transfer, around 95% nucleotide identity. Deviations from the 95% threshold (e.g., species showing lower intraspecies diversity) may be caused by ecological differentiation that imposes barriers to otherwise frequent gene transfer. While this hypothesis that clusters are driven by ecological differentiation coupled to recombination frequency (i.e., higher recombination within vs. between groups) is appealing, the supporting evidence remains anecdotal. The data needed to rigorously test the hypothesis toward advancing the species concept are also outlined.

Ancient Rapid Radiation Explains Most Conflicts Among Gene Trees and Well-Supported Phylogenomic Trees of Nostocalean Cyanobacteria

Prokaryotic genomes are often considered to be mosaics of genes that do not necessarily share the same evolutionary history due to widespread horizontal gene transfers (HGTs). Consequently, representing evolutionary relationships of prokaryotes as bifurcating trees has long been controversial. However, studies reporting conflicts among gene trees derived from phylogenomic data sets have shown that these conflicts can be the result of artifacts or evolutionary processes other than HGT, such as incomplete lineage sorting, low phylogenetic signal, and systematic errors due to substitution model misspecification. Here, we present the results of an extensive exploration of phylogenetic conflicts in the cyanobacterial order Nostocales, for which previous studies have inferred strongly supported conflicting relationships when using different concatenated phylogenomic data sets. We found that most of these conflicts are concentrated in deep clusters of short internodes of the Nostocales phylogeny, where the great majority of individual genes have low resolving power. We then inferred phylogenetic networks to detect HGT events while also accounting for incomplete lineage sorting. Our results indicate that most conflicts among gene trees are likely due to incomplete lineage sorting linked to an ancient rapid radiation, rather than to HGTs. Moreover, the short internodes of this radiation fit the expectations of the anomaly zone, i.e., a region of the tree parameter space where a species tree is discordant with its most likely gene tree. We demonstrated that concatenation of different sets of loci can recover up to 17 distinct and well-supported relationships within the putative anomaly zone of Nostocales, corresponding to the observed conflicts among well-supported trees based on concatenated data sets from previous studies. Our findings highlight the important role of rapid radiations as a potential cause of strongly conflicting phylogenetic relationships when using phylogenomic data sets of bacteria. We propose that polytomies may be the most appropriate phylogenetic representation of these rapid radiations that are part of anomaly zones, especially when all possible genomic markers have been considered to infer these phylogenies. [Anomaly zone; bacteria; horizontal gene transfer; incomplete lineage sorting; Nostocales; phylogenomic conflict; rapid radiation; Rhizonema.].

Leveraging comparative genomics to uncover alien genes in bacterial genomes

A significant challenge in bacterial genomics is to catalogue genes acquired through the evolutionary process of horizontal gene transfer (HGT). Both comparative genomics and sequence composition-based methods have often been invoked to quantify horizontally acquired genes in bacterial genomes. Comparative genomics methods rely on completely sequenced genomes and therefore the confidence in their predictions increases as the databases become more enriched in completely sequenced genomes. Recent developments including in microbial genome sequencing call for reassessment of alien genes based on information-rich resources currently available. We revisited the comparative genomics approach and developed a new algorithm for alien gene detection. Our algorithm compared favourably with the existing comparative genomics-based methods and is capable of detecting both recent and ancient transfers. It can be used as a standalone tool or in concert with other complementary algorithms for comprehensively cataloguing alien genes in bacterial genomes.

Reconstructing horizontal gene flow network to understand prokaryotic evolution

Horizontal gene transfer (HGT) is a major source of phenotypic innovation and a mechanism of niche adaptation in prokaryotes. Quantification of HGT is critical to decipher its myriad roles in microbial evolution and adaptation. Advances in genome sequencing and bioinformatics have augmented our ability to understand the microbial world, particularly the direct or indirect influence of HGT on diverse life forms. Methods for detecting HGT can be classified into phylogenetic-based and parametric or composition-based approaches. Here, we exploited the complementary strengths of both the approaches to construct a high confidence horizontal gene flow network. Our network is unique in its ability to detect the transfer of native genes of a genome to genomes from other taxa, thus establishing donor and recipient organisms (taxa), rather than through a post hoc analysis as is the practice with several other approaches. The scale-free horizontal gene flow network presented here provides new insights into modes of transfer for the exchange of genetic information and also illuminates differential gene flow across phyla.

Reconstructing horizontal gene flow network to understand prokaryotic evolution

Horizontal gene transfer (HGT) is a major source of phenotypic innovation and a mechanism of niche adaptation in prokaryotes. Quantification of HGT is critical to decipher its myriad roles in microbial evolution and adaptation. Advances in genome sequencing and bioinformatics have augmented our ability to understand the microbial world, particularly the direct or indirect influence of HGT on diverse life forms. Methods for detecting HGT can be classified into phylogenetic-based and parametric or composition-based approaches. Here, we exploited the complementary strengths of both the approaches to construct a high confidence horizontal gene flow network. Our network is unique in its ability to detect the transfer of native genes of a genome to genomes from other taxa, thus establishing donor and recipient organisms (taxa), rather than through a post hoc analysis as is the practice with several other approaches. The scale-free horizontal gene flow network presented here provides new insights into modes of transfer for the exchange of genetic information and also illuminates differential gene flow across phyla.

Genomic analysis reveals the role of integrative and conjugative elements in plant pathogenic bacteria

Background

ICEs are mobile genetic elements found integrated into bacterial chromosomes that can excise and be transferred to a new cell. They play an important role in horizontal gene transmission and carry accessory genes that may provide interesting phenotypes for the bacteria. Here, we seek to research the presence and the role of ICEs in 300 genomes of phytopathogenic bacteria with the greatest scientific and economic impact.

Results

Seventy-eight ICEs (45 distinct elements) were identified and characterized in chromosomes of Agrobacterium tumefaciens, Dickeya dadantii, and D. solani, Pectobacterium carotovorum and P. atrosepticum, Pseudomonas syringae, Ralstonia solanacearum Species Complex, and Xanthomonas campestris. Intriguingly, the co-occurrence of four ICEs was observed in some P. syringae strains. Moreover, we identified 31 novel elements, carrying 396 accessory genes with potential influence on virulence and fitness, such as genes coding for functions related to T3SS, cell wall degradation and resistance to heavy metals. We also present the analysis of previously reported data on the expression of cargo genes related to the virulence of P. atrosepticum ICEs, which evidences the role of these genes in the infection process of tobacco plants.

Conclusions

Altogether, this paper has highlighted the potential of ICEs to affect the pathogenicity and lifestyle of these phytopathogens and direct the spread of significant putative virulence genes in phytopathogenic bacteria.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13100-022-00275-1.

Bacteria grow swiftly and live thriftily

Bacteria have evolved numerous strategies to use resources efficiently. However, bacterial economies depend on both the physiological context of the organisms as well as their growth state - whether they are growing, non-growing or reinitiating growth. In this essay, we discuss some of the features that make bacteria efficient under these different conditions and during the transitions between them. We also highlight the many outstanding questions regarding the physiology of non-growing bacterial cells. Lastly, we examine how efficiency is apparent in both the mode and tempo of bacterial evolution.

The microbial ecology of Escherichia coli in the vertebrate gut

Escherichia coli has a rich history as biology's 'rock star', driving advances across many fields. In the wild, E. coli resides innocuously in the gut of humans and animals but is also a versatile pathogen commonly associated with intestinal and extraintestinal infections and antimicrobial resistance-including large foodborne outbreaks such as the one that swept across Europe in 2011, killing 54 individuals and causing approximately 4000 infections and 900 cases of haemolytic uraemic syndrome. Given that most E. coli are harmless gut colonizers, an important ecological question plaguing microbiologists is what makes E. coli an occasionally devastating pathogen? To address this question requires an enhanced understanding of the ecology of the organism as a commensal. Here, we review how our knowledge of the ecology and within-host diversity of this organism in the vertebrate gut has progressed in the 137 years since E. coli was first described. We also review current approaches to the study of within-host bacterial diversity. In closing, we discuss some of the outstanding questions yet to be addressed and prospects for future research.

Recombination resolves the cost of horizontal gene transfer in experimental populations of Helicobacter pylori

Horizontal gene transfer (HGT)—the transfer of DNA between lineages—is responsible for a large proportion of the genetic variation that contributes to evolution in microbial populations. While HGT can bring beneficial genetic innovation, the transfer of DNA from other species or strains can also have deleterious effects. In this study, we evolve populations of the bacteria Helicobacter pylori and use DNA sequencing to identify over 40,000 genetic variants transferred by HGT. We measure the cost of many of these and find that both strongly beneficial mutations and deleterious mutations are genetic variants transferred by natural transformation. Importantly, we also show how recombination that separates linked beneficial and deleterious mutations resolves the cost of HGT.

Horizontal gene transfer (HGT) is important for microbial evolution, yet we know little about the fitness effects and dynamics of horizontally transferred genetic variants. In this study, we evolve laboratory populations of Helicobacter pylori, which take up DNA from their environment by natural transformation, and measure the fitness effects of thousands of transferred genetic variants. We find that natural transformation increases the rate of adaptation but comes at the cost of significant genetic load. We show that this cost is circumvented by recombination, which increases the efficiency of selection by decoupling deleterious and beneficial genetic variants. Our results show that adaptation with HGT, pervasive in natural microbial populations, is shaped by a combination of selection, recombination, and genetic drift not accounted for in existing models of evolution.

MetaPop: a pipeline for macro- and microdiversity analyses and visualization of microbial and viral metagenome-derived populations

BACKGROUND

Microbes and their viruses are hidden engines driving Earth's ecosystems from the oceans and soils to humans and bioreactors. Though gene marker approaches can now be complemented by genome-resolved studies of inter-(macrodiversity) and intra-(microdiversity) population variation, analytical tools to do so remain scattered or under-developed.

RESULTS

Here, we introduce MetaPop, an open-source bioinformatic pipeline that provides a single interface to analyze and visualize microbial and viral community metagenomes at both the macro- and microdiversity levels. Macrodiversity estimates include population abundances and α- and β-diversity. Microdiversity calculations include identification of single nucleotide polymorphisms, novel codon-constrained linkage of SNPs, nucleotide diversity (π and θ), and selective pressures (pN/pS and Tajima's D) within and fixation indices (FST) between populations. MetaPop will also identify genes with distinct codon usage. Following rigorous validation, we applied MetaPop to the gut viromes of autistic children that underwent fecal microbiota transfers and their neurotypical peers. The macrodiversity results confirmed our prior findings for viral populations (microbial shotgun metagenomes were not available) that diversity did not significantly differ between autistic and neurotypical children. However, by also quantifying microdiversity, MetaPop revealed lower average viral nucleotide diversity (π) in autistic children. Analysis of the percentage of genomes detected under positive selection was also lower among autistic children, suggesting that higher viral π in neurotypical children may be beneficial because it allows populations to better "bet hedge" in changing environments. Further, comparisons of microdiversity pre- and post-FMT in autistic children revealed that the delivery FMT method (oral versus rectal) may influence viral activity and engraftment of microdiverse viral populations, with children who received their FMT rectally having higher microdiversity post-FMT. Overall, these results show that analyses at the macro level alone can miss important biological differences.

CONCLUSIONS

These findings suggest that standardized population and genetic variation analyses will be invaluable for maximizing biological inference, and MetaPop provides a convenient tool package to explore the dual impact of macro- and microdiversity across microbial communities. Video abstract.

Hidden diversity of double-stranded DNA phages in symbiotic Rhizobium species

In this study, we addressed the extent of diversification of phages associated with nitrogen-fixing symbiotic Rhizobium species. Despite the ecological and economic importance of the Rhizobium genus, little is known about the diversity of the associated phages. A thorough assessment of viral diversity requires investigating both lytic phages and prophages harboured in diverse Rhizobium genomes. Protein-sharing networks identified 56 viral clusters (VCs) among a set of 425 isolated phages and predicted prophages. The VCs formed by phages had more proteins in common and a higher degree of synteny, and they group together in clades in the associated phylogenetic tree. By contrast, the VCs of prophages showed significant genetic variation and gene loss, with selective pressure on the remaining genes. Some VCs were found in various Rhizobium species and geographical locations, suggesting that they have wide host ranges. Our results indicate that the VCs represent distinct taxonomic units, probably representing taxa equivalent to genera or even species. The finding of previously undescribed phage taxa indicates the need for further exploration of the diversity of phages associated with Rhizobium species. This article is part of the theme issue 'The secret lives of microbial mobile genetic elements'.

Genomic analysis and characterization of Scenedesmus glucoliberatum PABB004: An unconventional sugar-secreting green alga

AIMS

In this report, we present Scenedesmus glucoliberatum PABB004, a microalga that was isolated from an association with Paramecium bursaria with the potential for application in fermentative processes and co-culture schemes due to its advantageous high sugar secretion phenotype.

METHODS AND RESULTS

We sequenced, assembled and annotated the draft genome and transcriptome for this newly reported strain. The nuclear genome has an exceptionally high GC content of 78%. Our results revealed significant sugar accumulation over a range from 6.2 to 7.8 pH units. The predicted proteome was compared with other green algae that show different sugar secretion phenotypes aiming to help uncover their common features for simple sugar secretion and those unique to S. glucoliberatum PABB004.

CONCLUSIONS

The evolutionary history of this organism, inferred from its genomic traits, expands our current understanding of algal mutualistic relationships involving photosynthate exchanges. S. glucoliberatum PABB004 secreted ready-to-use fermentable sugars (glucose and maltose) directly to the extracellular media achieving concentrations greater than 2.7 g/L of free glucose and 1.2 g/L of maltose in batch cultures.

SIGNIFICANCE AND IMPACT OF THE STUDY

A draft genome is provided for a new member of an important class of green algae. Scenedesmus glucoliberatum PABB004 secretes high levels of simple sugars over a broad pH range.

Systematic Detection of Large-Scale Multigene Horizontal Transfer in Prokaryotes

Horizontal gene transfer (HGT) is central to prokaryotic evolution. However, little is known about the “scale” of individual HGT events. In this work, we introduce the first computational framework to help answer the following fundamental question: How often does more than one gene get horizontally transferred in a single HGT event? Our method, called HoMer, uses phylogenetic reconciliation to infer single-gene HGT events across a given set of species/strains, employs several techniques to account for inference error and uncertainty, combines that information with gene order information from extant genomes, and uses statistical analysis to identify candidate horizontal multigene transfers (HMGTs) in both extant and ancestral species/strains. HoMer is highly scalable and can be easily used to infer HMGTs across hundreds of genomes. We apply HoMer to a genome-scale data set of over 22,000 gene families from 103 Aeromonas genomes and identify a large number of plausible HMGTs of various scales at both small and large phylogenetic distances. Analysis of these HMGTs reveals interesting relationships between gene function, phylogenetic distance, and frequency of multigene transfer. Among other insights, we find that 1) the observed relative frequency of HMGT increases as divergence between genomes increases, 2) HMGTs often have conserved gene functions, and 3) rare genes are frequently acquired through HMGT. We also analyze in detail HMGTs involving the zonula occludens toxin and type III secretion systems. By enabling the systematic inference of HMGTs on a large scale, HoMer will facilitate a more accurate and more complete understanding of HGT and microbial evolution.

Detection of horizontal gene transfer in the genome of the choanoflagellate Salpingoeca rosetta

Horizontal gene transfer (HGT), the movement of heritable materials between distantly related organisms, is crucial in eukaryotic evolution. However, the scale of HGT in choanoflagellates, the closest unicellular relatives of metazoans, and its possible roles in the evolution of animal multicellularity remains unexplored. We identified at least 175 candidate HGTs in the genome of the colonial choanoflagellate Salpingoeca rosetta using sequence-based tests. The majority of these were orthologous to genes in bacterial and microalgal lineages, yet displayed genomic features consistent with the rest of the S. rosetta genome-evidence of ancient acquisition events. Putative functions include enzymes involved in amino acid and carbohydrate metabolism, cell signaling, and the synthesis of extracellular matrix components. Functions of candidate HGTs may have contributed to the ability of choanoflagellates to assimilate novel metabolites, thereby supporting adaptation, survival in diverse ecological niches, and response to external cues that are possibly critical in the evolution of multicellularity in choanoflagellates.

Rapid discovery of novel prophages using biological feature engineering and machine learning

Prophages are phages that are integrated into bacterial genomes and which are key to understanding many aspects of bacterial biology. Their extreme diversity means they are challenging to detect using sequence similarity, yet this remains the paradigm and thus many phages remain unidentified. We present a novel, fast and generalizing machine learning method based on feature space to facilitate novel prophage discovery. To validate the approach, we reanalyzed publicly available marine viromes and single-cell genomes using our feature-based approaches and found consistently more phages than were detected using current state-of-the-art tools while being notably faster. This demonstrates that our approach significantly enhances bacteriophage discovery and thus provides a new starting point for exploring new biologies.

A One Health Perspective for Defining and Deciphering Escherichia coli Pathogenic Potential in Multiple Hosts

E. coli is one of the most common species of bacteria colonizing humans and animals. The singularity of E. coli 's genus and species underestimates its multifaceted nature, which is represented by different strains, each with different combinations of distinct virulence factors. In fact, several E. coli pathotypes, or hybrid strains, may be associated with both subclinical infection and a range of clinical conditions, including enteric, urinary, and systemic infections. E. coli may also express DNA-damaging toxins that could impact cancer development. This review summarizes the different E. coli pathotypes in the context of their history, hosts, clinical signs, epidemiology, and control. The pathotypic characterization of E. coli in the context of disease in different animals, including humans, provides comparative and One Health perspectives that will guide future clinical and research investigations of E. coli infections.

Whole genome phylogenies reflect the distributions of recombination rates for many bacterial species

Although recombination is accepted to be common in bacteria, for many species robust phylogenies with well-resolved branches can be reconstructed from whole genome alignments of strains, and these are generally interpreted to reflect clonal relationships. Using new methods based on the statistics of single-nucleotide polymorphism (SNP) splits, we show that this interpretation is incorrect. For many species, each locus has recombined many times along its line of descent, and instead of many loci supporting a common phylogeny, the phylogeny changes many thousands of times along the genome alignment. Analysis of the patterns of allele sharing among strains shows that bacterial populations cannot be approximated as either clonal or freely recombining but are structured such that recombination rates between lineages vary over several orders of magnitude, with a unique pattern of rates for each lineage. Thus, rather than reflecting clonal ancestry, whole genome phylogenies reflect distributions of recombination rates.

Plasmids encode niche-specific traits in Lactobacillaceae

Species belonging to the family Lactobacillaceae are found in highly diverse environments and play an important role in fermented foods and probiotic products. Many of these species have been individually reported to harbour plasmids that encode important genes. In this study, we performed comparative genomic analysis of publicly available data for 512 plasmids from 282 strains represented by 51 species of this family and correlated the genomic features of plasmids with the ecological niches in which these species are found. Two-thirds of the species had at least one plasmid-harbouring strain. Plasmid abundance and GC content were significantly lower in vertebrate-adapted species as compared to nomadic and free-living species. Hierarchical clustering highlighted the distinct nature of plasmids from the nomadic and free-living species than those from the vertebrate-adapted species. EggNOG-assisted functional annotation revealed that genes associated with transposition, conjugation, DNA repair and recombination, exopolysaccharide production, metal ion transport, toxin–antitoxin system, and stress tolerance were significantly enriched on the plasmids of the nomadic and in some cases nomadic and free-living species. On the other hand, genes related to anaerobic metabolism, ABC transporters and the major facilitator superfamily were overrepresented on the plasmids of the vertebrate-adapted species. These genomic signatures correlate with the comparatively nutrient-depleted, stressful and dynamic environments of nomadic and free-living species and nutrient-rich and anaerobic environments of vertebrate-adapted species. Thus, these results indicate the contribution of the plasmids in the adaptation of lactobacilli to their respective habitats. This study also underlines the potential application of these plasmids in improving the technological and probiotic properties of lactic acid bacteria.

Deep phenotypic characterization of immunization-induced antibacterial IgG repertoires in mice using a single-antibody bioassay

Antibodies with antibacterial activity need to bind to the bacterial surface with affinity, specificity, and sufficient density to induce efficient elimination. To characterize the anti-bacterial antibody repertoire, we developed an in-droplet bioassay with single-antibody resolution. The assay not only allowed us to identify whether the secreted antibodies recognized a bacterial surface antigen, but also to estimate the apparent dissociation constant (KD app) of the interaction and the density of the recognized epitope on the bacteria. Herein, we found substantial differences within the KD app/epitope density profiles in mice immunized with various species of heat-killed bacteria. The experiments further revealed a high cross-reactivity of the secreted IgG repertoires, binding to even unrelated bacteria with high affinity. This application confirmed the ability to quantify the anti-bacterial antibody repertoire and the utility of the developed bioassay to study the interplay between bacteria and the humoral response.

Selection for ancient periodic motifs that do not impart DNA bending

A ~10-11 bp periodicity in dinucleotides imparting DNA bending, with shorter periods found in organisms with positively-supercoiled DNA and longer periods found in organisms with negatively-supercoiled DNA, was previously suggested to assist in DNA compaction. However, when measured with more robust methods, variation in the observed periods between organisms with different growth temperatures is not consistent with that hypothesis. We demonstrate that dinucleotide periodicity does not arise solely by mutational biases but is under selection. We found variation between genomes in both the period and the suite of dinucleotides that are periodic. Whereas organisms with similar growth temperatures have highly variable periods, differences in periods increase with phylogenetic distance between organisms. In addition, while the suites of dinucleotides under selection for periodicity become more dissimilar among more distantly-related organisms, there is a core set of dinucleotides that are strongly periodic among genomes in all domains of life. Notably, this core set of periodic motifs are not involved in DNA bending. These data indicate that dinucleotide periodicity is an ancient genomic architecture which may play a role in shaping the evolution of genes and genomes.

Biogeochemical Implications of N2O-Reducing Thermophilic Campylobacteria in Deep-Sea Vent Fields, and the Description of Nitratiruptor labii sp. nov

Nitrous oxide (N2O) is a potent greenhouse gas and has significantly increased in the atmosphere. Deep-sea hydrothermal fields are representative environments dominated by mesophilic to thermophilic members of the class Campylobacteria that possess clade II nosZ encoding nitrous oxide reductase. Here, we report a strain HRV44T representing the first thermophilic campylobacterium capable of growth by H2 oxidation coupled to N2O reduction. On the basis of physiological and genomic properties, it is proposed that strain HRV44T (=JCM 34002 = DSM 111345) represents a novel species of the genus Nitratiruptor, Nitratiruptor labii sp. nov. The comparison of the N2O consumption ability of strain HRV44T with those of additional Nitratiruptor and other campylobacterial strains revealed the highest level in strain HRV44T and suggests the N2O-respiring metabolism might be the common physiological trait for the genus Nitratiruptor. Our findings provide insights into contributions of thermophilic Campylobacteria to the N2O sink in deep-sea hydrothermal environments.

Evolutionary Mechanisms of Long-Term Genome Diversification Associated With Niche Partitioning in Marine Picocyanobacteria

Marine picocyanobacteria of the genera Prochlorococcus and Synechococcus are the most abundant photosynthetic organisms on Earth, an ecological success thought to be linked to the differential partitioning of distinct ecotypes into specific ecological niches. However, the underlying processes that governed the diversification of these microorganisms and the appearance of niche-related phenotypic traits are just starting to be elucidated. Here, by comparing 81 genomes, including 34 new Synechococcus, we explored the evolutionary processes that shaped the genomic diversity of picocyanobacteria. Time-calibration of a core-protein tree showed that gene gain/loss occurred at an unexpectedly low rate between the different lineages, with for instance 5.6 genes gained per million years (My) for the major Synechococcus lineage (sub-cluster 5.1), among which only 0.71/My have been fixed in the long term. Gene content comparisons revealed a number of candidates involved in nutrient adaptation, a large proportion of which are located in genomic islands shared between either closely or more distantly related strains, as identified using an original network construction approach. Interestingly, strains representative of the different ecotypes co-occurring in phosphorus-depleted waters (Synechococcus clades III, WPC1, and sub-cluster 5.3) were shown to display different adaptation strategies to this limitation. In contrast, we found few genes potentially involved in adaptation to temperature when comparing cold and warm thermotypes. Indeed, comparison of core protein sequences highlighted variants specific to cold thermotypes, notably involved in carotenoid biosynthesis and the oxidative stress response, revealing that long-term adaptation to thermal niches relies on amino acid substitutions rather than on gene content variation. Altogether, this study not only deciphers the respective roles of gene gains/losses and sequence variation but also uncovers numerous gene candidates likely involved in niche partitioning of two key members of the marine phytoplankton.

Potential evolutionary impact of integrative and conjugative elements (ICEs) and genomic islands in the Ralstonia solanacearum species complex

Ralstonia solanacearum, a soil-borne plant pathogen, encompasses a large number of strains known as R. solanacearum species complex (RSSC). Although it has been suggested that mobile genetic elements (MGEs) may play an important role in the RSSC genome, the evolutionary impact of these elements remains unknown. Here, we identified and analysed Integrative and Conjugative Elements (ICEs) and Genomic Islands (GIs) in the 121 genomes published for Ralstonia spp., including RSSC strains and three other non-plant pathogenic Ralstonia spp. Our results provided a dataset of 12 ICEs and 31 GIs distributed throughout Ralstonia spp. Four novel ICEs in RSSC were found. Some of these elements cover 5% of the host genome and carry accessory genes with a potential impact on the fitness and pathogenicity of RSSC. In addition, phylogenetic analysis revealed that these MGEs clustered to the same species, but there is evidence of strains from different countries that host the same element. Our results provide novel insight into the RSSC adaptation, opening new paths to a better understanding of how these elements affect this soil-borne plant pathogen.

Population genomics of Vibrionaceae isolated from an endangered oasis reveals local adaptation after an environmental perturbation

Background

In bacteria, pan-genomes are the result of an evolutionary “tug of war” between selection and horizontal gene transfer (HGT). High rates of HGT increase the genetic pool and the effective population size (N_e), resulting in open pan-genomes. In contrast, selective pressures can lead to local adaptation by purging the variation introduced by HGT and mutation, resulting in closed pan-genomes and clonal lineages. In this study, we explored both hypotheses, elucidating the pan-genome of Vibrionaceae isolates after a perturbation event in the endangered oasis of Cuatro Ciénegas Basin (CCB), Mexico, and looking for signals of adaptation to the environments in their genomes.

Results

We obtained 42 genomes of Vibrionaceae distributed in six lineages, two of them did not showed any close reference strain in databases. Five of the lineages showed closed pan-genomes and were associated to either water or sediment environment; their high N_e estimates suggest that these lineages are not from a recent origin. The only clade with an open pan-genome was found in both environments and was formed by ten genetic groups with low N_e, suggesting a recent origin. The recombination and mutation estimators (r/m) ranged from 0.005 to 2.725, which are similar to oceanic Vibrionaceae estimations. However, we identified 367 gene families with signals of positive selection, most of them found in the core genome; suggesting that despite recombination, natural selection moves the Vibrionaceae CCB lineages to local adaptation, purging the genomes and keeping closed pan-genome patterns. Moreover, we identify 598 SNPs associated with an unstructured environment; some of the genes associated with these SNPs were related to sodium transport.

Conclusions

Different lines of evidence suggest that the sampled Vibrionaceae, are part of the rare biosphere usually living under famine conditions. Two of these lineages were reported for the first time. Most Vibrionaceae lineages of CCB are adapted to their micro-habitats rather than to the sampled environments. This pattern of adaptation is concordant with the association of closed pan-genomes and local adaptation.

Determination of Plasmid pSN1216-29 Host Range and the Similarity in Oligonucleotide Composition Between Plasmid and Host Chromosomes

Plasmids are extrachromosomal DNA that can be horizontally transferred between different bacterial cells by conjugation. Horizontal gene transfer of plasmids can promote rapid evolution and adaptation of bacteria by imparting various traits involved in antibiotic resistance, virulence, and metabolism to their hosts. The host range of plasmids is an important feature for understanding how they spread in environmental microbial communities. Earlier bioinformatics studies have demonstrated that plasmids are likely to have similar oligonucleotide (k-mer) compositions to their host chromosomes and that evolutionary host ranges of plasmids could be predicted from this similarity. However, there are no complementary studies to assess the consistency between the predicted evolutionary host range and experimentally determined replication/transfer host range of a plasmid. In the present study, the replication/transfer host range of a model plasmid, pSN1216-29, exogenously isolated from cow manure as a newly discovered self-transmissible plasmid, was experimentally determined within microbial communities extracted from soil and cow manure. In silico prediction of evolutionary host range was performed with the pSN1216-29 using its oligonucleotide compositions independently. The results showed that oligonucleotide compositions of the plasmid pSN1216-29 had more similarities to those of hosts (transconjugants genera) than those of non-hosts (other genera). These findings can contribute to the understanding of how plasmids behave in microbial communities, and aid in the designing of appropriate plasmid vectors for different bacteria.

Codon usage trend in genes associated with obesity

Obesity is not only a social menace but also an economic burden as it reduces productivity and increases health care cost. We used bioinformatic tools to analyze the CUB of obesity associated genes and compared with housekeeping genes (control) to explore the similarities and differences between two data sets as no work was reported yet. The mean effective number of codons (ENC) in genes associated with obesity and housekeeping gene was 50.45 and 52.03 respectively, indicating low CUB. The relative synonymous codon usage (RSCU) suggested that codons namely CTG and GTG were over-represented in both obesity and housekeeping genes while under-represented codons were TCG, TTA, CTA, CCG, CAA, CGT, ATA, ACG, GTA and GCG in obesity genes and TCG, TTA, CCG, ATA, ACG, GTA, and GCG in housekeeping genes. t test analysis suggested that 11 codons namely TTA (Leu), TTG (Leu), CCG (Pro), CAC (His), CAA (Gln), CAG (Gln), CGT (Arg), AGA (Arg), ATA (Ile), ATT (Ile) and GCG (Ala) were significantly differed (p < 0.05 or p < 0.01) between obesity and housekeeping genes. Highly significant correlation was observed between GC12 and GC3 in obesity and housekeeping genes i.e. r = 0.580** and r = 0.498** (p < 0.01) respectively indicating the effect of directional mutation pressure present in all codon positions.

Early fate of exogenous promoters in E. coli

Gene gain by horizontal gene transfer is a major pathway of genome innovation in bacteria. The current view posits that acquired genes initially need to be silenced and that a bacterial chromatin protein, H-NS, plays a role in this silencing. However, we lack direct observation of the early fate of a horizontally transferred gene to prove this theory. We combine sequencing, flow cytometry and sorting, followed by microscopy to monitor gene expression and its variability after large-scale random insertions of a reporter gene in a population of Escherichia coli bacteria. We find that inserted promoters have a wide range of gene-expression variability related to their location. We find that high-expression clones carry insertions that are not correlated with H-NS binding. Conversely, binding of H-NS correlates with silencing. Finally, while most promoters show a common level of extrinsic noise, some insertions show higher noise levels. Analysis of these high-noise clones supports a scenario of switching due to transcriptional interference from divergent ribosomal promoters. Altogether, our findings point to evolutionary pathways where newly-acquired genes are not necessarily silenced, but may immediately explore a wide range of expression levels to probe the optimal ones.

Excessive Promoters as Silencers of Genes Horizontally Acquired by Escherichia coli

Horizontally acquired genes are usually transcriptionally inactive, although most of them are associated with genomic loci enriched with promoter-like sequences forming “promoter islands.” We hypothesized that lateral DNA transfer induces local mutagenesis, accumulating AT base pairs and creating promoter-like sequences, whose occupancy with RNA polymerase and a specific silencer H-NS suppresses the transcription of foreign genes. Error-prone mutagenesis was implemented for the “promoter island” of a foreign gene appY and the promoter region of an inherent gene dps. Derivatives with changed transcriptional activity were selected using a reporter plasmid pET28_eGFP. Only one cycle of mutagenesis with negative selection suppressed the activity of the main dps promoter to the background level due to a single substitution in its -10 element, while positive selection gave a sequence with improved -35 element, thus testifying feasibility of the approach. The same suppression for appY was achieved by three cycles, while eightfold transcription activation required nine iterations of mutagenesis. In both cases, the number of potential start points decreased resulting in an ordinary regulatory region with only one dominant promoter in the case of positive selection. Efficiency of H-NS binding remained virtually unchanged in all mutant constructs. Based on these findings we conclude that excessive promoters can adversely affect transcription by providing a platform for interference between several RNA polymerase molecules, which can act as a silencer at promoter-dense regions.

Herpesvirus-bacteria synergistic interaction in periodontitis

The etiopathogenesis of severe periodontitis includes herpesvirus-bacteria coinfection. This article evaluates the pathogenicity of herpesviruses (cytomegalovirus and Epstein-Barr virus) and periodontopathic bacteria (Aggregatibacter actinomycetemcomitans and Porphyromonas gingivalis) and coinfection of these infectious agents in the initiation and progression of periodontitis. Cytomegalovirus and A. actinomycetemcomitans/P. gingivalis exercise synergistic pathogenicity in the development of localized ("aggressive") juvenile periodontitis. Cytomegalovirus and Epstein-Barr virus are associated with P. gingivalis in adult types of periodontitis. Periodontal herpesviruses that enter the general circulation may also contribute to disease development in various organ systems. A 2-way interaction is likely to occur between periodontal herpesviruses and periodontopathic bacteria, with herpesviruses promoting bacterial upgrowth, and bacterial factors reactivating latent herpesviruses. Bacterial-induced gingivitis may facilitate herpesvirus colonization of the periodontium, and herpesvirus infections may impede the antibacterial host defense and alter periodontal cells to predispose for bacterial adherence and invasion. Herpesvirus-bacteria synergistic interactions, are likely to comprise an important pathogenic determinant of aggressive periodontitis. However, mechanistic investigations into the molecular and cellular interaction between periodontal herpesviruses and bacteria are still scarce. Herpesvirus-bacteria coinfection studies may yield significant new discoveries of pathogenic determinants, and drug and vaccine targets to minimize or prevent periodontitis and periodontitis-related systemic diseases.

In silico genome analysis of an acid mine drainage species, Acidiphilium multivorum, for potential commercial acetic acid production and biomining

The genome of Acidiphilium multivorum strain AIU 301, acidophilic, aerobic Gram-negative bacteria, was investigated for potential metabolic pathways associated with organic acid production and metal uptake. The genome was compared to other acidic mine drainage isolates, Acidiphilium cryptum JF-5 and Acidithiobacillus ferrooxidans ATCC 23270, as well as Acetobacter pasteurianus 386B, which ferments cocoa beans. Plasmids between two Acidiphilium spp. were compared, and only two of the sixteen plasmids were identified as potentially similar. Comparisons of the genome size to the number of protein coding sequences indicated that A. multivorum and A. cryptum follow the line of best fit unlike A. pasteurianus 386B, which suggests that it was improperly annotated in the database. Pathways between these four species were analyzed bioinformatically and are discussed here. A. multivorum AIU 301, shares pathways with A. pasteurianus 386B including aldehyde and alcohol dehydrogenase pathways, which are used in the generation of vinegar. Mercury reductase, arsenate reductase and sulfur utilization proteins were identified and discussed at length. The absence of sulfur utilization proteins from A. multivorum AIU 301 suggests that this species uses previously undefined pathways for sulfur acquisition. Bioinformatic examination revealed novel pathways that may benefit commercial fields including acetic acid production and biomining.

IslandCafe: Compositional Anomaly and Feature Enrichment Assessment for Delineation of Genomic Islands

One of the evolutionary forces driving bacterial genome evolution is the acquisition of clusters of genes through horizontal gene transfer (HGT). These genomic islands may confer adaptive advantages to the recipient bacteria, such as, the ability to thwart antibiotics, become virulent or hypervirulent, or acquire novel metabolic traits. Methods for detecting genomic islands either search for markers or features typical of islands or examine anomaly in oligonucleotide composition against the genome background. The former tends to underestimate, missing islands that have the markers either lost or degraded, while the latter tends to overestimate, due to their inability to discriminate compositional atypicality arising because of HGT from those that are a consequence of other biological factors. We propose here a framework that exploits the strengths of both these approaches while bypassing the pitfalls of either. Genomic islands lacking markers are identified by their association with genomic islands with markers. This was made possible by performing marker enrichment and phyletic pattern analyses within an integrated framework of recursive segmentation and clustering. The proposed method, IslandCafe, compared favorably with frequently used methods for genomic island detection on synthetic test datasets and on a test-set of known islands from 15 well-characterized bacterial species. Furthermore, IslandCafe identified novel islands with imprints of likely horizontal acquisition.

Current and Promising Approaches to Identify Horizontal Gene Transfer Events in Metagenomes

High-throughput shotgun metagenomics sequencing has enabled the profiling of myriad natural communities. These data are commonly used to identify gene families and pathways that were potentially gained or lost in an environment and which may be involved in microbial adaptation. Despite the widespread interest in these events, there are no established best practices for identifying gene gain and loss in metagenomics data. Horizontal gene transfer (HGT) represents several mechanisms of gene gain that are especially of interest in clinical microbiology due to the rapid spread of antibiotic resistance genes in natural communities. Several additional mechanisms of gene gain and loss, including gene duplication, gene loss-of-function events, and de novo gene birth are also important to consider in the context of metagenomes but have been less studied. This review is largely focused on detecting HGT in prokaryotic metagenomes, but methods for detecting these other mechanisms are first discussed. For this article to be self-contained, we provide a general background on HGT and the different possible signatures of this process. Lastly, we discuss how improved assembly of genomes from metagenomes would be the most straight-forward approach for improving the inference of gene gain and loss events. Several recent technological advances could help improve metagenome assemblies: long-read sequencing, determining the physical proximity of contigs, optical mapping of short sequences along chromosomes, and single-cell metagenomics. The benefits and limitations of these advances are discussed and open questions in this area are highlighted.

Long-term evolution of the natural isolate of Escherichia coli 536 in the mouse gut colonized after maternal transmission reveals convergence in the constitutive expression of the lactose operon

In vitro experimental evolution has taught us many lessons on the molecular bases of adaptation. To move towards more natural settings, evolution in the mice gut has been successfully performed. Yet, these experiments suffered from the use of laboratory strains as well as the use of axenic or streptomycin-treated mice to maintain the inoculated strains. To circumvent these limitations, we conducted a one-year experimental evolution in vivo using a natural isolate of E. coli, strain 536, in conditions mimicking as much as possible natural environment with mother-to-offspring microbiota transmission. Mice were then distributed in 24 independent cages and separated into two different diets: a regular one (chow diet, CD) and high-fat and high-sugar one (Western Diet, WD). Genome sequences revealed an early and rapid selection during the breastfeeding period that selected the constitutive expression of the well-characterized lactose operon. E. coli was lost significantly more in CD than WD; however, we could not detect any genomic signature of selection, nor any diet specificities during the later part of the experiments. The apparently neutral evolution presumably due to low population size maintained nevertheless at high frequency the early selected mutations affecting lactose regulation. The rapid loss of lactose operon regulation challenges the idea that plastic gene expression is both optimal and stable in the wild.

Alignment-free inference of hierarchical and reticulate phylogenomic relationships

We are amidst an ongoing flood of sequence data arising from the application of high-throughput technologies, and a concomitant fundamental revision in our understanding of how genomes evolve individually and within the biosphere. Workflows for phylogenomic inference must accommodate data that are not only much larger than before, but often more error prone and perhaps misassembled, or not assembled in the first place. Moreover, genomes of microbes, viruses and plasmids evolve not only by tree-like descent with modification but also by incorporating stretches of exogenous DNA. Thus, next-generation phylogenomics must address computational scalability while rethinking the nature of orthogroups, the alignment of multiple sequences and the inference and comparison of trees. New phylogenomic workflows have begun to take shape based on so-called alignment-free (AF) approaches. Here, we review the conceptual foundations of AF phylogenetics for the hierarchical (vertical) and reticulate (lateral) components of genome evolution, focusing on methods based on k-mers. We reflect on what seems to be successful, and on where further development is needed.

Prophage protein RacR activates lysozyme LysN, causing the growth defect of E. coli JM83

Prophage enriched the prokaryotic genome, and their transcriptional factors improved the protein expression network of the host. In this study, we uncovered a new prophage-prophage interaction in E. coli JM83. The Rac prophage protein RacR (GenBank accession no. AVI55875.1) directly activated the transcription of φ80dlacZΔM15 prophage lysozyme encoding gene 19 (GenBank accession no. ACB02445.1, renamed it lysN, lysozyme nineteen), resulting in the growth defect of JM83. This phenomenon also occurred in DH5α, but not in BL21(DE3) and MG1655 due to the genotype differences. However, deletion of lysN could not completely rescued JM83 from the growth arrest, indicating that RacR may regulate other related targets. In addition, passivation of RacR regulation was found in the late period of growth of JM83, and it was transmissible to daughter cells. Altogether, our study revealed part of RacR regulatory network, which suggested some advanced genetic strategies in bacteria.