rbrothers: R Package for Bayesian Multiple Change-Point Recombination Detection.

Phylogenetic recombination detection is a fundamental task in bioinformatics and evolutionary biology. Most of the computational tools developed to attack this important problem are not integrated into the growing suite of R packages for statistical analysis of molecular sequences. Here, we present an R package, rbrothers, that makes a Bayesian multiple change-point model, one of the most sophisticated model-based phylogenetic recombination tools, available to R users. Moreover, we equip the Bayesian change-point model with a set of pre- and post- processing routines that will broaden the application domain of this recombination detection framework. Specifically, we implement an algorithm that forms the set of input trees required by multiple change-point models. We also provide functionality for checking Markov chain Monte Carlo convergence and creating estimation result summaries and graphics. Using rbrothers, we perform a comparative analysis of two Salmonella enterica genes, fimA and fimH, that encode major and adhesive subunits of the type 1 fimbriae, respectively. We believe that rbrothers, available at R-Forge: http://evolmod.r-forge.r-project.org/, will allow researchers to incorporate recombination detection into phylogenetic workflows already implemented in R

Synonymous and Nonsynonymous Distances Help Untangle Convergent Evolution and Recombination

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.Comment: 21 pages, 8 figures, updated abstrac

ZPS: visualization of recent adaptive evolution of proteins

<p>Abstract</p> <p>Background</p> <p>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.</p> <p>Results</p> <p>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.</p> <p>Conclusions</p> <p>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).</p

Complete genome sequence of the epidemic and highly virulent CTX-M-15-producing H30-Rx subclone of Escherichia coli ST131

We report the complete genome sequence, including five complete plasmid sequences, of Escherichia coli ST131 isolate JJ1886. The isolate was obtained in 2007 in the United States from a patient with fatal urosepsis and belongs to the virulent, CTX-M-15-producing H30-Rx sublineage

Uropathogenic Escherichia coli Metabolite-Dependent Quiescence and Persistence May Explain Antibiotic Tolerance during Urinary Tract Infection

In the present study, it is shown that although Escherichia coli CFT073, a human uropathogenic (UPEC) strain, grows in liquid glucose M9 minimal medium, it fails to grow on glucose M9 minimal medium agar plates seeded with ≤106 CFU. The cells on glucose plates appear to be in a “quiescent” state that can be prevented by various combinations of lysine, methionine, and tyrosine. Moreover, the quiescent state is characteristic of ~80% of E. coli phylogenetic group B2 multilocus sequence type 73 strains, as well as 22.5% of randomly selected UPEC strains isolated from community-acquired urinary tract infections in Denmark. In addition, E. coli CFT073 quiescence is not limited to glucose but occurs on agar plates containing a number of other sugars and acetate as sole carbon sources. It is also shown that a number of E. coliCFT073 mini-Tn5 metabolic mutants (gnd, gdhA, pykF, sdhA, and zwf) are nonquiescent on glucose M9 minimal agar plates and that quiescence requires a complete oxidative tricarboxylic acid (TCA) cycle. In addition, evidence is presented that, although E. coli CFT073 quiescence and persistence in the presence of ampicillin are alike in that both require a complete oxidative TCA cycle and each can be prevented by amino acids, E. coli CFT073 quiescence occurs in the presence or absence of a functional rpoS gene, whereas maximal persistence requires a nonfunctional rpoS. Our results suggest that interventions targeting specific central metabolic pathways may mitigate UPEC infections by interfering with quiescence and persistence

Uropathogenic <i>Escherichia coli</i> metabolite-dependent quiescence and persistence may explain antibiotic tolerance during urinary tract infection

