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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Prevalence-dependent costs of parasite virulence.

Costs of parasitism are commonly measured by comparing the performance of infected groups of individuals to that of uninfected control groups. This measure potentially underestimates the cost of parasitism because it ignores indirect costs, which may result from the modification of the competitiveness of the hosts by the parasite. In this context, we used the host-parasite system consisting of the yellow fever mosquito Aedes aegypti and the microsporidian parasite Vavraia culicis to address this question: Do infected individuals exert a more or less intense intraspecific competition than uninfected individuals? Our experimental results show that, indeed, infected hosts incur a direct cost of parasitism: It takes them longer to become adults than uninfected individuals. They also incur an indirect cost, however, which is actually larger than the direct cost: When grown in competition with uninfected individuals they develop even slower. The consequence of this modification of competitiveness is that, in our system, the cost of parasitism is underestimated by the traditional measure. Moreover, because the indirect cost depends on the frequency of interactions between infected and uninfected individuals, our results suggest that the real cost of parasitism, i.e., virulence, is negatively correlated with the prevalence of the parasite. This link between prevalence and virulence may have dynamical consequences, such as reducing the invasion threshold of the parasite, and evolutionary consequences, such as creating a selection pressure maintaining the host's constitutive resistance to the parasite.

Sexual conflict and environmental change: trade-offs within and between the sexes during the evolution of desiccation resistance.

Intralocus sexual conflict occurs when males and females experience sex-specific selection on a shared genome. With several notable exceptions, intralocus sexual conflict has been investigated in constant environments to which the study organisms have had an opportunity to adapt. However, a change in the environment can result in differential or even opposing selection pressures on males and females, creating sexual conflict. We used experimental evolution to explore the interaction between intralocus sexual conflict, sexual dimorphism and environmental variation in Drosophila melanogaster. Six populations were selected for adult desiccation resistance (D), with six matched control populations maintained in parallel (C). After 46 generations, the D populations had increased in survival time under arid conditions by 68% and in body weight by 20% compared to the C populations. The increase in size was the result of both extended development and faster growth rate of D juveniles. Adaptation to the stress came at a cost in terms of preadult viability and female fecundity. Because males are innately less tolerant of desiccation stress, very few D males survived desiccation-selection; while potentially a windfall for survivors, these conditions mean that most males' fitness was determined posthumously. We conjectured that selection for early maturation and mating in males was in conflict with selection for survival and later reproduction in females. Consistent with this prediction, the sexes showed different patterns of age-specific desiccation resistance and resource acquisition, and there was a trend towards increasingly female-biased sexual size dimorphism. However, levels of desiccation resistance were unaffected, with D males and females increasing in parallel. Either there is a strong positive genetic correlation between the sexes that limits independent evolution of desiccation resistance, or fitness pay-offs from the strategy of riding out the stress bout are great enough to sustain concordant selection on the two sexes. We discuss the forces that mould fitness in males under a regimen where trade-offs between survival and reproduction may be considerable.

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.

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.

Parasitism increases and decreases the costs of insecticide resistance in mosquitoes.

Adaptations conferring resistance to xenobiotics (antibiotics, insecticides, herbicides, etc.) are often costly to the organism's fitness in the absence of the selecting agent. In such conditions, and unless other mutations compensate for the costs of resistance, sensitive individuals are expected to out-reproduce resistant individuals and drive resistance alleles to a low frequency, with the rate and magnitude of this decline being proportional to the costs of resistance. However, this evolutionary dynamic is open to modification by other sources of selection acting on the relative fitness of susceptible and resistant individuals. Here we show parasitism not only as a source of selection capable of modifying the costs of organophosphate insecticide resistance in mosquitoes, but also that qualitatively different interactions (increasing or decreasing the relative fitness of resistant individuals) occurred depending on the particular form of resistance involved. As estimates of the parasite's fitness also varied according to its host's form of resistance, our data illustrate the potential for epidemiological feedbacks to influence the strength and direction of selection acting on resistance mutations in untreated environments.

Virulence reaction norms across a food gradient.

Host-parasite interactions involve competition for nutritional resources between hosts and the parasites growing within them. Consuming part of a host's resources is one cause of a parasite's virulence, i.e. part of the fitness cost imposed on the host by the parasite. The influence of a host's nutritional conditions on the virulence of a parasite was experimentally tested using the mosquito Aedes aegypti and the microsporidian parasite Vavraia culicis. A condition-dependent expression of virulence was found and a positive relation between virulence and transmissibility was established. Spore production was positively influenced by host food availability, indicating that the parasite's within-host growth is limited by host condition. We also investigated how the fitness of each partner varied across the nutritional gradient and demonstrated that the sign of the correlation between host fitness and parasite fitness depended on the amount of nutritional resources available to the host.

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.