How and at what pace bacteria evolve when colonizing healthy hosts remains unclear. Here, by monitoring evolution for more than six thousand generations in the mouse gut, we show that the successful colonization of an invader Escherichia coli depends on the diversity of the existing microbiota and the presence of a closely related strain. Following colonization, two modes of evolution were observed: one in which diversifying selection leads to long-term coexistence of ecotypes and a second in which directional selection propels selective sweeps. These modes can be quantitatively distinguished by the statistics of mutation trajectories. In our experiments, diversifying selection was marked by the emergence of metabolic mutations, and directional selection by acquisition of prophages, which bring their own benefits and costs. In both modes, we observed parallel evolution, with mutation accumulation rates comparable to those typically observed in vitro on similar time scales. Our results show how rapid ecotype formation and phage domestication can be in the mammalian gut.

DNA mismatch repair (MMR), an evolutionarily conserved repair pathway shared by prokaryotic and eukaryotic species alike, influences molecular evolution by detecting and correcting mismatches, thereby protecting genetic fidelity, reducing the mutational load, and preventing lethality. Herein we conduct the first genome-wide evaluation of the alterations to the mutation rate and spectrum under impaired activity of the MutS alpha homolog, msh-2, in Caenorhabditis elegans male-female fog-2(lf) lines. We performed mutation accumulation (MA) under RNAi-induced knockdown of msh-2 for up to 50 generations, followed by next-generation sequencing of 19 MA lines and the ancestral control. msh-2 impairment in the male-female background substantially increased the frequency of nuclear base substitutions (similar to 23x) and small indels (similar to 328x) relative to wildtype hermaphrodites. However, we observed no increase in the mutation rates of mtDNA, and copy-number changes of single-copy genes. There was a marked increase in copy-number variation of rDNA genes under MMR impairment. In C. elegans, msh-2 repairs transitions more efficiently than transversions and increases the AT mutational bias relative to wildtype. The local sequence context, including sequence complexity, G + C-content, and flanking bases influenced the mutation rate. The X chromosome exhibited lower substitution and higher indel rates than autosomes, which can either result from sex-specific mutation rates or a nonrandom distribution of mutable sites between chromosomes. Provided the observed difference in mutational pattern is mostly due to MMR impairment, our results indicate that the specificity of MMR varies between taxa, and is more efficient in detecting and repairing small indels in eukaryotes relative to prokaryotes.

ABSTRACTTo ascertain the effect of relaxation of selection on global gene expression,Caenorhabditis elegansmutation accumulation (MA) lines were propagated under varying degrees of efficiency of selection determined by their different population sizes (N= 1, 10 and 100). Both the mutational variance (Vm), and the residual variance (Vr) were greatest in MA lines with the lowest efficiency of natural selection. The results suggest that gene expression is under strong balancing selection. Furthermore, mutations resulting in increased transcriptional noise or sensitivity to microenvironmental variation accumulate most under extreme genetic drift. In contrast, theVm/Vrratio was lowest in theN=1 lines. Chromatin domains associated with broad gene silencing and active transcription exhibited the greatest and the smallest increase in transcriptional variation, respectively. Furthermore, the preponderance of overexpressed genes was especially pronounced in mitochondrial respiration, stress response, and immune system pathways, especially in low fitnessN= 1 lines.

The ecology and environment of the microbes that inhabit the mammalian intestine undergoes several changes as the host ages. Here, we ask if the selection pressure experienced by a new strain colonizing the aging gut differs from that in the gut of young adults. Using experimental evolution in mice after a short antibiotic treatment, as a model for a common clinical situation, we show that a new colonizing E. coli strain rapidly adapts to the aging gut via both mutation accumulation and bacteriophage-mediated horizontal gene transfer (HGT). The pattern of evolution of E. coli in aging mice is characterized by a larger number of transposable element insertions and intergenic mutations compared to that in young mice, which is consistent with the gut of aging hosts harboring a stressful and iron limiting environment.

Experimental investigations into the rates and fitness effects of spontaneous mutations are fundamental to our understanding of the evolutionary process. To gain insights into the molecular and fitness consequences of spontaneous mutations, we conducted a mutation accumulation (MA) experiment at varying population sizes in the nematode Caenorhabditis elegans, evolving 35 lines in parallel for 409 generations at three population sizes (N = 1, 10, and 100 individuals). Here, we focus on nuclear SNPs and small insertion/deletions (indels) under minimal influence of selection, as well as their accrual rates in larger populations under greater selection efficacy. The spontaneous rates of base substitutions and small indels are 1.84 (95% C.I. ± 0.14) × 10-9 substitutions and 6.84 (95% C.I. ± 0.97) × 10-10 changes/site/generation, respectively. Small indels exhibit a deletion bias with deletions exceeding insertions by threefold. Notably, there was no correlation between the frequency of base substitutions, nonsynonymous substitutions, or small indels with population size. These results contrast with our previous analysis of mitochondrial DNA mutations and nuclear copy-number changes in these MA lines, and suggest that nuclear base substitutions and small indels are under less stringent purifying selection compared to the former mutational classes. A transition bias was observed in exons as was a near universal base substitution bias toward A/T. Strongly context-dependent base substitutions, where 5'-Ts and 3'-As increase the frequency of A/T → T/A transversions, especially at the boundaries of A or T homopolymeric runs, manifest as higher mutation rates in (i) introns and intergenic regions relative to exons, (ii) chromosomal cores vs. arms and tips, and (iii) germline-expressed genes.

