Genetic dominance governs the evolution and spread of mobile genetic elements in bacteria.

Mobile genetic elements (MGEs), such as plasmids, promote bacterial evolution through horizontal gene transfer (HGT). However, the rules governing the repertoire of traits encoded on MGEs remain unclear. In this study, we uncovered the central role of genetic dominance shaping genetic cargo in MGEs, using antibiotic resistance as a model system. MGEs are typically present in more than one copy per host bacterium, and as a consequence, genetic dominance favors the fixation of dominant mutations over recessive ones. In addition, genetic dominance also determines the phenotypic effects of horizontally acquired MGE-encoded genes, silencing recessive alleles if the recipient bacterium already carries a wild-type copy of the gene. The combination of these two effects governs the catalog of genes encoded on MGEs. Our results help to understand how MGEs evolve and spread, uncovering the neglected influence of genetic dominance on bacterial evolution. Moreover, our findings offer a framework to forecast the spread and evolvability of MGE-encoded genes, which encode traits of key human interest, such as virulence or antibiotic resistance.

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.

The genomic basis of adaptation to the fitness cost of rifampicin resistance in Pseudomonas aeruginosa.

Antibiotic resistance carries a fitness cost that must be overcome in order for resistance to persist over the long term. Compensatory mutations that recover the functional defects associated with resistance mutations have been argued to play a key role in overcoming the cost of resistance, but compensatory mutations are expected to be rare relative to generally beneficial mutations that increase fitness, irrespective of antibiotic resistance. Given this asymmetry, population genetics theory predicts that populations should adapt by compensatory mutations when the cost of resistance is large, whereas generally beneficial mutations should drive adaptation when the cost of resistance is small. We tested this prediction by determining the genomic mechanisms underpinning adaptation to antibiotic-free conditions in populations of the pathogenic bacterium Pseudomonas aeruginosa that carry costly antibiotic resistance mutations. Whole-genome sequencing revealed that populations founded by high-cost rifampicin-resistant mutants adapted via compensatory mutations in three genes of the RNA polymerase core enzyme, whereas populations founded by low-cost mutants adapted by generally beneficial mutations, predominantly in the quorum-sensing transcriptional regulator gene lasR. Even though the importance of compensatory evolution in maintaining resistance has been widely recognized, our study shows that the roles of general adaptation in maintaining resistance should not be underestimated and highlights the need to understand how selection at other sites in the genome influences the dynamics of resistance alleles in clinical settings.

Here's to the losers: evolvable residents accelerate the evolution of high-fitness invaders.

Recent work has shown that evolvability plays a key role in determining the long-term population dynamics of asexual clones. However, simple considerations suggest that the evolvability of a focal lineage of bacteria should also be influenced by the evolvability of its competitors. First, evolvable competitors should accelerate evolution by impeding the fixation of the focal lineage through a clonal interference-like mechanism. Second, evolvable competitors should increase the strength of selection by rapidly degrading the environment, increasing selection for adaptive mutations. Here we tested these ideas by allowing a high-fitness clone of the bacterium Pseudomonas aeruginosa to invade populations of two low-fitness resident clones that differ in their evolvability. Both competition from mutations in the resident lineage and environmental degradation lead to faster adaptation in the invader through fixing single mutations with a greater fitness advantage. The results suggest that competition from mutations in both the successful invader and the unsuccessful resident shapes the adaptive trajectory of the invader through both direct competition and indirect environmental effects. Therefore, to predict evolutionary outcomes, it will be necessary to consider the evolvability of all members of the community and the effects of adaptation on the quality of the environment. This is particularly relevant to mixed microbial communities where lineages differ in their adaptive potential, a common feature of chronic infections.

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.

Emergence of novel domains in proteins.

BACKGROUND: Proteins 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. RESULTS: To 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. CONCLUSIONS: We 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.

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.

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.

Structure and age jointly influence rates of protein evolution.

What factors determine a protein's rate of evolution are actively debated. Especially unclear is the relative role of intrinsic factors of present-day proteins versus historical factors such as protein age. Here we study the interplay of structural properties and evolutionary age, as determinants of protein evolutionary rate. We use a large set of one-to-one orthologs between human and mouse proteins, with mapped PDB structures. We report that previously observed structural correlations also hold within each age group - including relationships between solvent accessibility, designabililty, and evolutionary rates. However, age also plays a crucial role: age modulates the relationship between solvent accessibility and rate. Additionally, younger proteins, despite being less designable, tend to evolve faster than older proteins. We show that previously reported relationships between age and rate cannot be explained by structural biases among age groups. Finally, we introduce a knowledge-based potential function to study the stability of proteins through large-scale computation. We find that older proteins are more stable for their native structure, and more robust to mutations, than younger ones. Our results underscore that several determinants, both intrinsic and historical, can interact to determine rates of protein evolution.

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.

