High frequency of hotspot mutations in core genes of Escherichia coli due to short-term positive selection

Unveiling intraspecific diversity and evolutionary dynamics of the foodborne pathogen Bacillus paranthracis through high-quality pan-genome analysis

Understanding the evolutionary dynamics of foodborne pathogens throughout host-associated habitats is of utmost importance. Bacterial pan-genomes, as dynamic entities, are strongly influenced by ecological lifestyles. As a phenotypically diverse species in the Bacillus cereus group, Bacillus paranthracis is recognized as an emerging foodborne pathogen and a probiotic simultaneously. This poorly understood species is a suitable study model for adaptive pan-genome evolution. In this study, we determined the biogeographic distribution, abundance, genetic diversity, and genotypic profiles of key genetic elements of B. paranthracis. Metagenomic read recruitment analyses demonstrated that B. paranthracis members are globally distributed and abundant in host-associated habitats. A high-quality pan-genome of B. paranthracis was subsequently constructed to analyze the evolutionary dynamics involved in ecological adaptation comprehensively. The open pan-genome indicated a flexible gene repertoire with extensive genetic diversity. Significant divergences in the phylogenetic relationships, functional enrichment, and degree of selective pressure between the different components demonstrated different evolutionary dynamics between the core and accessory genomes driven by ecological forces. Purifying selection and gene loss are the main signatures of evolutionary dynamics in B. paranthracis pan-genome. The plasticity of the accessory genome is characterized by horizontal gene transfer (HGT), massive gene losses, and weak purifying or positive selection, which might contribute to niche-specific adaptation. In contrast, although the core genome dominantly undergoes purifying selection, its association with HGT and positively selected mutations indicates its potential role in ecological diversification. Furthermore, host fitness-related dynamics are characterized by the loss of secondary metabolite biosynthesis gene clusters (BGCs) and CAZyme-encoding genes and the acquisition of antimicrobial resistance (AMR) and virulence genes via HGT. This study offers a case study of pan-genome evolution to investigate the ecological adaptations reflected by biogeographical characteristics, thereby advancing the understanding of intraspecific diversity and evolutionary dynamics of foodborne pathogens.

Pan-genome insights into adaptive evolution of bacterial symbionts in mixed host-microbe symbioses represented by human gut microbiota Bacteroides cellulosilyticus

Animal hosts harbor diverse assemblages of microbial symbionts that play crucial roles in the host's lifestyle. The link between microbial symbiosis and host development remains poorly understood. In particular, little is known about the adaptive evolution of gut bacteria in host-microbe symbioses. Recently, symbiotic relationships have been categorized as open, closed, or mixed, reflecting their modes of inter-host transmission and resulting in distinct genomic features. Members of the genus Bacteroides are the most abundant human gut microbiota and possess both probiotic and pathogenic potential, providing an excellent model for studying pan-genome evolution in symbiotic systems. Here, we determined the complete genome of an novel clinical strain PL2022, which was isolated from a blood sample and performed pan-genome analyses on a representative set of Bacteroides cellulosilyticus strains to quantify the influence of the symbiotic relationship on the evolutionary dynamics. B. cellulosilyticus exhibited correlated genomic features with both open and closed symbioses, suggesting a mixed symbiosis. An open pan-genome is characterized by abundant accessory gene families, potential horizontal gene transfer (HGT), and diverse mobile genetic elements (MGEs), indicating an innovative gene pool, mainly associated with genomic islands and plasmids. However, massive parallel gene loss, weak purifying selection, and accumulation of positively selected mutations were the main drivers of genome reduction in B. cellulosilyticus. Metagenomic read recruitment analyses showed that B. cellulosilyticus members are globally distributed and active in human gut habitats, in line with predominant vertical transmission in the human gut. However, existence and/or high abundance were also detected in non-intestinal tissues, other animal hosts, and non-host environments, indicating occasional horizontal transmission to new niches, thereby creating arenas for the acquisition of novel genes. This case study of adaptive evolution under a mixed host-microbe symbiosis advances our understanding of symbiotic pan-genome evolution. Our results highlight the complexity of genetic evolution in this unusual intestinal symbiont.

Low-frequency and rare genetic variants associated with rheumatoid arthritis risk

Rheumatoid arthritis (RA) has an estimated heritability of nearly 50%, which is particularly high in seropositive RA. HLA alleles account for a large proportion of this heritability, in addition to many common single-nucleotide polymorphisms with smaller individual effects. Low-frequency and rare variants, such as those captured by next-generation sequencing, can also have a large role in heritability in some individuals. Rare variant discovery has informed the development of drugs such as inhibitors of PCSK9 and Janus kinases. Some 34 low-frequency and rare variants are currently associated with RA risk. One variant (19:10352442G>C in TYK2) was identified in five separate studies, and might therefore represent a promising therapeutic target. Following a set of best practices in future studies, including studying diverse populations, using large sample sizes, validating RA and serostatus, replicating findings, adjusting for other variants and performing functional assessment, could help to ensure the relevance of identified variants. Exciting opportunities are now on the horizon for genetics in RA, including larger datasets and consortia, whole-genome sequencing and direct applications of findings in the management, and especially treatment, of RA.

Gene as a dynamical notion: An extensive and integrative vision. Redefining the gene concept, from traditional to genic-interaction, as a new dynamical version

The current concept of gene has been very useful during the 20th and 21st centuries. However, recent advances in molecular biology and bioinformatics, which have further diversified the functional and adaptive profile of genetic information and its integration with cell physiology and environmental response, have contributed to focusing on additional new gene properties besides the traditional definition. Considering the inherent complexity of gene expression, whose adaptive objective must be referred to the Tortoise-Hare model, in which two tendencies converge, one focused on rapid adaptation to achieve survival, and the other that prevents an over-adaptation effect. In this context, a revision of the gene concept must be made, which must include these new mechanisms and approaches. In this paper, we propose a new conception of the idea of a gene that moves from a static and defined version of hereditary information to a dynamic idea that preponderates gene interaction (circumscribed to that established between protein-protein, protein-nucleic acid, and nucleic acid-nucleic acid) and the selection it exerts, as the irreducible element that works in a coordinated way in a genomic regulatory network (GRN).

Mutational hotspots lead to robust but suboptimal adaptive outcomes in certain environments

The observed mutational spectrum of adaptive outcomes can be constrained by many factors. For example, mutational biases can narrow the observed spectrum by increasing the rate of mutation at isolated sites in the genome. In contrast, complex environments can shift the observed spectrum by defining fitness consequences of mutational routes. We investigate the impact of different nutrient environments on the evolution of motility in Pseudomonas fluorescens Pf0-2x (an engineered non-motile derivative of Pf0-1) in the presence and absence of a strong mutational hotspot. Previous work has shown that this mutational hotspot can be built and broken via six silent mutations, which provide rapid access to a mutation that rescues swimming motility and confers the strongest swimming phenotype in specific environments. Here, we evolved a hotspot and non-hotspot variant strain of Pf0-2x for motility under nutrient-rich (LB) and nutrient-limiting (M9) environmental conditions. We observed the hotspot strain consistently evolved faster across all environmental conditions and its mutational spectrum was robust to environmental differences. However, the non-hotspot strain had a distinct mutational spectrum that changed depending on the nutrient environment. Interestingly, while alternative adaptive mutations in nutrient-rich environments were equal to, or less effective than, the hotspot mutation, the majority of these mutations in nutrient-limited conditions produced superior swimmers. Our competition experiments mirrored these findings, underscoring the role of environment in defining both the mutational spectrum and the associated phenotype strength. This indicates that while mutational hotspots working in concert with natural selection can speed up access to robust adaptive mutations (which can provide a competitive advantage in evolving populations), they can limit exploration of the mutational landscape, restricting access to potentially stronger phenotypes in specific environments.

Commensal Fitness Advantage May Contribute to the Global Dissemination of Multidrug-Resistant Lineages of Bacteria-The Case of Uropathogenic E. coli

It is widely accepted that favorable fitness in commensal colonization is one of the prime facilitators of clonal dissemination in bacteria. The question arises as to what kind of fitness advantage may be wielded by uropathogenic strains of the two predominant fluoroquinolone- and multidrug-resistant clonal groups of E. coli-ST131-H30 and ST1193, which has permitted their unprecedented pandemic-like global expansion in the last few decades. The colonization-associated genes' content, carriage of low-cost plasmids, and integrons with weak promoters could certainly contribute to the fitness of the pandemic groups, although those genetic factors are common among other clonal groups as well. Also, ST131-H30 and ST1193 strains harbor fluoroquinolone-resistance conferring mutations targeting serine residues in DNA gyrase (GyrA-S83) and topoisomerase IV (ParC-S80) that, in those clonal backgrounds, might result in a commensal fitness benefit, i.e., beyond the antibiotic resistance per se. This fitness gain might have contributed not only to the widespread dissemination of these major clones in the healthcare setting but also to their long-term colonization of healthy individuals and, thus, circulation in the community, even in a low or no fluoroquinolone use environment. This evolutionary shift affecting commensal E. coli, initiated by mutations co-favorable in both antibiotics-treated patients and healthy individuals warrants more in-depth studies to monitor further changes in the epidemiological situation and develop effective measures to reduce the antibiotic resistance spread.

