Genetic drift, selection and the evolution of the mutation rate.

As one of the few cellular traits that can be quantified across the tree of life, DNA-replication fidelity provides an excellent platform for understanding fundamental evolutionary processes. Furthermore, because mutation is the ultimate source of all genetic variation, clarifying why mutation rates vary is crucial for understanding all areas of biology. A potentially revealing hypothesis for mutation-rate evolution is that natural selection primarily operates to improve replication fidelity, with the ultimate limits to what can be achieved set by the power of random genetic drift. This drift-barrier hypothesis is consistent with comparative measures of mutation rates, provides a simple explanation for the existence of error-prone polymerases and yields a formal counter-argument to the view that selection fine-tunes gene-specific mutation rates.

Asymmetric Context-Dependent Mutation Patterns Revealed through Mutation-Accumulation Experiments.

Despite the general assumption that site-specific mutation rates are independent of the local sequence context, a growing body of evidence suggests otherwise. To further examine context-dependent patterns of mutation, we amassed 5,645 spontaneous mutations in wild- type (WT) and mismatch-repair deficient (MMR(-)) mutation-accumulation (MA) lines of the gram-positive model organism Bacillus subtilis. We then analyzed>7,500 spontaneous base-substitution mutations across B. subtilis, Escherichia coli, and Mesoplasma florum WT and MMR(-) MA lines, finding a context-dependent mutation pattern that is asymmetric around the origin of replication. Different neighboring nucleotides can alter site-specific mutation rates by as much as 75-fold, with sites neighboring G:C base pairs or dimers involving alternating pyrimidine-purine and purine-pyrimidine nucleotides having significantly elevated mutation rates. The influence of context-dependent mutation on genome architecture is strongest in M. florum, consistent with the reduced efficiency of selection in organisms with low effective population size. If not properly accounted for, the disparities arising from patterns of context-dependent mutation can significantly influence interpretations of positive and purifying selection.

Hybridization and the Origin of Contagious Asexuality in Daphnia pulex.

Hybridization plays a potentially important role in the origin of obligate parthenogenesis (OP) in many organisms. However, it remains controversial whether hybridization directly triggers the transition from sexual reproduction to obligate asexuality or a hybrid genetic background enables asexual species to persist. Furthermore, we know little about the specific genetic elements from the divergent, yet still hybridizing lineages responsible for this transition and how these elements are further spread to create other OP lineages. In this study, we address these questions in Daphnia pulex, where cyclically parthenogenetic (CP) and OP lineages coexist. Ancestry estimates and whole-genome association mapping using 32 OP isolates suggest that a complex hybridization history between the parental species D. pulex and D. pulicaria is responsible for the introgression of a set of 647 D. pulicaria single nucleotide polymorphism alleles that show perfect association with OP. Crossing experiments using males of OP lineages and females of CP lineages strongly support a polygenic basis for OP. Single-sperm analyses show that although normal meiotic recombination occurs in the production of haploid sperm by males of OP lineages, a significant proportion of such sperm are polyploid, suggesting that the spread of asexual elements through these males (i.e., contagious asexuality) is much less efficient than previously envisioned. Although the current Daphnia genome annotation does not provide mechanistic insight into the nature of the asexuality-associated alleles, these alleles should be considered as candidates for future investigations on the genetic underpinnings of OP.

Population-genomic insights into the evolutionary origin and fate of obligately asexual Daphnia pulex.

