Macarena Toll-Riera

Positive selection and compensatory adaptation interact to stabilize non-transmissible plasmids

Plasmids are important drivers of bacterial evolution, but it is challenging to understand how plasmids persist over the long term because plasmid carriage is costly. Classical models predict that horizontal transfer is necessary for plasmid persistence, but recent work shows that almost half of plasmids are non-transmissible. Here we use a combination of mathematical modelling and experimental evolution to investigate how a costly, non-transmissible plasmid, pNUK73, can be maintained in populations of Pseudomonas aeruginosa. Compensatory adaptation increases plasmid stability by eliminating the cost of plasmid carriage. However, positive selection for plasmid-encoded antibiotic resistance is required to maintain the plasmid by offsetting reductions in plasmid frequency due to segregational loss. Crucially, we show that compensatory adaptation and positive selection reinforce each other’s effects. Our study provides a new understanding of how plasmids persist in bacterial populations, and it helps to explain why resistance can be maintained after antibiotic use is stopped.

Natural selection drives the accumulation of amino acid tandem repeats in human proteins

Amino acid tandem repeats are found in a large number of eukaryotic proteins. They are often encoded by trinucleotide repeats and exhibit high intra- and interspecies size variability due to the high mutation rate associated with replication slippage. The extent to which natural selection is important in shaping amino acid repeat evolution is a matter of debate. On one hand, their high frequency may simply reflect their high probability of expansion by slippage, and they could essentially evolve in a neutral manner. On the other hand, there is experimental evidence that changes in repeat size can influence protein-protein interactions, transcriptional activity, or protein subcellular localization, indicating that repeats could be functionally relevant and thus shaped by selection. To gauge the relative contribution of neutral and selective forces in amino acid repeat evolution, we have performed a comparative analysis of amino acid repeat conservation in a large set of orthologous proteins from 12 vertebrate species. As a neutral model of repeat evolution we have used sequences with the same DNA triplet composition as the coding sequences--and thus expected to be subject to the same mutational forces--but located in syntenic noncoding genomic regions. The results strongly indicate that selection has played a more important role than previously suspected in amino acid tandem repeat evolution, by increasing the repeat retention rate and by modulating repeat size. The data obtained in this study have allowed us to identify a set of 92 repeats that are postulated to play important functional roles due to their strong selective signature, including five cases with direct experimental evidence.

Interactions between horizontally acquired genes create a fitness cost in Pseudomonas aeruginosa.

Horizontal gene transfer (HGT) plays a key role in bacterial evolution, especially with respect to antibiotic resistance. Fitness costs associated with mobile genetic elements (MGEs) are thought to constrain HGT, but our understanding of these costs remains fragmentary, making it difficult to predict the success of HGT events. Here we use the interaction between P. aeruginosa and a costly plasmid (pNUK73) to investigate the molecular basis of the cost of HGT. Using RNA-Seq, we show that the acquisition of pNUK73 results in a profound alteration of the transcriptional profile of chromosomal genes. Mutations that inactivate two genes encoded on chromosomally integrated MGEs recover these fitness costs and transcriptional changes by decreasing the expression of the pNUK73 replication gene. Our study demonstrates that interactions between MGEs can compromise bacterial fitness via altered gene expression, and we argue that conflicts between mobile elements impose a general constraint on evolution by HGT.

Role of low-complexity sequences in the formation of novel protein coding sequences

Low-complexity sequences are extremely abundant in eukaryotic proteins for reasons that remain unclear. One hypothesis is that they contribute to the formation of novel coding sequences, facilitating the generation of novel protein functions. Here, we test this hypothesis by examining the content of low-complexity sequences in proteins of different age. We show that recently emerged proteins contain more low-complexity sequences than older proteins and that these sequences often form functional domains. These data are consistent with the idea that low-complexity sequences may play a key role in the emergence of novel genes.

Lineage-Specific Variation in Intensity of Natural Selection in Mammals

