Hypermutator emergence in experimental Escherichia coli populations is stress-type dependent.

Genotypes exhibiting an increased mutation rate, called hypermutators, can propagate in microbial populations because they can have an advantage due to the higher supply of beneficial mutations needed for adaptation. Although this is a frequently observed phenomenon in natural and laboratory populations, little is known about the influence of parameters such as the degree of maladaptation, stress intensity, and the genetic architecture for adaptation on the emergence of hypermutators. To address this knowledge gap, we measured the emergence of hypermutators over ~1,000 generations in experimental Escherichia coli populations exposed to different levels of osmotic or antibiotic stress. Our stress types were chosen based on the assumption that the genetic architecture for adaptation differs between them. Indeed, we show that the size of the genetic basis for adaptation is larger for osmotic stress compared to antibiotic stress. During our experiment, we observed an increased emergence of hypermutators in populations exposed to osmotic stress but not in those exposed to antibiotic stress, indicating that hypermutator emergence rates are stress type dependent. These results support our hypothesis that hypermutator emergence is linked to the size of the genetic basis for adaptation. In addition, we identified other parameters that covaried with stress type (stress level and IS transposition rates) that might have contributed to an increased hypermutator provision and selection. Our results provide a first comparison of hypermutator emergence rates under varying stress conditions and point towards complex interactions of multiple stress-related factors on the evolution of mutation rates.

Ecology, more than antibiotics consumption, is the major predictor for the global distribution of aminoglycoside-modifying enzymes.

Antibiotic consumption and its abuses have been historically and repeatedly pointed out as the major driver of antibiotic resistance emergence and propagation. However, several examples show that resistance may persist despite substantial reductions in antibiotic use, and that other factors are at stake. Here, we study the temporal, spatial, and ecological distribution patterns of aminoglycoside resistance, by screening more than 160,000 publicly available genomes for 27 clusters of genes encoding aminoglycoside-modifying enzymes (AME genes). We find that AME genes display a very ubiquitous pattern: about 25% of sequenced bacteria carry AME genes. These bacteria were sequenced from all the continents (except Antarctica) and terrestrial biomes, and belong to a wide number of phyla. By focusing on European countries between 1997 and 2018, we show that aminoglycoside consumption has little impact on the prevalence of AME-gene-carrying bacteria, whereas most variation in prevalence is observed among biomes. We further analyze the resemblance of resistome compositions across biomes: soil, wildlife, and human samples appear to be central to understand the exchanges of AME genes between different ecological contexts. Together, these results support the idea that interventional strategies based on reducing antibiotic use should be complemented by a stronger control of exchanges, especially between ecosystems.

Following the Dynamics of Structural Variants in Experimentally Evolved Populations.

Structural variants (SVs) (i.e., deletions, insertions, duplications, and inversions) are now known to play an important role in phenotypic variation, and consequently in processes such as disease determination or adaptation to a new environment. However, single-nucleotide variants receive much more attention than SVs, probably because they are easier to detect, and their phenotypic effects are easier to predict. The development of short- and long-read deep sequencing technologies have strongly improved the detection of SVs, but the quantification of their frequency from pooled sequencing (poolseq) data is still technically complex and expensive. Here, we present a rather simple and inexpensive method, which allows researchers to follow the dynamics of SV allele frequency. As an example of application, we follow the frequency of an insertion sequence (IS) insertion in experimental evolution populations of bacteria. This method is based on the design of triplets of primers around the structural variant borders, such that the amplicons produced by amplification of the wild-type (WT) and derived alleles differ in size by at least 5%, and that their amplification efficiency is similar. The quantity of each amplicon is then determined by parallel capillary electrophoresis and normalized to a calibration curve. This method can be easily extended to the quantification of the frequency of other structural variants (deletions, duplications, and inversions) and to pool-seq approaches of natural populations, including within-patient pathogen populations.

Evolutionary responses to codon usage of horizontally transferred genes in Pseudomonas aeruginosa: gene retention, amelioration and compensatory evolution.