Despite much theoretical work, the molecular-genetic causes and evolutionary consequences of asexuality remain largely undetermined. Asexual animal species are rare, evolutionarily short-lived, and thought to suffer mutational meltdown as a result of lack of recombination. Whole-genome analysis of 11 sexual and 11 asexual genotypes of Daphnia pulex indicates that current asexual lineages are in fact very young, exhibit no signs of purifying selection against accumulating mutations, and have extremely high rates of gene conversion and deletion. The reconstruction of chromosomal haplotypes in regions containing SNP markers associated with asexuality (chromosomes VIII and IX) indicates that introgression from a sister species, Daphnia pulicaria, underlies the origin of the asexual phenotype. Silent-site divergence of the shared chromosomal haplotypes of asexuals indicates that the spread of asexuality is as recent as 1,250 y, although the origin of the meiosis-suppressing element or elements could be substantially older. In addition, using previous estimates of the gene conversion rate from Daphnia mutation accumulation lines, we are able to age each asexual lineage. Although asexual lineages originate from wide crosses that introduce elevated individual heterozygosities on clone foundation, they also appear to be constrained by the inbreeding-like effect of loss of heterozygosity that accrues as gene conversion and hemizygous deletion expose preexisting recessive deleterious alleles of asexuals, limiting their evolutionary longevity. Our study implies that the buildup of newly introduced deleterious mutations (i.e., Muller's ratchet) may not be the dominant force imperiling nonrecombining populations of D. pulex, as previously proposed.

Polynucleobacter necessarius, a model for genome reduction in both free-living and symbiotic bacteria.

We present the complete genomic sequence of the essential symbiont Polynucleobacter necessarius (Betaproteobacteria), which is a valuable case study for several reasons. First, it is hosted by a ciliated protist, Euplotes; bacterial symbionts of ciliates are still poorly known because of a lack of extensive molecular data. Second, the single species P. necessarius contains both symbiotic and free-living strains, allowing for a comparison between closely related organisms with different ecologies. Third, free-living P. necessarius strains are exceptional by themselves because of their small genome size, reduced metabolic flexibility, and high worldwide abundance in freshwater systems. We provide a comparative analysis of P. necessarius metabolism and explore the peculiar features of a genome reduction that occurred on an already streamlined genome. We compare this unusual system with current hypotheses for genome erosion in symbionts and free-living bacteria, propose modifications to the presently accepted model, and discuss the potential consequences of translesion DNA polymerase loss.

Drift-barrier hypothesis and mutation-rate evolution.

Mutation dictates the tempo and mode of evolution, and like all traits, the mutation rate is subject to evolutionary modification. Here, we report refined estimates of the mutation rate for a prokaryote with an exceptionally small genome and for a unicellular eukaryote with a large genome. Combined with prior results, these estimates provide the basis for a potentially unifying explanation for the wide range in mutation rates that exists among organisms. Natural selection appears to reduce the mutation rate of a species to a level that scales negatively with both the effective population size (N(e)), which imposes a drift barrier to the evolution of molecular refinements, and the genomic content of coding DNA, which is proportional to the target size for deleterious mutations. As a consequence of an expansion in genome size, some microbial eukaryotes with large N(e) appear to have evolved mutation rates that are lower than those known to occur in prokaryotes, but multicellular eukaryotes have experienced elevations in the genome-wide deleterious mutation rate because of substantial reductions in N(e).

Estimating Seven Coefficients of Pairwise Relatedness Using Population-Genomic Data.

Population structure can be described by genotypic-correlation coefficients between groups of individuals, the most basic of which are the pairwise relatedness coefficients between any two individuals. There are nine pairwise relatedness coefficients in the most general model, and we show that these can be reduced to seven coefficients for biallelic loci. Although all nine coefficients can be estimated from pedigrees, six coefficients have been beyond empirical reach. We provide a numerical optimization procedure that estimates all seven reduced coefficients from population-genomic data. Simulations show that the procedure is nearly unbiased, even at 3× coverage, and errors in five of the seven coefficients are statistically uncorrelated. The remaining two coefficients have a negative correlation of errors, but their sum provides an unbiased assessment of the overall correlation of heterozygosity between two individuals. Application of these new methods to four populations of the freshwater crustacean Daphnia pulex reveal the occurrence of half siblings in our samples, as well as a number of identical individuals that are likely obligately asexual clone mates. Statistically significant negative estimates of these pairwise relatedness coefficients, including inbreeding coefficients that were typically negative, underscore the difficulties that arise when interpreting genotypic correlations as estimations of the probability that alleles are identical by descent.

A New Reference Genome Assembly for the Microcrustacean Daphnia pulex.

Comparing genomes of closely related genotypes from populations with distinct demographic histories can help reveal the impact of effective population size on genome evolution. For this purpose, we present a high quality genome assembly of Daphnia pulex (PA42), and compare this with the first sequenced genome of this species (TCO), which was derived from an isolate from a population with >90% reduction in nucleotide diversity. PA42 has numerous similarities to TCO at the gene level, with an average amino acid sequence identity of 98.8 and >60% of orthologous proteins identical. Nonetheless, there is a highly elevated number of genes in the TCO genome annotation, with ∼7000 excess genes appearing to be false positives. This view is supported by the high GC content, lack of introns, and short length of these suspicious gene annotations. Consistent with the view that reduced effective population size can facilitate the accumulation of slightly deleterious genomic features, we observe more proliferation of transposable elements (TEs) and a higher frequency of gained introns in the TCO genome.

Evolution of the Insertion-Deletion Mutation Rate Across the Tree of Life.

Mutations are the ultimate source of variation used for evolutionary adaptation, while also being predominantly deleterious and a source of genetic disorders. Understanding the rate of insertion-deletion mutations (indels) is essential to understanding evolutionary processes, especially in coding regions, where such mutations can disrupt production of essential proteins. Using direct estimates of indel rates from 14 phylogenetically diverse eukaryotic and bacterial species, along with measures of standing variation in such species, we obtain results that imply an inverse relationship of mutation rate and effective population size. These results, which corroborate earlier observations on the base-substitution mutation rate, appear most compatible with the hypothesis that natural selection reduces mutation rates per effective genome to the point at which the power of random genetic drift (approximated by the inverse of effective population size) becomes overwhelming. Given the substantial differences in DNA metabolism pathways that give rise to these two types of mutations, this consistency of results raises the possibility that refinement of other molecular and cellular traits may be inversely related to species-specific levels of random genetic drift.

A Male-Specific Genetic Map of the Microcrustacean Daphnia pulex Based on Single-Sperm Whole-Genome Sequencing.

Genetic linkage maps are critical for assembling draft genomes to a meaningful chromosome level and for deciphering the genomic underpinnings of biological traits. The estimates of recombination rates derived from genetic maps also play an important role in understanding multiple aspects of genomic evolution such as nucleotide substitution patterns and accumulation of deleterious mutations. In this study, we developed a high-throughput experimental approach that combines fluorescence-activated cell sorting, whole-genome amplification, and short-read sequencing to construct a genetic map using single-sperm cells. Furthermore, a computational algorithm was developed to analyze single-sperm whole-genome sequencing data for map construction. These methods allowed us to rapidly build a male-specific genetic map for the freshwater microcrustacean Daphnia pulex, which shows significant improvements compared to a previous map. With a total of mapped 1672 haplotype blocks and an average intermarker distance of 0.87 cM, this map spans a total genetic distance of 1451 Kosambi cM and comprises 90% of the resolved regions in the current Daphnia reference assembly. The map also reveals the mistaken mapping of seven scaffolds in the reference assembly onto chromosome II by a previous microsatellite map based on F2 crosses. Our approach can be easily applied to many other organisms and holds great promise for unveiling the intragenomic and intraspecific variation in the recombination rates.

Mutation rate, spectrum, topology, and context-dependency in the DNA mismatch repair-deficient Pseudomonas fluorescens ATCC948.

High levels of genetic diversity exist among natural isolates of the bacterium Pseudomonas fluorescens, and are especially elevated around the replication terminus of the genome, where strain-specific genes are found. In an effort to understand the role of genetic variation in the evolution of Pseudomonas, we analyzed 31,106 base substitutions from 45 mutation accumulation lines of P. fluorescens ATCC948, naturally deficient for mismatch repair, yielding a base-substitution mutation rate of 2.34 × 10(-8) per site per generation (SE: 0.01 × 10(-8)) and a small-insertion-deletion mutation rate of 1.65 × 10(-9) per site per generation (SE: 0.03 × 10(-9)). We find that the spectrum of mutations in prophage regions, which often contain virulence factors and antibiotic resistance, is highly similar to that in the intergenic regions of the host genome. Our results show that the mutation rate varies around the chromosome, with the lowest mutation rate found near the origin of replication. Consistent with observations from other studies, we find that site-specific mutation rates are heavily influenced by the immediately flanking nucleotides, indicating that mutations are context dependent.