Life history and deleterious mutation rate coevolution

The cost of germline maintenance gives rise to a trade-off between lowering the deleterious mutation rate and investing in life history functions. Therefore, life history and the mutation rate coevolve, but this coevolution is not well understood. We develop a mathematical model to analyse the evolution of resource allocation traits, which simultaneously affect life history and the deleterious mutation rate. First, we show that the invasion fitness of such resource allocation traits can be approximated by the basic reproductive number of the least-loaded class; the expected lifetime production of offspring without deleterious mutations born to individuals without deleterious mutations. Second, we apply the model to investigate (i) the coevolution of reproductive effort and germline maintenance and (ii) the coevolution of age-at-maturity and germline maintenance. This analysis provides two resource allocation predictions when exposure to environmental mutagens is higher. First, selection favours higher allocation to germline maintenance, even if it comes at the expense of life history functions, and leads to a shift in allocation towards reproduction rather than survival. Second, life histories tend to be faster, characterised by individuals with shorter lifespans and smaller body sizes at maturity. Our results suggest that mutation accumulation via the cost of germline maintenance can be a major force shaping life-history traits.

de novo structural variations of Escherichia coli detected by Nanopore long-read sequencing

Spontaneous mutations power evolution, while large-scale structural variations (SVs) remain poorly studied, primarily because of the lack of long-read sequencing techniques and powerful analytical tools. Here, we explore the SVs of Escherichia coli by running 67 wild-type (WT) and 37 MMR-deficient (ΔmutS) mutation accumulation lines, each experiencing more than 4000 cell divisions, by applying Nanopore long-read sequencing and Illumina PE150 sequencing, and verifying the results by Sanger sequencing. In addition to precisely repeating previous mutation rates of base-pair substitutions and indels mutation rates, we do find significant improvement in insertion and deletion detection using long-read sequencing. The long-read sequencing and corresponding softwares can particularly detect bacterial SVs in both simulated and real datasets with high accuracy. These lead to SV rates of 2.77 × 10-4 (WT) and 5.26 × 10-4 (MMR-deficient) per cell division per genome, which is comparable with previous reports. This study provides the SV rates of E. coli by applying long-read sequencing and SV-detection programs, revealing a broader and more accurate picture of spontaneous mutations in bacteria.

Sex without crossing over in the yeast Saccharomycodes ludwigii

BACKGROUND

Intermixing of genomes through meiotic reassortment and recombination of homologous chromosomes is a unifying theme of sexual reproduction in eukaryotic organisms and is considered crucial for their adaptive evolution. Previous studies of the budding yeast species Saccharomycodes ludwigii suggested that meiotic crossing over might be absent from its sexual life cycle, which is predominated by fertilization within the meiotic tetrad.

RESULTS

We demonstrate that recombination is extremely suppressed during meiosis in Sd. ludwigii. DNA double-strand break formation by the conserved transesterase Spo11, processing and repair involving interhomolog interactions are required for normal meiosis but do not lead to crossing over. Although the species has retained an intact meiotic gene repertoire, genetic and population analyses suggest the exceptionally rare occurrence of meiotic crossovers in its genome. A strong AT bias of spontaneous mutations and the absence of recombination are likely responsible for its unusually low genomic GC level.

CONCLUSIONS

Sd. ludwigii has followed a unique evolutionary trajectory that possibly derives fitness benefits from the combination of frequent mating between products of the same meiotic event with the extreme suppression of meiotic recombination. This life style ensures preservation of heterozygosity throughout its genome and may enable the species to adapt to its environment and survive with only minimal levels of rare meiotic recombination. We propose Sd. ludwigii as an excellent natural forum for the study of genome evolution and recombination rates.

Mutational property of newly identified mutagen l-glutamic acid γ-hydrazide in Escherichia coli