ABSTRACT In the present study, it is shown that although Escherichia coli CFT073, a human uropathogenic (UPEC) strain, grows in liquid glucose M9 minimal medium, it fails to grow on glucose M9 minimal medium agar plates seeded with ≤106 CFU. The cells on glucose plates appear to be in a “quiescent” state that can be prevented by various combinations of lysine, methionine, and tyrosine. Moreover, the quiescent state is characteristic of ~80% of E. coli phylogenetic group B2 multilocus sequence type 73 strains, as well as 22.5% of randomly selected UPEC strains isolated from community-acquired urinary tract infections in Denmark. In addition, E. coli CFT073 quiescence is not limited to glucose but occurs on agar plates containing a number of other sugars and acetate as sole carbon sources. It is also shown that a number of E. coli CFT073 mini-Tn5 metabolic mutants (gnd, gdhA, pykF, sdhA, and zwf) are nonquiescent on glucose M9 minimal agar plates and that quiescence requires a complete oxidative tricarboxylic acid (TCA) cycle. In addition, evidence is presented that, although E. coli CFT073 quiescence and persistence in the presence of ampicillin are alike in that both require a complete oxidative TCA cycle and each can be prevented by amino acids, E. coli CFT073 quiescence occurs in the presence or absence of a functional rpoS gene, whereas maximal persistence requires a nonfunctional rpoS. Our results suggest that interventions targeting specific central metabolic pathways may mitigate UPEC infections by interfering with quiescence and persistence. IMPORTANCE Recurrent urinary tract infections (UTIs) affect 10 to 40% of women. In up to 77% of those cases, the recurrent infections are caused by the same uropathogenic E. coli (UPEC) strain that caused the initial infection. Upon infection of urothelial transitional cells in the bladder, UPEC appear to enter a nongrowing quiescent intracellular state that is thought to serve as a reservoir responsible for recurrent UTIs. Here, we report that many UPEC strains enter a quiescent state when ≤106 CFU are seeded on glucose M9 minimal medium agar plates and show that mutations in several genes involved in central carbon metabolism prevent quiescence, as well as persistence, possibly identifying metabolic pathways involved in UPEC quiescence and persistence in vivo

The epidemic of extended-spectrum-beta-lactamase-producing Escherichia coli ST131 is driven by a single highly pathogenic subclone, H30-Rx

The Escherichia coli sequence type 131 (ST131) clone is notorious for extraintestinal infections, fluoroquinolone resistance, and extended-spectrum beta-lactamase (ESBL) production, attributable to a CTX-M-15-encoding mobile element. Here, we applied pulsed-field gel electrophoresis (PFGE) and whole-genome sequencing to reconstruct the evolutionary history of the ST131 clone. PFGE-based cluster analyses suggested that both fluoroquinolone resistance and ESBL production had been acquired by multiple ST131 sublineages through independent genetic events. In contrast, the more robust whole-genome-sequence-based phylogenomic analysis revealed that fluoroquinolone resistance was confined almost entirely to a single, rapidly expanding ST131 subclone, designated H30-R. Strikingly, 91% of the CTX-M-15-producing isolates also belonged to a single, well-defined clade nested within H30-R, which was named H30-Rx due to its more extensive resistance. Despite its tight clonal relationship with H30Rx, the CTX-M-15 mobile element was inserted variably in plasmid and chromosomal locations within the H30-Rx genome. Screening of a large collection of recent clinical E. coli isolates both confirmed the global clonal expansion of H30-Rx and revealed its disproportionate association with sepsis (relative risk, 7.5; P < 0.001). Together, these results suggest that the high prevalence of CTX-M-15 production among ST131 isolates is due primarily to the expansion of a single, highly virulent subclone, H30-Rx. IMPORTANCE We applied an advanced genomic approach to study the recent evolutionary history of one of the most important Escherichia coli strains in circulation today. This strain, called sequence type 131 (ST131), causes multidrug-resistant bladder, kidney, and bloodstream infections around the world. The rising prevalence of antibiotic resistance in E. coli is making these infections more difficult to treat and is leading to increased mortality. Past studies suggested that many different ST131 strains gained resistance to extended-spectrum cephalosporins independently. In contrast, our research indicates that most extended-spectrum-cephalosporin-resistant ST131 strains belong to a single highly pathogenic subclone, called H30-Rx. The clonal nature of H30-Rx may provide opportunities for vaccine or transmission prevention-based control strategies, which could gain importance as H30-Rx and other extraintestinal pathogenic E. coli subclones become resistant to our best antibiotics