Significance Copy-number variants are ubiquitous in nature, yet their immediate functional consequences are obscure. We conducted a spontaneous mutation accumulation experiment at varying sizes in Caenorhabditis elegans , thereby enabling the simultaneous investigation of the mutational input and strength of selection on the evolution of copy-number changes. Whole-genome sequencing reveals the highest genome-wide rate of gene duplication for any species thus far. Our transcriptome analysis further demonstrates that gene duplication frequently results in a greater than two-fold change in transcription. Despite the adaptive role of duplication as the primary source of novel genes, we find duplications and deletions of highly transcribed genes to be more detrimental to fitness and evidence for selection against increase in transcript abundance.

Mitochondrial genomes of metazoans, given their elevated rates of evolution, have served as pivotal markers for phylogeographic studies and recent phylogenetic events. In order to determine the dynamics of spontaneous mitochondrial mutations in small populations in the absence and presence of selection, we evolved mutation accumulation (MA) lines of Caenorhabditis elegans in parallel over 409 consecutive generations at three varying population sizes of N = 1, 10, and 100 hermaphrodites. The N =1 populations should have a minimal influence of natural selection to provide the spontaneous mutation rate and the expected rate of neutral evolution, whereas larger population sizes should experience increasing intensity of selection. New mutations were identified by Illumina paired-end sequencing of 86 mtDNA genomes across 35 experimental lines and compared with published genomes of natural isolates. The spontaneous mitochondrial mutation rate was estimated at 1.05 × 10-7/site/generation. A strong G/C→A/T mutational bias was observed in both the MA lines and the natural isolates. This suggests that the low G + C content at synonymous sites is the product of mutation bias rather than selection as previously proposed. The mitochondrial effective population size per worm generation was estimated to be 62. Although it was previously concluded that heteroplasmy was rare in C. elegans, the vast majority of mutations in this study were heteroplasmic despite an experimental regime exceeding 400 generations. The frequencies of frameshift and nonsynonymous mutations were negatively correlated with population size, which suggests their deleterious effects on fitness and a potent role for selection in their eradication.

Deoxyribonucleoside kinases (dNKs) carry out the rate-determining step in the nucleoside salvage pathway within all domains of life where the pathway is present, and, hence, are an indication on whether or not a species/genus retains the ability to salvage deoxyribonucleosides. Here, a phylogenetic tree is constructed for the thymidine kinase 2-like dNK gene family in metazoa. Each enzyme class (deoxycytidine, deoxyguanosine, and deoxythymidine kinases, as well as the multisubstrate dNKs) falls into a monophyletic clade. However, in vertebrates, dCK contains an apparent duplication with one paralog lost in mammals, and a number of crustacean genomes (like Caligus rogercresseyi and Lepeophtheirus salmonis) unexpectedly contain not only the multisubstrate dNKs, related to Drosophila multisubstrate dNK, but also a TK2-like kinase. Additionally, crustaceans (Daphnia, Caligus, and Lepeophtheirus) and some insects (Tribolium, Danaus, Pediculus, and Acyrthosiphon) contain several multisubstrate dNK-like enzymes which group paraphyletically within the arthropod clade. This might suggest that the multisubstrate dNKs underwent multiple rounds of duplications with differential retention of duplicate copies between insect families and more complete retention within some crustaceans and insects. Genomes of several basal animalia contain more than one dNK-like sequence, some of which group outside the remaining eukaryotes (both plants and animals) and/or with bacterial dNKs. Within the vertebrates, the mammalian genomes do not contain the second dCK, while birds, fish, and amphibians do retain it. Phasianidae (chicken and turkey) have lost dGK, while it has been retained in other bird lineages, like zebra finch. Reconstruction of the ancestral sequence between the multisubstrate arthropod dNKs and the TK2 clade of vertebrates followed by homology modeling and discrete molecular dynamics calculations on this sequence were performed to examine the evolutionary path which led to the two different enzyme classes. The structural models showed that the carboxyl terminus of the ancestral sequence is more helical than dNK, in common with TK2, although any implications of this for enzyme specificity will require biochemical validation. Finally, rate-shift and conservation-shift analysis between clades with different specificities uncovered candidate residues outside the active site pocket which may have contributed to differentiation in substrate specificity between enzyme clades.