We previously found that an l-glutamine analog l-glutamic acid γ-hydrazide has high mutagenic activity through the high-throughput laboratory evolution of Escherichia coli. In this study, mutagenicity and mutational property of l-glutamic acid γ-hydrazide were examined by the Ames test and mutation accumulation experiments using E. coli. The Ames test revealed that l-glutamic acid γ-hydrazide showed higher mutagenic activity without metabolic activation than known mutagens 2-aminoanthracene, and cobalt(II) acetate tetrahydrate. This result indicates that l-glutamic acid γ-hydrazide does not require metabolic activation for mutagenic activity in E. coli. Mutation accumulation experiments and whole-genome sequencing analysis revealed the number and spectrum of the accumulated mutations with or without l-glutamic acid γ-hydrazide. In the presence of l-glutamic acid γ-hydrazide, MDS42 strain accumulated 392.3 ± 116.2 point mutations during 30 passages corresponding to 777 generations, while MDS42 strain accumulated 1.5 ± 2.5 point mutations without l-glutamic acid γ-hydrazide during 50 passages corresponding to 1341 generations. The mutational spectrum of l-glutamic acid γ-hydrazide was G/C to A/T transition (82.2 ± 4.3 %) and A/T to G/C transition (17.4 ± 4.3 %). These results indicated that l-glutamic acid γ-hydrazide has a strong mutagenic activity.

Genome-wide DNA mutations in Arabidopsis plants after multigenerational exposure to high temperatures

BACKGROUND

Elevated temperatures can cause physiological, biochemical, and molecular responses in plants that can greatly affect their growth and development. Mutations are the most fundamental force driving biological evolution. However, how long-term elevations in temperature influence the accumulation of mutations in plants remains unknown.

RESULTS

Multigenerational exposure of Arabidopsis MA (mutation accumulation) lines and MA populations to extreme heat and moderate warming results in significantly increased mutation rates in single-nucleotide variants (SNVs) and small indels. We observe distinctive mutational spectra under extreme and moderately elevated temperatures, with significant increases in transition and transversion frequencies. Mutation occurs more frequently in intergenic regions, coding regions, and transposable elements in plants grown under elevated temperatures. At elevated temperatures, more mutations accumulate in genes associated with defense responses, DNA repair, and signaling. Notably, the distribution patterns of mutations among all progeny differ between MA populations and MA lines, suggesting that stronger selection effects occurred in populations. Methylation is observed more frequently at mutation sites, indicating its contribution to the mutation process at elevated temperatures. Mutations occurring within the same genome under elevated temperatures are significantly biased toward low gene density regions, special trinucleotides, tandem repeats, and adjacent simple repeats. Additionally, mutations found in all progeny overlap significantly with genetic variations reported in 1001 Genomes, suggesting non-uniform distribution of de novo mutations through the genome.

CONCLUSION

Collectively, our results suggest that elevated temperatures can accelerate the accumulation, and alter the molecular profiles, of DNA mutations in plants, thus providing significant insight into how environmental temperatures fuel plant evolution.

The insect-killing bacterium Photorhabdus luminescens has the lowest mutation rate among bacteria

Mutation is a primary source of genetic variation that is used to power evolution. Many studies, however, have shown that most mutations are deleterious and, as a result, extremely low mutation rates might be beneficial for survival. Using a mutation accumulation experiment, an unbiased method for mutation study, we found an extremely low base-substitution mutation rate of 5.94 × 10^-11 per nucleotide site per cell division (95% Poisson confidence intervals: 4.65 × 10^-11, 7.48 × 10^-11) and indel mutation rate of 8.25 × 10^-12 per site per cell division (95% confidence intervals: 3.96 × 10^-12, 1.52 × 10^-11) in the bacterium Photorhabdus luminescens ATCC29999. The mutations are strongly A/T-biased with a mutation bias of 10.28 in the A/T direction. It has been hypothesized that the ability for selection to lower mutation rates is inversely proportional to the effective population size (drift-barrier hypothesis) and we found that the effective population size of this bacterium is significantly greater than most other bacteria. This finding further decreases the lower-bounds of bacterial mutation rates and provides evidence that extreme levels of replication fidelity can evolve within organisms that maintain large effective population sizes.

