Mathematical Models of Plasmid Population Dynamics.

With plasmid-mediated antibiotic resistance thriving and threatening to become a serious public health problem, it is paramount to increase our understanding of the forces that enable the spread and maintenance of drug resistance genes encoded in mobile genetic elements. The relevance of plasmids as vehicles for the dissemination of antibiotic resistance genes, in addition to the extensive use of plasmid-derived vectors for biotechnological and industrial purposes, has promoted the in-depth study of the molecular mechanisms controlling multiple aspects of a plasmids' life cycle. This body of experimental work has been paralleled by the development of a wealth of mathematical models aimed at understanding the interplay between transmission, replication, and segregation, as well as their consequences in the ecological and evolutionary dynamics of plasmid-bearing bacterial populations. In this review, we discuss theoretical models of plasmid dynamics that span from the molecular mechanisms of plasmid partition and copy-number control occurring at a cellular level, to their consequences in the population dynamics of complex microbial communities. We conclude by discussing future directions for this exciting research topic.

Respiratory Heme A-Containing Oxidases Originated in the Ancestors of Iron-Oxidizing Bacteria.

Respiration is a major trait shaping the biology of many environments. Cytochrome oxidase containing heme A (COX) is a common terminal oxidase in aerobic bacteria and is the only one in mammalian mitochondria. The synthesis of heme A is catalyzed by heme A synthase (CtaA/Cox15), an enzyme that most likely coevolved with COX. The evolutionary origin of COX in bacteria has remained unknown. Using extensive sequence and phylogenetic analysis, we show that the ancestral type of heme A synthases is present in iron-oxidizing Proteobacteria such as Acidithiobacillus spp. These bacteria also contain a deep branching form of the major COX subunit (COX1) and an ancestral variant of CtaG, a protein that is specifically required for COX biogenesis. Our work thus suggests that the ancestors of extant iron-oxidizers were the first to evolve COX. Consistent with this conclusion, acidophilic iron-oxidizing prokaryotes lived on emerged land around the time for which there is the earliest geochemical evidence of aerobic respiration on earth. Hence, ecological niches of iron oxidation have apparently promoted the evolution of aerobic respiration.

An object-oriented framework for evolutionary pangenome analysis.

Pangenome analysis is fundamental to explore molecular evolution occurring in bacterial populations. Here, we introduce Pagoo, an R framework that enables straightforward handling of pangenome data. The encapsulated nature of Pagoo allows the storage of complex molecular and phenotypic information using an object-oriented approach. This facilitates to go back and forward to the data using a single programming environment and saving any stage of analysis (including the raw data) in a single file, making it sharable and reproducible. Pagoo provides tools to query, subset, compare, visualize, and perform statistical analyses, in concert with other microbial genomics packages available in the R ecosystem. As working examples, we used 1,000 Escherichia coli genomes to show that Pagoo is scalable, and a global dataset of Campylobacter fetus genomes to identify evolutionary patterns and genomic markers of host-adaptation in this pathogen.

Evolution of bacteria seen through their essential genes: the case of Pseudomonas aeruginosa and Azotobacter vinelandii.

Pseudomonas aeruginosa is a metabolically versatile bacterium and also an important opportunistic pathogen. It has a remarkable genomic structure since the genetic information encoding its pathogenicity-related traits belongs to its core-genome while both environmental and clinical isolates are part of the same population with a highly conserved genomic sequence. Unexpectedly, considering the high level of sequence identity and homologue gene number shared between different P. aeruginosa isolates, the presence of specific essential genes of the two type strains PAO1 and PA14 has been reported to be highly variable. Here we report the detailed bioinformatics analysis of the essential genes of P. aeruginosa PAO1 and PA14 that have been previously experimentally identified and show that the reported gene variability was owed to sequencing and annotation inconsistencies, but that in fact they are highly conserved. This bioinformatics analysis led us to the definition of 348 P. aeruginosa general essential genes. In addition we show that 342 of these 348 essential genes are conserved in Azotobacter vinelandii, a nitrogen-fixing, cyst-forming, soil bacterium. These results support the hypothesis of A. vinelandii having a polyphyletic origin with a Pseudomonads genomic backbone, and are a challenge to the accepted theory of bacterial evolution.

Diversification of Vibrio anguillarum Driven by the Bacteriophage CHOED.