Reduced polymorphism of Plasmodium vivax early transcribed membrane protein (PvETRAMP) 11.2

BACKGROUND

ETRAMP11.2 (PVX_003565) is a well-characterized protein with antigenic potential. It is considered to be a serological marker for diagnostic tools, and it has been suggested as a potential vaccine candidate. Despite its immunological relevance, the polymorphism of the P. vivax ETRAMP11.2 gene (pvetramp11.2) remains undefined. The genetic variability of an antigen may limit the effectiveness of its application as a serological surveillance tool and in vaccine development and, therefore, the aim of this study was to investigate the genetic diversity of pvetramp11.2 in parasite populations from Amazonian regions and worldwide. We also evaluated amino acid polymorphism on predicted B-cell epitopes. The low variability of the sequence encoding PvETRAMP11.2 protein suggests that it would be a suitable marker in prospective serodiagnostic assays for surveillance strategies or in vaccine design against P. vivax malaria.

METHODS

The pvetramp11.2 of P. vivax isolates collected from Brazil (n = 68) and Peru (n = 36) were sequenced and analyzed to assess nucleotide polymorphisms, allele distributions, population differentiation, genetic diversity and signature of selection. In addition, sequences (n = 104) of seven populations from different geographical regions were retrieved from the PlasmoDB database and included in the analysis to study the worldwide allele distribution. Potential linear B-cell epitopes and their polymorphisms were also explored.

RESULTS

The multiple alignments of 208 pvetramp11.2 sequences revealed a low polymorphism and a marked geographical variation in allele diversity. Seven polymorphic sites and 11 alleles were identified. All of the alleles were detected in isolates from the Latin American region and five alleles were detected in isolates from the Southeast Asia/Papua New Guinea (SEA/PNG) region. Three alleles were shared by all Latin American populations (H1, H6 and H7). The H1 allele (reference allele from Salvador-1 strain), which was absent in the SEA/PNG populations, was the most represented allele in populations from Brazil (54%) and was also detected at high frequencies in populations from all other Latin America countries (range: 13.0% to 33.3%). The H2 allele was the major allele in SEA/PNG populations, but was poorly represented in Latin America populations (only in Brazil: 7.3%). Plasmodium vivax populations from Latin America showed a marked inter-population genetic differentiation (fixation index [Fst]) in contrast to SEA/PNG populations. Codon bias measures (effective number of codons [ENC] and Codon bias index [CBI]) indicated preferential use of synonymous codons, suggesting selective pressure at the translation level. Only three amino acid substitutions, located in the C-terminus, were detected. Linear B-cell epitope mapping predicted two epitopes in the Sal-1 PvETRAMP11.2 protein, one of which was fully conserved in all of the parasite populations analyzed.

CONCLUSIONS

We provide an overview of the allele distribution and genetic differentiation of ETRAMP11.2 antigen in P. vivax populations from different endemic areas of the world. The reduced polymorphism and the high degree of protein conservation supports the application of PvETRAMP11.2 protein as a reliable antigen for application in serological assays or vaccine design. Our findings provide useful information that can be used to inform future study designs.

Pan-Chromosome and Comparative Analysis of Agrobacterium fabrum Reveal Important Traits Concerning the Genetic Diversity, Evolutionary Dynamics, and Niche Adaptation of the Species

Agrobacterium fabrum has been critical for the development of plant genetic engineering and agricultural biotechnology due to its ability to transform eukaryotic cells. However, the gene composition, evolutionary dynamics, and niche adaptation of this species is still unknown. Therefore, we established a comparative genomic analysis based on a pan-chromosome data set to evaluate the genetic diversity of A. fabrum. Here, 25 A. fabrum genomes were selected for analysis by core genome phylogeny combined with the average nucleotide identity (ANI), amino acid identity (AAI), and in silico DNA-DNA hybridization (DDH) values. An open pan-genome of A. fabrum exhibits genetic diversity with variable accessorial genes as evidenced by a consensus pan-genome of 12 representative genomes. The genomic plasticity of A. fabrum is apparent in its putative sequences for mobile genetic elements (MGEs), limited horizontal gene transfer barriers, and potentially horizontally transferred genes. The evolutionary constraints and functional enrichment in the pan-chromosome were measured by the Clusters of Orthologous Groups (COG) categories using eggNOG-mapper software, and the nonsynonymous/synonymous rate ratio (dN/dS) was determined using HYPHY software. Comparative analysis revealed significant differences in the functional enrichment and the degree of purifying selection between the core genome and non-core genome. We demonstrate that the core gene families undergo stronger purifying selection but have a significant bias to contain one or more positively selected sites. Furthermore, although they shared similar genetic diversity, we observed significant differences between chromosome 1 (Chr I) and the chromid in their functional features and evolutionary constraints. We demonstrate that putative genetic elements responsible for plant infection, ecological adaptation, and speciation represent the core genome, highlighting their importance in the adaptation of A. fabrum to plant-related niches. Our pan-chromosome analysis of A. fabrum provides comprehensive insights into the genetic properties, evolutionary patterns, and niche adaptation of the species. IMPORTANCE Agrobacterium spp. live in diverse plant-associated niches such as soil, the rhizosphere, and vegetation, which are challenged by multiple stressors such as diverse energy sources, plant defenses, and microbial competition. They have evolved the ability to utilize diverse resources, escape plant defenses, and defeat competitors. However, the underlying genetic diversity and evolutionary dynamics of Agrobacterium spp. remain unexplored. We examined the phylogeny and pan-genome of A. fabrum to define intraspecies evolutionary relationships. Our results indicate an open pan-genome and numerous MGEs and horizontally transferred genes among A. fabrum genomes, reflecting the flexibility of the chromosomes and the potential for genetic exchange. Furthermore, we observed significant differences in the functional features and evolutionary constraints between the core and accessory genomes and between Chr I and the chromid, respectively.

Pitfalls of genotyping microbial communities with rapidly growing genome collections

Detecting genetic variants in metagenomic data is a priority for understanding the evolution, ecology, and functional characteristics of microbial communities. Many tools that perform this metagenotyping rely on aligning reads of unknown origin to a database of sequences from many species before calling variants. In this synthesis, we investigate how databases of increasingly diverse and closely related species have pushed the limits of current alignment algorithms, thereby degrading the performance of metagenotyping tools. We identify multi-mapping reads as a prevalent source of errors and illustrate a trade-off between retaining correct alignments versus limiting incorrect alignments, many of which map reads to the wrong species. Then we evaluate several actionable mitigation strategies and review emerging methods showing promise to further improve metagenotyping in response to the rapid growth in genome collections. Our results have implications beyond metagenotyping to the many tools in microbial genomics that depend upon accurate read mapping.

The synergistic activity of SBC3 in combination with Ebselen against Escherichia coli infection

Escherichia coli ranks as the number one clinical isolate in the past years in China according to The China Antimicrobial Surveillance Network (CHINET), and its multidrug-resistant (MDR) pathogenic strains account for over 160 million cases of dysentery and one million deaths per year. Here, our work demonstrates that E. coli is highly sensitive to the synergistic combination of SBC3 [1,3-Dibenzyl-4,5-diphenyl-imidazol-2-ylidene silver (I) acetate] and Ebselen, which shows no synergistic toxicity on mammalian cells. The proposed mechanism for the synergistic antibacterial effect of SBC3 in combination with Ebselen is based on directly inhibiting E. coli thioredoxin reductase and rapidly depleting glutathione, resulting in the increase of reactive oxygen species that cause bacterial cell death. Furthermore, the bactericidal efficacy of SBC3 in combination with Ebselen has been confirmed in mild and acute peritonitis mice. In addition, the five most difficult to treat Gram-negative bacteria (including E. coli, Acinetobacter baumannii, Enterobacter cloacae, Klebsiella pneumoniae, and Pseudomonas aeruginosa) are also highly sensitive to a synergistic combination of SBC3 and Ebselen. Thus, SBC3 in combination with Ebselen has potential as a treatment for clinically important Gram-negative bacterial infections.

Genomics and pathotypes of the many faces of Escherichia coli