Prokaryote genome evolution is characterized by the frequent gain of genes through horizontal gene transfer (HGT). For a gene, being horizontally transferred can represent a strong change in its genomic and physiological context. If the codon usage of a transferred gene deviates from that of the receiving organism, the fitness benefits it provides can be reduced due to a mismatch with the expression machinery. Consequently, transferred genes with a deviating codon usage can be selected against or elicit evolutionary responses that enhance their integration, such as gene amelioration and compensatory evolution. Within bacterial species, the extent and relative importance of these different mechanisms has never been considered altogether. In this study, a phylogeny-based method was used to investigate the occurrence of these different evolutionary responses in Pseudomonas aeruginosa. Selection on codon usage of genes acquired through HGT was observed over evolutionary time, with the overall codon usage converging towards that of the core genome. Gene amelioration, through the accumulation of synonymous mutations after HGT, did not seem to systematically affect transferred genes. This pattern therefore seemed to be mainly driven by selective retention of transferred genes with an initial codon usage similar to that of the core genes. Additionally, variation in the copy number of tRNA genes was often associated with the acquisition of genes for which the observed variation could enhance their expression. This provides evidence that compensatory evolution might be an important mechanism for the integration of horizontally transferred genes.

PlasForest: a homology-based random forest classifier for plasmid detection in genomic datasets.

BACKGROUND: Plasmids are mobile genetic elements that often carry accessory genes, and are vectors for horizontal transfer between bacterial genomes. Plasmid detection in large genomic datasets is crucial to analyze their spread and quantify their role in bacteria adaptation and particularly in antibiotic resistance propagation. Bioinformatics methods have been developed to detect plasmids. However, they suffer from low sensitivity (i.e., most plasmids remain undetected) or low precision (i.e., these methods identify chromosomes as plasmids), and are overall not adapted to identify plasmids in whole genomes that are not fully assembled (contigs and scaffolds). RESULTS: We developed PlasForest, a homology-based random forest classifier identifying bacterial plasmid sequences in partially assembled genomes. Without knowing the taxonomical origin of the samples, PlasForest identifies contigs as plasmids or chromosomes with a F1 score of 0.950. Notably, it can detect 77.4% of plasmid contigs below 1 kb with 2.8% of false positives and 99.9% of plasmid contigs over 50 kb with 2.2% of false positives. CONCLUSIONS: PlasForest outperforms other currently available tools on genomic datasets by being both sensitive and precise. The performance of PlasForest on metagenomic assemblies are currently well below those of other k-mer-based methods, and we discuss how homology-based approaches could improve plasmid detection in such datasets.

Read between the Lines: Diversity of Nontranslational Selection Pressures on Local Codon Usage.

Protein coding genes can contain specific motifs within their nucleotide sequence that function as a signal for various biological pathways. The presence of such sequence motifs within a gene can have beneficial or detrimental effects on the phenotype and fitness of an organism, and this can lead to the enrichment or avoidance of this sequence motif. The degeneracy of the genetic code allows for the existence of alternative synonymous sequences that exclude or include these motifs, while keeping the encoded amino acid sequence intact. This implies that locally, there can be a selective pressure for preferentially using a codon over its synonymous alternative in order to avoid or enrich a specific sequence motif. This selective pressure could-in addition to mutation, drift and selection for translation efficiency and accuracy-contribute to shape the codon usage bias. In this review, we discuss patterns of avoidance of (or enrichment for) the various biological signals contained in specific nucleotide sequence motifs: transcription and translation initiation and termination signals, mRNA maturation signals, and antiviral immune system targets. Experimental data on the phenotypic or fitness effects of synonymous mutations in these sequence motifs confirm that they can be targets of local selection pressures on codon usage. We also formulate the hypothesis that transposable elements could have a similar impact on codon usage through their preferred integration sequences. Overall, selection on codon usage appears to be a combination of a global selection pressure imposed by the translation machinery, and a patchwork of local selection pressures related to biological signals contained in specific sequence motifs.

Evolutionary Changes after Translational Challenges Imposed by Horizontal Gene Transfer.

