Impact of mutation rate on the adaptation of gut bacteria

To study the role of mutator bacteria in the evolution of bacterial populations, we followed the impact of the mutation rate of Escherichia coli strains in the colonisation of the gut of axenic mice and the evolution of the mutation rate of bacterial populations living in the gut. We show that mutator bacteria have an advantage during the colonization. This adaptive advantage comes from their ability to generate adaptive mutations faster than wild type strains, mutations that allow their maintenance in the ecosystem. However, while mutator bacteria are becoming specialised to the environment they are living in, they accumulate mutations that may be deleterious or lethal in secondary environments. By following the evolution of the mutation rate of bacterial populations living in the gut of mice receiving antibiotics, we show that this therapy selects not only for antibiotic resistant mutants but also for mutator alleles that enhance mutation rates and are responsible for the appearance of the resistance. The costs of a high mutation rate, due to the accumulation of mutations, is seen in environments where changes are recurrent. In an ever-changing situation where every change is new, mutator bacteria might help the evolution of bacterial populations.     

A trade-off between oxidative stress resistance and DNA repair plays a role in the evolution of elevated mutation rates in bacteria

The dominant paradigm for the evolution of mutator alleles in bacterial populations is that they spread by indirect selection for linked beneficial mutations when bacteria are poorly adapted. In this paper, we challenge the ubiquity of this paradigm by demonstrating that a clinically important stressor, hydrogen peroxide, generates direct selection for an elevated mutation rate in the pathogenic bacterium Pseudomonas aeruginosa as a consequence of a trade-off between the fidelity of DNA repair and hydrogen peroxide resistance. We demonstrate that the biochemical mechanism underlying this trade-off in the case of mutS is the elevated secretion of catalase by the mutator strain. Our results provide, to our knowledge, the first experimental evidence that direct selection can favour mutator alleles in bacterial populations, and pave the way for future studies to understand how mutation and DNA repair are linked to stress responses and how this affects the evolution of bacterial mutation rates.     

Enrichment and elimination of mutY mutators in Escherichia coli populations

The kinetics of mutator sweeps was followed in two independent populations of Escherichia coli grown for up to 350 generations in glucose-limited continuous culture. A rapid elevation of mutation rates was observed in both populations within 120-150 generations, as was apparent from major increases in the proportion of the populations with unselected mutations in fhuA. The increase in mutation rates was due to sweeps by mutY mutators. In both cultures, the enrichment of mutators resulted from hitchhiking with identified beneficial mutations increasing fitness under glucose limitation; mutY hitchhiked with mgl mutations in one culture and ptsG in the other. In both cases, mutators were enriched to constitute close to 100% of the population before a periodic selection event reduced the frequency of unselected mutations and mutators in the cultures. The high proportion of mutators persisted for 150 generations in one population but began to be eliminated within 50 generations in the other. The persistence of mutator, as well as experimental data showing that mutY bacteria were as fit as near-isogenic mutY(+) bacteria in competition experiments, suggest that mutator load by deleterious mutations did not explain the rapidly diminishing proportion of mutators in the populations. The nonmutators sweeping out mutators were also unlikely to have arisen by reversion or antimutator mutations; the mutY mutations were major deletions in each case and the bacteria sweeping out mutators contained intact mutY. By following mgl allele frequencies in one population, we discovered that mutators were outcompeted by bacteria that had rare mgl mutations previously as well as additional beneficial mutation(s). The pattern of appearance of mutY, but not its elimination, conforms to current models of mutator sweeps in bacterial populations. A mutator with a narrow mutational spectrum like mutY may be lost if the requirement for beneficial mutations is for changes other than GC --> TA transversions. Alternatively, epistatic interactions between mutator mutation and beneficial mutations need to be postulated to explain mutator elimination.     

Fitness evolution and the rise of mutator alleles in experimental Escherichia coli populations

We studied the evolution of high mutation rates and the evolution of fitness in three experimental populations of Escherichia coli adapting to a glucose-limited environment. We identified the mutations responsible for the high mutation rates and show that their rate of substitution in all three populations was too rapid to be accounted for simply by genetic drift. In two of the populations, large gains in fitness relative to the ancestor occurred as the mutator alleles rose to fixation, strongly supporting the conclusion that mutator alleles fixed by hitchhiking with beneficial mutations at other loci. In one population, no significant gain in fitness relative to the ancestor occurred in the population as a whole while the mutator allele rose to fixation, but a substantial and significant gain in fitness occurred in the mutator subpopulation as the mutator neared fixation. The spread of the mutator allele from rarity to fixation took >1000 generations in each population. We show that simultaneous adaptive gains in both the mutator and wild-type subpopulations (clonal interference) retarded the mutator fixation in at least one of the populations. We found little evidence that the evolution of high mutation rates accelerated adaptation in these populations.     