The molecular clock hypothesis states that protein-coding genes evolve at an approximately constant rate. However, this is only expected to be true as long as the function and the tertiary structure of the molecule remain unaltered. An important implication of this statement is that significant deviations in the rate of evolution of a gene with respect to the species clock are likely to reflect functional and/or structural alterations. Here, we present a method to identify such deviations and apply it to a data set of 2,929 high-quality coding sequence alignments corresponding to one-to-one orthologous genes from six mammalian species--human, macaque, mouse, rat, cow, and dog. Deviated branches are defined as those that present significant alterations in both the rate of nonsynonymous substitutions (dN) and the selective pressure (dN/dS). Strikingly, we find that as many as 24.5% of the genes show branch-specific deviations in dN and dN/dS, though this is a relatively well-conserved set of genes. Around half of these genes show branch-specific acceleration of evolutionary rates. Positive selection (PS) tests based on divergence data only identify 17.7% of the accelerated branches. Failure to identify PS in accelerated branches with an excess of radical amino acid replacements suggests that these tests are conservative. Interestingly, genes with accelerated branches are significantly enriched in neural proteins, indicating that this type of protein might play a more important role than previously thought in species diversification, although they are generally not detected by PS tests. We discuss in detail several examples of genes that show lineage-specific evolutionary rate acceleration and are involved in synaptic transmission, chemosensory perception, and ubiquitination.

Fitness Is Strongly Influenced by Rare Mutations of Large Effect in a Microbial Mutation Accumulation Experiment

Our understanding of the evolutionary consequences of mutation relies heavily on estimates of the rate and fitness effect of spontaneous mutations generated by mutation accumulation (MA) experiments. We performed a classic MA experiment in which frequent sampling of MA lines was combined with whole genome resequencing to develop a high-resolution picture of the effect of spontaneous mutations in a hypermutator (ΔmutS) strain of the bacterium Pseudomonas aeruginosa. After ∼644 generations of mutation accumulation, MA lines had accumulated an average of 118 mutations, and we found that average fitness across all lines decayed linearly over time. Detailed analyses of the dynamics of fitness change in individual lines revealed that a large fraction of the total decay in fitness (42.3%) was attributable to the fixation of rare, highly deleterious mutations (comprising only 0.5% of fixed mutations). Furthermore, we found that at least 0.64% of mutations were beneficial and probably fixed due to positive selection. The majority of mutations that fixed (82.4%) were base substitutions and we failed to find any signatures of selection on nonsynonymous or intergenic mutations. Short indels made up a much smaller fraction of the mutations that were fixed (17.4%), but we found evidence of strong selection against indels that caused frameshift mutations in coding regions. These results help to quantify the amount of natural selection present in microbial MA experiments and demonstrate that changes in fitness are strongly influenced by rare mutations of large effect.

Emergence of novel domains in proteins

BackgroundProteins are composed of a combination of discrete, well-defined, sequence domains, associated with specific functions that have arisen at different times during evolutionary history. The emergence of novel domains is related to protein functional diversification and adaptation. But currently little is known about how novel domains arise and how they subsequently evolve.ResultsTo gain insights into the impact of recently emerged domains in protein evolution we have identified all human young protein domains that have emerged in approximately the past 550 million years. We have classified them into vertebrate-specific and mammalian-specific groups, and compared them to older domains. We have found 426 different annotated young domains, totalling 995 domain occurrences, which represent about 12.3% of all human domains. We have observed that 61.3% of them arose in newly formed genes, while the remaining 38.7% are found combined with older domains, and have very likely emerged in the context of a previously existing protein. Young domains are preferentially located at the N-terminus of the protein, indicating that, at least in vertebrates, novel functional sequences often emerge there. Furthermore, young domains show significantly higher non-synonymous to synonymous substitution rates than older domains using human and mouse orthologous sequence comparisons. This is also true when we compare young and old domains located in the same protein, suggesting that recently arisen domains tend to evolve in a less constrained manner than older domains.ConclusionsWe conclude that proteins tend to gain domains over time, becoming progressively longer. We show that many proteins are made of domains of different age, and that the fastest evolving parts correspond to the domains that have been acquired more recently.