Bacteriophages are an important factor in bacterial evolution. Some reports suggest that lytic bacteriophages can select for resistant mutant strains with reduced virulence. The present study explores the role of the CHOED bacteriophage in the diversification and virulence of its host Vibrio anguillarum. Nine phage-resistant strains were analyzed for their phenotype and different virulence factors, showing alterations in their fitness, motility, biofilm formation, lipopolysaccharide profiles and/or protease activity. Seven of the nine phage-resistant strains showed virulence reduction in a Sparus aurata larvae model. However, this is not generalized since two of the resistant strains show equal virulence compared with the parental strain. The genomic analysis of representative resistant strains displayed that the majority of the mutations are specific for each isolate, affecting genes related to lipopolysaccharide biosynthesis, quorum sensing, motility, toxin and membrane transport. The observed mutations were coherent with the phenotypic and virulence differences observed. These results suggest that the CHOED phage acts as a selective pressure on V. anguillarum, allowing proliferation of resistant strains with different genotypes, phenotypes and degrees of virulence, contributing to bacterial diversification.

Oxygen Reductases in Alphaproteobacterial Genomes: Physiological Evolution From Low to High Oxygen Environments.

Oxygen reducing terminal oxidases differ with respect to their subunit composition, heme groups, operon structure, and affinity for O2. Six families of terminal oxidases are currently recognized, all of which occur in alphaproteobacterial genomes, two of which are also present in mitochondria. Many alphaproteobacteria encode several different terminal oxidases, likely reflecting ecological versatility with respect to oxygen levels. Terminal oxidase evolution likely started with the advent of O2 roughly 2.4 billion years ago and terminal oxidases diversified in the Proterozoic, during which oxygen levels remained low, around the Pasteur point (ca. 2 µM O2). Among the alphaproteobacterial genomes surveyed, those from members of the Rhodospirillaceae reveal the greatest diversity in oxygen reductases. Some harbor all six terminal oxidase types, in addition to many soluble enzymes typical of anaerobic fermentations in mitochondria and hydrogenosomes of eukaryotes. Recent data have it that O2 levels increased to current values (21% v/v or ca. 250 µM) only about 430 million years ago. Ecological adaptation brought forth different lineages of alphaproteobacteria and different lineages of eukaryotes that have undergone evolutionary specialization to high oxygen, low oxygen, and anaerobic habitats. Some have remained facultative anaerobes that are able to generate ATP with or without the help of oxygen and represent physiological links to the ancient proteobacterial lineage at the origin of mitochondria and eukaryotes. Our analysis reveals that the genomes of alphaproteobacteria appear to retain signatures of ancient transitions in aerobic metabolism, findings that are relevant to mitochondrial evolution in eukaryotes as well.

Genomic Epidemiology of NDM-1-Encoding Plasmids in Latin American Clinical Isolates Reveals Insights into the Evolution of Multidrug Resistance.

Bacteria that produce the broad-spectrum Carbapenem antibiotic New Delhi Metallo-ß-lactamase (NDM) place a burden on health care systems worldwide, due to the limited treatment options for infections caused by them and the rapid global spread of this antibiotic resistance mechanism. Although it is believed that the associated resistance gene blaNDM-1 originated in Acinetobacter spp., the role of Enterobacteriaceae in its dissemination remains unclear. In this study, we used whole genome sequencing to investigate the dissemination dynamics of blaNDM-1-positive plasmids in a set of 21 clinical NDM-1-positive isolates from Colombia and Mexico (Providencia rettgeri, Klebsiella pneumoniae, and Acinetobacter baumannii) as well as six representative NDM-1-positive Escherichia coli transconjugants. Additionally, the plasmids from three representative P. rettgeri isolates were sequenced by PacBio sequencing and finished. Our results demonstrate the presence of previously reported plasmids from K. pneumoniae and A. baumannii in different genetic backgrounds and geographically distant locations in Colombia. Three new previously unclassified plasmids were also identified in P. rettgeri from Colombia and Mexico, plus an interesting genetic link between NDM-1-positive P. rettgeri from distant geographic locations (Canada, Mexico, Colombia, and Israel) without any reported epidemiological links was discovered. Finally, we detected a relationship between plasmids present in P. rettgeri and plasmids from A. baumannii and K. pneumoniae. Overall, our findings suggest a Russian doll model for the dissemination of blaNDM-1 in Latin America, with P. rettgeri playing a central role in this process, and reveal new insights into the evolution and dissemination of plasmids carrying such antibiotic resistance genes.