Oscillations in continuous culture populations of Streptococcus pneumoniae: population dynamics and the evolution of clonal suicide

Agents that kill or induce suicide in the organisms that produce them or other individuals of the same genotype are intriguing puzzles for ecologists and evolutionary biologists. When those organisms are pathogenic bacteria, these suicidal toxins have the added appeal as candidates for the development of narrow spectrum antibiotics to kill the pathogens that produce them. We show that when clinical as well as laboratory strains of Streptococcus pneumoniae are maintained in continuous culture (chemostats), their densities oscillate by as much as five orders of magnitude with an apparently constant period. This dynamic, which is unanticipated for single clones of bacteria in chemostats, can be attributed to population-wide die-offs and recoveries. Using a combination of mathematical models and experiments with S. pneumoniae, we present evidence that these die-offs can be attributed to the autocatalytic production of a toxin that lyses or induces autolysis in members of the clone that produces it. This toxin, which our evidence indicates is a protein, appears to be novel; S. pneumoniae genetic constructs knocked out for lytA and other genes coding for known candidates for this agent oscillate in chemostat culture. Since this toxin lyses different strains of S. pneumoniae as well as other closely related species of Streptococcus, we propose that its ecological role is as an allelopathic agent. Using a mathematical model, we explore the conditions under which toxins that kill members of the same clone that produces them can prevent established populations from invasion by different strains of the same or other species. We postulate that the production of the toxin observed here as well as other bacteria-produced toxins that kill members of the same genotype, ‘clonal suicide’, evolved and are maintained to prevent colonization of established populations by different strains of the same and closely related species.     

The evolution of bacterial mutation rates under simultaneous selection by interspecific and social parasitism

Many bacterial populations harbour substantial numbers of hypermutable bacteria, in spite of hypermutation being associated with deleterious mutations. One reason for the persistence of hypermutators is the provision of novel mutations, enabling rapid adaptation to continually changing environments, for example coevolving virulent parasites. However, hypermutation also increases the rate at which intraspecific parasites (social cheats) are generated. Interspecific and intraspecific parasitism are therefore likely to impose conflicting selection pressure on mutation rate. Here, we combine theory and experiments to investigate how simultaneous selection from inter- and intraspecific parasitism affects the evolution of bacterial mutation rates in the plant-colonizing bacterium Pseudomonas fluorescens. Both our theoretical and experimental results suggest that phage presence increases and selection for public goods cooperation (the production of iron-scavenging siderophores) decreases selection for mutator bacteria. Moreover, phages imposed a much greater growth cost than social cheating, and when both selection pressures were imposed simultaneously, selection for cooperation did not affect mutation rate evolution. Given the ubiquity of infectious phages in the natural environment and clinical infections, our results suggest that phages are likely to be more important than social interactions in determining mutation rate evolution.     

Low Selection Pressure Aids the Evolution of Cooperative Ribozyme Mutations in Cells

Understanding the evolution of functional RNA molecules is important for our molecular understanding of biology. Here we tested experimentally how two evolutionary parameters, selection pressure and recombination, influenced the evolution of an evolving RNA population. This was done using four parallel evolution experiments that employed low or gradually increasing selection pressure, and recombination events either at the end or dispersed throughout the evolution. As model system, a trans-splicing group I intron ribozyme was evolved in Escherichia coli cells over 12 rounds of selection and amplification, including mutagenesis and recombination. The low selection pressure resulted in higher efficiency of the evolved ribozyme populations, whereas differences in recombination did not have a strong effect. Five mutations were responsible for the highest efficiency. The first mutation swept quickly through all four evolving populations, whereas the remaining four mutations accumulated later and more efficiently under low selection pressure. To determine why low selection pressure aided this evolution, all evolutionary intermediates between the wild type and the 5-mutation variant were constructed, and their activities at three different selection pressures were determined. The resulting fitness profiles showed a high cooperativity among the four late mutations, which can explain why high selection pressure led to inefficient evolution. These results show experimentally how low selection pressure can benefit the evolution of cooperative mutations in functional RNAs.     