Genes acquired by horizontal gene transfer (HGT) may provide the recipient organism with potentially new functions, but proper expression level and integration of the transferred genes in the novel environment are not granted. Notably, transferred genes can differ from the receiving genome in codon usage preferences, leading to impaired translation and reduced functionality. Here, we characterize the genomic and proteomic changes undergone during experimental evolution of Escherichia coli after HGT of three synonymous versions, presenting very different codon usage preference, of an antibiotic resistance gene. The experimental evolution was conducted with and without the corresponding antibiotic and the mutational patterns and proteomic profiles after 1,000 generations largely depend on the experimental growth conditions (e.g., mutations in antibiotic off-target genes), and on the synonymous gene version transferred (e.g., mutations in genes responsive to translational stress). The transfer of an exogenous gene extensively modifies the whole proteome, and these proteomic changes are different for the different version of the transferred gene. Additionally, we identified conspicuous changes in global regulators and in intermediate metabolism, confirmed the evolutionary ratchet generated by mutations in DNA repair genes and highlighted the plasticity of bacterial genomes accumulating large and occasionally transient duplications. Our results support a central role of HGT in fuelling evolution as a powerful mechanism promoting rapid, often dramatic genotypic and phenotypic changes. The profound reshaping of the pre-existing geno/phenotype allows the recipient bacteria to explore new ways of functioning, far beyond the mere acquisition of a novel function.

Assessing parallel gene histories in viral genomes.

BACKGROUND: The increasing abundance of sequence data has exacerbated a long known problem: gene trees and species trees for the same terminal taxa are often incongruent. Indeed, genes within a genome have not all followed the same evolutionary path due to events such as incomplete lineage sorting, horizontal gene transfer, gene duplication and deletion, or recombination. Considering conflicts between gene trees as an obstacle, numerous methods have been developed to deal with these incongruences and to reconstruct consensus evolutionary histories of species despite the heterogeneity in the history of their genes. However, inconsistencies can also be seen as a source of information about the specific evolutionary processes that have shaped genomes. RESULTS: The goal of the approach here proposed is to exploit this conflicting information: we have compiled eleven variables describing phylogenetic relationships and evolutionary pressures and submitted them to dimensionality reduction techniques to identify genes with similar evolutionary histories. To illustrate the applicability of the method, we have chosen two viral datasets, namely papillomaviruses and Turnip mosaic virus (TuMV) isolates, largely dissimilar in genome, evolutionary distance and biology. Our method pinpoints viral genes with common evolutionary patterns. In the case of papillomaviruses, gene clusters match well our knowledge on viral biology and life cycle, illustrating the potential of our approach. For the less known TuMV, our results trigger new hypotheses about viral evolution and gene interaction. CONCLUSIONS: The approach here presented allows turning phylogenetic inconsistencies into evolutionary information, detecting gene assemblies with similar histories, and could be a powerful tool for comparative pathogenomics.

Cancer, Warts, or Asymptomatic Infections: Clinical Presentation Matches Codon Usage Preferences in Human Papillomaviruses.

Viruses rely completely on the hosts' machinery for translation of viral transcripts. However, for most viruses infecting humans, codon usage preferences (CUPrefs) do not match those of the host. Human papillomaviruses (HPVs) are a showcase to tackle this paradox: they present a large genotypic diversity and a broad range of phenotypic presentations, from asymptomatic infections to productive lesions and cancer. By applying phylogenetic inference and dimensionality reduction methods, we demonstrate first that genes in HPVs are poorly adapted to the average human CUPrefs, the only exception being capsid genes in viruses causing productive lesions. Phylogenetic relationships between HPVs explained only a small proportion of CUPrefs variation. Instead, the most important explanatory factor for viral CUPrefs was infection phenotype, as orthologous genes in viruses with similar clinical presentation displayed similar CUPrefs. Moreover, viral genes with similar spatiotemporal expression patterns also showed similar CUPrefs. Our results suggest that CUPrefs in HPVs reflect either variations in the mutation bias or differential selection pressures depending on the clinical presentation and expression timing. We propose that poor viral CUPrefs may be central to a trade-off between strong viral gene expression and the potential for eliciting protective immune response.

Emerging viruses: why they are not jacks of all trades?

