Muller's Ratchet Research Articles

Evidence of a putative CO2 delivery system to the chromatophore in the photosynthetic amoeba Paulinella

AbstractThe photosynthetic amoeba, Paulinella provides a recent (ca. 120 Mya) example of primary plastid endosymbiosis. Given the extensive data demonstrating host lineage‐driven endosymbiont integration, we analysed nuclear genome and transcriptome data to investigate mechanisms that may have evolved in Paulinella micropora KR01 (hereinafter, KR01) to maintain photosynthetic function in the novel organelle, the chromatophore. The chromatophore is of α‐cyanobacterial provenance and has undergone massive gene loss due to Muller's ratchet, but still retains genes that encode the ancestral α‐carboxysome and the shell carbonic anhydrase, two critical components of the biophysical CO2 concentrating mechanism (CCM) in cyanobacteria. We identified KR01 nuclear genes potentially involved in the CCM that arose via duplication and divergence and are upregulated in response to high light and downregulated under elevated CO2. We speculate that these genes may comprise a novel CO2 delivery system (i.e., a biochemical CCM) to promote the turnover of the RuBisCO carboxylation reaction and counteract photorespiration. We posit that KR01 has an inefficient photorespiratory system that cannot fully recycle the C2 product of RuBisCO oxygenation back to the Calvin‐Benson cycle. Nonetheless, both these systems appear to be sufficient to allow Paulinella to persist in environments dominated by faster‐growing phototrophs.

A dynamical limit to evolutionary adaptation

Natural selection makes evolutionary adaptation possible even if the overwhelming majority of new mutations are deleterious. However, in rapidly evolving populations where numerous linked mutations occur and segregate simultaneously, clonal interference and genetic hitchhiking can limit the efficiency of selection, allowing deleterious mutations to accumulate over time. This can in principle overwhelm the fitness increases provided by beneficial mutations, leading to an overall fitness decline. Here, we analyze the conditions under which evolution will tend to drive populations to higher versus lower fitness. Our analysis focuses on quantifying the boundary between these two regimes, as a function of parameters such as population size, mutation rates, and selection pressures. This boundary represents a state in which adaptation is precisely balanced by Muller's ratchet, and we show that it can be characterized by rapid molecular evolution without any net fitness change. Finally, we consider the implications of global fitness-mediated epistasis and find that under some circumstances, this can drive populations toward the boundary state, which can thus represent a long-term evolutionary attractor.

Parthenogenomics: Insights on mutation rates and nucleotide diversity in parthenogenetic Panagrolaimus nematodes.

Asexual reproduction is assumed to lead to the accumulation of deleterious mutations, and reduced heterozygosity due to the absence of recombination. Panagrolaimid nematode species display different modes of reproduction. Sexual reproduction with distinct males and females, asexual reproduction through parthenogenesis in the genus Panagrolaimus, and hermaphroditism in Propanagrolaimus. Here, we compared genomic features of free-living nematodes in populations and species isolated from geographically distant regions to study diversity, and genome-wide differentiation under different modes of reproduction. We firstly estimated genome-wide spontaneous mutation rates in a triploid parthenogenetic Panagrolaimus, and a diploid hermaphroditic Propanagrolaimus via long-term mutation accumulation lines. Secondly, we calculated population genetic parameters including nucleotide diversity, and fixation index (F ST) between populations of asexually and sexually reproducing nematodes. Thirdly, we used phylogenetic network methods on sexually and asexually reproducing Panagrolaimus populations to understand evolutionary relationships between them. The estimated mutation rate was slightly lower for the asexual population, as expected for taxa with this reproductive mode. Natural polyploid asexual populations revealed higher nucleotide diversity. Despite their common ancestor, a gene network revealed a high level of genetic differentiation among asexual populations. The elevated heterozygosity found in the triploid parthenogens could be explained by the third genome copy. Given their tendentially lower mutation rates it can be hypothesized that this is part of the mechanism to evade Muller's ratchet. Our findings in parthenogenetic triploid nematode populations seem to challenge common expectations of evolution under asexuality.

Stochastic drift in discrete waves of nonlocally interacting particles.

