The Conserved nhaAR Operon Is Drastically Divergent between B2 and Non-B2 Escherichia coli and Is Involved in Extra-Intestinal Virulence

The Escherichia coli species is divided in phylogenetic groups that differ in their virulence and commensal distribution. Strains belonging to the B2 group are involved in extra-intestinal pathologies but also appear to be more prevalent as commensals among human occidental populations. To investigate the genetic specificities of B2 sub-group, we used 128 sequenced genomes and identified genes of the core genome that showed marked difference between B2 and non-B2 genomes. We focused on the gene and its surrounding region with the strongest divergence between B2 and non-B2, the antiporter gene nhaA. This gene is part of the nhaAR operon, which is in the core genome but flanked by mobile regions, and is involved in growth at high pH and high sodium concentrations. Consistently, we found that a panel of non-B2 strains grew faster than B2 at high pH and high sodium concentrations. However, we could not identify differences in expression of the nhaAR operon using fluorescence reporter plasmids. Furthermore, the operon deletion had no differential impact between B2 and non-B2 strains, and did not result in a fitness modification in a murine model of gut colonization. Nevertheless, sequence analysis and experiments in a murine model of septicemia revealed that recombination in nhaA among B2 strains was observed in strains with low virulence. Finally, nhaA and nhaAR operon deletions drastically decreased virulence in one B2 strain. This effect of nhaAR deletion appeared to be stronger than deletion of all pathogenicity islands. Thus, a population genetic approach allowed us to identify an operon in the core genome without strong effect in commensalism but with an important role in extra-intestinal virulence, a landmark of the B2 strains

Coevolutionary Landscape Inference and the Context-Dependence of Mutations in Beta-Lactamase TEM-1

The quantitative characterization of mutational landscapes is a task of outstanding importance in evolutionary and medical biology: It is, for example, of central importance for our understanding of the phenotypic effect of mutations related to disease and antibiotic drug resistance. Here we develop a novel inference scheme for mutational landscapes, which is based on the statistical analysis of large alignments of homologs of the protein of interest. Our method is able to capture epistatic couplings between residues, and therefore to assess the dependence of mutational effects on the sequence context where they appear. Compared with recent large-scale mutagenesis data of the beta-lactamase TEM-1, a protein providing resistance against beta-lactam antibiotics, our method leads to an increase of about 40% in explicative power as compared with approaches neglecting epistasis. We find that the informative sequence context extends to residues at native distances of about 20 Å from the mutated site, reaching thus far beyond residues in direct physical contact

The Error and Repair Catastrophes: A Two-Dimensional Phase Diagram in the Quasispecies Model

This paper develops a two gene, single fitness peak model for determining the equilibrium distribution of genotypes in a unicellular population which is capable of genetic damage repair. The first gene, denoted by $\sigma_{via}$, yields a viable organism with first order growth rate constant $k > 1$ if it is equal to some target ``master'' sequence $\sigma_{via, 0}$. The second gene, denoted by $\sigma_{rep}$, yields an organism capable of genetic repair if it is equal to some target ``master'' sequence $\sigma_{rep, 0}$. This model is analytically solvable in the limit of infinite sequence length, and gives an equilibrium distribution which depends on \mu \equiv L\eps , the product of sequence length and per base pair replication error probability, and \eps_r , the probability of repair failure per base pair. The equilibrium distribution is shown to exist in one of three possible ``phases.'' In the first phase, the population is localized about the viability and repairing master sequences. As \eps_r exceeds the fraction of deleterious mutations, the population undergoes a ``repair'' catastrophe, in which the equilibrium distribution is still localized about the viability master sequence, but is spread ergodically over the sequence subspace defined by the repair gene. Below the repair catastrophe, the distribution undergoes the error catastrophe when $\mu$ exceeds \ln k/\eps_r , while above the repair catastrophe, the distribution undergoes the error catastrophe when $\mu$ exceeds $\ln k/f_{del}$, where $f_{del}$ denotes the fraction of deleterious mutations.Comment: 14 pages, 3 figures. Submitted to Physical Review