In order to limit the impact of the recent pandemics ignited by viral host jumps, it is necessary to better understand the ecological and evolutionary factors influencing the early steps of emergence and the interactions between them. Antagonistic pleiotropy, that is, the negative fitness effect in the primary host of mutations allowing the infection of and the multiplication in a new host, has long been thought to be the main limitation to the evolution of generalist viruses and thus to emergence. However, the accumulation of experimental examples contradicting the hypothesis of antagonistic pleiotropy has highlighted the importance of other factors such as the epistasis between mutations increasing the adaptation to a new host. Epistasis is pervasive in viruses, affects the shape of the adaptive landscape and consequently the accessibility of evolutionary pathways. Finally, recent studies have gone steps further in the complexity of viral fitness determinism and stressed the potential importance of the epistatic pleiotropy and of the impact of host living conditions.

Social interactions modulate the virulence of avian malaria infection.

There is an increasing understanding of the context-dependent nature of parasite virulence. Variation in parasite virulence can occur when infected individuals compete with conspecifics that vary in infection status; virulence may be higher when competing with uninfected competitors. In vertebrates with social hierarchies, we propose that these competition-mediated costs of infection may also vary with social status. Dominant individuals have greater competitive ability than competing subordinates, and consequently may pay a lower prevalence-mediated cost of infection. In this study we investigated whether costs of malarial infection were affected by the occurrence of the parasite in competitors and social status in domestic canaries (Serinus canaria). We predicted that infected subordinates competing with non-infected dominants would pay higher costs than infected subordinates competing with infected dominants. We also predicted that these occurrence-mediated costs of infection would be ameliorated in infected dominant birds. We found that social status and the occurrence of parasites in competitors significantly interacted to change haematocrit in infected birds. Namely, subordinate and dominant infected birds differed in haematocrit depending on the infection status of their competitors. However, in contrast to our prediction, dominants fared better with infected subordinates, whereas subordinates fared better with uninfected dominants. Moreover, we found additional effects of parasite occurrence on mortality in canaries. Ultimately, we provide evidence for costs of parasitism mediated by social rank and the occurrence of parasites in competitors in a vertebrate species. This has important implications for our understanding of the evolutionary processes that shape parasite virulence and group living.

Genotypic but not phenotypic historical contingency revealed by viral experimental evolution.

BACKGROUND: The importance of historical contingency in determining the potential of viral populations to evolve has been largely unappreciated. Identifying the constraints imposed by past adaptations is, however, of importance for understanding many questions in evolutionary biology, such as the evolution of host usage dynamics by multi-host viruses or the emergence of escape mutants that persist in the absence of antiviral treatments. To address this issue, we undertook an experimental approach in which sixty lineages of Tobacco etch potyvirus that differ in their past evolutionary history and degree of adaptation to Nicotiana tabacum were allowed to adapt to this host for 15 rounds of within host multiplication and transfer. We thereafter evaluated the degree of adaptation to the new host as well as to the original ones and characterized the consensus sequence of each lineage. RESULTS: We found that past evolutionary history did not determine the phenotypic outcome of this common host evolution phase, and that the signal of local adaptation to past hosts had largely disappeared. By contrast, evolutionary history left footprints at the genotypic level, since the majority of host-specific mutations present at the beginning of this experiment were retained in the end-point populations and may have affected which new mutations were consequently fixed. This resulted in further divergence between the sequences despite a shared selective environment. CONCLUSIONS: The present experiment reinforces the idea that the answer to the question "How important is historical contingency in evolution?" strongly depends on the level of integration of the traits studied. A strong historical contingency was found for TEV genotype, whereas a weak effect of on phenotypic evolution was revealed. In an applied context, our results imply that viruses are not easily trapped into suboptimal phenotypes and that (re)emergence is not evolutionarily constrained.

Multihost experimental evolution of a plant RNA virus reveals local adaptation and host-specific mutations.

For multihost pathogens, adaptation to multiple hosts has important implications for both applied and basic research. At the applied level, it is one of the main factors determining the probability and the severity of emerging disease outbreaks. At the basic level, it is thought to be a key mechanism for the maintenance of genetic diversity both in host and pathogen species. Using Tobacco etch potyvirus (TEV) and four natural hosts, we have designed an evolution experiment whose strength and novelty are the use of complex multicellular host organism as hosts and a high level of replication of different evolutionary histories and lineages. A pattern of local adaptation, characterized by a higher infectivity and virulence on host(s) encountered during the experimental evolution was found. Local adaptation only had a cost in terms of performance on other hosts in some cases. We could not verify the existence of a cost for generalists, as expected to arise from antagonistic pleiotropy and other genetic mechanisms generating a fitness trade-off between hosts. This observation confirms that this classical theoretical prediction lacks empirical support. We discuss the reasons for this discrepancy between theory and experiment in the light of our results. The analysis of full genome consensus sequences of the evolved lineages established that all mutations shared between lineages were host specific. A low degree of parallel evolution was observed, possibly reflecting the various adaptive pathways available for TEV in each host. Altogether, these results reveal a strong adaptive potential of TEV to new hosts without severe evolutionary constraints.