Cloning and expression of the Escherichia coli glgC gene from a mutant containing an ADPglucose pyrophosphorylase with altered allosteric properties.

A mutant strain of Escherichia coli K-12, designated 618, accumulates glycogen at a faster rate than wild-type strain 356. The mutation affects the ADPglucose pyrophosphorylase regulatory properties (N. Creuzat-Sigal, M. Latil-Damotte, J. Cattaneo, and J. Puig, p. 647-680, in R. Piras and H. G. Pontis, ed., Biochemistry of the Glycocide Linkage, 1972). The enzyme is less dependent on the activator, fructose 1,6 bis-phosphate for activity and is less sensitive to inhibition by the inhibitor, 5'-AMP. The structural gene, glgC, for this allosteric mutant enzyme was cloned into the bacterial plasmid pBR322 by inserting the chromosomal DNA at the PstI site. The glycogen biosynthetic genes were selected by cotransformation of the neighboring asd gene into an E. coli mutant also defective in branching enzyme (glgB) activity. Two recombinant plasmids, pEBL1 and pEBL3, that had PstI chromosomal DNA inserts containing glgC and glgB were isolated. Branching enzyme and ADPglucose pyrophosphorylase activities were increased 240- and 40-fold, respectively, in the asd glgB mutant, E. coli K-12 6281. The E. coli K-12 618 mutant glgC gene product was characterized after transformation of an E. coli B ADPglucose pyrophosphorylase mutant with the recombinant plasmid pEBL3. The kinetic properties of the cloned ADPglucose pyrophosphorylase were similar to those of the E. coli K-12 618 enzyme. The inserted DNA in pEBL1 was arranged in opposite orientation to that in pEBL3.     

Spontaneously arising mutL mutators in evolving Escherichia coli populations are the result of changes in repeat length

Over the course of thousands of generations of growth in a glucose-limited environment, 3 of 12 experimental populations of Escherichia coli spontaneously and independently evolved greatly increased mutation rates. In two of the populations, the mutations responsible for this increased mutation rate lie in the same region of the mismatch repair gene mutL. In this region, a 6-bp repeat is present in three copies in the gene of the wild-type ancestor of the experimental populations but is present in four copies in one of the experimental populations and two copies in the other. These in-frame mutations either add or delete the amino acid sequence LA in the MutL protein. We determined that the replacement of the wild-type sequence with either of these mutations was sufficient to increase the mutation rate of the wild-type strain to a level comparable to that of the mutator strains. Complementation of strains bearing the mutator mutations with wild-type copies of either mutL or the mismatch repair gene uvrD rescued the wild-type mutation rate. The position of the mutator mutations-in the region of MutL known as the ATP lid-suggests a possible deficiency in MutL's ATPase activity as the cause of the mutator phenotype. The similarity of the two mutator mutations (despite the independent evolutionary histories of the populations that gave rise to them) leads to a discussion of the potential adaptive role of DNA repeats.     

Altered regulation of the OmpF porin by Fis in Escherichia coli during an evolution experiment and between B and K-12 strains

The phenotypic plasticity of global regulatory networks provides bacteria with rapid acclimation to a wide range of environmental conditions, while genetic changes in those networks provide additional flexibility as bacteria evolve across long time scales. We previously identified mutations in the global regulator-encoding gene fis that enhanced organismal fitness during a long-term evolution experiment with Escherichia coli. To gain insight into the effects of these mutations, we produced two-dimensional protein gels with strains carrying different fis alleles, including a beneficial evolved allele and one with an in-frame deletion. We found that Fis controls the expression of the major porin-encoding gene ompF in the E. coli B-derived ancestral strain used in the evolution experiment, a relationship that has not been described before. We further showed that this regulatory connection evolved over two different time scales, perhaps explaining why it was not observed before. On the longer time scale, we showed that this regulation of ompF by Fis is absent from the more widely studied K-12 strain and thus is specific to the B strain. On a shorter time scale, this regulatory linkage was lost during 20,000 generations of experimental evolution of the B strain. Finally, we mapped the Fis binding sites in the ompF regulatory region, and we present a hypothetical model of ompF expression that includes its other known regulators.     

