Impacts of mutation effects and population size on mutation rate in asexual populations: a simulation study

Effects of urban-induced mutations on ecology, evolution and health

Increasing evidence suggests that urbanization is associated with higher mutation rates, which can affect the health and evolution of organisms that inhabit cities. Elevated pollution levels in urban areas can induce DNA damage, leading to de novo mutations. Studies on mutations induced by urban pollution are most prevalent in humans and microorganisms, whereas studies of non-human eukaryotes are rare, even though increased mutation rates have the potential to affect organisms and their populations in contemporary time. Our Perspective explores how higher mutation rates in urban environments could impact the fitness, ecology and evolution of populations. Most mutations will be neutral or deleterious, and higher mutation rates associated with elevated pollution in urban populations can increase the risk of cancer in humans and potentially other species. We highlight the potential for urban-driven increased deleterious mutational loads in some organisms, which could lead to a decline in population growth of a wide diversity of organisms. Although beneficial mutations are expected to be rare, we argue that higher mutation rates in urban areas could influence adaptive evolution, especially in organisms with short generation times. Finally, we explore avenues for future research to better understand the effects of urban-induced mutations on the fitness, ecology and evolution of city-dwelling organisms.

Improvement of Akanthomyces lecanii resistance to tebuconazole through UV-C radiation and selective pressure on microbial evolution and growth arenas

Tebuconazole (TEB) is a fungicide widely used in agriculture; however, its constant application has increased the emergence of resistant plant pathogenic fungal strains and reduced the effectiveness of fungi as biological control agents; for instance, the entomopathogenic and hyperparasitic fungus Akanthomyces lecanii, suitable for simultaneous biological control of insect pest and plant pathogenic fungi, is highly sensitive to fungicides. We carried out the induction of resistance to TEB in two wild type strains of A. lecanii by UV radiation and selective pressure in increasing fungicide gradients using a modified Microbial Evolution and Growth Arena (MEGA), to produce A. lecanii strains that can be used as biological control agent in the presence of tebuconazole. Nine UV-induced and three naturally adapted A. lecanii strains were resistant to TEB at the agriculturally recommended dose, and three irradiated strains were resistant to TEB concentration ten times higher; moreover, growth, sporulation rates, production of hydrolytic enzymes, and virulence against the hemipteran Coccus viridis, a major pest of coffee crops, were not affected in the TEB-resistant strains. These A. lecanii TEB-resistant strains would have a greater opportunity to develop and to establish themselves in fields where the fungicide is present and can be used in a combined biological-chemical strategy to improve insect and plant pathogenic fungal control in agriculture. Also, the selective pressure through modified MEGA plate methodology can be used for the adaptation of entomopathogenic filamentous fungi to withstand other chemical or abiotic stresses that limits its effectiveness for pest control.

Experimental and In Silico Analysis of TEM β-Lactamase Adaptive Evolution

Multiple mutations often have non-additive (epistatic) phenotypic effects. Epistasis is of fundamental biological relevance but is not well understood mechanistically. Adaptive evolution, i.e., the evolution of new biochemical activities, is rich in epistatic interactions. To better understand the principles underlying epistasis during genetic adaptation, we studied the evolution of TEM-1 β-lactamase variants exhibiting cefotaxime resistance. We report the collection of a library of 487 observed evolutionary trajectories for TEM-1 and determine the epistasis status based on cefotaxime resistance phenotype for 206 combinations of 2-3 TEM-1 mutations involving 17 positions under adaptive selective pressure. Gain-of-function (GOF) mutations are gatekeepers for adaptation. To see if GOF phenotypes can be inferred based solely on sequence data, we calculated the enrichment of GOF mutations in the different categories of epistatic pairs. Our results suggest that this is possible because GOF mutations are particularly enriched in sign and reciprocal sign epistasis, which leave a major imprint on the sequence space accessible to evolution. We also used FoldX to explore the relationship between thermodynamic stability and epistasis. We found that mutations in observed evolutionary trajectories tend to destabilize the folded structure of the protein, albeit their cumulative effects are consistently below the protein's free energy of folding. The destabilizing effect is stronger for epistatic pairs, suggesting that modest or local alterations in folding stability can modulate catalysis. Finally, we report a significant relationship between epistasis and the degree to which two protein positions are structurally and dynamically coupled, even in the absence of ligand.