Streptomycin and nalidixic acid elevate the spontaneous genome-wide mutation rate in Escherichia coli

Since antibiotic resistance is a growing public health problem worldwide, it is important to understand how antibiotics and spontaneous mutations cooperate and shape the genome-wide mutation rate and spectrum. Here, we quantitatively evaluate genome-wide mutational profiles of Escherichia coli after long-term subinhibitory exposure to a broad-spectrum (streptomycin) and a narrow-spectrum antibiotic (nalidixic acid), using a mutation accumulation design combined with whole-genome resequencing of replicate lines as a mutagenicity test. We determined that, while the genome-wide mutation rate is slightly higher in the streptomycin-treated lines compared to the control lines, there is a significant increase in the nalidixic acid-treated lines. Our findings suggest that both broad and narrow-spectrum antibiotics may elevate the mutation rates in E. coli, but mechanisms of action may affect the consequence, thus contribute to accelerating the rate of adaptation and conferring antibiotic resistance.

Metal exposure causes rDNA copy number to fluctuate in mutation accumulation lines of Daphnia pulex

Ribosomal (r)DNA is a highly dynamic, conserved, multigene family whose sequence homogeneity is thought to be maintained by intra- and interchromosomal recombination, which are capable of changing rDNA copy number. It is generally not known how environmental stress such as sublethal exposure to environmentally relevant concentrations of metals impacts rDNA copy number. To determine how chronic metal exposure affects rDNA, we measured copy number of the 18S rRNA gene in 355 copper and nickel-exposed samples and 132 metal-free samples derived from 325 mutation accumulation (MA) lines of two genetically distinct Daphnia pulex lineages. The MA lines were sampled at four time points over 100+ generations of clonal propagation. The copy number of rDNA was also measured in 15 individuals sampled from a metal-free non-MA control population established from the same progenitor as one of the MA lineages. We found that mean rDNA copy number fluctuated across lines exposed to metals with a tendency to decrease over time. In contrast, mean rDNA copy number in the metal-free control lines and the non-MA population remained stable over time. It is generally accepted that extreme rDNA loss results in the loss of organism fitness. Thus, fluctuations in rDNA copy number, including losses, could affect the long-term viability of natural populations of Daphnia in metal-contaminated habitats.

New complexities of SOS-induced "untargeted" mutagenesis in Escherichia coli as revealed by mutation accumulation and whole-genome sequencing

When its DNA is damaged, Escherichia coli induces the SOS response, which consists of about 40 genes that encode activities to repair or tolerate the damage. Certain alleles of the major SOS-control genes, recA and lexA, cause constitutive expression of the response, resulting in an increase in spontaneous mutations. These mutations, historically called "untargeted", have been the subject of many previous studies. Here we re-examine SOS-induced mutagenesis using mutation accumulation followed by whole-genome sequencing (MA/WGS), which allows a detailed picture of the types of mutations induced as well as their sequence-specificity. Our results confirm previous findings that SOS expression specifically induces transversion base-pair substitutions, with rates averaging about 60-fold above wild-type levels. Surprisingly, the rates of G:C to C:G transversions, normally an extremely rare mutation, were induced an average of 160-fold above wild-type levels. The SOS-induced transversion showed strong sequence specificity, the most extreme of which was the G:C to C:G transversions, 60% of which occurred at the middle base of 5'GGC3'+5'GCC3' sites, although these sites represent only 8% of the G:C base pairs in the genome. SOS-induced transversions were also DNA strand-biased, occurring, on average, 2- to 4- times more often when the purine was on the leading-strand template and the pyrimidine on the lagging-strand template than in the opposite orientation. However, the strand bias was also sequence specific, and even of reverse orientation at some sites. By eliminating constraints on the mutations that can be recovered, the MA/WGS protocol revealed new complexities of SOS "untargeted" mutations.

Revamping the evolutionary theories of aging