A limit on the evolutionary rescue of an Antarctic bacterium from rising temperatures.

Climate change is gradual, but it can also cause brief extreme heat waves that can exceed the upper thermal limit of any one organism. To study the evolutionary potential of upper thermal tolerance, we evolved the cold-adapted Antarctic bacterium Pseudoalteromonas haloplanktis to survive at 30°C, beyond its ancestral thermal limit. This high-temperature adaptation occurred rapidly and in multiple populations. It involved genomic changes that occurred in a highly parallel fashion and mitigated the effects of protein misfolding. However, it also confronted a physiological limit, because populations failed to grow beyond 30°C. Our experiments aimed to facilitate evolutionary rescue by using a small organism with large populations living at temperatures several degrees below their upper thermal limit. Larger organisms with smaller populations and living at temperatures closer to their upper thermal tolerances are even more likely to go extinct during extreme heat waves.

Staphylococcal phages and pathogenicity islands drive plasmid evolution.

Conjugation has classically been considered the main mechanism driving plasmid transfer in nature. Yet bacteria frequently carry so-called non-transmissible plasmids, raising questions about how these plasmids spread. Interestingly, the size of many mobilisable and non-transmissible plasmids coincides with the average size of phages (~40 kb) or that of a family of pathogenicity islands, the phage-inducible chromosomal islands (PICIs, ~11 kb). Here, we show that phages and PICIs from Staphylococcus aureus can mediate intra- and inter-species plasmid transfer via generalised transduction, potentially contributing to non-transmissible plasmid spread in nature. Further, staphylococcal PICIs enhance plasmid packaging efficiency, and phages and PICIs exert selective pressures on plasmids via the physical capacity of their capsids, explaining the bimodal size distribution observed for non-conjugative plasmids. Our results highlight that transducing agents (phages, PICIs) have important roles in bacterial plasmid evolution and, potentially, in antimicrobial resistance transmission.

Origin of primate orphan genes: a comparative genomics approach.

Genomes contain a large number of genes that do not have recognizable homologues in other species and that are likely to be involved in important species-specific adaptive processes. The origin of many such "orphan" genes remains unknown. Here we present the first systematic study of the characteristics and mechanisms of formation of primate-specific orphan genes. We determine that codon usage values for most orphan genes fall within the bulk of the codon usage distribution of bona fide human proteins, supporting their current protein-coding annotation. We also show that primate orphan genes display distinctive features in relation to genes of wider phylogenetic distribution: higher tissue specificity, more rapid evolution, and shorter peptide size. We estimate that around 24% are highly divergent members of mammalian protein families. Interestingly, around 53% of the orphan genes contain sequences derived from transposable elements (TEs) and are mostly located in primate-specific genomic regions. This indicates frequent recruitment of TEs as part of novel genes. Finally, we also obtain evidence that a small fraction of primate orphan genes, around 5.5%, might have originated de novo from mammalian noncoding genomic regions.

Evolution of primate orphan proteins.

Genomes contain a large number of genes that do not have recognizable homologues in other species. These genes, found in only one or a few closely related species, are known as orphan genes. Their limited distribution implies that many of them are probably involved in lineage-specific adaptive processes. One important question that has remained elusive to date is how orphan genes originate. It has been proposed that they might have arisen by gene duplication followed by a period of very rapid sequence divergence, which would have erased any traces of similarity to other evolutionarily related genes. However, this explanation does not seem plausible for genes lacking homologues in very closely related species. In the present article, we review recent efforts to identify the mechanisms of formation of primate orphan genes. These studies reveal an unexpected important role of transposable elements in the formation of novel protein-coding genes in the genomes of primates.

New insights on Pseudoalteromonas haloplanktis TAC125 genome organization and benchmarks of genome assembly applications using next and third generation sequencing technologies.