The evolution of low mutation rates in experimental mutator populations of Saccharomyces cerevisiae

Mutation is the source of both beneficial adaptive variation and deleterious genetic load, fueling the opposing selective forces than shape mutation rate evolution. This dichotomy is well illustrated by the evolution of the mutator phenotype, a genome-wide 10- to 100-fold increase in mutation rate. This phenotype has often been observed in clonally expanding populations exposed to novel or frequently changing conditions. Although studies of both experimental and natural populations have shed light on the evolutionary forces that lead to the spread of the mutator allele through a population, significant gaps in our understanding of mutator evolution remain. Here we use an experimental evolution approach to investigate the conditions required for the evolution of a reduction in mutation rate and the mechanisms by which populations tolerate the accumulation of deleterious mutations. We find that after ～6,700 generations, four out of eight experimental mutator lines had evolved a decreased mutation rate. We provide evidence that the accumulation of deleterious mutations leads to selection for reduced mutation rate clones in populations of mutators. Finally, we test the long-term consequences of the mutator phenotype, finding that mutator lines follow different evolutionary trajectories, some of which lead to drug resistance.     

The Molecular and Genetic Basis of Repeatable Coevolution between Escherichia coli and Bacteriophage T3 in a Laboratory Microcosm

The objective of this study was to determine the genomic changes that underlie coevolution between Escherichia coli B and bacteriophage T3 when grown together in a laboratory microcosm. We also sought to evaluate the repeatability of their evolution by studying replicate coevolution experiments inoculated with the same ancestral strains. We performed the coevolution experiments by growing Escherichia coli B and the lytic bacteriophage T3 in seven parallel continuous culture devices (chemostats) for 30 days. In each of the chemostats, we observed three rounds of coevolution. First, bacteria evolved resistance to infection by the ancestral phage. Then, a new phage type evolved that was capable of infecting the resistant bacteria as well as the sensitive bacterial ancestor. Finally, we observed second-order resistant bacteria evolve that were resistant to infection by both phage types. To identify the genetic changes underlying coevolution, we isolated first- and second-order resistant bacteria as well as a host-range mutant phage from each chemostat and sequenced their genomes. We found that first-order resistant bacteria consistently evolved resistance to phage via mutations in the gene, waaG, which codes for a glucosyltransferase required for assembly of the bacterial lipopolysaccharide (LPS). Phage also showed repeatable evolution, with each chemostat producing host-range mutant phage with mutations in the phage tail fiber gene T3p48 which binds to the bacterial LPS during adsorption. Two second-order resistant bacteria evolved via mutations in different genes involved in the phage interaction. Although a wide range of mutations occurred in the bacterial waaG gene, mutations in the phage tail fiber were restricted to a single codon, and several phage showed convergent evolution at the nucleotide level. These results are consistent with previous studies in other systems that have documented repeatable evolution in bacteria at the level of pathways or genes and repeatable evolution in viruses at the nucleotide level. Our data are also consistent with the expectation that adaptation via loss-of-function mutations is less constrained than adaptation via gain-of-function mutations.     

Rates of DNA sequence evolution in experimental populations of Escherichia coli during 20,000 generations

We examined rates of DNA sequence evolution in 12 populations of Escherichia coli propagated in a glucose minimal medium for 20,000 generations. Previous work saw mutations mediated by mobile elements in these populations, but the extent of other genomic changes was not investigated. Four of the populations evolved defects in DNA repair and became mutators. Some 500 bp was sequenced in each of 36 genes for 50 clones, including 2 ancestral variants, 2 clones from each population at generation 10,000, and 2 from each at generation 20,000. Ten mutations were found in total, all point mutations including mostly synonymous substitutions and nonsynonymous polymorphisms; all 10 were found in mutator populations. We compared the observed sequence evolution to predictions based on different scenarios. The number of synonymous substitutions is lower than predicted from measured mutation rates in E. coli, but the number is higher than rates based on comparing E. coli and Salmonella genomes. Extrapolating to the entire genome, these data predict about 250 synonymous substitutions on average per mutator population, but only about 3 synonymous substitutions per nonmutator population, during 20,000 generations. These data illustrate the challenge of finding sequence variation among bacterial isolates that share such a recent ancestor. However, this limited variation also provides a useful baseline for research aimed at finding the beneficial substitutions in these populations.     