Epidemic spreading under mutually independent intra- and inter-host pathogen evolution

The dynamics of epidemic spreading is often reduced to the single control parameter R₀ (reproduction-rate), whose value, above or below unity, determines the state of the contagion. If, however, the pathogen evolves as it spreads, R₀ may change over time, potentially leading to a mutation-driven spread, in which an initially sub-pandemic pathogen undergoes a breakthrough mutation. To predict the boundaries of this pandemic phase, we introduce here a modeling framework to couple the inter-host network spreading patterns with the intra-host evolutionary dynamics. We find that even in the extreme case when these two process are driven by mutually independent selection forces, mutations can still fundamentally alter the pandemic phase-diagram. The pandemic transitions, we show, are now shaped, not just by R₀, but also by the balance between the epidemic and the evolutionary timescales. If mutations are too slow, the pathogen prevalence decays prior to the appearance of a critical mutation. On the other hand, if mutations are too rapid, the pathogen evolution becomes volatile and, once again, it fails to spread. Between these two extremes, however, we identify a broad range of conditions in which an initially sub-pandemic pathogen can breakthrough to gain widespread prevalence.

Population Genetic Analysis of the Threatened Plant Leavenworthia exigua var. laciniata (Brassicaceae) Reveals Virtually No Genetic Diversity and a Unique Mating System

Leavenworthia (Brassicaceae) has served as a model group for investigating the evolution of mating systems in plants, yet several Leavenworthia species remain understudied. One such taxon is Leavenworthia exigua var. laciniata, one of three varieties of L. exigua, a winter-annual plant endemic to the central United States. Because L. exigua var. laciniata occupies a narrow geographic range and is experiencing major habitat loss, it was recently listed as threatened; however, little is known about its genetic diversity and implications for conservation. We conducted a range-wide population genetic study of L. exigua var. laciniata and L. exigua var. exigua to understand: (1) levels of genetic diversity within and among populations, (2) whether L. exigua var. laciniata is genetically distinct from L. exigua var. exigua, and (3) implications for conservation. L. exigua var. laciniata showed identical genotypes at all 16 microsatellite loci across most of its range, fixed heterozygosity at some loci, and significant heterozygote excesses, consistent with a lack of recombination associated with an asexual mating system, which has not been documented previously in Leavenworthia. Because L. exigua var. laciniata is an annual and the same genotype occurs across multiple populations, asexuality may be caused by apomixis, asexual reproduction via seed. In contrast, most populations of L. exigua var. exigua demonstrated population genetic patterns consistent with a self-compatible mating system. Because L. exigua var. laciniata is morphologically, geographically, and genetically distinct, it should be recognized as an evolutionarily significant unit for conservation. We recommend maintaining large population sizes to conserve evolutionary potential in L. exigua var. laciniata, as the likelihood that facultative sexual reproduction may occur may be greater in larger populations. Additional research in L. exigua var. laciniata is needed to confirm the occurrence of asexuality and apomixis, clarify its reproductive isolation from other taxa, and to understand whether it exhibits residual sexual reproduction, epigenetic variation, or phenotypic plasticity to help it persist in response to environmental variation. In the future, L. exigua var. laciniata may serve as an important model in which to investigate the conservation of threatened plant species with little genetic variation in a changing climate.

High Polymorphism in the Dmrt2a Gene Is Incompletely Sex-Linked in Spotted Scat, Scatophagus argus