Radical lifespan disparities exist in the animal kingdom. While the ocean quahog can survive for half a millennium, the mayfly survives for less than 48 h. The evolutionary theories of aging seek to explain why such stark longevity differences exist and why a deleterious process like aging evolved. The classical mutation accumulation, antagonistic pleiotropy, and disposable soma theories predict that increased extrinsic mortality should select for the evolution of shorter lifespans and vice versa. Most experimental and comparative field studies conform to this prediction. Indeed, animals with extreme longevity (e.g., Greenland shark, bowhead whale, giant tortoise, vestimentiferan tubeworms) typically experience minimal predation. However, data from guppies, nematodes, and computational models show that increased extrinsic mortality can sometimes lead to longer evolved lifespans. The existence of theoretically immortal animals that experience extrinsic mortality - like planarian flatworms, panther worms, and hydra - further challenges classical assumptions. Octopuses pose another puzzle by exhibiting short lifespans and an uncanny intelligence, the latter of which is often associated with longevity and reduced extrinsic mortality. The evolutionary response to extrinsic mortality is likely dependent on multiple interacting factors in the organism, population, and ecology, including food availability, population density, reproductive cost, age-mortality interactions, and the mortality source.

A little bit of sex prevents mutation accumulation even in apomictic polyploid plants

Background

In the absence of sex and recombination, genomes are expected to accumulate deleterious mutations via an irreversible process known as Muller’s ratchet, especially in the case of polyploidy. In contrast, no genome-wide mutation accumulation was detected in a transcriptome of facultative apomictic, hexaploid plants of the Ranunculus auricomus complex. We hypothesize that mutations cannot accumulate in flowering plants with facultative sexuality because sexual and asexual development concurrently occurs within the same generation. We assume a strong effect of purging selection on reduced gametophytes in the sexual developmental pathway because previously masked recessive deleterious mutations would be exposed to selection.

Results

We test this hypothesis by modeling mutation elimination using apomictic hexaploid plants of the R. auricomus complex. To estimate mean recombination rates, the mean number of recombinants per generation was calculated by genotyping three F1 progeny arrays with six microsatellite markers and character incompatibility analyses. We estimated the strength of purging selection in gametophytes by calculating abortion rates of sexual versus apomictic development at the female gametophyte, seed and offspring stage. Accordingly, we applied three selection coefficients by considering effects of purging selection against mutations on (1) male and female gametophytes in the sexual pathway (additive, s = 1.000), (2) female gametophytes only (s = 0.520), and (3) on adult plants only (sporophytes, s = 0.212). We implemented recombination rates into a mathematical model considering the three different selection coefficients, and a genomic mutation rate calculated from genome size of our plants and plant-specific mutation rates. We revealed a mean of 6.05% recombinants per generation. This recombination rate eliminates mutations after 138, 204 or 246 generations, depending on the respective selection coefficients (s = 1.000, 0.520, and 0.212).

Conclusions

Our results confirm that the empirically observed frequencies of facultative recombination suffice to prevent accumulation of deleterious mutations via Muller’s ratchet even in a polyploid genome. The efficiency of selection is in flowering plants strongly increased by acting on the haplontic (reduced) gametophyte stage.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1495-z) contains supplementary material, which is available to authorized users.

Running out of developmental program and selfish anti-aging: a new hypothesis explaining the aging process in primates

Of the three complementary theories of aging, two (mutation accumulation and antagonistic pleiotropy) were formulated over fifty years ago before the introduction of molecular biology, and the third (disposable soma) is over thirty years old. Despite rigorous research in the past fifty years, none have gained substantial experimental support. Here, I review these theories and introduce a new hypothesis called the selfish anti-aging (SAA). Aging happens because natural selection is indifferent to the organism's life beyond reproduction; however, many mammalian species acquired anti-aging genes, which are providing instructions following completion of developmental, ontogeny, program. Such instructor-genes might be responsible for the elongation of lifespans of primates as a byproduct of parental care program. According to the SAA hypothesis, the animal models used in aging research could be divided into three groups, based on the degree of perceived presence and action of instructor-genes in each group. This new hypothesis is grounded in evolutionary theory and describes the unique primate aging process.