Nucleoside phosphorylases are essential for the salvage and catabolism of nucleotides in bacteria and other organisms, and members of this enzyme superfamily have been of interest for the development of antimicrobial and cancer therapies. The nucleotide phosphorylase superfamily 1 encompasses a number of different enzymes which share a general superfold and catalytic mechanism, while they differ in the nature of the nucleophiles used and in the nature of characteristic active site residues. Recently, one subfamily, the uridine phosphorylases, has been subdivided into two types which differ with respect to the mechanism of transition state stabilization, as dictated by differences in critical amino acid residues. Little is known about the phylogenetic distribution and relationship of the two different types, as well as the relationship to other NP-1 superfamily members. Here comparative genomic analysis illustrates that UP-1s and UP-2s fall into monophyletic groups and are biased with respect to species representation. UP-1 evolved in Gram negative bacteria, while Gram positive species tend to predominantly contain UP-2. PNP (a sister clade to all UPs) contains both Gram positive and Gram negative species. The findings imply that the nucleoside phosphorylase superfamily 1 evolved through a series of three important duplications, leading to the separate, monophyletic enzyme families, coupled to individual lateral transfer events. Extensive horizontal transfer explains the occurrence of unexpected uridine phosphorylases in some genomes. This study provides a basis for understanding the evolution of uridine and purine nucleoside phosphorylases with respect to DNA/RNA metabolism and with potential utility in the design of antimicrobial and anti-tumor drugs.

Deoxyribonucleoside kinases (dNKs) are important to DNA metabolism, especially in environments where nucleosides are freely available to be absorbed and used for the salvage pathway. Little has previously been known about the complement of dNKs in different bacterial genomes. However, it was believed that Gram-negative bacteria had a single dNK, while Gram-positive bacteria possessed several. An analysis of 992 fully sequenced bacterial genomes, including both Gram-positive and Gram-negative organisms, was conducted to investigate the phylogenetic relationship of all TK1-like and non-TK1-like dNKs. It was illustrated that both gene families evolved through a number of duplications and horizontal gene transfers, leading to the presence of multiple dNKs in different types of bacteria. The findings of this study provide a backbone for further studies into the evolution of the interplay between the de novo and salvage pathways in DNA synthesis with respect to environmental availability of deoxyribonucleosides and metabolic processes generating the provisions of different dNTPs.

Deoxyribonucleoside kinases (dNKs) are essential in the mammalian cell but their 'importance' in bacteria, especially aquatic ones, is less clear. We studied two aquatic bacteria, Gram-negative Flavobacterium psychrophilum JIP02/86 and Polaribacter sp. MED152, for their ability to salvage deoxyribonucleosides (dNs). Both had a Gram-positive-type thymidine kinase (TK1), which could phosphorylate thymidine, and one non-TK1 dNK, which could efficiently phosphorylate deoxyadenosine and slightly also deoxycytosine. Surprisingly, the four tested dNKs could not phosphorylate deoxyguanosine, and apparently, these two bacteria are missing this activity. When tens of available aquatic bacteria genomes were examined for the presence of dNKs, a majority had at least a TK1-like gene, but several lacked any dNKs. Apparently, among aquatic bacteria, the role of the dN salvage varies.

Gene duplication is an important process in the functional divergence of genes and genomes. Several processes have been described that lead to duplicate gene retention over different timescales after both smaller-scale events and whole-genome duplication, including neofunctionalization, subfunctionalization, and dosage balance. Two common modes of duplicate gene loss include nonfunctionalization and loss due to population dynamics (failed fixation). Previous work has characterized expectations of duplicate gene retention under the neofunctionalization and subfunctionalization models. Here, that work is extended to dosage balance using simulations. A general model for duplicate gene loss/retention is then presented that is capable of fitting expectations under the different models, is defined at t = 0, and decays to an orthologous asymptotic rate rather than zero, based upon a modified Weibull hazard function. The model in a maximum likelihood framework shows the property of identifiability, recovering the evolutionary mechanism and parameters of simulation. This model is also capable of recovering the evolutionary mechanism of simulation from data generated using an unrelated network population genetic model. Lastly, the general model is applied as part of a mixture model to recent gene duplicates from the Oikopleura dioica genome, suggesting that neofunctionalization may be an important process leading to duplicate gene retention in that organism.

Genomes of animals as different as sponges and humans show conservation of global architecture. Here we show that multiple genomic features including transposon diversity, developmental gene repertoire, physical gene order, and intron-exon organization are shattered in the tunicate Oikopleura, belonging to the sister group of vertebrates and retaining chordate morphology. Ancestral architecture of animal genomes can be deeply modified and may therefore be largely nonadaptive. This rapidly evolving animal lineage thus offers unique perspectives on the level of genome plasticity. It also illuminates issues as fundamental as the mechanisms of intron gain.