In this paper, we investigate a generalized model of N particles undergoing second-order nonlocal interactions on a lattice. Our results have applications across many research areas, including the modeling of migration, information dynamics, and Muller's ratchet-the irreversible accumulation of deleterious mutations in an evolving population. Strikingly, numerical simulations of the model are observed to deviate significantly from its mean-field approximation even for large population sizes. We show that the disagreement between deterministic and stochastic solutions stems from finite-size effects that change the propagation speed and cause the position of the wave to fluctuate. These effects are shown to decay anomalously as (lnN)^{-2} and (lnN)^{-3}, respectively-much slower than the usual N^{-1/2} factor. Our results suggest that the accumulation of deleterious mutations in a Muller's ratchet and the loss of awareness in a population may occur much faster than predicted by the corresponding deterministic models. The general applicability of our model suggests that this unexpected scaling could be important in a wide range of real-world applications.

Differences in homologous recombination and maintenance of heteropolyploidy between Haloferax volcanii and Haloferax mediterranei.

Polyploidy, the phenomenon of having more than one copy of the genome in an organism, is common among haloarchaea. While providing short-term benefits for DNA repair, polyploidy is generally regarded as an "evolutionary trap" that by the notion of the Muller's ratchet will inevitably conclude in the species' decline or even extinction due to a gradual reduction in fitness. In most reported cases of polyploidy in archaea, the genetic state of the organism is considered as homoploidy i.e. all copies of the genome are identical. Here we demonstrate that while this is indeed the prevalent genetic status in the halophilic archaeon Haloferax volcanii, its close relative H. mediterranei maintains a prolonged heteroploidy state in a nonselective environment once a second allele is introduced. Moreover, a strong genetic linkage was observed between two distant loci in H. mediterranei indicating a low rate of homologous recombination while almost no such linkage was shown in H. volcanii indicating a high rate of recombination in the latter species. We suggest that H. volcanii escapes Muller's ratchet by means of an effective chromosome-equalizing gene-conversion mechanism facilitated by highly active homologous recombination, whereas H. mediterranei must elude the ratchet via a different, yet to be elucidated mechanism.

Repeat sequences limit the effectiveness of lateral gene transfer and favored the evolution of meiotic sex in early eukaryotes

The transition from prokaryotic lateral gene transfer to eukaryotic meiotic sex is poorly understood. Phylogenetic evidence suggests that it was tightly linked to eukaryogenesis, which involved an unprecedented rise in both genome size and the density of genetic repeats. Expansion of genome size raised the severity of Muller's ratchet, while limiting the effectiveness of lateral gene transfer (LGT) at purging deleterious mutations. In principle, an increase in recombination length combined with higher rates of LGT could solve this problem. Here, we show using a computational model that this solution fails in the presence of genetic repeats prevalent in early eukaryotes. The model demonstrates that dispersed repeat sequences allow ectopic recombination, which leads to the loss of genetic information and curtails the capacity of LGT to prevent mutation accumulation. Increasing recombination length in the presence of repeat sequences exacerbates the problem. Mutational decay can only be resisted with homology along extended sequences of DNA. We conclude that the transition to homologous pairing along linear chromosomes was a key innovation in meiotic sex, which was instrumental in the expansion of eukaryotic genomes and morphological complexity.

The extinction time under mutational meltdown driven by high mutation rates.

Mutational meltdown describes an eco‐evolutionary process in which the accumulation of deleterious mutations causes a fitness decline that eventually leads to the extinction of a population. Possible applications of this concept include medical treatment of RNA virus infections based on mutagenic drugs that increase the mutation rate of the pathogen. To determine the usefulness and expected success of such an antiviral treatment, estimates of the expected time to mutational meltdown are necessary. Here, we compute the extinction time of a population under high mutation rates, using both analytical approaches and stochastic simulations. Extinction is the result of three consecutive processes: (a) initial accumulation of deleterious mutations due to the increased mutation pressure; (b) consecutive loss of the fittest haplotype due to Muller's ratchet; (c) rapid population decline toward extinction. We find accurate analytical results for the mean extinction time, which show that the deleterious mutation rate has the strongest effect on the extinction time. We confirm that intermediate‐sized deleterious selection coefficients minimize the extinction time. Finally, our simulations show that the variation in extinction time, given a set of parameters, is surprisingly small.