DNA polymerase V-dependent mutator activity in an SOS-induced Escherichia coli strain with a temperature-sensitive DNA polymerase III.

The temperature-sensitive DNA polymerase III (Pol III) encoded by the dnaE486 allele confers a spontaneous mutator activity in SOS-induced bacteria that is largely dependent upon DNA polymerase V (Pol V), encoded by umuD, C. This mutator activity is influenced by the defective proof-reading sub-unit of Pol III encoded by the dnaQ905 (mutD5) allele arguing that Pol V is most likely fixing mutations arising from mismatched primer termini produced by Pol III(486). The size of the dnaQ effect is, however, modest leaving open the possibility that Pol V may be responsible for some of the mutator effect by engaging in bursts of processive activity.     

Mismatch repair genes of Streptococcus pneumoniae: HexA confers a mutator phenotype in Escherichia coli by negative complementation.

DNA repair systems able to correct base pair mismatches within newly replicated DNA or within heteroduplex molecules produced during recombination are widespread among living organisms. Evidence that such generalized mismatch repair systems evolved from a common ancestor is particularly strong for two of them, the Hex system of the gram-positive Streptococcus pneumoniae and the Mut system of the gram-negative Escherichia coli and Salmonella typhimurium. The homology existing between HexA and MutS and between HexB and MutL prompted us to investigate the effect of expressing hex genes in E. coli. Complementation of mutS or mutL mutations, which confer a mutator phenotype, was assayed by introducing on a multicopy plasmid the hexA and hexB genes, under the control of an inducible promoter, either individually or together in E. coli strains. No decrease in mutation rate was conferred by either hexA or hexB gene expression. However, a negative complementation effect was observed in wild-type E. coli cells: expression of hexA resulted in a typical Mut- mutator phenotype. hexB gene expression did not increase the mutation rate either individually or in conjunction with hexA. Since expression of hexA did not affect the mutation rate in mutS mutant cells and the hexA-induced mutator effect was recA independent, it is concluded that this effect results from inhibition of the Mut system. We suggest that HexA, like its homolog MutS, binds to mismatches resulting from replication errors, but in doing so it protects them from repair by the Mut system. In agreement with this hypothesis, an increase in mutS gene copy number abolished the hexA-induced mutator phenotype. HexA protein could prevent repair either by being unable to interact with Mut proteins or by producing nonfunctional repair complexes.     

Hypermutability impedes cooperation in pathogenic bacteria

When the supply of beneficial mutations limits adaptation, bacterial mutator alleles can reach high frequencies by hitchhiking with advantageous mutations. However, when populations are well adapted to their environments, the increased rate of deleterious mutations makes hypermutability selectively disadvantageous. Here, we consider a further cost of hypermutability: its potential to break down cooperation (group-beneficial behavior that is costly to the individual). This probably occurs for three reasons. First, an increased rate at which 'cheating' genotypes are generated; second, an increased probability of producing efficient cheats; and third, a decrease in relatedness (not addressed in the present study). We used Pseudomonas aeruginosa's production of extracellular iron-scavenging molecules, siderophores, to determine if cheating evolved more readily in mutator populations. Siderophore production is costly to individual bacteria but benefits all nearby cells. Siderophore-deficient cheats therefore have a selective advantage within populations. We observed the de novo evolution and subsequent increase in frequency of siderophore cheats within both wild-type and mutator populations for 200 generations. Cheats appeared and increased in frequency more rapidly in mutator populations. The presence of cheats was costly to the group, as shown by a negative correlation between cheat frequency and population density.     

Activation of the cryptic PhnE permease promotes rapid adaptive evolution in a population of Escherichia coli K-12 starved for phosphate

Escherichia coli K-12 suffers acetic acid stress during prolonged incubation in glucose minimal medium containing a limiting concentration of inorganic phosphate (0.1 mM P(i)), which decreases the number of viable cells from 6 × 10(8) to ≤10 CFU/ml between days 6 and 14 of incubation. Here we show that following two serial transfers into P(i)-limiting medium, evolved mutants survived prolonged incubation (≈10(7) CFU/ml on day 14 of incubation). The evolved strains that overtook the populations were generally PhnE(+), whereas the ancestral K-12 strain carries an inactive phnE allele, which prevents the transport of phosphonates. The switching in phnE occurred with a high frequency as a result of the deletion of an 8-bp repeated sequence. In a mixed culture starved for P(i) that contained the K-12 ancestral strain in majority, evolved strains grew through PhnE-dependent scavenging of probably organic phosphate esters (not phosphonates or P(i)) released by E. coli K-12 between days 1 and 3, before acetic acid excreted by E. coli K-12 reached toxic levels. The growth yield of phnE(+) strains in mixed culture was dramatically enhanced by mutations that affect glucose metabolism, such as an rpoS mutation inactivating the alternative sigma factor RpoS. The long-term viability of evolved populations was generally higher when the ancestral strain carried an inactive rather than an active phnE allele, which indicates that cross-feeding of phosphorylated products as a result of the phnE polymorphism may be essential for the spread of mutants which eventually help populations to survive under P(i) starvation conditions.     