Pseudoalteromonas haloplanktis TAC125 is among the most commonly studied bacteria adapted to cold environments. Aside from its ecological relevance, P. haloplanktis has a potential use for biotechnological applications. Due to its importance, we decided to take advantage of next generation sequencing (Illumina) and third generation sequencing (PacBio and Oxford Nanopore) technologies to resequence its genome. The availability of a reference genome, obtained using whole genome shotgun sequencing, allowed us to study and compare the results obtained by the different technologies and draw useful conclusions for future de novo genome assembly projects. We found that assembly polishing using Illumina reads is needed to achieve a consensus accuracy over 99.9% when using Oxford Nanopore sequencing, but not in PacBio sequencing. However, the dependency of consensus accuracy on coverage is lower in Oxford Nanopore than in PacBio, suggesting that a cost-effective solution might be the use of low coverage Oxford Nanopore sequencing together with Illumina reads. Despite the differences in consensus accuracy, all sequencing technologies revealed the presence of a large plasmid, pMEGA, which was undiscovered until now. Among the most interesting features of pMEGA is the presence of a putative error-prone polymerase regulated through the SOS response. Aside from the characterization of the newly discovered plasmid, we confirmed the sequence of the small plasmid pMtBL and uncovered the presence of a potential partitioning system. Crucially, this study shows that the combination of next and third generation sequencing technologies give us an unprecedented opportunity to characterize our bacterial model organisms at a very detailed level.

Integrative analysis of fitness and metabolic effects of plasmids in Pseudomonas aeruginosa PAO1.

Horizontal gene transfer (HGT) mediated by the spread of plasmids fuels evolution in prokaryotes. Although plasmids provide bacteria with new adaptive genes, they also produce physiological alterations that often translate into a reduction in bacterial fitness. The fitness costs associated with plasmids represent an important limit to plasmid maintenance in bacterial communities, but their molecular origins remain largely unknown. In this work, we combine phenomics, transcriptomics and metabolomics to study the fitness effects produced by a collection of diverse plasmids in the opportunistic pathogen Pseudomonas aeruginosa PAO1. Using this approach, we scan the physiological changes imposed by plasmids and test the generality of some main mechanisms that have been proposed to explain the cost of HGT, including increased biosynthetic burden, reduced translational efficiency, and impaired chromosomal replication. Our results suggest that the fitness effects of plasmids have a complex origin, since none of these mechanisms could individually provide a general explanation for the cost of plasmid carriage. Interestingly, our results also showed that plasmids alter the expression of a common set of metabolic genes in PAO1, and produce convergent changes in host cell metabolism. These surprising results suggest that there is a common metabolic response to plasmids in P. aeruginosa PAO1.

Mistranslation can enhance fitness through purging of deleterious mutations.

Phenotypic mutations are amino acid changes caused by mistranslation. How phenotypic mutations affect the adaptive evolution of new protein functions is unknown. Here we evolve the antibiotic resistance protein TEM-1 towards resistance on the antibiotic cefotaxime in an Escherichia coli strain with a high mistranslation rate. TEM-1 populations evolved in such strains endow host cells with a general growth advantage, not only on cefotaxime but also on several other antibiotics that ancestral TEM-1 had been unable to deactivate. High-throughput sequencing of TEM-1 populations shows that this advantage is associated with a lower incidence of weakly deleterious genotypic mutations. Our observations show that mistranslation is not just a source of noise that delays adaptive evolution. It could even facilitate adaptive evolution by exacerbating the effects of deleterious mutations and leading to their more efficient purging. The ubiquity of mistranslation and its effects render mistranslation an important factor in adaptive protein evolution.

Epistatic interactions between ancestral genotype and beneficial mutations shape evolvability in Pseudomonas aeruginosa.

The idea that interactions between mutations influence adaptation by driving populations to low and high fitness peaks on adaptive landscapes is deeply ingrained in evolutionary theory. Here, we investigate the impact of epistasis on evolvability by challenging populations of two Pseudomonas aeruginosa clones bearing different initial mutations (in rpoB conferring rifampicin resistance, and the type IV pili gene network) to adaptation to a medium containing l-serine as the sole carbon source. Despite being initially indistinguishable in fitness, populations founded by the two ancestral genotypes reached different fitness following 300 generations of evolution. Genome sequencing revealed that the difference could not be explained by acquiring mutations in different targets of selection; the majority of clones from both ancestors converged on one of the following two strategies: (1) acquiring mutations in either PA2449 (gcsR, an l-serine-metabolism RpoN enhancer binding protein) or (2) protease genes. Additionally, populations from both ancestors converged on loss-of-function mutations in the type IV pili gene network, either due to ancestral or acquired mutations. No compensatory or reversion mutations were observed in RNA polymerase (RNAP) genes, in spite of the large fitness costs typically associated with mutations in rpoB. Although current theory points to sign epistasis as the dominant constraint on evolvability, these results suggest that the role of magnitude epistasis in constraining evolvability may be underappreciated. The contribution of magnitude epistasis is likely to be greatest under the biologically relevant mutation supply rates that make back mutations probabilistically unlikely.

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.