Escherichia coli is the most researched microbial organism in the world. Its varied impact on human health, consisting of commensalism, gastrointestinal disease, or extraintestinal pathologies, has generated a separation of the species into at least eleven pathotypes (also known as pathovars). These are broadly split into two groups, intestinal pathogenic E. coli (InPEC) and extraintestinal pathogenic E. coli (ExPEC). However, components of E. coli's infinite open accessory genome are horizontally transferred with substantial frequency, creating pathogenic hybrid strains that defy a clear pathotype designation. Here, we take a birds-eye view of the E. coli species, characterizing it from historical, clinical, and genetic perspectives. We examine the wide spectrum of human disease caused by E. coli, the genome content of the bacterium, and its propensity to acquire, exchange, and maintain antibiotic resistance genes and virulence traits. Our portrayal of the species also discusses elements that have shaped its overall population structure and summarizes the current state of vaccine development targeted at the most frequent E. coli pathovars. In our conclusions, we advocate streamlining efforts for clinical reporting of ExPEC, and emphasize the pathogenic potential that exists throughout the entire species.

Distinct evolutionary trajectories in the Escherichia coli pangenome occur within sequence types

The Escherichia coli species contains a diverse set of sequence types and there remain important questions regarding differences in genetic content within this population that need to be addressed. Pangenomes are useful vehicles for studying gene content within sequence types. Here, we analyse 21 E. coli sequence type pangenomes using comparative pangenomics to identify variance in both pangenome structure and content. We present functional breakdowns of sequence type core genomes and identify sequence types that are enriched in metabolism, transcription and cell membrane biogenesis genes. We also uncover metabolism genes that have variable core classification, depending on which allele is present. Our comparative pangenomics approach allows for detailed exploration of sequence type pangenomes within the context of the species. We show that ongoing gene gain and loss in the E. coli pangenome is sequence type-specific, which may be a consequence of distinct sequence type-specific evolutionary drivers.

Signatures of selection in core and accessory genomes indicate different ecological drivers of diversification among Bacillus cereus clades

Bacterial clades are often ecologically distinct, despite extensive horizontal gene transfer (HGT). How selection works on different parts of bacterial pan-genomes to drive and maintain the emergence of clades is unclear. Focusing on the three largest clades in the diverse and well-studied Bacillus cereus sensu lato group, we identified clade-specific core genes (present in all clade members) and then used clade-specific allelic diversity to identify genes under purifying and diversifying selection. Clade-specific accessory genes (present in a subset of strains within a clade) were characterized as being under selection using presence/absence in specific clades. Gene ontology analyses of genes under selection revealed that different gene functions were enriched in different clades. Furthermore, some gene functions were enriched only amongst clade-specific core or accessory genomes. Genes under purifying selection were often clade-specific, while genes under diversifying selection showed signs of frequent HGT. These patterns are consistent with different selection pressures acting on both the core and the accessory genomes of different clades and can lead to ecological divergence in both cases. Examining variation in allelic diversity allows us to uncover genes under clade-specific selection, allowing ready identification of strains and their ecological niche.

A critical evaluation of Mycobacterium bovis pangenomics, with reference to its utility in outbreak investigation

The increased accessibility of next generation sequencing has allowed enough genomes from a given bacterial species to be sequenced to describe the distribution of genes in the pangenome, without limiting analyses to genes present in reference strains. Although some taxa have thousands of whole genome sequences available on public databases, most genomes were sequenced with short read technology, resulting in incomplete assemblies. Studying pangenomes could lead to important insights into adaptation, pathogenicity, or molecular epidemiology, however given the known information loss inherent in analyzing contig-level assemblies, these inferences may be biased or inaccurate. In this study we describe the pangenome of a clonally evolving pathogen, Mycobacterium bovis , and examine the utility of gene content variation in M. bovis outbreak investigation. We constructed the M. bovis pangenome using 1463 de novo assembled genomes. We tested the assumption of strict clonal evolution by studying evidence of recombination in core genes and analyzing the distribution of accessory genes among core monophyletic groups. To determine if gene content variation could be utilized in outbreak investigation, we carefully examined accessory genes detected in a well described M. bovis outbreak in Minnesota. We found significant errors in accessory gene classification. After accounting for these errors, we show that M. bovis has a much smaller accessory genome than previously described and provide evidence supporting ongoing clonal evolution and a closed pangenome, with little gene content variation generated over outbreaks. We also identified frameshift mutations in multiple genes, including a mutation in glpK , which has recently been associated with antibiotic tolerance in Mycobacterium tuberculosis . A pangenomic approach enables a more comprehensive analysis of genome dynamics than is possible with reference-based approaches; however, without critical evaluation of accessory gene content, inferences of transmission patterns employing these loci could be misguided.

Distinct Potentially Adaptive Accumulation of Truncation Mutations in Salmonella enterica serovar Typhi and Salmonella enterica serovar Paratyphi A

Gene inactivation through the accumulation of truncation (or premature stop codon) mutations is a common mode of evolution in bacteria. It is frequently believed to result from reductive evolutionary processes allowing purging of superfluous traits. However, several works have demonstrated that, similar to the occurrences of inactivating nonsynonymous (i.e., amino acid replacement) mutations under positive selection pressures, truncation mutations can also be adaptive where specific traits deleterious in particular environmental conditions need to be inactivated for survival. Here, we performed a comparative analysis of genome-wide accumulation of truncation mutations in Salmonella enterica serovar Typhi and Salmonella enterica serovar Paratyphi A. Considering the known convergent evolutionary trajectories in these two serovars, we expected a strong overlap of truncated genes in S. Typhi and S. Paratyphi A, emerging through either reductive or adaptive dynamics. However, we detected a distinct set of core truncated genes encoding different overrepresented functional clusters in each serovar. In 54% and 28% truncated genes in S. Typhi and S. Paratyphi A, respectively, inactivating mutations were acquired by only different subsets of isolates, instead of all isolates analyzed for that serovar. Importantly, 62% truncated genes (P < 0.001) in S. Typhi and S. Paratyphi A were also targeted by convergent amino acid mutations in different serovars, suggesting those genes to be under selection pressures. Our findings indicate significant presence of potentially adaptive truncation mutations in conjunction with the ones emerging due to reductive evolution. Further experimental and large-scale bioinformatic studies are necessary to better explore the impact of such adaptive footprints of truncation mutations in the evolution of bacterial virulence. IMPORTANCE Detecting the adaptive mutations leading to gene inactivation or loss of function is crucial for understanding their contribution in the evolution of bacterial virulence and antibiotic resistance. Such inactivating mutations, apart from being of nonsynonymous (i.e., amino acid replacement) nature, can also be truncation mutations, abruptly trimming the length of encoded proteins. Importantly, the notion of reductive evolutionary dynamics is primarily accepted toward the accumulation of truncation mutations. However, our case study on S. Typhi and S. Paratyphi A, two human-restricted systemically invasive pathogens exerting similar clinical manifestations, indicated that a significant proportion of truncation mutations emerge from positive selection pressures. The candidate genes from our study will enable directed functional assays for deciphering the adaptive role of truncation mutations in S. Typhi and S. Paratyphi A pathogenesis. Also, our genome-level analytical approach will pave the way to understand the contribution of truncation mutations in the adaptive evolution of other bacterial pathogens.

Fast and accurate metagenotyping of the human gut microbiome with GT-Pro

Single nucleotide polymorphisms (SNPs) in metagenomics are used to quantify population structure, track strains and identify genetic determinants of microbial phenotypes. However, existing alignment-based approaches for metagenomic SNP detection require high-performance computing and enough read coverage to distinguish SNPs from sequencing errors. To address these issues, we developed the GenoTyper for Prokaryotes (GT-Pro), a suite of methods to catalog SNPs from genomes and use unique k-mers to rapidly genotype these SNPs from metagenomes. Compared to methods that use read alignment, GT-Pro is more accurate and two orders of magnitude faster. Using high-quality genomes, we constructed a catalog of 104 million SNPs in 909 human gut species and used unique k-mers targeting this catalog to characterize the global population structure of gut microbes from 7,459 samples. GT-Pro enables fast and memory-efficient metagenotyping of millions of SNPs on a personal computer.

Assessment of selection pressure exerted on genes from complete pangenomes helps to improve the accuracy in the prediction of new genes

Bacterial genomes are massively sequenced, and they provide valuable data to better know the complete set of genes of a species. The analysis of thousands of bacterial strains can identify both shared genes and those appearing only in the pathogenic ones. Current computational gene finders facilitate this task but often miss some existing genes. However, the present availability of different genomes from the same species is useful to estimate the selective pressure applied on genes of complete pangenomes. It may assist in evaluating gene predictions either by checking the certainty of a new gene or annotating it as a gene under positive selection. Here, we estimated the selective pressure of 19 271 genes that are part of the pangenome of the human opportunistic pathogen Acinetobacter baumannii and found that most genes in this bacterium are subject to negative selection. However, 23% of them showed values compatible with positive selection. These latter were mainly uncharacterized proteins or genes required to evade the host defence system including genes related to resistance and virulence whose changes may be favoured to acquire new functions. Finally, we evaluated the utility of measuring selection pressure in the detection of sequencing errors and the validation of gene prediction.

