Search

We analyzed the within-household evolution of two household-associated Escherichia coli strains from pandemic clonal group ST131- H 30, using isolates recovered from five individuals within two families, each of which had a distinct strain. Family 1's strain was represented by a urine isolate from the index patient (older sister) with recurrent cystitis and a blood isolate from her younger sister with fatal urosepsis. Family 2's strain was represented by a urine isolate from the index patient (father) with pyelonephritis and renal abscesses, blood and kidney drainage isolates from the daughter with emphysematous pyelonephritis, and urine and fecal isolates from the mother with cystitis. Collectively, the several variants of each family's strain had accumulated a total of 8 (family 1) and 39 (family 2) point mutations; no two isolates were identical. Of the 47 total mutations, 36 resulted in amino acid changes or truncation of coded proteins. Fourteen such mutations (39%) targeted genes encoding transcriptional regulators, and 9 (25%) involved DNA-binding transcription factors (TFs), which significantly exceeded the relative contribution of TF genes to the isolates' genomes (∼6%). At least one-half of the transcriptional regulator mutations were inactivating, based on phenotypic and/or transcriptional analysis. In particular, inactivating mutations in the global regulator LrhA (repressor of type 1 fimbriae and flagella) occurred in the blood isolates from both households and increased the virulence of E. coli strains in a murine sepsis model. The results indicate that E. coli undergoes adaptive evolution between and/or within hosts, generating subpopulations with distinctive phenotypes and virulence potential. IMPORTANCE The clonal evolution of bacterial strains associated with interhost transmission is poorly understood. We characterized the genome sequences of clonal descendants of two Escherichia coli strains, recovered at different time points from multiple individuals within two households who had different types of urinary tract infection. We found evidence that the E. coli strains underwent extensive mutational diversification between and within these individuals, driven disproportionately by inactivation of transcriptional regulators. In urosepsis isolates, the mutations observed in the global regulator LrhA increased bacterial virulence in a murine sepsis model. Our findings help in understanding the adaptive dynamics and strategies of E. coli during short-term natural evolution., (Copyright © 2017 American Society for Microbiology.)

Horizontal acquisition of novel chromosomal genes is considered to be a key process in the evolution of bacterial pathogens. However, the identification of gene presence or absence could be hindered by the inconsistencies in bacterial genome annotations. Here, we performed a cross-annotation of omnipresent core and mosaic accessory genes in the chromosome of Salmonella enterica serovar Typhimurium across a total of 20 fully assembled genomes deposited into GenBank. Cross-annotation resulted in a 32% increase in the number of core genes and a 3-fold decrease in the number of genes identified as mosaic genes (i.e., genes present in some strains only) by the original annotation. Of the remaining noncore genes, the vast majority were prophage genes, and 255 of the nonphage genes were actually of core origin but lost in some strains upon the emergence of the S Typhimurium serovar, suggesting that the chromosomal portion of the S Typhimurium genome acquired a very limited number of novel genes other than prophages. Only horizontally acquired nonphage genes related to bacterial fitness or virulence were found in four recently sequenced isolates, all located on three different genomic islands that harbor multidrug resistance determinants. Thus, the extensive use of antimicrobials could be the main selection force behind the new fitness gene acquisition and the emergence of novel Salmonella pathotypes., Importance: Significant discrepancies in the annotations of bacterial genomes could mislead the conclusions about evolutionary origin of chromosomal genes, as we demonstrate here via a cross-annotation-based analysis of Salmonella Typhimurium genomes from GenBank. We conclude that despite being able to infect a broad range of vertebrate hosts, the genomic diversity of S Typhimurium strains is almost exclusively limited to gene loss and the transfer of prophage DNA. Only nonphage chromosomal genes acquired after the emergence of the serovar are linked to the genomic islands harboring multidrug resistance factors. Since the fitness factors could lead to increased virulence, this poses an important research question: could overuse or misuse of antimicrobials act as selection forces for the emergence of more pathogenic strains of Salmonella?, (Copyright © 2016, American Society for Microbiology. All Rights Reserved.)