Rapid cell division of Staphylococcus aureus during colonization of the human nose

BACKGROUND

Staphylococcus aureus is an important opportunistic pathogen and a commensal bacterium, thriving in the nasal cavities of 20% of the human population. Little is known about the dynamics of asymptomatic colonization and the occasional transition to infectious disease.

RESULTS

In this study, we inferred that S. aureus cells replicate every one to three hours on average while colonizing the human nose, based on two independent lines of genomic evidence. First, we collected nasal swab samples from human subjects, extracted and sequenced metagenomic DNA, and analyzed the distribution of sequencing coverage along the staphylococcal chromosome. Calibration of this data by comparison to a laboratory culture enabled measuring S. aureus cell division rates in nasal samples. Second, we applied mutation accumulation experiments paired with genome sequencing to measure spontaneous mutation rates at a genome scale. Relating these mutation rates to annual evolutionary rates confirmed that nasal S. aureus continuously pass several thousand cell divisions per year when averaged over large, globally distributed populations and over many years, corresponding to generation times of less than two hours.

CONCLUSIONS

The cell division rates we determined were higher than the fastest documented rates during fulminant disease progression (in a mouse model of systemic infection) and much higher than those previously measured in expectorated sputum from cystic fibrosis patients. This paper supplies absolute in-vivo generation times for an important bacterial commensal, indicating that colonization of the human upper respiratory tract is characterized by a highly dynamic equilibrium between bacterial growth and removal.

Cooperation between non-essential DNA polymerases contributes to genome stability in Saccharomyces cerevisiae

DNA polymerases influence genome stability through their involvement in DNA replication, response to DNA damage, and DNA repair processes. Saccharomyces cerevisiae possess four non-essential DNA polymerases, Pol λ, Pol η, Pol ζ, and Rev1, which have varying roles in genome stability. In order to assess the contribution of the non-essential DNA polymerases in genome stability, we analyzed the pol4Δ rev1Δ rev3Δ rad30Δ quadruple mutant in microhomology mediated repair, due to recent studies linking some of these DNA polymerases to this repair pathway. Our results suggest that the length and quality of microhomology influence both the overall efficiency of repair and the involvement of DNA polymerases. Furthermore, the non-essential DNA polymerases demonstrate overlapping and redundant functions when repairing double-strand breaks using short microhomologies containing mismatches. Then, we examined genome-wide mutation accumulation in the pol4Δ rev1Δ rev3Δ rad30Δ quadruple mutant compared to wild type cells. We found a significant decrease in the overall rate of mutation accumulation in the quadruple mutant cells compared to wildtype, but an increase in frameshift mutations and a shift towards transversion base-substitution with a preference for G:C to T:A or C:G. Thus, the non-essential DNA polymerases have an impact on the nature of the mutational spectrum. The sequence and functional homology shared between human and S. cerevisiae non-essential DNA polymerases suggest these DNA polymerases may have a similar role in human cells.

Temperature responses of mutation rate and mutational spectrum in an Escherichia coli strain and the correlation with metabolic rate

Background

Temperature is a major determinant of spontaneous mutation, but the precise mode, and the underlying mechanisms, of the temperature influences remain less clear. Here we used a mutation accumulation approach combined with whole-genome sequencing to investigate the temperature dependence of spontaneous mutation in an Escherichia coli strain. Experiments were performed under aerobic conditions at 25, 28 and 37 °C, three temperatures that were non-stressful for the bacterium but caused significantly different bacterial growth rates.

Results

Mutation rate did not differ between 25 and 28 °C, but was higher at 37 °C. Detailed analyses of the molecular spectrum of mutations were performed; and a particularly interesting finding is that higher temperature led to a bias of mutation to coding, relative to noncoding, DNA. Furthermore, the temperature response of mutation rate was extremely similar to that of metabolic rate, consistent with an idea that metabolic rate predicts mutation rate.

Conclusions