Interplay between pleiotropy and secondary selection determines rise and fall of mutators in stress response

Dramatic rise of mutators has been found to accompany adaptation of bacteria in response to many kinds of stress. Two views on the evolutionary origin of this phenomenon emerged: the pleiotropic hypothesis positing that it is a byproduct of environmental stress or other specific stress response mechanisms and the second order selection which states that mutators hitchhike to fixation with unrelated beneficial alleles. Conventional population genetics models could not fully resolve this controversy because they are based on certain assumptions about fitness landscape. Here we address this problem using a microscopic multiscale model, which couples physically realistic molecular descriptions of proteins and their interactions with population genetics of carrier organisms without assuming any a priori fitness landscape. We found that both pleiotropy and second order selection play a crucial role at different stages of adaptation: the supply of mutators is provided through destabilization of error correction complexes or fluctuations of production levels of prototypic mismatch repair proteins (pleiotropic effects), while rise and fixation of mutators occur when there is a sufficient supply of beneficial mutations in replication-controlling genes. This general mechanism assures a robust and reliable adaptation of organisms to unforeseen challenges. This study highlights physical principles underlying physical biological mechanisms of stress response and adaptation

Monotonicity of Fitness Landscapes and Mutation Rate Control

A common view in evolutionary biology is that mutation rates are minimised. However, studies in combinatorial optimisation and search have shown a clear advantage of using variable mutation rates as a control parameter to optimise the performance of evolutionary algorithms. Much biological theory in this area is based on Ronald Fisher's work, who used Euclidean geometry to study the relation between mutation size and expected fitness of the offspring in infinite phenotypic spaces. Here we reconsider this theory based on the alternative geometry of discrete and finite spaces of DNA sequences. First, we consider the geometric case of fitness being isomorphic to distance from an optimum, and show how problems of optimal mutation rate control can be solved exactly or approximately depending on additional constraints of the problem. Then we consider the general case of fitness communicating only partial information about the distance. We define weak monotonicity of fitness landscapes and prove that this property holds in all landscapes that are continuous and open at the optimum. This theoretical result motivates our hypothesis that optimal mutation rate functions in such landscapes will increase when fitness decreases in some neighbourhood of an optimum, resembling the control functions derived in the geometric case. We test this hypothesis experimentally by analysing approximately optimal mutation rate control functions in 115 complete landscapes of binding scores between DNA sequences and transcription factors. Our findings support the hypothesis and find that the increase of mutation rate is more rapid in landscapes that are less monotonic (more rugged). We discuss the relevance of these findings to living organisms

CRISPR/Cas9 recombineering-mediated deep mutational scanning of essential genes in Escherichia coli

Deep mutational scanning can provide significant insights into the function of essential genes in bacteria. Here, we developed a high-throughput method for mutating essential genes of Escherichia coli in their native genetic context. We used Cas9-mediated recombineering to introduce a library of mutations, created by error-prone PCR, within a gene fragment on the genome using a single gRNA pre-validated for high efficiency. Tracking mutation frequency through deep sequencing revealed biases in the position and the number of the introduced mutations. We overcame these biases by increasing the homology arm length and blocking mismatch repair to achieve a mutation efficiency of 85% for non-essential genes and 55% for essential genes. These experiments also improved our understanding of poorly characterized recombineering process using dsDNA donors with single nucleotide changes. Finally, we applied our technology to target rpoB, the beta subunit of RNA polymerase, to study resistance against rifampicin. In a single experiment, we validate multiple biochemical and clinical observations made in the previous decades and provide insights into resistance compensation with the study of double mutants.</p

Differential Dynamics of Transposable Elements during Long-Term Diploidization of Nicotiana Section Repandae (Solanaceae) Allopolyploid Genomes

PubMed ID: 23185607This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited

Multidrug-resistant Escherichia coli from canine urinary tract infections tend to have commensal phylotypes, lower prevalence of virulence determinants and ampC-replicons