When estimating a phylogeny from a multiple sequence alignment, researchers often assume the absence of recombination. However, if recombination is present, then tree estimation and all downstream analyses will be impacted, because different segments of the sequence alignment support different phylogenies. Similarly, convergent selective pressures at the molecular level can also lead to phylogenetic tree incongruence across the sequence alignment. Current methods for detection of phylogenetic incongruence are not equipped to distinguish between these two different mechanisms and assume that the incongruence is a result of recombination or other horizontal transfer of genetic information. We propose a new recombination detection method that can make this distinction, based on synonymous codon substitution distances. Although some power is lost by discarding the information contained in the nonsynonymous substitutions, our new method has lower false positive probabilities than the comparable recombination detection method when the phylogenetic incongruence signal is due to convergent evolution. We apply our method to three empirical examples, where we analyze: (1) sequences from a transmission network of the human immunodeficiency virus, (2) tlpB gene sequences from a geographically diverse set of 38 Helicobacter pylori strains, and (3) hepatitis C virus sequences sampled longitudinally from one patient.

Analysis of genetic differences (gene presence/absence and nucleotide polymorphisms) among strains of a bacterial species is crucial to understanding molecular mechanisms of bacterial pathogenesis and selecting targets for novel antibacterial therapeutics. However, lack of genome-wide association studies on large and epidemiologically well-defined strain collections from the same species makes it difficult to identify the genes under positive selection and define adaptive polymorphisms in those genes. To address this need and to overcome existing limitations, we propose to create a "microbial variome"--a species-specific resource database of genomic variations based on molecular evolutionary analysis. Here, we present prototype variome databases of Escherichia coli and Salmonella enterica subspecies enterica (http://depts.washington.edu/sokurel/variome, last accessed March 26, 2013). The prototypes currently include the point mutations data of core protein-coding genes from completely sequenced genomes of 22 E. coli and 17 S. enterica strains. These publicly available databases allow for single- and multiple-field sorting, filtering, and searching of the gene variability data and the potential adaptive significance. Such resource databases would immensely help experimental research, clinical diagnostics, epidemiology, and environmental control of human pathogens.

Unlabelled: Salmonella and Escherichia coli mannose-binding type 1 fimbriae exhibit highly similar receptor specificities, morphologies, and mechanisms of assembly but are nonorthologous in nature, i.e., not closely related evolutionarily. Their operons differ in chromosomal location, gene arrangement, and regulatory components. In the current study, we performed a comparative genetic and structural analysis of the major structural subunit, FimA, from Salmonella and E. coli and found that FimA pilins undergo diverse evolutionary adaptation in the different species. Whereas the E. coli fimA locus is characterized by high allelic diversity, frequent intragenic recombination, and horizontal movement, Salmonella fimA shows structural diversity that is more than 5-fold lower without strong evidence of gene shuffling or homologous recombination. In contrast to Salmonella FimA, the amino acid substitutions in the E. coli pilin heavily target the protein regions that are predicted to be exposed on the external surface of fimbriae. Altogether, our results suggest that E. coli, but not Salmonella, type 1 fimbriae display a high level of structural diversity consistent with a strong selection for antigenic variation under immune pressure. Thus, type 1 fimbriae in these closely related bacterial species appear to function in distinctly different physiological environments., Importance: E. coli and Salmonella are enteric bacteria that are closely related from an evolutionary perspective. They are both notorious human pathogens, though with somewhat distinct ecologies and virulence mechanisms. Type 1 fimbriae are rod-shaped surface appendages found in most E. coli and Salmonella isolates. In both species, they mediate bacterial adhesion to mannose receptors on host cells and share essentially the same morphology and assembly mechanisms. Here we show that despite the strong resemblances in function and structure, they are exposed to very different natural selection environments. Sequence analysis indicates that E. coli, but not Salmonella, fimbriae are subjected to strong immune pressure, resulting in a high level of major fimbrial protein gene shuffling and interbacterial transfer. Thus, evolutionary analysis tools can provide evidence of divergent physiological roles of functionally similar traits in different bacterial species.

One of the strongest signals of adaptive molecular evolution of proteins is the occurrence of convergent hot spot mutations: repeated changes in the same amino acid positions. We performed a comparative genome-wide analysis of mutation-driven evolution of core (omnipresent) genes in 17 strains of Salmonella enterica subspecies I and 22 strains of Escherichia coli. More than 20% of core genes in both Salmonella and E. coli accumulated hot spot mutations, with a predominance of identical changes having recent evolutionary origin. There is a significant overlap in the functional categories of the adaptively evolving genes in both species, although mostly via separate molecular mechanisms. As a strong evidence of the link between adaptive mutations and virulence in Salmonella, two human-restricted serovars, Typhi and Paratyphi A, shared the highest number of genes with serovar-specific hot spot mutations. Many of the core genes affected by Typhi/Paratyphi A-specific mutations have known virulence functions. For each species, a list of nonrecombinant core genes (and the hot spot mutations therein) under positive selection is provided.

Class 5 fimbriae of enterotoxigenic Escherichia coli (ETEC) comprise eight serologically discrete colonization factors that mediate small intestinal adhesion. Their differentiation has been attributed to the pressure imposed by host adaptive immunity. We sequenced the major pilin and minor adhesin subunit genes of a geographically diverse population of ETEC elaborating CFA/I (n = 31), CS17 (n = 20), and CS2 (n = 18) and elucidated the functional effect of microevolutionary processes. Between the fimbrial types, the pairwise nucleotide diversity for the pilin or adhesin genes ranged from 35-43%. Within each fimbrial type, there were 17 non-synonymous and 1 synonymous point mutations among all pilin or adhesin gene copies, implying that each fimbrial type was acquired by ETEC strains very recently, consistent with a recent origin of this E. coli pathotype. The 17 non-synonymous allelic differences occurred in the CFA/I pilin gene cfaB (two changes) and adhesin gene cfaE (three changes), and CS17 adhesin gene csbD (12 changes). All but one amino acid change in the adhesins clustered around the predicted ligand-binding pocket. Functionally, these changes conferred an increase in cell adhesion in a flow chamber assay. In contrast, the two mutations in the non-adhesive CfaB subunit localized to the intersubunit interface and significantly reduced fimbrial adhesion in this assay. In conclusion, naturally occurring mutations in the ETEC adhesive and non-adhesive subunits altered function, were acquired under positive selection, and are predicted to impact bacteria-host interactions.

In the rabbit model of syphilis, infection phenotypes associated with the Nichols and Chicago strains of Treponema pallidum (T. pallidum), though similar, are not identical. Between these strains, significant differences are found in expression of, and antibody responses to some candidate virulence factors, suggesting the existence of functional genetic differences between isolates. The Chicago strain genome was therefore sequenced and compared to the Nichols genome, available since 1998. Initial comparative analysis suggested the presence of 44 single nucleotide polymorphisms (SNPs), 103 small (≤3 nucleotides) indels, and 1 large (1204 bp) insertion in the Chicago genome with respect to the Nichols genome. To confirm the above findings, Sanger sequencing was performed on most loci carrying differences using DNA from Chicago and the Nichols strain used in the original T. pallidum genome project. A majority of the previously identified differences were found to be due to errors in the published Nichols genome, while the accuracy of the Chicago genome was confirmed. However, 20 SNPs were confirmed between the two genomes, and 16 (80.0%) were found in coding regions, with all being of non-synonymous nature, strongly indicating action of positive selection. Sequencing of 16 genomic loci harboring SNPs in 12 additional T. pallidum strains, (SS14, Bal 3, Bal 7, Bal 9, Sea 81-3, Sea 81-8, Sea 86-1, Sea 87-1, Mexico A, UW231B, UW236B, and UW249C), was used to identify "Chicago-" or "Nichols -specific" differences. All but one of the 16 SNPs were "Nichols-specific", with Chicago having identical sequences at these positions to almost all of the additional strains examined. These mutations could reflect differential adaptation of the Nichols strain to the rabbit host or pathoadaptive mutations acquired during human infection. Our findings indicate that SNPs among T. pallidum strains emerge under positive selection and, therefore, are likely to be functional in nature.

Type IV pili contribute to virulence in Vibrio vulnificus, the bacterium responsible for the majority of fatal seafood-related infections. Here, we performed within- and between-species evolutionary analysis of the gene that encodes the major structural subunit of the pilus, pilA, by comparing it with pilD and gyrB, the genes encoding the type IV prepilin peptidase and beta subunit of DNA gyrase, respectively. Although the diversity in pilD and gyrB is similar to each other and likely to have accumulated after speciation of V. vulnificus, pilA is several times more diverse at both nonsynonymous and synonymous levels. Also, in contrast to pilD and gyrB, there are virtually unrestricted and highly localized horizontal movements of pilA alleles between the major phylogenetic groups of V. vulnificus. The frequent movement of pilA involves homologous recombination of the entire gene with no evidence for intragenic recombination between the alleles. We propose that pilA allelic diversity and horizontal movement is maintained in the population by both diversifying and frequency-dependent selection most likely to escape shellfish innate immunity defense or lytic phages. Other possibilities leading to such selection dynamics of V. vulnificus pilA could involve adaptation to diverse host populations or within-host compartments, or natural DNA uptake and transformation. We show that the history of nucleotide diversification in pilA predates V. vulnificus speciation and this diversification started at or before the time of the last common ancestor for V. vulnificus, Vibrio parahaemolyticus, and Vibrio cholerae. At the same time, it appears that within the various pilA groups of V. vulnificus, there is no positive selection for structural mutations and consequently no evidence for source-sink selection. In contrast, pilD has accumulated a number of apparently adaptive mutations in the regions encoding the membrane-spanning portions of the prepilin peptidase, possibly affecting fimbrial expression and/or function, and is being subjected to source-sink selection dynamics.

Background: Detection of adaptive amino acid changes in proteins under recent short-term selection is of great interest for researchers studying microevolutionary processes in microbial pathogens or any other biological species. However, independent occurrence of such point mutations within genetically diverse haplotypes makes it difficult to detect the selection footprint by using traditional molecular evolutionary analyses. The recently developed Zonal Phylogeny (ZP) has been shown to be a useful analytic tool for identifying the footprints of short-term positive selection. ZP separates protein-encoding genes into evolutionarily long-term (with silent diversity) and short-term (without silent diversity) categories, or zones, followed by statistical analysis to detect signs of positive selection in the short-term zone. However, successful broad application of ZP for analysis of large haplotype datasets requires automation of the relatively labor-intensive computational process., Results: Here we present Zonal Phylogeny Software (ZPS), an application that describes the distribution of single nucleotide polymorphisms (SNPs) of synonymous (silent) and non-synonymous (replacement) nature along branches of the DNA tree for any given protein-coding gene locus. Based on this information, ZPS separates the protein variant haplotypes with silent variability (Primary zone) from those that have recently evolved from the Primary zone variants by amino acid changes (External zone). Further comparative analysis of mutational hot-spot frequencies and haplotype diversity between the two zones allows determination of whether the External zone haplotypes emerged under positive selection., Conclusions: As a visualization tool, ZPS depicts the protein tree in a DNA tree, indicating the most parsimonious numbers of synonymous and non-synonymous changes along the branches of a maximum-likelihood based DNA tree, along with information on homoplasy, reversion and structural mutation hot-spots. Through zonal differentiation, ZPS allows detection of recent adaptive evolution via selection of advantageous structural mutations, even when the advantage conferred by such mutations is relatively short-term (as in the case of "source-sink" evolutionary dynamics, which may represent a major mode of virulence evolution in microbes).

Spread of biological species from primary into novel habitats leads to within-species adaptive niche differentiation and is commonly driven by acquisition of point mutations in individual genes that increase fitness in the alternative environment. However, finding footprints of adaptive niche differentiation in specific genes remains a challenge. Here we describe a novel method to analyze the footprint of pathogenicity-adaptive, or pathoadaptive, mutations in the Escherichia coli gene encoding FimH-the major, mannose-sensitive adhesin. Analysis of distribution of mutations across the nodes and branches of the FimH phylogenetic network shows (1) zonal separation of evolutionary primary structural variants of FimH and recently derived ones, (2) dramatic differences in the ratio of synonymous and nonsynonymous changes between nodes from different zones, (3) evidence for replacement hot-spots in the FimH protein, (4) differential zonal distribution of FimH variants from commensal and uropathogenic E. coli, and (5) pathoadaptive functional changes in FimH brought by the mutations. The selective footprint in fimH indicates that the pathoadaptive niche differentiation of E. coli is either in its initial stages or undergoing an evolutionary "source/sink" dynamic.