Population genetics of polymorphism and divergence in rapidly evolving populations.

In rapidly evolving populations, numerous beneficial and deleterious mutations can arise and segregate within a population at the same time. In this regime, evolutionary dynamics cannot be analyzed using traditional population genetic approaches that assume that sites evolve independently. Instead, the dynamics of many loci must be analyzed simultaneously. Recent work has made progress by first analyzing the fitness variation within a population, and then studying how individual lineages interact with this traveling fitness wave. However, these "traveling wave" models have previously been restricted to extreme cases where selection on individual mutations is either much faster or much slower than the typical coalescent timescale Tc. In this work, we show how the traveling wave framework can be extended to intermediate regimes in which the scaled fitness effects of mutations (Tcs) are neither large nor small compared to one. This enables us to describe the dynamics of populations subject to a wide range of fitness effects, and in particular, in cases where it is not immediately clear which mutations are most important in shaping the dynamics and statistics of genetic diversity. We use this approach to derive new expressions for the fixation probabilities and site frequency spectra of mutations as a function of their scaled fitness effects, along with related results for the coalescent timescale Tc and the rate of adaptation or Muller's ratchet. We find that competition between linked mutations can have a dramatic impact on the proportions of neutral and selected polymorphisms, which is not simply summarized by the scaled selection coefficient Tcs. We conclude by discussing the implications of these results for population genetic inferences.

Mutation accumulation in inbreeding populations under evolution of the selfing rate.

It is theoretically established that self-fertilization can facilitate mutation accumulation, thus increasing extinction risk. However, in previous studies, selfing rates are often set as fixed parameters, but in natural systems, evolution of selfing rates and deleterious mutations may mutually affect each other. I carried out simulations to investigate the dynamics of selfing rates and mutation accumulation, by allowing deleterious mutations to coevolve with alleles that modify the selfing rate (selfing modifiers). I found that selfing rates will often fluctuate over time, due to successive invasion of alleles that increase selfing and outcrossing. Since mutation fixation is mainly caused by Muller's ratchet, its rate is sensitive to the change of the selfing rate mutations will accumulate in a punctuated pattern. The dynamics are influenced by several factors, such as recombination and the selfing rate effects of selfing modifier loci. Also, such temporal variation produces variation of selfing rates and mutation accumulation rates between multiple conspecific populations, which can increase the average fitness across populations. As factors, such as the genomic mutation rate of deleterious mutations, can simultaneously influence the selfing rate and mutation fixation, effects of these factors on mutation accumulation rates can be complicated and non-monotonic.

Ribozyme Mutagenic Evolution: Mechanisms of Survival.

Primeval populations replicating at high error rates required a mechanism to overcome the accumulation of mutations and information deterioration. Known strategies to overcome mutation pressures include RNA processivity, epistasis, selection, and quasispecies. We investigated the mechanism by which small molecular ribozyme populations can survive under high error rates by propagating several lineages under different mutagen concentrations. We found that every population that evolved without mutagen went extinct, while those subjected to mutagenic evolution survived. To understand how they survived, we characterized the evolved genotypic diversity, the formation of genotype-genotype interaction networks, the fitness of the most common mutants for each enzymatic step, and changes in population size along the course of evolution. We found that the elevated mutation rate was necessary for the populations to survive in the novel environment, in which all the steps of the metabolism worked to promote the survival of even less catalytically efficient ligases. Besides, an increase in population size and the mutational coupling of genotypes in close-knit networks, which helped maintain or recover lost genotypes making their disappearance transient, prevented Muller's ratchet and extinction.

Spike protein mutational landscape in India during the complete lockdown phase: Could Muller's ratchet be a future game-changer for COVID-19?