Unlike mammals and birds, many fishes have young sex chromosomes, providing excellent models to study sex chromosome differentiation at early stages. Previous studies showed that spotted scat possesses an XX-XY sex determination system. The X has a complete Dmrt3 copy (termed normal) and a truncated copy of Dmrt1 (called Dmrt1b), while the Y has the opposite (normal Dmrt1, which is male-specific, and a truncated Dmrt3 called Dmrt3△-Y). Dmrt1 is the candidate sex determination gene, while the differentiation of other sex-linked genes remains unknown. The spotted scat has proven to be a good model to study the evolution of sex chromosomes in vertebrates. Herein, we sequenced a neighbor gene of this family, Dmrt2, positioned farther from Dmrt1 and closer to Dmrt3 in the spotted scat, and analyzed its sequence variation and expression profiles. The physical locations of the three genes span across an estimated size of >40 kb. The open reading frames of Dmrt2a and its paralog Dmrt2b are 1578 bp and 1311 bp, encoding peptides of 525 and 436 amino acid residues, respectively. Dmrt2a is positioned close to Dmrt3 but farther from Dmrt1 on the same chromosome, while Dmrt2b is not. Sequence analysis revealed several mutations; insertions, and deletions (indels) on Dmrt2a non-coding regions and single-nucleotide polymorphisms (SNPs) on the Dmrt2a transcript. These indels and SNPs are sex-linked and showed high male heterogeneity but do not affect gene translation. The markers designed to span the mutation sites tested on four different populations showed varied concordance with the genetic sexes. Dmrt2a is transcribed solely in the gonads and gills, while Dmrt2b exists in the gonads, hypothalamus, gills, heart, and spleen. The Dmrt2a and Dmrt2b transcripts are profoundly expressed in the male gonads. Analyses of the transcriptome data from five other fish species (Hainan medaka (Oryzias curvinotus), silver sillago (Sillago sihama), Nile tilapia (Oreochromis niloticus), Hong Kong catfish (Clarias fuscus), and spot-fin porcupine fish (Diodon hystrix)) revealed testes-biased expression of Dmrt1 in all, similar to spotted scat. Additionally, the expression of Dmrt2a is higher in the testes than the ovaries in spotted scat and Hainan medaka. The Dmrt2a transcript was not altered in the coding regions as found in Dmrt1 and Dmrt3 in spotted scat. This could be due to the functional importance of Dmrt2a in development. Another possibility is that because Dmrt2a is positioned farther from Dmrt1 and the chromosome is still young, meaning it is only a matter of time before it differentiates. This study undeniably will aid in understanding the functional divergence of the sex-linked genes in fish.

Inferring the Nature of Missing Heritability in Human Traits Using Data from the GWAS Catalog

Thousands of genes responsible for many diseases and other common traits in humans have been detected by Genome Wide Association Studies (GWAS) in the last decade. However, candidate causal variants found so far usually explain only a small fraction of the heritability estimated by family data. The most common explanation for this observation is that the missing heritability corresponds to variants, either rare or common, with very small effect, which pass undetected due to a lack of statistical power. We carried out a meta-analysis using data from the NHGRI-EBI GWAS Catalog in order to explore the observed distribution of locus effects for a set of 42 complex traits and to quantify their contribution to narrow-sense heritability. With the data at hand, we were able to predict the expected distribution of locus effects for 16 traits and diseases, their expected contribution to heritability, and the missing number of loci yet to be discovered to fully explain the familial heritability estimates. Our results indicate that, for 6 out of the 16 traits, the additive contribution of a great number of loci is unable to explain the familial (broad-sense) heritability, suggesting that the gap between GWAS and familial estimates of heritability may not ever be closed for these traits. In contrast, for the other 10 traits, the additive contribution of hundreds or thousands of loci yet to be found could potentially explain the familial heritability estimates, if this were the case. Computer simulations are used to illustrate the possible contribution from nonadditive genetic effects to the gap between GWAS and familial estimates of heritability.

Human-Mediated Gene Flow Contributes to Metapopulation Genetic Structure of the Pathogenic Fungus Alternaria alternata from Potato

Metapopulation structure generated by recurrent extinctions and recolonizations plays an important role in the evolution of species but is rarely considered in agricultural systems. In this study, generation and mechanism of metapopulation structure were investigated by microsatellite assaying 725 isolates of Alternaria alternata sampled from potato hosts at 16 locations across China. We found a single major cluster, no isolate-geography associations and no bottlenecks in the A. alternata isolates, suggesting a metapopulation genetic structure of the pathogen. We also found weak isolation-by-distance, lower among than within cropping region population differentiation, concordant moving directions of potato products and net gene flow and the highest gene diversity in the region with the most potato imports. These results indicate that in addition to natural dispersal, human-mediated gene flow also contributes to the generation and dynamics of the metapopulation genetic structure of A. alternata in China. Metapopulation structure increases the adaptive capacity of the plant pathogen as a result of enhanced genetic variation and reduced population fragmentation. Consequently, rigid quarantine regulations may be required to reduce population connectivity and the evolutionary potential of A. alternata and other pathogens with a similar population dynamics for a sustainable plant disease management.

Genomic Instability in Cancer: Teetering on the Limit of Tolerance