Temperature affects mutation rate and the types of mutation supply, both being crucial for the opportunity of natural selection. Our results help understand how temperature drives evolutionary speed of organisms and thus the global patterns of biodiversity. This study also lend support to the metabolic theory of ecology for linking metabolic rate and molecular evolution rate.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1252-8) contains supplementary material, which is available to authorized users.

Insertion sequence elements-mediated structural variations in bacterial genomes

Mobile genetic elements (MGEs) impact the evolution and stability of their host genomes. Insertion sequence (IS) elements are the most common MGEs in bacterial genomes and play a crucial role in mediating large-scale variations in bacterial genomes. It is understood that IS elements and MGEs in general coexist in a dynamical equilibrium with their respective hosts. Current studies indicate that the spontaneous movement of IS elements does not follow a constant rate in different bacterial genomes. However, due to the paucity and sparsity of the data, these observations are yet to be conclusive. In this paper, we conducted a comparative analysis of the IS-mediated genome structural variations in ten mutation accumulation (MA) experiments across eight strains of five bacterial species containing IS elements, including four strains of the E. coli. We used GRASPER algorithm, a denovo structural variation (SV) identification algorithm designed to detect SVs involving repetitive sequences in the genome. We observed highly diverse rates of IS insertions and IS-mediated recombinations across different bacterial species as well as across different strains of the same bacterial species. We also observed different rates of the elements from the same IS family in different bacterial genomes, suggesting that the distinction in rates might not be due to the different composition of IS elements across bacterial genomes.

Increase in viability due to the accumulation of X chromosome mutations in Drosophila melanogaster males

To increase our understanding of the role of new X-chromosome mutations in adaptive evolution, single-X Drosophila melanogaster males were mated with attached-X chromosome females, allowing the male X chromosome to accumulate mutations over 28 generations. Contrary to our hypothesis that male viability would decrease over time, due to the accumulation and expression of X-linked recessive deleterious mutations in hemizygous males, viability significantly increased. This increase may be attributed to germinal selection and to new X-linked beneficial or compensatory mutations, possibly supporting the faster-X hypothesis.

Muver, a computational framework for accurately calling accumulated mutations

BACKGROUND

Identification of mutations from next-generation sequencing data typically requires a balance between sensitivity and accuracy. This is particularly true of DNA insertions and deletions (indels), that can impart significant phenotypic consequences on cells but are harder to call than substitution mutations from whole genome mutation accumulation experiments. To overcome these difficulties, we present muver, a computational framework that integrates established bioinformatics tools with novel analytical methods to generate mutation calls with the extremely low false positive rates and high sensitivity required for accurate mutation rate determination and comparison.

RESULTS

Muver uses statistical comparison of ancestral and descendant allelic frequencies to identify variant loci and assigns genotypes with models that include per-sample assessments of sequencing errors by mutation type and repeat context. Muver identifies maximally parsimonious mutation pathways that connect these genotypes, differentiating potential allelic conversion events and delineating ambiguities in mutation location, type, and size. Benchmarking with a human gold standard father-son pair demonstrates muver's sensitivity and low false positive rates. In DNA mismatch repair (MMR) deficient Saccharomyces cerevisiae, muver detects multi-base deletions in homopolymers longer than the replicative polymerase footprint at rates greater than predicted for sequential single-base deletions, implying a novel multi-repeat-unit slippage mechanism.

CONCLUSIONS

Benchmarking results demonstrate the high accuracy and sensitivity achieved with muver, particularly for indels, relative to available tools. Applied to an MMR-deficient Saccharomyces cerevisiae system, muver mutation calls facilitate mechanistic insights into DNA replication fidelity.

Population-based statistical inference for temporal sequence of somatic mutations in cancer genomes

BACKGROUND

It is well recognized that accumulation of somatic mutations in cancer genomes plays a role in carcinogenesis; however, the temporal sequence and evolutionary relationship of somatic mutations remain largely unknown.

METHODS

In this study, we built a population-based statistical framework to infer the temporal sequence of acquisition of somatic mutations. Using the model, we analyzed the mutation profiles of 1954 tumor specimens across eight tumor types.