The dire need of effective preventive measures and treatment approaches against SARS-CoV-2 virus, causing COVID-19 pandemic, calls for an in-depth understanding of its evolutionary dynamics with attention to specific geographic locations, since lockdown and social distancing to prevent the virus spread could lead to distinct localized dynamics of virus evolution within and between countries owing to different environmental and host-specific selection pressures. To decipher any correlation between SARS-CoV-2 evolution and its epidemiology in India, we studied the mutational diversity of spike glycoprotein, the key player for the attachment, fusion and entry of virus to the host cell. For this, we analyzed the sequences of 630 Indian isolates as available in GISAID database till June 07, 2020 (during the time-period before the start of Unlock 1.0 in India on and from June 08, 2020), and detected the spike protein variants to emerge from two major ancestors – Wuhan-Hu-1/2019 and its D614G variant. Average stability of the docked spike protein – host receptor (S-R) complexes for these variants correlated strongly (R2 = 0.96) with the fatality rates across Indian states. However, while more than half of the variants were found unique to India, 67% of all variants showed lower stability of S-R complex than the respective ancestral variants, indicating a possible fitness loss in recently emerged variants, despite a continuous increase in mutation rate. These results conform to the sharply declining fatality rate countrywide (>7-fold during April 11 – June 28, 2020). Altogether, while we propose the potential of S-R complex stability to track disease severity, we urge an immediate need to explore if SARS-CoV-2 is approaching mutational meltdown in India.

Two linked loci under mutation-selection balance and Muller’s ratchet

We report complete analysis of a deterministic model of deleterious mutations and negative selection against them at two haploid loci without recombination. As long as mutation is a weaker force than selection, mutant alleles remain rare at the only stable equilibrium, and otherwise, a variety of dynamics are possible. If the mutation-free genotype is absent, generally the only stable equilibrium is the one that corresponds to fixation of the mutant allele at the locus where it is less deleterious. This suggests that fixation of a deleterious allele that follows a click of the Muller's ratchet is governed by natural selection, instead of random drift.

Avoiding organelle mutational meltdown across eukaryotes with or without a germline bottleneck.

Mitochondrial DNA (mtDNA) and plastid DNA (ptDNA) encode vital bioenergetic apparatus, and mutations in these organelle DNA (oDNA) molecules can be devastating. In the germline of several animals, a genetic "bottleneck" increases cell-to-cell variance in mtDNA heteroplasmy, allowing purifying selection to act to maintain low proportions of mutant mtDNA. However, most eukaryotes do not sequester a germline early in development, and even the animal bottleneck remains poorly understood. How then do eukaryotic organelles avoid Muller's ratchet-the gradual buildup of deleterious oDNA mutations? Here, we construct a comprehensive and predictive genetic model, quantitatively describing how different mechanisms segregate and decrease oDNA damage across eukaryotes. We apply this comprehensive theory to characterise the animal bottleneck with recent single-cell observations in diverse mouse models. Further, we show that gene conversion is a particularly powerful mechanism to increase beneficial cell-to-cell variance without depleting oDNA copy number, explaining the benefit of observed oDNA recombination in diverse organisms which do not sequester animal-like germlines (for example, sponges, corals, fungi, and plants). Genomic, transcriptomic, and structural datasets across eukaryotes support this mechanism for generating beneficial variance without a germline bottleneck. This framework explains puzzling oDNA differences across taxa, suggesting how Muller's ratchet is avoided in different eukaryotes.

The Novel Genetic Background of Infectious Bursal Disease Virus Strains Emerging from the Action of Positive Selection.

The circulation in Europe of novel reassortant strains of infectious bursal disease virus (IBDV), containing a unique genetic background composition, represents a serious problem for animal health. Since the emergence of this novel IBDV mosaic was first described in Poland, this scenario has become particularly attractive to uncover the evolutionary forces driving the genetic diversity of IBDV populations. This study additionally addressed the phenotypic characterization of these emergent strains, as well as the main features affecting the viral fitness during the competition process of IBDV lineages in the field. Our results showed how different evolutionary mechanisms modulate the genetic diversity of co-existent IBDV lineages, leading to the error catastrophe effect, Muller ratchet effect, or prevalence, depending on their genetic compositions. We also determined that the action of the positive selection pressure, depending on the genomic segment on which it is acting, can drive two main phenotypes for IBDV: immune-escaping strains from the selection on segment A or strains with functional advantages from the selection on segment B. This last group seems to possess an increased fitness landscape in the viral quasispecies composition, presenting better adaptability to dissimilar environmental conditions and likely becoming the dominant population. The reassortant strains also exhibited a lower mortality rate compared with the well-known vvIBDV strains, which can facilitate their spreading.