The epidemic of extended-spectrum-β-lactamase-producing Escherichia coli ST131 is driven by a single highly pathogenic subclone, H30-Rx

The Escherichia coli sequence type 131 (ST131) clone is notorious for extraintestinal infections, fluoroquinolone resistance, and extended-spectrum beta-lactamase (ESBL) production, attributable to a CTX-M-15-encoding mobile element. Here, we applied pulsed-field gel electrophoresis (PFGE) and whole-genome sequencing to reconstruct the evolutionary history of the ST131 clone. PFGE-based cluster analyses suggested that both fluoroquinolone resistance and ESBL production had been acquired by multiple ST131 sublineages through independent genetic events. In contrast, the more robust whole-genome-sequence-based phylogenomic analysis revealed that fluoroquinolone resistance was confined almost entirely to a single, rapidly expanding ST131 subclone, designated H30-R. Strikingly, 91% of the CTX-M-15-producing isolates also belonged to a single, well-defined clade nested within H30-R, which was namedH30-Rx due to its more extensive resistance. Despite its tight clonal relationship with H30Rx, the CTX-M-15 mobile element was inserted variably in plasmid and chromosomal locations within the H30-Rx genome. Screening of a large collection of recent clinical E. coli isolates both confirmed the global clonal expansion of H30-Rx and revealed its disproportionate association with sepsis (relative risk, 7.5; P \u3c 0.001). Together, these results suggest that the high prevalence of CTX-M-15 production among ST131 isolates is due primarily to the expansion of a single, highly virulent subclone, H30-Rx

Footprint of Positive Selection in Treponema pallidum subsp. pallidum Genome Sequences Suggests Adaptive Microevolution of the Syphilis Pathogen

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 nonsynonymous 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

Evolution of Salmonella enterica Virulence via Point Mutations in the Fimbrial Adhesin

Whereas the majority of pathogenic Salmonella serovars are capable of infecting many different animal species, typically producing a self-limited gastroenteritis, serovars with narrow host-specificity exhibit increased virulence and their infections frequently result in fatal systemic diseases. In our study, a genetic and functional analysis of the mannose-specific type 1 fimbrial adhesin FimH from a variety of serovars of Salmonella enterica revealed that specific mutant variants of FimH are common in host-adapted (systemically invasive) serovars. We have found that while the low-binding shear-dependent phenotype of the adhesin is preserved in broad host-range (usually systemically non-invasive) Salmonella, the majority of host-adapted serovars express FimH variants with one of two alternative phenotypes: a significantly increased binding to mannose (as in S. Typhi, S. Paratyphi C, S. Dublin and some isolates of S. Choleraesuis), or complete loss of the mannose-binding activity (as in S. Paratyphi B, S. Choleraesuis and S. Gallinarum). The functional diversification of FimH in host-adapted Salmonella results from recently acquired structural mutations. Many of the mutations are of a convergent nature indicative of strong positive selection. The high-binding phenotype of FimH that leads to increased bacterial adhesiveness to and invasiveness of epithelial cells and macrophages usually precedes acquisition of the non-binding phenotype. Collectively these observations suggest that activation or inactivation of mannose-specific adhesive properties in different systemically invasive serovars of Salmonella reflects their dynamic trajectories of adaptation to a life style in specific hosts. In conclusion, our study demonstrates that point mutations are the target of positive selection and, in addition to horizontal gene transfer and genome degradation events, can contribute to the differential pathoadaptive evolution of Salmonella