Sperm competition and mate harm unresponsive to male-limited selection in Drosophila: an evolving genetic architecture under domestication.

Earlier research by W.R. Rice showed that experimentally limiting gene expression to males in Drosophila melanogaster leads to the rapid evolution of higher fitness. Using a similar male-limited (ML) selection protocol, we confirmed that result and showed that eliminating intralocus sexual conflict results in a comprehensive remodeling of the sexually dimorphic phenotype. However, despite starting from laboratory-evolved descendants of the same founder population used in earlier work, we found no evidence for the increased performance in sperm competition or increased postmating harm to females previously demonstrated. We employed females with both ancestral population genotypes and those of the special "clone generator" females used in ML selection. Despite strong differences in sperm storage or usage patterns between these females, there was no detectable adaptation by males to the specific female stock used in the selection protocol. The lack of evolution of postcopulatory traits suggests either that requisite genetic variation was eliminated by long-term domestication of the base population, or that complex male-by-male-by-female interactions made these traits unavailable to selection. The different evolutionary outcomes produced by two very similar experiments done at different time points underscores the potential for cryptic adaptation in the laboratory to qualitatively affect inferences made using quantitative genetic methodologies.

Virus infection suppresses Nicotiana benthamiana adaptive phenotypic plasticity.

Competition and parasitism are two important selective forces that shape life-histories, migration rates and population dynamics. Recently, it has been shown in various pathosystems that parasites can modify intraspecific competition, thus generating an indirect cost of parasitism. Here, we investigated if this phenomenon was present in a plant-potyvirus system using two viruses of different virulence (Tobacco etch virus and Turnip mosaic virus). Moreover, we asked if parasitism interacted with the shade avoidance syndrome, the plant-specific phenotypic plasticity in response to intraspecific competition. Our results indicate that the modification of intraspecific competition by parasitism is not present in the Nicotiana benthamiana--potyvirus system and suggests that this phenomenon is not universal but depends on the peculiarities of each pathosystem. However, whereas the healthy N. benthamiana presented a clear shade avoidance syndrome, this phenotypic plasticity totally disappeared when the plants were infected with TEV and TuMV, very likely resulting in a fitness loss and being another form of indirect cost of parasitism. This result suggests that the suppression or the alteration of adaptive phenotypic plasticity might be a component of virulence that is often overlooked.

The evolutionary genetics of emerging plant RNA viruses.

Over the years, agriculture across the world has been compromised by a succession of devastating epidemics caused by new viruses that spilled over from reservoir species or by new variants of classic viruses that acquired new virulence factors or changed their epidemiological patterns. Viral emergence is usually associated with ecological change or with agronomical practices bringing together reservoirs and crop species. The complete picture is, however, much more complex, and results from an evolutionary process in which the main players are ecological factors, viruses' genetic plasticity, and host factors required for virus replication, all mixed with a good measure of stochasticity. The present review puts emergence of plant RNA viruses into the framework of evolutionary genetics, stressing that viral emergence begins with a stochastic process that involves the transmission of a preexisting viral strain into a new host species, followed by adaptation to the new host.

Male-limited evolution suggests no extant intralocus sexual conflict over the sexually dimorphic cuticular hydrocarbons of Drosophila melanogaster.