Sex and recombination purge the genome of deleterious alleles: An Individual Based Modeling Approach

The prevalence of sexual reproduction in most animal species despite its considerable costs such as useless males, energy spent on mating, the cost of meiosis and genome dilution remains a puzzle in evolutionary theory. One prominent single factor attempt to solve this persistent puzzle is the claim that sexual reproduction is instrumental in eliminating deleterious alleles from the species genome by the mechanism of recombination. There are three major versions of the deleterious allele hypothesis: First, the mutational deterministic hypothesis (MDH), which rests on the assumption of negative epistasis, predicts that recombination will help to purge the species genome of deleterious alleles by breaking apart linkages between these alleles. The assumption is that the joint negative effects of linked deleterious alleles is sometimes greater than the effects of the alleles considered separately. Second, there is the hypothesis that sexual reproduction speeds up purifying (negative) selection, which purges the genome of deleterious alleles. Alleles that are less deleterious than the wild type are naturally selected. These alleles, attained via recombination, are sometimes ‘leaky’ mutations giving rise to reduced functionality of attendant proteins. This hypothesis does not necessarily rest on the assumption of negative epistasis, which some argue is relatively rare in nature (Kouyos, Silander and Bonhoeffer (2012)) and which arguably could be seen as a virtue of the purifying selection hypothesis vs. the MDH. Third, Muller's ratchet hypothesis predicts that recombination will help to prevent the buildup of deleterious mutations by the mechanism of recombination. In this study, we focus primarily on testing the purifying selection hypothesis. We performed an individual-based model computer simulation using the program EcoSim to test this hypothesis. The experimental runs for sexual reproduction, asexual reproduction and facultative reproduction involved introducing a deleterious allele into the genome, which exacts an intermediate-level energy penalty on individuals. It was found that whereas on average, deleteriousness consistently declined over 18,000 time-steps due to recombination in sexual reproduction, deleteriousness did not decline for asexual and facultative runs. These results corroborate the hypothesis that recombination due to sexual reproduction helps to eliminate deleterious alleles from the genome through the selection of reduced function mutations.

No evidence for accumulation of deleterious mutations and fitness degradation in clonal fish hybrids: Abandoning sex without regrets.

Despite its inherent costs, sexual reproduction is ubiquitous in nature, and the mechanisms to protect it from a competitive displacement by asexuality remain unclear. Popular mutation‐based explanations, like the Muller's ratchet and the Kondrashov's hatchet, assume that purifying selection may not halt the accumulation of deleterious mutations in the nonrecombining genomes, ultimately leading to their degeneration. However, empirical evidence is scarce and it remains particularly unclear whether mutational degradation proceeds fast enough to ensure the decay of clonal organisms and to prevent them from outcompeting their sexual counterparts. To test this hypothesis, we jointly analysed the exome sequences and the fitness‐related phenotypic traits of the sexually reproducing fish species and their clonal hybrids, whose evolutionary ages ranged from F1 generations to 300 ky. As expected, mutations tended to accumulate in the clonal genomes in a time‐dependent manner. However, contrary to the predictions, we found no trend towards increased nonsynonymity of mutations acquired by clones, nor higher radicality of their amino acid substitutions. Moreover, there was no evidence for fitness degeneration in the old clones compared with that in the younger ones. In summary, although an efficacy of purifying selection may still be reduced in the asexual genomes, our data indicate that its efficiency is not drastically decreased. Even the oldest investigated clone was found to be too young to suffer fitness consequences from a mutation accumulation. This suggests that mechanisms other than mutation accumulation may be needed to explain the competitive advantage of sex in the short term.