The Transcriptional Regulators of the CRP Family Regulate Different Essential Bacterial Functions and Can Be Inherited Vertically and Horizontally.

One of the best-studied transcriptional regulatory proteins in bacteria is the Escherichia coli catabolite repressor protein (CRP) that when complexed with 3'-5'-cyclic AMP (cAMP) changes its conformation and interacts with specific DNA-sequences. CRP DNA-binding can result in positive or negative regulation of gene expression depending on the position of its interaction with respect to RNA polymerase binding site. The aim of this work is to review the biological role and phylogenetic relations that some members of the CRP family of transcriptional regulators (also known as cAMP receptor protein family) have in different bacterial species. This work is not intended to give an exhaustive revision of bacterial CRP-orthologs, but to provide examples of the role that these proteins play in the expression of genes that are fundamental for the life style of some bacterial species. We highlight the conservation of their structural characteristics and of their binding to conserved-DNA sequences, in contrast to their very diverse repertoire of gene activation. CRP activates a wide variety of fundamental genes for the biological characteristic of each bacterial species, which in several instances form part of their core-genome (defined as the gene sequences present in all members of a bacterial species). We present evidence that support the fact that some of the transcriptional regulators that belong to the CRP family in different bacterial species, and some of the genes that are regulated by them, can be inherited by horizontal gene transfer. These data are discussed in the framework of bacterial evolution models.

Genome Analysis of Structure-Function Relationships in Respiratory Complex I, an Ancient Bioenergetic Enzyme.

Respiratory complex I (NADH:ubiquinone oxidoreductase) is a ubiquitous bioenergetic enzyme formed by over 40 subunits in eukaryotes and a minimum of 11 subunits in bacteria. Recently, crystal structures have greatly advanced our knowledge of complex I but have not clarified the details of its reaction with ubiquinone (Q). This reaction is essential for bioenergy production and takes place in a large cavity embedded within a conserved module that is homologous to the catalytic core of Ni-Fe hydrogenases. However, how a hydrogenase core has evolved into the protonmotive Q reductase module of complex I has remained unclear. This work has exploited the abundant genomic information that is currently available to deduce structure-function relationships in complex I that indicate the evolutionary steps of Q reactivity and its adaptation to natural Q substrates. The results provide answers to fundamental questions regarding various aspects of complex I reaction with Q and help re-defining the old concept that this reaction may involve two Q or inhibitor sites. The re-definition leads to a simplified classification of the plethora of complex I inhibitors while throwing a new light on the evolution of the enzyme function.

PHYLOGENETIC ESTIMATION OF PLASMID EXCHANGE IN BACTERIA.

The existence of differential horizontal gene transfer may be assessed by comparing the phylogenetic trees derived from two different genes. We use this concept to estimate quantitatively the amount of plasmid exchange that has occurred in a bacterial population. By means of computer simulations we studied the effect of gene transfer on the topological distortion between two phylogenetic trees: one obtained from an euchromosomal gene and another from a plasmid-borne sequence, which may be subjected to horizontal transfer. The basic assumptions of our simulations were (a) that plasmid exchange had occurred recently (after the last population split); and (b) that either the amount of chromosomal horizontal exchange was negligible or that it was only a fraction of the amount of plasmid exchange in which case we will be estimating relative amounts of plasmid transfer. We found that the topological difference between two such trees is a function of the number of plasmid exchange events that have occurred. It can be explained by a logistic model that relates the average distortion index between two trees (dT ) to the number of transfer events (x). The behavior remains the same under different conditions that were tested (symmetry of the topology, number of taxa in the tree, effect of reconstruction errors, mutation after plasmid transfer). We have also tried our method on empirical data from the literature and estimated the amount of gene transfer that may have occurred among Sym plasmids in agricultural field populations of Rhizobium leguminosarum biovar phaseoli. We found that between 15.77 to 29.98% of all genetic types in these populations have been either the source or the target of a plasmid transfer event. When the comparisons were made among trees derived exclusively from plasmid probes this value dropped to 2.00%. Phylogenetic trees derived from symbiotic and nonsymbiotic sequences were also used to infer the number of gene transfer events among 11 isolates from R. galegae. The estimated number of transfer events of symbiotic sequences was 10.515 (although we do not know out of how many genetic types). We concluded that intraspecific transfer of symbiotic sequences is widespread in these two species of the genus Rhizobium.