AbstractMultidrug-resistant Escherichia coli is an emerging clinical challenge in domestic species. Treatment options in many cases are limited. This study characterized MDR E. coli isolates from urinary tract infections in dogs, collected between 2002 and 2011. Isolates were evaluated in terms of β-lactamase production, phylogenetic group, ST type, replicon type and virulence marker profile. Comparisons were made with antibiotic susceptible isolates also collected from dogs with urinary tract infections. AmpC β-lactamase was produced in 67% of the MDR isolates (12/18). Of these, 8 could be specifically attributed to the CMY-2 gene. None of the isolates tested in either group expressed ESBLs. Phylo-group distribution was as expected in the susceptible isolates, with an over representation of the pathogenic B2 phylo-group (67%). In contrast, the phylogenetic background for the MDR group was mixed, with representation of commensal phylo-groups A and B1. The B2 phylo-group represented the smallest proportion (A, B1, B2 or D was 28%, 22%, 11% and 33%, respectively). Virulence marker profiles, evaluated using Identibac® microarray, discriminated between the two groups. Marker sequences for a core panel of virulence determinants were identified in most of the susceptible isolates, but not in most of the MDR isolates. These findings indicate that for MDR isolates, plasmid-mediated AmpC is an important resistance mechanism, and while still capable of causing clinical disease, there is evidence for a shift towards phylogenetic groups of reduced inferred virulence potential. There was no evidence of zoonotic potential in either the susceptible or MDR urinary tract isolates in this study

The genome-wide dynamics of purging during selfing in maize

Self-fertilization (also known as selfing) is an important reproductive strategy in plants and a widely applied tool for plant genetics and plant breeding. Selfing can lead to inbreeding depression by uncovering recessive deleterious variants, unless these variants are purged by selection. Here we investigated the dynamics of purging in a set of eleven maize lines that were selfed for six generations. We show that heterozygous, putatively deleterious single nucleotide polymorphisms are preferentially lost from the genome during selfing. Deleterious single nucleotide polymorphisms were lost more rapidly in regions of high recombination, presumably because recombination increases the efficacy of selection by uncoupling linked variants. Overall, heterozygosity decreased more slowly than expected, by an estimated 35% to 40% per generation instead of the expected 50%, perhaps reflecting pervasive associative overdominance. Finally, three lines exhibited marked decreases in genome size due to the purging of transposable elements. Genome loss was more likely to occur for lineages that began with larger genomes with more transposable elements and chromosomal knobs. These three lines purged an average of 398 Mb from their genomes, an amount equivalent to three Arabidopsis thaliana genomes per lineage, in only a few generations

DNA Barcoding the Geometrid Fauna of Bavaria (Lepidoptera): Successes, Surprises, and Questions

BACKGROUND: The State of Bavaria is involved in a research program that will lead to the construction of a DNA barcode library for all animal species within its territorial boundaries. The present study provides a comprehensive DNA barcode library for the Geometridae, one of the most diverse of insect families. METHODOLOGY/PRINCIPAL FINDINGS: This study reports DNA barcodes for 400 Bavarian geometrid species, 98 per cent of the known fauna, and approximately one per cent of all Bavarian animal species. Although 98.5% of these species possess diagnostic barcode sequences in Bavaria, records from neighbouring countries suggest that species-level resolution may be compromised in up to 3.5% of cases. All taxa which apparently share barcodes are discussed in detail. One case of modest divergence (1.4%) revealed a species overlooked by the current taxonomic system: Eupithecia goossensiata Mabille, 1869 stat.n. is raised from synonymy with Eupithecia absinthiata (Clerck, 1759) to species rank. Deep intraspecific sequence divergences (>2%) were detected in 20 traditionally recognized species. CONCLUSIONS/SIGNIFICANCE: The study emphasizes the effectiveness of DNA barcoding as a tool for monitoring biodiversity. Open access is provided to a data set that includes records for 1,395 geometrid specimens (331 species) from Bavaria, with 69 additional species from neighbouring regions. Taxa with deep intraspecific sequence divergences are undergoing more detailed analysis to ascertain if they represent cases of cryptic diversity