Cancer genomic instability contributes to the phenomenon of intratumoral genetic heterogeneity, provides the genetic diversity required for natural selection, and enables the extensive phenotypic diversity that is frequently observed among patients. Genomic instability has previously been associated with poor prognosis. However, we have evidence that for solid tumors of epithelial origin, extreme levels of genomic instability, where more than 75% of the genome is subject to somatic copy number alterations, are associated with a potentially better prognosis compared with intermediate levels under this threshold. This has been observed in clonal subpopulations of larger size, especially when genomic instability is shared among a limited number of clones. We hypothesize that cancers with extreme levels of genomic instability may be teetering on the brink of a threshold where so much of their genome is adversely altered that cells rarely replicate successfully. Another possibility is that tumors with high levels of genomic instability are more immunogenic than other cancers with a less extensive burden of genetic aberrations. Regardless of the exact mechanism, but hinging on our ability to quantify how a tumor's burden of genetic aberrations is distributed among coexisting clones, genomic instability has important therapeutic implications. Herein, we explore the possibility that a high genomic instability could be the basis for a tumor's sensitivity to DNA-damaging therapies. We primarily focus on studies of epithelial-derived solid tumors. Cancer Res; 77(9); 2179-85. ©2017 AACR.

Effect of mismatch repair on the mutation rate of bacteriophage ϕX174

Viral mutation rates vary widely in nature, yet the mechanistic and evolutionary determinants of this variability remain unclear. Small DNA viruses mutate orders of magnitude faster than their hosts despite using host-encoded polymerases for replication, which suggests these viruses may avoid post-replicative repair. Supporting this, the genome of bacteriophage ϕX174 is completely devoid of GATC sequence motifs, which are required for methyl-directed mismatch repair in Escherichia coli. Here, we show that restoration of the randomly expected number of GATC sites leads to an eightfold reduction in the rate of spontaneous mutation of the phage, without severely impairing its replicative capacity over the short term. However, the efficacy of mismatch repair in the presence of GATC sites is limited by inefficient methylation of the viral DNA. Therefore, both GATC avoidance and DNA under-methylation elevate the mutation rate of the phage relative to that of the host. We also found that the effects of GATC sites on the phage mutation rate vary extensively depending on their specific location within the phage genome. Finally, the mutation rate reduction afforded by GATC sites is fully reverted under stress conditions, which up-regulate repair pathways and expression of error-prone host polymerases such as heat and treatment with the base analog 5-fluorouracil, suggesting that access to repair renders the phage sensitive to stress-induced mutagenesis.

Delayed lysis confers resistance to the nucleoside analogue 5-fluorouracil and alleviates mutation accumulation in the single-stranded DNA bacteriophage ϕX174

Rates of spontaneous mutation determine viral fitness and adaptability. In RNA viruses, treatment with mutagenic nucleoside analogues selects for polymerase variants with increased fidelity, showing that viral mutation rates can be adjusted in response to imposed selective pressures. However, this type of resistance is not possible in viruses that do not encode their own polymerases, such as single-stranded DNA viruses. We previously showed that serial passaging of bacteriophage ϕX174 in the presence of the nucleoside analogue 5-fluorouracil (5-FU) favored substitutions in the lysis protein E (P. Domingo-Calap, M. Pereira-Gomez, and R. Sanjuán, J. Virol. 86:: 9640-9646, 2012, doi:10.1128/JVI.00613-12). Here, we found that approximately half (6/12) of the amino acid replacements in the N-terminal region of this protein led to delayed lysis, and two of these changes (V2A and D8A) also conferred partial resistance to 5-FU. By delaying lysis, the V2A and D8A substitutions allowed the virus to increase the burst size per cell in the presence of 5-FU. Furthermore, these substitutions tended to alleviate drug-induced mutagenesis by reducing the number of rounds of copying required for population growth, revealing a new mechanism of resistance. This form of mutation rate regulation may also be utilized by other viruses whose replication mode is similar to that of bacteriophage ϕX174.

IMPORTANCE