Accumulation of trehalose by Escherichia coli K-12 at high osmotic pressure depends on the presence of amber suppressors

When grown at high osmotic pressure, some strains of Escherichia coli K-12 synthesized substantial levels of free sugar and accumulated proline if it was present in the growth medium. The sugar was identified as trehalose by chemical reactivity, gas-liquid chromatography, and nuclear magnetic resonance spectroscopy. Strains of E. coli K-12 could be divided into two major classes with respect to osmoregulation. Those of class A showed a large increase in trehalose levels with increasing medium osmolarity and also accumulated proline from the medium, whereas those in class B showed no accumulation of trehalose or proline. Most class A strains carried suppressor mutations which arose during their derivation from the wild type, whereas the osmodefective strains of class B were suppressor free. When amber suppressor mutations at the supD, supE, or supF loci were introduced into such sup0 osmodefective strains, they became osmotolerant and gained the ability to accumulate trehalose in response to elevated medium osmolarity. It appears that the original K-12 strain of E. coli carries an amber mutation in a gene affecting osmoregulation. Mutants lacking ADP-glucose synthetase (glgC) accumulated trehalose normally, whereas mutants lacking UDP-glucose synthetase (galU) did not make trehalose and grew poorly in medium of high osmolarity. Trehalose synthesis was repressed by exogenous glycine betaine but not by proline.     

Grazing protozoa and the evolution of the Escherichia coli O157:H7 Shiga toxin-encoding prophage

Humans play little role in the epidemiology of Escherichia coli O157:H7, a commensal bacterium of cattle. Why then does E. coli O157:H7 code for virulence determinants, like the Shiga toxins (Stxs), responsible for the morbidity and mortality of colonized humans? One possibility is that the virulence of these bacteria to humans is coincidental and these virulence factors evolved for and are maintained for other roles they play in the ecology of these bacteria. Here, we test the hypothesis that the carriage of the Stx-encoding prophage of E. coli O157:H7 increases the rate of survival of E. coli in the presence of grazing protozoa, Tetrahymena pyriformis. In the presence but not the absence of Tetrahymena, the carriage of the Stx-encoding prophage considerably augments the fitness of E. coli K-12 as well as clinical isolates of E. coli O157 by increasing the rate of survival of the bacteria in the food vacuoles of these ciliates. Grazing protozoa in the environment or natural host are likely to play a significant role in the ecology and maintenance of the Stx-encoding prophage of E. coli O157:H7 and may well contribute to the evolution of the virulence of these bacteria to colonize humans.     

Sexual ornamentation reflects antibacterial activity of ejaculates in mallards

Bacteria present in ejaculates can impair sperm function and reduce male reproductive success. Thus, selection should favour the evolution of antimicrobial defences to limit the detrimental effects of sperm-associated bacteria. Additionally, current hypotheses suggest that ornamental traits may signal information about the infection status of an individual or the ability of an individual to resist bacterial-induced sperm damage. However, despite the evolutionary implications of ejaculate antimicrobials, and the putative importance of pathogens for the evolution of male ornamentation, tests of these hypotheses are lacking. We examined the antibacterial activity of semen from mallard ducks (Anas platyrhynchos) and tested whether the bactericidal capacity of semen was associated with bill coloration, a sexually selected trait. We show that mallard semen exhibits significant antibacterial activity, as measured by the in vitro capacity to kill Escherichia coli and Staphylococcus aureus. Furthermore, we demonstrate that males with more colourful bills have semen with superior bacterial-killing ability. These results suggest that females could use male phenotypic traits to avoid sexually transmitted pathogens and acquire partners whose sperm suffer less bacteria-induced damage.