SARS‐CoV‐2 variants: Relevance for symptom granularity, epidemiology, immunity (herd, vaccines), virus origin and containment?

SummaryThe origin of the SARS‐CoV‐2 virus remains enigmatic. It is likely to be a continuum resulting from inevitable mutations and recombination events. These genetic changes keep developing in the present epidemic. Mutations tending to deplete the genome in its cytosine content will progressively lead to attenuation as a consequence of Muller's ratchet, but this is counteracted by recombination when different mutants co‐infect the same host, in particular, in clusters of infection. Monitoring as a function of time the genome sequences in closely related cases is critical to anticipate the future of SARS‐CoV‐2 and hence of COVID‐19.

The Chromatin Extrusion Phenomenon in Amoeba proteus Cell Cycle.

Amoeba proteus is possibly the best known of all unicellular eukaryotes. At the same time, several quintessential issues of its biology, including some aspects of the cell cycle, remain unsolved. Here, we show that this obligate agamic amoebae and related species have a special type of cyclic polyploidy. Their nucleus has an euploid status only for a small fraction of the cell cycle, during metaphase and telophase. The rest of the time it has an aneuploid status, which is a consequence of polyploidization. Extrusion of "excess" chromatin from the nucleus in late interphase and during prophase results in depolyploidization. Such a strategy of life cycle in unicellular eukaryotes is thought to be the main mechanism of "resetting" the Muller's ratchet and is a satisfactory alternative to the meiotic recombination for agamic protists.

Genetic and epigenetic Muller's ratchet as a mechanism of frailty and morbidity during aging: a demographic genetic model.

Mutation accumulation has been proposed as a cause of senescence. During this process, age-related genetic and epigenetic mutations steadily accumulate. Cascading deleterious effects of mutations might initiate a steady "accumulation of deficits" in cells, despite the existence of repair mechanisms, leading to cellular senescence and functional decline of tissues and organs, which ultimately manifest as frailty and disease. Here, we investigate several of these aspects in differentiating cell populations through modeling and simulation using the Moran birth-death (demographic) process, under several scenarios of mutation accumulation. Deleterious mutations seem to rapidly accumulate particularly early in the course of life, during which the rate of cell division is high, thereby exerting a greater effect on subsequent cellular senescence. Our results are compatible with the principle of the Muller's ratchet taking place in asexually reproducing organisms. The ratchet speed in a given tissue depends on the size of the cell population, mutation rate and the impact of such mutations on cell phenotypes. It varies substantially among cells in different tissues and organs due to heterogeneity in relation to cell and organ-specific demographic features. Ratchet accelerates particularly after middle age, resulting in a synergistic fitness decay at the level of cell populations. We extend Fisher's average excess concept and rank order scale to interpret differential phenotypic effects of the increase of the mutation load among cell populations within a given tissue. We postulate that classical evolutionary genetic models can explain, at least in part, the origins of frailty, subclinical conditions, morbidity and the health consequences of senescence.

Consequences of recombination for the evolution of the mating type locus in Chlamydomonas reinhardtii.

Recombination suppression in sex chromosomes and mating type loci can lead to degeneration as a result of reduced selection efficacy and Muller's ratchet effects. However, genetic exchange in the form of noncrossover gene conversions may still take place within crossover-suppressed regions. Recent work has found evidence that gene conversion may explain the low degrees of allelic differentiation in the dimorphic mating-type locus (MT) of the isogamous alga Chlamydomonas reinhardtii. However, no one has tested whether gene conversion is sufficient to avoid the degeneration of functional sequence within MT. Here, we calculate degree of linkage disequilibrium (LD) across MT as a proxy for recombination rate and investigate its relationship to patterns of population genetic variation and the efficacy of selection in the region. We find that degree of LD predicts selection efficacy across MT, and that purifying selection is stronger in shared genes than in MT-limited genes to the point of being equivalent to that of autosomal genes. We argue that while crossover suppression is needed in the mating-type loci of many isogamous systems, these loci are less likely to experience selection to differentiate further. Thus, recombination can act in these regions and prevent degeneration caused by Hill-Robertson effects.