Genome Plasticity of African Swine Fever Virus: Implications for Diagnostics and Live-Attenuated Vaccines

African swine fever (ASF) is a highly contagious transboundary viral hemorrhagic disease of domestic and wild pigs presenting a significant threat to the global swine industry. Following its introduction in Caucasus, Georgia, in 2007, the genome of the genotype II of African swine fever virus (ASFV) strain Georgia-07 and its derivatives accumulated significant mutations, resulting in the emergence of genetic variants within short epidemiological timescales as it spreads and infects different hosts in diverse ecosystems, causing outbreaks in Europe, South Asia, South East Asia and the Caribbean. This suggests that ASFV, with a comparatively large and complex DNA genome, is susceptible to genetic mutations including deletions and that although the virus is environmentally stable, it is genetically unstable. This has implications for the development of vaccines and diagnostic tests for disease detection and surveillance. Analysis of the ASFV genome revealed recombination hotspots, which in double-stranded DNA (dsDNA) viruses represent key drivers of genetic diversity. The ability of pox virus, a dsDNA virus with a genome complexity similar to ASFV, regaining virulence following the deletion of a virulence gene via gene amplification, coupled with the recent emergence and spread of live-attenuated ASFV vaccine strains causing disease and death in pigs in China, raise legitimate concerns around the use of live-attenuated ASFV vaccines in non-endemic regions to control the potential introduction. Further research into the risk of using live-attenuated ASFV in non-endemic regions is highly needed.

Comparative phylogenetic analysis of SARS-CoV-2 spike protein-possibility effect on virus spillover

Coronavirus disease 2019 has developed into a dramatic pandemic with tremendous global impact. The receptor-binding motif (RBM) region of the causative virus, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), binds to host angiotensin-converting enzyme 2 (ACE2) receptors for infection. As ACE2 receptors are highly conserved within vertebrate species, SARS-CoV-2 can infect significant animal species as well as human populations. An analysis of SARS-CoV-2 genotypes isolated from human and significant animal species was conducted to compare and identify mutation and adaptation patterns across different animal species. The phylogenetic data revealed seven distinct phylogenetic clades with no significant relationship between the clades and geographical locations. A high rate of variation within SARS-CoV-2 mink isolates implies that mink populations were infected before human populations. Positions of most single-nucleotide polymorphisms (SNPs) within the spike (S) protein of SARS-CoV-2 genotypes from the different hosts are mostly accumulated in the RBM region and highlight the pronounced accumulation of variants with mutations in the RBM region in comparison with other variants. These SNPs play a crucial role in viral transmission and pathogenicity and are keys in identifying other animal species as potential intermediate hosts of SARS-CoV-2. The possible roles in the emergence of new viral strains and the possible implications of these changes, in compromising vaccine effectiveness, deserve urgent considerations.

Detection of novel sequence types and zoonotic transmission potentiality among strains of Shiga toxigenic Escherichia coli (STEC) from dairy calves, animal handlers and associated environments

Shiga toxigenic Escherichia coli (STEC) is one of the most important food-borne zoonotic bacterial pathogens responsible for causing gastrointestinal infections, haemorrhagic colitis and haemolytic uremic syndrome. The present study was aimed to isolate and characterize STEC from neonatal dairy calves, animal handlers and their surrounding environment and to establish the genetic relationship among isolates by multilocus sequence typing (MLST). A total number of 115 samples were collected and processed for the isolation of E. coli. The occurrence rate of E. coli was 92.2% (106/115), of which, 18 were typed as STEC. Antibacterial susceptibility analysis revealed 11 (61.1%) strains as multiple drug-resistant (MDR). MLST analysis has delineated 16 sequence types (STs) including nine novel STs. Among STs, ST58 dominated with three strains and was recovered from the environment and neonatal calves. Strains from neonatal calves and humans showed genetic relatedness with significant bootstrap support values indicative of zoonotic transmission potentiality. Analysis of 211 global isolates belonging to 61 STs indicated predominant STs (ST 21, ST 33 and ST 3416) that can be either host-specific (ST 33 and ST 3416) or can be shared among human and bovine hosts (ST 21). The MLST analysis indicates genetic relatedness among isolates and the results predispose inter-host transmission and zoonotic spread.

Protein-Coding Genes of Helicobacter pylori Predominantly Present Purifying Selection though Many Membrane Proteins Suffer from Selection Pressure: A Proposal to Analyze Bacterial Pangenomes

The current availability of complete genome sequences has allowed knowing that bacterial genomes can bear genes not present in the genome of all the strains from a specific species. So, the genes shared by all the strains comprise the core of the species, but the pangenome can be much greater and usually includes genes appearing in one only strain. Once the pangenome of a species is estimated, other studies can be undertaken to generate new knowledge, such as the study of the evolutionary selection for protein-coding genes. Most of the genes of a pangenome are expected to be subject to purifying selection that assures the conservation of function, especially those in the core group. However, some genes can be subject to selection pressure, such as genes involved in virulence that need to escape to the host immune system, which is more common in the accessory group of the pangenome. We analyzed 180 strains of Helicobacter pylori, a bacterium that colonizes the gastric mucosa of half the world population and presents a low number of genes (around 1500 in a strain and 3000 in the pangenome). After the estimation of the pangenome, the evolutionary selection for each gene has been calculated, and we found that 85% of them are subject to purifying selection and the remaining genes present some grade of selection pressure. As expected, the latter group is enriched with genes encoding for membrane proteins putatively involved in interaction to host tissues. In addition, this group also presents a high number of uncharacterized genes and genes encoding for putative spurious proteins. It suggests that they could be false positives from the gene finders used for identifying them. All these results propose that this kind of analyses can be useful to validate gene predictions and functionally characterize proteins in complete genomes.

Genomic analysis of trimethoprim-resistant extraintestinal pathogenic Escherichia coli and recurrent urinary tract infections