Sequence shortening in the rodent ancestor

Insertions and deletions (indels), together with nucleotide substitutions, are major drivers of sequence evolution. An excess of deletions over insertions in genomic sequences-the so-called deletional bias-has been reported in a wide range of species, including mammals. However, this bias has not been found in the coding sequences of some mammalian species, such as human and mouse. To determine the strength of the deletional bias in mammals, and the influence of mutation and selection, we have quantified indels in both neutrally evolving noncoding sequences and protein-coding sequences, in six mammalian branches: human, macaque, ancestral primate, mouse, rat, and ancestral rodent. The results obtained with an improved algorithm for the placement of insertions in multiple alignments, Prank(+F), indicate that contrary to previous results, the only mammalian branch with a strong deletional bias is the rodent ancestral branch. We estimate that such a bias has resulted in an ~2.5% sequence loss of mammalian syntenic region in the ancestor of the mouse and rat. Further, a comparison of coding and noncoding sequences shows that negative selection is acting more strongly against mutations generating amino acid insertions than against mutations resulting in amino acid deletions. The strength of selection against indels is found to be higher in the rodent branches than in the primate branches, consistent with the larger effective population sizes of the rodents.

Sequencing of plasmids pAMBL1 and pAMBL2 from Pseudomonas aeruginosa reveals a blaVIM-1 amplification causing high-level carbapenem resistance

BACKGROUND Carbapenemases are a major concern for the treatment of infectious diseases caused by Gram-negative bacteria. Although plasmids are responsible for the spread of resistance genes among these pathogens, there is limited information on the nature of the mobile genetic elements carrying carbapenemases in Pseudomonas aeruginosa. METHODS We combined data from two different next-generation sequencing platforms, Illumina HiSeq2000 and PacBio RSII, to obtain the complete nucleotide sequences of two blaVIM-1-carrying plasmids (pAMBL1 and pAMBL2) isolated from P. aeruginosa clinical isolates. RESULTS Plasmid pAMBL1 has 26 440 bp and carries a RepA_C family replication protein. pAMBL1 is similar to plasmids pNOR-2000 and pKLC102 from P. aeruginosa and pAX22 from Achromobacter xylosoxidans, which also carry VIM-type carbapenemases. pAMBL2 is a 24 133 bp plasmid with a replication protein that belongs to the Rep_3 family. It shows a high degree of homology with a fragment of the blaVIM-1-bearing plasmid pPC9 from Pseudomonas putida. Plasmid pAMBL2 carries three copies of the blaVIM-1 cassette in an In70 class 1 integron conferring, unlike pAMBL1, high-level resistance to carbapenems. CONCLUSIONS We present two new plasmids coding for VIM-1 carbapenemase from P. aeruginosa and report that the presence of three copies of blaVIM-1 in pAMBL2 produces high-level resistance to carbapenems.

The Genomic Basis of Evolutionary Innovation in Pseudomonas aeruginosa.

Novel traits play a key role in evolution, but their origins remain poorly understood. Here we address this problem by using experimental evolution to study bacterial innovation in real time. We allowed 380 populations of Pseudomonas aeruginosa to adapt to 95 different carbon sources that challenged bacteria with either evolving novel metabolic traits or optimizing existing traits. Whole genome sequencing of more than 80 clones revealed profound differences in the genetic basis of innovation and optimization. Innovation was associated with the rapid acquisition of mutations in genes involved in transcription and metabolism. Mutations in pre-existing duplicate genes in the P. aeruginosa genome were common during innovation, but not optimization. These duplicate genes may have been acquired by P. aeruginosa due to either spontaneous gene amplification or horizontal gene transfer. High throughput phenotype assays revealed that novelty was associated with increased pleiotropic costs that are likely to constrain innovation. However, mutations in duplicate genes with close homologs in the P. aeruginosa genome were associated with low pleiotropic costs compared to mutations in duplicate genes with distant homologs in the P. aeruginosa genome, suggesting that functional redundancy between duplicates facilitates innovation by buffering pleiotropic costs.