RESULTS

As a result, we identified tumor type-specific directed networks composed of 2-15 cancer-related genes (nodes) and their mutational orders (edges). The most common ancestors identified in pairwise comparison of somatic mutations were TP53 mutations in breast, head/neck, and lung cancers. The known relationship of KRAS to TP53 mutations in colorectal cancers was identified, as well as potential ancestors of TP53 mutation such as NOTCH1, EGFR, and PTEN mutations in head/neck, lung and endometrial cancers, respectively. We also identified apoptosis-related genes enriched with ancestor mutations in lung cancers and a relationship between APC hotspot mutations and TP53 mutations in colorectal cancers.

CONCLUSION

While evolutionary analysis of cancers has focused on clonal versus subclonal mutations identified in individual genomes, our analysis aims to further discriminate ancestor versus descendant mutations in population-scale mutation profiles that may help select cancer drivers with clinical relevance.

Mutation independently affects reproductive traits and dauer larvae development in mutation accumulation lines of Caenorhabditis elegans

Developmental decisions are important in organismal fitness. For the nematode Caenorhabditis elegans, which is naturally found in the ephemeral food patches formed by rotting plant material, correctly committing to dauer or non-dauer larval development is key to genotype survival. To investigate the link between reproductive traits, which will determine how populations grow, and dauer larvae formation, we have analysed these traits in mutation accumulation lines of C. elegans. We find that reproductive traits of individual worms-the total number of progeny and the timing of progeny production-are highly correlated with the population size observed in growing populations. In contrast, we find no relationship between reproduction traits and the number of dauer larvae observed in growing populations. We also do not observe a mutational bias in dauer larvae formation. These results indicate that the control of dauer larvae formation is distinct from the control of reproduction and that differences in dauer larvae formation can evolve rapidly.

The role of cell location and spatial gradients in the evolutionary dynamics of colon and intestinal crypts

BACKGROUND

Colon and intestinal crypts serve as an important model system for adult stem cell proliferation and differentiation. We develop a spatial stochastic model to study the rate of somatic evolution in a normal crypt, focusing on the production of two-hit mutants that inactivate a tumor suppressor gene. We investigate the effect of cell division pattern along the crypt on mutant production, assuming that the division rate of each cell depends on its location.

RESULTS

We find that higher probability of division at the bottom of the crypt, where the stem cells are located, leads to a higher rate of double-hit mutant production. The optimal case for delaying mutations occurs when most of the cell divisions happen at the top of the crypt. We further consider an optimization problem where the "evolutionary" penalty for double-hit mutant generation is complemented with a "functional" penalty that assures that fully differentiated cells at the top of the crypt cannot divide.

CONCLUSION

The trade-off between the two types of objectives leads to the selection of an intermediate division pattern, where the cells in the middle of the crypt divide with the highest rate. This matches the pattern of cell divisions obtained experimentally in murine crypts.

REVIEWERS

This article was reviewed by David Axelrod (nominated by an Editorial Board member, Marek Kimmel), Yang Kuang and Anna Marciniak-Czochra. For the full reviews, please go to the Reviewers' comments section.

The repeatability of genome-wide mutation rate and spectrum estimates

Over the last decade, mutation studies have grown in popularity due to the affordability and accessibility of whole genome sequencing. As the number of species in which spontaneous mutation has been directly estimated approaches 20 across two domains of life, questions arise over the repeatability of results in such experiments. Five species were identified in which duplicate mutation studies have been performed. Across these studies the difference in estimated spontaneous mutation rate is at most, weakly significant (p < 0.01). However, a highly significant (p < 10(-5)), threefold difference in the rate of insertions/deletions (indels) exists between two recent studies in Schizosaccharomyces pombe. Upon investigation of the ancestral genome sequence for both studies, a possible anti-mutator allele was identified. The observed variation in indel rate may imply that the use of indel markers, such as microsatellites, for the investigation of genetic diversity within and among populations may be inappropriate because of the assumption of uniform mutation rate within a species.