Urinary tract infections (UTIs) are the most common bacterial infections requiring medical attention and a leading justification for antibiotic prescription. Trimethoprim is prescribed empirically for uncomplicated cases. UTIs are primarily caused by extraintestinal pathogenic Escherichia coli (ExPEC) and ExPEC strains play a central role in disseminating antimicrobial-resistance genes worldwide. Here, we describe the whole-genome sequences of trimethoprim-resistant ExPEC and/or ExPEC from recurrent UTIs (67 in total) from patients attending a regional Australian hospital from 2006 to 2008. Twenty-three sequence types (STs) were observed, with ST131 predominating (28 %), then ST69 and ST73 (both 7 %). Co-occurrence of trimethoprim-resistance genes with genes conferring resistance to extended-spectrum β-lactams, heavy metals and quaternary ammonium ions was a feature of the ExPEC described here. Seven trimethoprim-resistance genes were identified, most commonly dfrA17 (38 %) and dfrA12 (18 %). An uncommon dfrB4 variant was also observed. Two blaCTX-M variants were identified - blaCTX-M-15 (16 %) and blaCTX-M-14 (10 %). The former was always associated with dfrA12, the latter with dfrA17, and all blaCTX-M genes co-occurred with chromate-resistance gene chrA. Eighteen class 1 integron structures were characterized, and chrA featured in eight structures; dfrA genes featured in seventeen. ST131 H30Rx isolates possessed distinct antimicrobial gene profiles comprising aac(3)-IIa, aac(6)-Ib-cr, aph(3')-Ia, aadA2, blaCTX-M-15, blaOXA-1 and dfrA12. The most common virulence-associated genes (VAGs) were fimH, fyuA, irp2 and sitA (all 91 %). Virulence profile clustering showed ST131 H30 isolates carried similar VAGs to ST73, ST405, ST550 and ST1193 isolates. The sole ST131 H27 isolate carried molecular predictors of enteroaggregative E. coli/ExPEC hybrid strains (aatA, aggR, fyuA). Seven isolates (10 %) carried VAGs suggesting ColV plasmid carriage. Finally, SNP analysis of serial UTI patients experiencing worsening sequelae demonstrated a high proportion of point mutations in virulence factors.

Genomic Diversity and Hotspot Mutations in 30,983 SARS-CoV-2 Genomes: Moving Toward a Universal Vaccine for the "Confined Virus"?

The COVID-19 pandemic has been ongoing since its onset in late November 2019 in Wuhan, China. Understanding and monitoring the genetic evolution of the virus, its geographical characteristics, and its stability are particularly important for controlling the spread of the disease and especially for the development of a universal vaccine covering all circulating strains. From this perspective, we analyzed 30,983 complete SARS-CoV-2 genomes from 79 countries located in the six continents and collected from 24 December 2019, to 13 May 2020, according to the GISAID database. Our analysis revealed the presence of 3206 variant sites, with a uniform distribution of mutation types in different geographic areas. Remarkably, a low frequency of recurrent mutations has been observed; only 169 mutations (5.27%) had a prevalence greater than 1% of genomes. Nevertheless, fourteen non-synonymous hotspot mutations (>10%) have been identified at different locations along the viral genome; eight in ORF1ab polyprotein (in nsp2, nsp3, transmembrane domain, RdRp, helicase, exonuclease, and endoribonuclease), three in nucleocapsid protein, and one in each of three proteins: Spike, ORF3a, and ORF8. Moreover, 36 non-synonymous mutations were identified in the receptor-binding domain (RBD) of the spike protein with a low prevalence (<1%) across all genomes, of which only four could potentially enhance the binding of the SARS-CoV-2 spike protein to the human ACE2 receptor. These results along with intra-genomic divergence of SARS-CoV-2 could indicate that unlike the influenza virus or HIV viruses, SARS-CoV-2 has a low mutation rate which makes the development of an effective global vaccine very likely.

Adaptation to host-specific bacterial pathogens drive rapid evolution of novel PhoP/PhoQ regulation pathway modulating the virulence

The presence of the PhoP-PhoQ system is usually different in various bacterial groups, suggesting that PhoP can control the expression of different genes in species. However, little is known about the evolution of the PhoP-PhoQ system among bacterial pathogens. Here, we study the evolution of PhoP and PhoQ regulation in 15 species of Enterobacteriaceae family. We have determined that the regulatory objectives adopted by PhoP and PhoQ are mainly different, due to the result of horizontal gene transfer events and even the change in the genetic content between closely related species. We have compared many possibilities tests (M1 vs. M2 and M7 with M8) to determine the positive selection. Estimating parameters at M1 and M2, with positive selection in M2 of the two proteins. The proportions of positive selection sites significant with ω = 4.53076 for PhoP and ω = 4.21041 PhQ. M8 was significant for PhoP and PhQ proteins. To further confirm the positive selection results, we used the Selecton server to confer positive selection on individual sites using the Mechanistic-Empirical Combination model, and we noticed that several sites had been identified under selection pressure during the evolution. There was a strong indication for the positive selection in bacterial genes of PhoP and PhoQ showed the results. By the use of REL and IFEL, the positive selection for PhoP was detected 14 and 11 sites respectively at different codon positions. The positively selected sites of amino acids such as Arginine, Alanine, Lysine, and Leucine are more important for the production of signals. Our results suggest that the positive selection of PhoP-PhoQ genes in host adaptation during evolution raises an intriguing possibility causes subtle variations in actions of PhoP-PhoQ and also increases the opportunities that cause modification in protein structure for the evolution of increasing pathogenicity in bacterial pathogens.

Prevalence of hypermutator isolates of Achromobacter spp. from cystic fibrosis patients

Bacteria colonising the lungs of cystic fibrosis (CF) patients encounter high selective pressures. Hypermutation facilitates adaptation to fluctuating environments, and hypermutator strains are frequently isolated from CF patients. We investigated the prevalence of hypermutator isolates of Achromobacter spp. among patients affiliated with the CF Centre in Aarhus, Denmark. By exposure to rifampicin, the mutation frequency was determined for 90 isolates of Achromobacter spp. cultured from 42 CF patients; 20 infections were categorised as chronic, 22 as intermittent. The genetic mechanisms of hypermutation were examined by comparing DNA repair gene sequences from hypermutator and normomutator isolates. Achromobacter spp. cultured from 11 patients were categorised as hypermutators, and this phenotype was exclusively associated with chronic infections. Isolates of the Danish epidemic strain (DES) of Achromobacter ruhlandii cultured from patients from both Danish CF centres showed elevated mutation frequencies. The hypermutator state of Achromobacter spp. was most commonly associated with nonsynonymous mutations in the DNA mismatch repair gene mutS; a single clone had developed a substitution in the S-adenosyl-L-methionine-dependent methyltransferase putatively involved in DNA repair mechanisms, but not previously linked to the hypermutator phenotype. Hypermutation is prevalent among clinical isolates of Achromobacter spp. and could be a key determinant for the extraordinary adaptation and persistence of DES.

Inactivation of a Mismatch-Repair System Diversifies Genotypic Landscape of Escherichia coli During Adaptive Laboratory Evolution

Adaptive laboratory evolution (ALE) is used to find causal mutations that underlie improved strain performance under the applied selection pressure. ALE studies have revealed that mutator populations tend to outcompete their non-mutator counterparts following the evolutionary trajectory. Among them, mutS-inactivated mutator cells, characterize d by a dysfunctional methyl-mismatch repair system, are frequently found in ALE experiments. Here, we examined mutS inactivation as an approach to facilitate ALE of Escherichia coli. The wild-type E. coli MG1655 and mutS knock-out derivative (ΔmutS) were evolved in parallel for 800 generations on lactate or glycerol minimal media in a serial-transfer experiment. Whole-genome re-sequencing of each lineage at 100-generation intervals revealed that (1) mutations emerge rapidly in the ΔmutS compared to in the wild-type strain; (2) mutations were more than fourfold higher in the ΔmutS strain at the end-point populations compared to the wild-type strain; and (3) a significant number of random mutations accumulated in the ΔmutS strains. We then measured the fitness of the end-point populations on an array of non-adaptive carbon sources. Interestingly, collateral fitness increases on non-adaptive carbon sources were more pronounced in the ΔmutS strains than the parental strain. Fitness measurement of single mutants revealed that the collateral fitness increase seen in the mutator lineages can be attributed to a pool of random mutations. Together, this study demonstrates that short-term mutator ALE extensively expands possible genotype space, resulting in versatile bacteria with elevated fitness levels across various carbon sources.

Mutational convergence acts as a major player in adaptive parallel evolution of Shigella spp

Shigella spp., emerging from multiple origins of Escherichia coli, poses a significant health threat as a causative agent of bacillary dysentery. While multiple serotypes of four different species have evolved via independent lineages, Shigella spp. are designated as a single pathotype, primarily because of their common mode of pathogenesis. Convergent horizontal transfer events have so far been attributed to the commonalities in the evolution of virulence across diverse lineages. However, the role of mutational convergence in such parallel evolution is not yet well understood. Here we have carried out a genome-wide analysis of Shigella strains from all four species to detect the core genes (i.e. the ones present in all analyzed strains) acquiring convergent mutations of evolutionarily recent origin. Simulation studies show non-neutral accumulation of these convergent mutations across species, suggesting their adaptive role in the evolution of Shigella virulence. S. dysenteriae strain 197, representing highly virulent type 1 (Sd1) clone, carries excessively high number of core genes with recent convergent mutations compared to other analyzed strains. We propose that this high frequency of adaptive convergence in S. dysenteriae strain 197 could be linked to recent re-emergence of the Sd1 clone and its increased resistance to antimicrobials.

Whole-genome sequencing of bloodstream Staphylococcus aureus isolates does not distinguish bacteraemia from endocarditis

Most Staphylococcus aureus isolates can cause invasive disease given the right circumstances, but it is unknown if some isolates are more likely to cause severe infections than others. S. aureus bloodstream isolates from 120 patients with definite infective endocarditis and 121 with S. aureus bacteraemia without infective endocarditis underwent whole-genome sequencing. Genome-wide association analysis was performed using a variety of bioinformatics approaches including SNP analysis, accessory genome analysis and k-mer based analysis. Core and accessory genome analyses found no association with either of the two clinical groups. In this study, the genome sequences of S. aureus bloodstream isolates did not discriminate between bacteraemia and infective endocarditis. Based on our study and the current literature, it is not convincing that a specific S. aureus genotype is clearly associated to infective endocarditis in patients with S. aureus bacteraemia.

Recombination-independent rapid convergent evolution of the gastric pathogen Helicobacter pylori

BACKGROUND

Helicobacter pylori is a human stomach pathogen, naturally-competent for DNA uptake, and prone to homologous recombination. Extensive homoplasy (i.e., phylogenetically-unlinked identical variations) observed in H. pylori genes is considered a hallmark of such recombination. However, H. pylori also exhibits a high mutation rate. The relative adaptive role of homologous recombination and mutation in species diversity is a highly-debated issue in biology. Recombination results in homoplasy. While convergent mutation can also account for homoplasy, its contribution is thought to be minor. We demonstrate here that, contrary to dogma, convergent mutation is a key contributor to Helicobacter pylori homoplasy, potentially driven by adaptive evolution of proteins.

RESULTS

Our present genome-wide analysis shows that homoplastic nonsynonymous (amino acid replacement) changes are not typically accompanied by homoplastic synonymous (silent) variations. Moreover, the majority of the codon positions with homoplastic nonsynonymous changes also contain different (i.e. non-homoplastic) nonsynonymous changes arising from mutation only. This indicates that, to a considerable extent, nonsynonymous homoplasy is due to convergent mutations. High mutation rate or limited availability of evolvable sites cannot explain this excessive convergence, as suggested by our simulation studies. Rather, the genes with convergent mutations are overrepresented in distinct functional categories, suggesting possible selective responses to conditions such as distinct micro-niches in single hosts, and to differences in host genotype, physiology, habitat and diet.

CONCLUSIONS

We propose that mutational convergence is a key player in H. pylori's adaptation and extraordinary persistence in human hosts. High frequency of mutational convergence could be due to saturation of evolvable sites capable of responding to selection pressures, while the number of mutable residues is far from saturation. We anticipate a similar scenario of mutational vs. recombinational genome dynamics or plasticity for other naturally competent microbes where strong positive selection could favor frequent convergent mutations in adaptive protein evolution.

Genes significantly associated with lineage II food isolates of Listeria monocytogenes

Background

Listeria monocytogenes is a widespread foodborne pathogen that can cause listeriosis, a potentially fatal infection. L. monocytogenes is subdivided into four phylogenetic lineages, with the highest incidence of listeriosis occurring within lineage I followed by lineage II. Strains of L. monocytogenes differ in their phenotypic characteristics, including virulence. However, the genetic bases for these observed differences are not well understood, and current efforts to monitor L. monocytogenes in food consider all strains to be equally virulent. We use a comparative genomics approach to identify genes and single nucleotide polymorphisms (SNPs) in 174 clinical and food isolates of L. monocytogenes that potentially contribute to virulence or the capacity to adapt to food environments.

Results

No SNPs are significantly associated with food or clinical isolates. No genes are significantly associated with food or clinical isolates from lineage I, but eight genes consisting of multiple homologues are associated with lineage II food isolates. These include three genes which encode hypothetical proteins, the cadmium resistance genes cadA and cadC, the multi-drug resistance gene ebrB, a quaternary ammonium compound resistance gene qac, and a regulatory gene. All eight genes are plasmid-borne, and most closed L. monocytogenes plasmids carry at least five of the genes (24/27). In addition, plasmids are more frequently associated with lineage II food isolates than with lineage II clinical isolates.

Conclusions

We identify eight genes that are significantly associated with food isolates in lineage II. Interestingly, the eight genes are virtually absent in lineage II outbreak isolates, are composed of homologues which show a nonrandom distribution among lineage I serotypes, and the sequences are highly conserved across 27 closed Listeria plasmids. The functions of these genes should be explored further and will contribute to our understanding of how L. monocytogenes adapts to the host and food environments. Moreover, these genes may also be useful as markers for risk assessment models of either pathogenicity or the ability to proliferate in food and the food processing environment.

Electronic supplementary material

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

Environmental pleiotropy and demographic history direct adaptation under antibiotic selection

Evolutionary rescue following environmental change requires mutations permitting population growth in the new environment. If change is severe enough to prevent most of the population reproducing, rescue becomes reliant on mutations already present. If change is sustained, the fitness effects in both environments, and how they are associated—termed ‘environmental pleiotropy’—may determine which alleles are ultimately favoured. A population’s demographic history—its size over time—influences the variation present. Although demographic history is known to affect the probability of evolutionary rescue, how it interacts with environmental pleiotropy during severe and sustained environmental change remains unexplored. Here, we demonstrate how these factors interact during antibiotic resistance evolution, a key example of evolutionary rescue fuelled by pre-existing mutations with pleiotropic fitness effects. We combine published data with novel simulations to characterise environmental pleiotropy and its effects on resistance evolution under different demographic histories. Comparisons among resistance alleles typically revealed no correlation for fitness—i.e., neutral pleiotropy—above and below the sensitive strain’s minimum inhibitory concentration. Resistance allele frequency following experimental evolution showed opposing correlations with their fitness effects in the presence and absence of antibiotic. Simulations demonstrated that effects of environmental pleiotropy on allele frequencies depended on demographic history. At the population level, the major influence of environmental pleiotropy was on mean fitness, rather than the probability of evolutionary rescue or diversity. Our work suggests that determining both environmental pleiotropy and demographic history is critical for predicting resistance evolution, and we discuss the practicalities of this during in vivo evolution.

Hotspot mutations and ColE1 plasmids contribute to the fitness of Salmonella Heidelberg in poultry litter

Salmonella enterica subsp. enterica serovar Heidelberg (S. Heidelberg) is a clinically-important serovar linked to food-borne illness, and commonly isolated from poultry. Investigations of a large, multistate outbreak in the USA in 2013 identified poultry litter (PL) as an important extra-intestinal environment that may have selected for specific S. Heidelberg strains. Poultry litter is a mixture of bedding materials and chicken excreta that contains chicken gastrointestinal (GI) bacteria, undigested feed, feathers, and other materials of chicken origin. In this study, we performed a series of controlled laboratory experiments which assessed the microevolution of two S. Heidelberg strains (SH-2813 and SH-116) in PL previously used to raise 3 flocks of broiler chickens. The strains are closely related at the chromosome level, differing from the reference genome by 109 and 89 single nucleotide polymorphisms/InDels, respectively. Whole genome sequencing was performed on 86 isolates recovered after 0, 1, 7 and 14 days of microevolution in PL. Only strains carrying an IncX1 (37kb), 2 ColE1 (4 and 6kb) and 1 ColpVC (2kb) plasmids survived more than 7 days in PL. Competition experiments showed that carriage of these plasmids was associated with increased fitness. This increased fitness was associated with an increased copy number of IncX1 and ColE1 plasmids. Further, all Col plasmid-bearing strains had hotspot mutations in 37 loci on the chromosome and in 3 loci on the IncX1 plasmid. Additionally, we observed a decrease in susceptibility to tobramycin, kanamycin, gentamicin, neomycin and fosfomycin for Col plasmid-bearing strains. Our study demonstrates how positive selection from poultry litter can change the evolutionary path of S. Heidelberg.

Mechanical interactions in bacterial colonies and the surfing probability of beneficial mutations

Bacterial conglomerates such as biofilms and microcolonies are ubiquitous in nature and play an important role in industry and medicine. In contrast to well-mixed cultures routinely used in microbial research, bacteria in a microcolony interact mechanically with one another and with the substrate to which they are attached. Here, we use a computer model of a microbial colony of rod-shaped cells to investigate how physical interactions between cells determine their motion in the colony and how this affects biological evolution. We show that the probability that a faster-growing mutant ‘surfs’ at the colony's frontier and creates a macroscopic sector depends on physical properties of cells (shape, elasticity and friction). Although all these factors contribute to the surfing probability in seemingly different ways, their effects can be summarized by two summary statistics that characterize the front roughness and cell alignment. Our predictions are confirmed by experiments in which we measure the surfing probability for colonies of different front roughness. Our results show that physical interactions between bacterial cells play an important role in biological evolution of new traits, and suggest that these interactions may be relevant to processes such as de novo evolution of antibiotic resistance.

Recurrent Reverse Evolution Maintains Polymorphism after Strong Bottlenecks in Commensal Gut Bacteria

The evolution of new strains within the gut ecosystem is poorly understood. We used a natural but controlled system to follow the emergence of intraspecies diversity of commensal Escherichia coli, during three rounds of adaptation to the mouse gut (∼1,300 generations). We previously showed that, in the first round, a strongly beneficial phenotype (loss-of-function for galactitol consumption; gat-negative) spread to >90% frequency in all colonized mice. Here, we show that this loss-of-function is repeatedly reversed when a gat-negative clone colonizes new mice. The regain of function occurs via compensatory mutation and reversion, the latter leaving no trace of past adaptation. We further show that loss-of-function adaptive mutants reevolve, after colonization with an evolved gat-positive clone. Thus, even under strong bottlenecks a regime of strong-mutation-strong-selection dominates adaptation. Coupling experiments and modeling, we establish that reverse evolution recurrently generates two coexisting phenotypes within the microbiota that can or not consume galactitol (gat-positive and gat-negative, respectively). Although the abundance of the dominant strain, the gat-negative, depends on the microbiota composition, gat-positive abundance is independent of the microbiota composition and can be precisely manipulated by supplementing the diet with galactitol. These results show that a specific diet is able to change the abundance of specific strains. Importantly, we find polymorphism for these phenotypes in indigenous Enterobacteria of mice and man. Our results demonstrate that natural selection can greatly overwhelm genetic drift at structuring the strain diversity of gut commensals and that competition for limiting resources may be a key mechanism for maintaining polymorphism in the gut.

Pathoadaptive Mutations in Uropathogenic Escherichia coli

Uropathogenic Escherichia coli (UPEC) are opportunistic human pathogens that primarily circulate as part of commensal intestinal microbiota. Though they have the ability to survive and proliferate in various urinary tract compartments, the urinary tract is a transient, occasional habitat for UPEC. Because of this, most of the UPEC traits have originally evolved to serve in intestinal colonization and transmission. Some of these bacterial traits serve as virulence factors - they are critical to or assist in survival of UPEC as pathogens, and the structure and/or function may be specialized for the infection. Other traits could serve as anti-virulence factors - they represent liability in the urinary tract and are under selection to be lost or inactivated during the infection. Inactivation, variation, or other changes of the bacterial genes that increase the pathogen's fitness during the infection are called pathoadaptive mutations. This chapter describes examples of pathoadaptive mutations in UPEC and provides rationale for their further in-depth study.

Core Genes Evolve Rapidly in the Long-Term Evolution Experiment with Escherichia coli

Bacteria can evolve rapidly under positive selection owing to their vast numbers, allowing their genes to diversify by adapting to different environments. We asked whether the same genes that evolve rapidly in the long-term evolution experiment (LTEE) with Escherichia coli have also diversified extensively in nature. To make this comparison, we identified ∼2000 core genes shared among 60 E. coli strains. During the LTEE, core genes accumulated significantly more nonsynonymous mutations than flexible (i.e., noncore) genes. Furthermore, core genes under positive selection in the LTEE are more conserved in nature than the average core gene. In some cases, adaptive mutations appear to modify protein functions, rather than merely knocking them out. The LTEE conditions are novel for E. coli, at least in relation to its evolutionary history in nature. The constancy and simplicity of the environment likely favor the complete loss of some unused functions and the fine-tuning of others.

Population Phylogenomics of Extraintestinal Pathogenic Escherichia coli

The emergence of genomics over the last 10 years has provided new insights into the evolution and virulence of extraintestinal Escherichia coli. By combining population genetics and phylogenetic approaches to analyze whole-genome sequences, it became possible to link genomic features to specific phenotypes, such as the ability to cause urinary tract infections. An E. coli chromosome can vary extensively in length, ranging from 4.3 to 6.2 Mb, encoding 4,084 to 6,453 proteins. This huge diversity is structured as a set of less than 2,000 genes (core genome) that are conserved between all the strains and a set of variable genes. Based on the core genome, the history of the species can be reliably reconstructed, revealing the recent emergence of phylogenetic groups A and B1 and the more ancient groups B2, F, and D. Urovirulence is most often observed in B2/F/D group strains and is a multigenic process involving numerous combinations of genes and specific alleles with epistatic interactions, all leading down multiple evolutionary paths. The genes involved mainly code for adhesins, toxins, iron capture systems, and protectins, as well as metabolic pathways and mutation-rate-control systems. However, the barrier between commensal and uropathogenic E. coli strains is difficult to draw as the factors that are responsible for virulence have probably also been selected to allow survival of E. coli as a commensal in the intestinal tract. Genomic studies have also demonstrated that infections are not the result of a unique and stable isolate, but rather often involve several isolates with variable levels of diversity that dynamically changes over time.

Adaptation of Escherichia coli traversing from the faecal environment to the urinary tract

The majority of extraintestinal pathogenic Escherichia coli (ExPEC) causing urinary tract infections (UTI) are found in the patient's own gut flora, but only limited knowledge is available on the potential adaptation that may occur in the bacteria in order to traverse the perineum and successfully infect the urinary tract. Here, matching pairs of faecal and UTI isolates from 42 patients were compared pairwise using in-depth whole-genome sequencing to investigate whether genetic changes were evident for successful colonization in these two different environments. The identified non-synonymous mutations (0-12 substitutions in each pair) were primarily associated to genes encoding virulence factors and nutrient metabolism; and indications of parallel evolution were observed in genes encoding the major phase-variable protein antigen 43, a toxin/antitoxin locus and haemolysin B. No differences in virulence potential were observed in a mouse UTI model for five matching faecal and UTI isolates with or without mutations in antigen 43 and haemolysin B. Variations in plasmid content were observed in only four of the 42 pairs. Although, we observed mutations in known UTI virulence genes for a few pairs, the majority showed no detectable differences with respect to mutations or mobilome when compared to their faecal counterpart. The results show that UPECs are successful in colonizing both the bladder and gut without adaptation.

Mutation-Driven Divergence and Convergence Indicate Adaptive Evolution of the Intracellular Human-Restricted Pathogen, Bartonella bacilliformis

Among all species of Bartonella, human-restricted Bartonella bacilliformis is the most virulent but harbors one of the most reduced genomes. Carrión’s disease, the infection caused by B. bacilliformis, has been afflicting poor rural populations for centuries in the high-altitude valleys of the South American Andes, where the pathogen’s distribution is probably restricted by its sand fly vector’s range. Importantly, Carrión’s disease satisfies the criteria set by the World Health Organization for a disease amenable to elimination. However, to date, there are no genome-level studies to identify potential footprints of B. bacilliformis (patho)adaptation. Our comparative genomic approach demonstrates that the evolution of this intracellular pathogen is shaped predominantly via mutation. Analysis of strains having publicly-available genomes shows high mutational divergence of core genes leading to multiple sub-species. We infer that the sub-speciation event might have happened recently where a possible adaptive divergence was accelerated by intermediate emergence of a mutator phenotype. Also, within a sub-species the pathogen shows inter-clonal adaptive evolution evidenced by non-neutral accumulation of convergent amino acid mutations. A total of 67 non-recombinant core genes (over-representing functional categories like DNA repair, glucose metabolic process, ATP-binding and ligase) were identified as candidates evolving via adaptive mutational convergence. Such convergence, both at the level of genes and their encoded functions, indicates evolution of B. bacilliformis clones along common adaptive routes, while there was little diversity within a single clone.

How host-restriction, intracellularity and genome reduction interplay to exert or maintain virulence is poorly characterized. The fact that B. bacilliformis is the most pathogenic Bartonella and has a highly reduced genome makes it an attractive model to gain insights into (patho)adaptive evolution of intracellular pathogens. Also, B. bacilliformis is known to lack many virulence genes present in other Bartonella, indicating unique strategies of (patho)adaptation. Our study reveals a prevalent nature of mutational force in B. bacilliformis evolution with two distinct outcomes: (a) mutational divergence leading to sub-speciation, possibly recently, via accelerated accumulation and fixation of favorable mutations mediated by a mutator phenotype; and (b) mutational convergence between clones of a sub-species exhibiting shared functional trajectories of adaptive evolution. Our findings highlight positions accumulating adaptive mutations in candidate genes, offering future functional studies to elucidate B. bacilliformis virulence evolution, and of broad application to intracellular pathogens with a reduced gene repertoire.

The Non-Fimbriate Phenotype Is Predominant among Salmonella enterica Serovar Choleraesuis from Swine and Those Non-Fimbriate Strains Possess Distinct Amino Acid Variations in FimH

Although most Salmonella serovars are able to infect a range of animal hosts, some have acquired the ability to cause systemic infections of specific hosts. For example, Salmonella enterica serovar Choleraesuis is primarily associated with systemic infection in swine. Adherence to host epithelial cells is considered a prerequisite for initial infection, and fimbrial appendages on the outer membrane of the bacteria are implicated in this process. Although type 1 fimbriae encoded by the fim gene cluster are commonly found in Salmonella serovars, it is not known whether S. Choleraesuis produces this fimbrial type and if and how fimbriae are involved in pathogenesis. In the present study, we demonstrated that only four out of 120 S. Choleraesuis isolates from pigs with salmonellosis produced type 1 fimbriae as assayed by the yeast agglutination test and electron microscopy. One of the 116 non-type 1 fimbria-producing isolates was transformed with plasmids carrying different fim genes from S. Typhimurium LB5010, a type 1 fimbria-producing strain. Our results indicate that non-type 1 fimbria-producing S. Choleraesuis required only an intact fimH to regain the ability to produce fimbrial appendages. Sequence comparison revealed six amino acid variations between the FimH of the non-type 1 fimbria-producing S. Choleraesuis isolates and those of the type 1 fimbria-producing S. Choleraesuis isolates. S. Choleraesuis that produced type 1 fimbriae contained FimH with an amino acid sequence identical to that of S. Typhimurium LB5010. Site-directed mutagenesis leading to the replacement of the non-conserved residues revealed that a change from glycine to valine at position of 63 (G63V) resulted in a non-type 1 fimbria-producing S. Choleraesuis being able to express type 1 fimbriae on its outer membrane. It is possible that this particular amino acid change prevents this polypeptide from proper interaction with other Fim subunits required for assembly of an intact type 1 fimbrial shaft in S. Choleraesuis; however, it remains to be determined if and how the absence of type 1 fimbriae production is related to the systemic infection of the swine host by S. Choleraesuis.

Practical Approaches for Detecting Selection in Microbial Genomes

Microbial genome evolution is shaped by a variety of selective pressures. Understanding how these processes occur can help to address important problems in microbiology by explaining observed differences in phenotypes, including virulence and resistance to antibiotics. Greater access to whole-genome sequencing provides microbiologists with the opportunity to perform large-scale analyses of selection in novel settings, such as within individual hosts. This tutorial aims to guide researchers through the fundamentals underpinning popular methods for measuring selection in pathogens. These methods are transferable to a wide variety of organisms, and the exercises provided are designed for researchers with any level of programming experience.

Trade-Offs of Escherichia coli Adaptation to an Intracellular Lifestyle in Macrophages

The bacterium Escherichia coli exhibits remarkable genomic and phenotypic variation, with some pathogenic strains having evolved to survive and even replicate in the harsh intra-macrophage environment. The rate and effects of mutations that can cause pathoadaptation are key determinants of the pace at which E. coli can colonize such niches and become pathogenic. We used experimental evolution to determine the speed and evolutionary paths undertaken by a commensal strain of E. coli when adapting to intracellular life. We estimated the acquisition of pathoadaptive mutations at a rate of 10⁻⁶ per genome per generation, resulting in the fixation of more virulent strains in less than a hundred generations. Whole genome sequencing of independently evolved clones showed that the main targets of intracellular adaptation involved loss of function mutations in genes implicated in the assembly of the lipopolysaccharide core, iron metabolism and di- and tri-peptide transport, namely rfaI, fhuA and tppB, respectively. We found a substantial amount of antagonistic pleiotropy in evolved populations, as well as metabolic trade-offs, commonly found in intracellular bacteria with reduced genome sizes. Overall, the low levels of clonal interference detected indicate that the first steps of the transition of a commensal E. coli into intracellular pathogens are dominated by a few pathoadaptive mutations with very strong effects.

PanCoreGen - Profiling, detecting, annotating protein-coding genes in microbial genomes

A large amount of genomic data, especially from multiple isolates of a single species, has opened new vistas for microbial genomics analysis. Analyzing the pan-genome (i.e. the sum of genetic repertoire) of microbial species is crucial in understanding the dynamics of molecular evolution, where virulence evolution is of major interest. Here we present PanCoreGen - a standalone application for pan- and core-genomic profiling of microbial protein-coding genes. PanCoreGen overcomes key limitations of the existing pan-genomic analysis tools, and develops an integrated annotation-structure for a species-specific pan-genomic profile. It provides important new features for annotating draft genomes/contigs and detecting unidentified genes in annotated genomes. It also generates user-defined group-specific datasets within the pan-genome. Interestingly, analyzing an example-set of Salmonella genomes, we detect potential footprints of adaptive convergence of horizontally transferred genes in two human-restricted pathogenic serovars - Typhi and Paratyphi A. Overall, PanCoreGen represents a state-of-the-art tool for microbial phylogenomics and pathogenomics study.

Faecal Escherichia coli isolates show potential to cause endogenous infection in patients admitted to the ICU in a tertiary care hospital

Nosocomial infections are acquired during hospital treatment or in a hospital environment. One such infecting agent, Escherichia coli, harbours many virulence genes that enable it to become pathogenic, causing damage to the host. The mechanism of the E. coli virulence factors provenance to cause infection in host environments is not clearly elucidated. We investigated the virulence and pathogenicity of E. coli affected by the host environment. For this, blood (n = 78) and faecal (n = 83) E. coli isolates were collected from patients with and without sepsis, respectively, who had been admitted to the intensive care unit. The E. coli genomic DNA was isolated; the phylogenetic grouping was conducted by triplex PCR. The occurrence of nine virulence genes among the all the isolates was confirmed by gene-specific PCR. The prevalence of E. coli in blood isolates was more in phylogenetic groups B2 and D compared to groups A and B1. However, in faecal isolates, there was no significant difference. The prevalence of adhesin and toxin (papG, sfa, afa, cnf1, hlyA) genes was higher in blood compared to faecal E. coli isolates. However, the prevalence of aer, traT and PAI was similar as well as higher among both of these groups. These observations indicate a role of external environment (hospital setting) on host susceptibility (development of infection) in the faecal E. coli isolates, thereby making the patient prone to a sepsis condition.

Measurably evolving pathogens in the genomic era

Current sequencing technologies have created unprecedented opportunities for studying microbial populations. For pathogens with comparatively low per-site mutation rates, such as DNA viruses and bacteria, whole-genome sequencing can reveal the accumulation of novel genetic variation between population samples taken at different times. The concept of 'measurably evolving populations' and related analytical approaches have provided powerful insights for fast-evolving RNA viruses, but their application to other pathogens is still in its infancy. We argue that previous distinctions between slow- and fast-evolving pathogens become blurred once evolution is assessed at a genome-wide scale, and we highlight important analytical challenges to be overcome to infer pathogen population dynamics from genomic data.

Bacterial Strain Diversity Within Wounds

Significance: Rare bacterial taxa (taxa of low relative frequency) are numerous and ubiquitous in virtually any sample-including wound samples. In addition, even the high-frequency genera and species contain multiple strains. These strains, individually, are each only a small fraction of the total bacterial population. Against the view that wounds contain relatively few kinds of bacteria, this newly recognized diversity implies a relatively high rate of migration into the wound and the potential for diversification during infection. Understanding the biological and medical importance of these numerous taxa is an important new element of wound microbiology. Recent Advances: Only recently have these numerous strains been discovered; the technology to detect, identify, and characterize them is still in its infancy. Multiple strains of both gram-negative and gram-positive bacteria have been found in a single wound. In the few cases studied, the distribution of the bacteria suggests microhabitats and biological interactions. Critical Issues: The distribution of the strains, their phenotypic diversity, and their interactions are still largely uncharacterized. The technologies to investigate this level of genomic detail are still developing and have not been largely deployed to investigate wounds. Future Directions: As advanced metagenomics, single-cell genomics, and advanced microscopy develop, the study of wound microbiology will better address the complex interplay of numerous individually rare strains with both the host and each other.

Frequency and fitness consequences of bacteriophage φ6 host range mutations

Viruses readily mutate and gain the ability to infect novel hosts, but few data are available regarding the number of possible host range-expanding mutations allowing infection of any given novel host, and the fitness consequences of these mutations on original and novel hosts. To gain insight into the process of host range expansion, we isolated and sequenced 69 independent mutants of the dsRNA bacteriophage Φ6 able to infect the novel host, Pseudomonas pseudoalcaligenes. In total, we found at least 17 unique suites of mutations among these 69 mutants. We assayed fitness for 13 of 17 mutant genotypes on P. pseudoalcaligenes and the standard laboratory host, P. phaseolicola. Mutants exhibited significantly lower fitnesses on P. pseudoalcaligenes compared to P. phaseolicola. Furthermore, 12 of the 13 assayed mutants showed reduced fitness on P. phaseolicola compared to wildtype Φ6, confirming the prevalence of antagonistic pleiotropy during host range expansion. Further experiments revealed that the mechanistic basis of these fitness differences was likely variation in host attachment ability. In addition, using computational protein modeling, we show that host-range expanding mutations occurred in hotspots on the surface of the phage's host attachment protein opposite a putative hydrophobic anchoring domain.

Emergence of a new epidemic Neisseria meningitidis serogroup A Clone in the African meningitis belt: high-resolution picture of genomic changes that mediate immune evasion

In the African “meningitis belt,” outbreaks of meningococcal meningitis occur in cycles, representing a model for the role of host-pathogen interactions in epidemic processes. The periodicity of the epidemics is not well understood, nor is it currently possible to predict them. In our longitudinal colonization and disease surveys, we have observed waves of clonal replacement with the same serogroup, suggesting that immunity to noncapsular antigens plays a significant role in natural herd immunity. Here, through comparative genomic analysis of 100 meningococcal isolates, we provide a high-resolution view of the evolutionary changes that occurred during clonal replacement of a hypervirulent meningococcal clone (ST-7) by a descendant clone (ST-2859). We show that the majority of genetic changes are due to homologous recombination of laterally acquired DNA, with more than 20% of these events involving acquisition of DNA from other species. Signals of adaptation to evade herd immunity were indicated by genomic hot spots of recombination. Most striking is the high frequency of changes involving the pgl locus, which determines the glycosylation patterns of major protein antigens. High-frequency changes were also observed for genes involved in the regulation of pilus expression and the synthesis of Maf3 adhesins, highlighting the importance of these surface features in host-pathogen interaction and immune evasion.

While established meningococcal capsule polysaccharide vaccines are protective through the induction of anticapsular antibodies, findings of our longitudinal studies in the African meningitis belt have indicated that immunity to noncapsular antigens plays a significant role in natural herd immunity. Our results show that meningococci evade herd immunity through the rapid homologous replacement of just a few key genomic loci that affect noncapsular cell surface components. Identification of recombination hot spots thus represents an eminent approach to gain insight into targets of protective natural immune responses. Moreover, our results highlight the role of the dynamics of the protein glycosylation repertoire in immune evasion by Neisseria meningitidis. These results have major implications for the design of next-generation protein-based subunit vaccines.