Sexually dimorphic traits are likely to have evolved through sexually antagonistic selection. However, recent empirical data suggest that intralocus sexual conflict often persists, even when traits have diverged between males and females. This implies that evolved dimorphism is often incomplete in resolving intralocus conflict, providing a mechanism for the maintenance of genetic variance in fitness-related traits. We used experimental evolution in Drosophila melanogaster to directly test for ongoing conflict over a suite of sexually dimorphic cuticular hydrocarbons (CHCs) that are likely targets of sex-specific selection. Using a set of experimental populations in which the transmission of genetic material had been restricted to males for 82 generations, we show that CHCs did not evolve, providing experimental evidence for the absence of current intralocus sexual conflict over these traits. The absence of ongoing conflict could indicate that CHCs have never been the target of sexually antagonistic selection, although this would require the existing dimorphism to have evolved via completely sexlinked mutations or as a result of former, but now absent, pleiotropic effects of the underlying loci on another trait under sexually antagonistic selection. An alternative interpretation, and which we believe to be more likely, is that the extensive CHC sexual dimorphism is the result of past intralocus sexual conflict that has been fully resolved, implying that these traits have evolved genetic independence between the sexes and that genetic variation in them is therefore maintained by alternative mechanisms. This latter interpretation is consistent with the known roles of CHCs in sexual communication in this species and with previous studies suggesting the genetic independence of CHCs between males and females. Nevertheless, direct estimates of sexually antagonistic selection will be important to fully resolve these alternatives.

Evolution in regulatory regions rapidly compensates the cost of nonoptimal codon usage.

The redundant genetic code contains synonymous codons, whose relative frequencies vary among species. Nonoptimal codon usage lowers gene translation efficiency, potentially leading to a fitness cost. This is particularly relevant for horizontal gene transfer, common among bacteria and a key player in antibiotic resistance propagation. By mimicking the horizontal transfer of an antibiotic resistance gene, we established that a nonoptimal codon usage renders Escherichia coli 10-20 times more sensitive to the antibiotic. After 350 generations of experimental evolution under antibiotic selection pressure, this cost was compensated through both in cis changes in the gene promoter and in trans changes in the host bacterial genome, without introducing mutations in the coding sequence of the resistance gene. Further, we have found experimental evidence for convergent molecular adaptive evolution. The high fitness cost of nonoptimal codon usage remains a minor obstacle to gene fixation upon horizontal transfer. Our results highlight the importance of rapid evolution of regulatory mechanisms in the adaptation to new environmental and genetic situations.

How does breeding system variation modulate sexual antagonism?

The study of sexually antagonistic (SA) traits remains largely limited to dioecious (separate sex), mobile animals. However, the occurrence of sexual conflict is restricted neither by breeding system (the mode of sexual reproduction, e.g. dioecy or hermaphroditism) nor by sessility. Here, we synthesize how variation in breeding system can affect the evolution and expression of intra- and inter-locus sexual conflicts in plants and animals. We predict that, in hermaphrodites, SA traits will (i) display lower levels of polymorphism; (ii) respond more quickly to selection; and (iii) involve unique forms of interlocus conflict over sex allocation, mating roles and selfing rates. Explicit modelling and empirical tests in a broader range of breeding systems are necessary to obtain a general understanding of the evolution of SA traits.

Reproductive behaviour evolves rapidly when intralocus sexual conflict is removed.

BACKGROUND: Intralocus sexual conflict can inhibit the evolution of each sex towards its own fitness optimum. In a previous study, we confirmed this prediction through the experimental removal of female selection pressures in Drosophila melanogaster, achieved by limiting the expression of all major chromosomes to males. Compared to the control populations (C(1-4)) where the genomes are exposed to selection in both sexes, the populations with male-limited genomes (ML(1-4)) showed rapid increases in male fitness, whereas the fitness of females expressing ML-evolved chromosomes decreased. METHODOLOGY/PRINCIPAL FINDINGS: Here we examine the behavioural phenotype underlying this sexual antagonism. We show that males expressing the ML genomes have a reduced courtship level but acquire the same number of matings. On the other hand, our data suggest that females expressing the ML genomes had reduced attractiveness, stimulating a lower rate of courtship from males. Moreover, females expressing ML genomes tend to display reduced yeast-feeding behaviour, which is probably linked to the reduction of their fecundity. CONCLUSION/SIGNIFICANCE: These results suggest that reproductive behaviour is shaped by opposing selection on males and females, and that loci influencing attractiveness and foraging were polymorphic for alleles with sexually antagonistic expression patterns prior to ML selection. Hence, intralocus sexual conflict appears to play a role in the evolution of a wide range of fitness-related traits and may be a powerful mechanism for the maintenance of genetic variation in fitness.