Many viruses display high rates of spontaneous mutations due to defects in proofreading or postreplicative repair, allowing them to rapidly adapt to changing environments. Viral mutation rates may have been optimized to achieve high adaptability without incurring an excessive genetic load. Supporting this, RNA viruses subjected to chemical mutagenesis treatments have been shown to evolve higher-fidelity polymerases. However, many viruses cannot modulate replication fidelity because they do not encode their own polymerase. Here, we show a new mechanism for regulating viral mutation rates. We found that, under mutagenic conditions, the single-stranded bacteriophage ϕX174 evolved delayed lysis, and that this allowed the virus to increase the amount of progeny produced per cell. As a result, the viral population was amplified in fewer infection cycles, reducing the chances for mutation appearance.

The advantage of arriving first: characteristic times in finite size populations of error-prone replicators

We study the evolution of a finite size population formed by mutationally isolated lineages of error-prone replicators in a two-peak fitness landscape. Computer simulations are performed to gain a stochastic description of the system dynamics. More specifically, for different population sizes, we compute the probability of each lineage being selected in terms of their mutation rates and the amplification factors of the fittest phenotypes. We interpret the results as the compromise between the characteristic time a lineage takes to reach its fittest phenotype by crossing the neutral valley and the selective value of the sequences that form the lineages. A main conclusion is drawn: for finite population sizes, the survival probability of the lineage that arrives first to the fittest phenotype rises significantly.

Correlation between mutation rate and genome size in riboviruses: mutation rate of bacteriophage Qβ

Genome sizes and mutation rates covary across all domains of life. In unicellular organisms and DNA viruses, they show an inverse relationship known as Drake’s rule. However, it is still unclear whether a similar relationship exists between genome sizes and mutation rates in RNA genomes. Coronaviruses, the RNA viruses with the largest genomes (∼30 kb), encode a proofreading 3′ exonuclease that allows them to increase replication fidelity. However, it is unknown whether, conversely, the RNA viruses with the smallest genomes tend to show particularly high mutation rates. To test this, we measured the mutation rate of bacteriophage Qβ, a 4.2-kb levivirus. Amber reversion-based Luria–Delbrück fluctuation tests combined with mutant sequencing gave an estimate of 1.4 × 10⁻⁴ substitutions per nucleotide per round of copying, the highest mutation rate reported for any virus using this method. This estimate was confirmed using a direct plaque sequencing approach and after reanalysis of previously published estimates for this phage. Comparison with other riboviruses (all RNA viruses except retroviruses) provided statistical support for a negative correlation between mutation rates and genome sizes. We suggest that the mutation rates of RNA viruses might be optimized for maximal adaptability and that the value of this optimum may in turn depend inversely on genome size.

Viral mutation rates: modelling the roles of within-host viral dynamics and the trade-off between replication fidelity and speed

Many viruses, particularly RNA viruses, mutate at a very high rate per genome per replication. One possible explanation is that high mutation rates are selected to meet the challenge of fluctuating environments, including the host immune response. Alternatively, recent studies argue that viruses evolve under a trade-off between replication speed and fidelity such that fast replication is selected, and, along with it, high mutation rates. Here, in addition to these factors, we consider the role of viral life-history properties: namely, the within-host dynamics of viruses resulting from their interaction with the host. We develop mathematical models incorporating factors occurring within and between hosts, including deleterious and advantageous mutations, host death owing to virulence and clearance of viruses by the host. Beneficial mutations confer both a within-host and a transmission advantage. First, we find that advantageous mutations have only a weak effect on the optimal genomic mutation rate. Second, viral life-history properties have a large effect on the mutation rate. Third, when the speed-fidelity trade-off is included, there can be two locally optimal mutation rates. Our analysis provides a way to consider how life-history properties combine with biochemical trade-offs to shape mutation rates.

The influence of deleterious mutations on adaptation in asexual populations

We study the dynamics of adaptation in asexual populations that undergo both beneficial and deleterious mutations. In particular, how the deleterious mutations affect the fixation of beneficial mutations was investigated. Using extensive Monte Carlo simulations, we find that in the “strong-selection weak mutation (SSWM)” regime or in the “clonal interference (CI)” regime, deleterious mutations rarely influence the distribution of “selection coefficients of the fixed mutations (SCFM)”; while in the “multiple mutations” regime, the accumulation of deleterious mutations would lead to a decrease in fitness significantly. We conclude that the effects of deleterious mutations on adaptation depend largely on the supply of beneficial mutations. And interestingly, the lowest adaptation rate occurs for a moderate value of selection coefficient of deleterious mutations.