Gene Importance Fails to Predict Mutation Rates in Major Microbes
Biology

Gene Importance Fails to Predict Mutation Rates in Major Microbes

Rigorous new genomic re-analysis challenges adaptationist theories and exposes widespread statistical artefacts in mutation research.

By Priya Sharma
Published:
Email this Article
Gene Importance Fails to Predict Mutation Rates in Major Microbes
Mutation accumulation data merging and mutation profiles. (A) Breakdown of mutation accumulation data combined from WT strains generated by Foster and colleagues [64] and Wei and colleagues [65]. Mutations are assembled as those occurring in the whole-genome sequence and then filtered to include only those in the coding sequence (CDS). The proportion of CDS mutations contributed by each data set to the full CDS set are Foster:Wei 0.65:0.35. For Foster data, mutations are tracked back to E. coli genome NC_000913.2 (with 4321 genes and 4073876 total CDS bp), for Wei data NC_000913.3 (with 4308 genes and 4005730 total CDS bp). See S1 Table for mutation annotations per gene. (B–D) Distribution of trinucleotide mutation differences for each pairwise codon site comparison. Shapiro-Wilk normality test results displayed on the plot, as well as those of a paired Wilcoxon Signed-Rank test. Trinucleotides were formed such that the middle focal base is the mutating one and can occur at within-codo Sofia Radrizzani et al. / PLoS biology

Eukaryotic and bacterial genomes display a surprising lack of targeted mutation shielding, according to a sweeping re-analysis of genetic data that challenges prevailing adaptationist theories.

For decades, evolutionary biologists have debated whether natural selection actively shapes mutation rates across different genes within an organism. Prominent hypotheses have argued that vital genes, such as those essential for survival or responsible for high levels of protein production, are shielded from mutations by being positioned on specific DNA strands or genomic coordinates. Conversely, genes linked to rapid adaptation have been theorized to possess elevated mutation rates to help organisms weather environmental stress.

However, a comprehensive investigation conducted by Sofia Radrizzani and Laurence D. Hurst from the University of Bath upends these assumptions. Their exhaustive study, published in PLOS Biology, demonstrates that a gene’s functional importance does not predict its mutation rate.

Methodological Pitfalls and the Gene Length Illusion

To uncover the true drivers of mutational variation, the researchers examined amassed mutation accumulation data from wild-type Escherichia coli and the budding yeast Saccharomyces cerevisiae. Previous studies supporting adaptive mutation models often relied on proxies like single-nucleotide polymorphisms or data from mismatch repair knockout strains, which are laboratory strains lacking key DNA spelling-error correction mechanisms.

Fig 2: Average mutation density per sequence feature. Mutation density determined as mutation count per gene normalised by gene length. Features include gene length (A), GC content (B), deviation of CpG content from that expected by mononucleotide content (C), relative position from origin (D), protein abundance from PaxDB [78] (E), and median transcript abundance from a compendium of RNA-seq data [79] (F), strand orientation relative to the origin of replication (G, scored as 0 for lagging (−) strand, 1 for leading (+) strand), gene essentiality classification from the PEC database [77] (H, scored as 2 for non-essential and 1 for essential). For all panels, genes binned according to ranging sequence feature in order to produce similarly sized bins (number of genes in each bin is reported on top of bars). Error bars represent standard error of the mean (SEM). For panels G and H, a Mann–Whitney U test was performed and the U statistic with respective P-value displayed in the plot. For panels A
Average mutation density per sequence feature. Mutation density determined as mutation count per gene normalised by gene length. Features include gene length (A), GC content (B), deviation of CpG content from that expected by mononucleotide content (C), relative position from origin (D), protein abundance from PaxDB [78] (E), and median transcript abundance from a compendium of RNA-seq data [79] (F), strand orientation relative to the origin of replication (G, scored as 0 for lagging (−) strand, 1 for leading (+) strand), gene essentiality classification from the PEC database [77] (H, scored as 2 for non-essential and 1 for essential). For all panels, genes binned according to ranging sequence feature in order to produce similarly sized bins (number of genes in each bin is reported on top of bars). Error bars represent standard error of the mean (SEM). For panels G and H, a Mann–Whitney U test was performed and the U statistic with respective P-value displayed in the plot. For panels A Credit: Sofia Radrizzani et al. / PLoS biology

Radrizzani and Hurst discovered that these older datasets frequently introduced statistical artefacts. Most notably, standard correlation methods previously indicated a strong positive relationship between gene length and mutational density, meaning longer genes appeared to accumulate more mutations per base pair. Through rigorous Monte Carlo simulations, the authors proved this relationship is an artifactual byproduct driven entirely by the high abundance of short genes that happen to harbor zero mutations.

When correcting for these statistical anomalies and controlling for genomic covariates such as GC content and replication timing, the alluring signals of adaptive mutation vanished entirely.

Fig 3: Lack of effect of gene size on mutation density. A. Mutation density by gene length for observed and simulated data (partial Kendall correlations). The partial tau values found for each of the 10,000 simulated mutation densities are shown in the blue histograms and their mean by the black line, whilst the value from the observed mutation density by the orange line. The Z score displayed in top right represents the deviation in observed to the simulated null calculated by subtracting the mean simulated value from the observed value and dividing it by the standard deviation of the simulated values. B. Average mutation density per gene length controlling for total gene sequence. To assess mutational density across gene sizes, genes were ranked by length and divided into bins of approximately constant total sequence length such that for a total sequence of N bases and M bins, each bin contained approximately N/M bases, resulting in fewer genes in bins with larger gene sizes. Mutational den
Lack of effect of gene size on mutation density. A. Mutation density by gene length for observed and simulated data (partial Kendall correlations). The partial tau values found for each of the 10,000 simulated mutation densities are shown in the blue histograms and their mean by the black line, whilst the value from the observed mutation density by the orange line. The Z score displayed in top right represents the deviation in observed to the simulated null calculated by subtracting the mean simulated value from the observed value and dividing it by the standard deviation of the simulated values. B. Average mutation density per gene length controlling for total gene sequence. To assess mutational density across gene sizes, genes were ranked by length and divided into bins of approximately constant total sequence length such that for a total sequence of N bases and M bins, each bin contained approximately N/M bases, resulting in fewer genes in bins with larger gene sizes. Mutational den Credit: Sofia Radrizzani et al. / PLoS biology

What Actually Drives Mutational Density Across Genomes

Across both bacterial and yeast genomes, the researchers found no evidence that essential genes or genes encoding abundant proteins experience lower mutation rates. Furthermore, analysis of bacterial strand orientation relative to the origin of replication revealed no systematic differences in mutation density between leading and lagging DNA strands.

Instead, the study highlights that transcription itself exerts a minor antimutagenic effect in E. coli, where actively transcribed regions display slightly lower mutation densities. This phenomenon likely stems from transcription-coupled DNA repair mechanisms rather than evolutionary selection protecting important functions. In yeast, replication timing remains the dominant structural predictor of mutation frequency, with later-replicating genomic regions accumulating more mutations.

Ultimately, the researchers concluded that any potential mutation rate variations linked to gene importance are so minuscule that they explain less than half a percent of total genomic variation. These findings strongly reaffirm the classic evolutionary paradigm that spontaneous mutation operates as a fundamentally directionless force.

The research was published in PLoS biology on November 1, 2026.

Fact Checked

This article has been fact checked for accuracy, with information verified against reputable sources. Learn more about us and our editorial process.

Last reviewed on .

Article history

  • Latest version

Reference(s)

  1. Radrizzani, Sofia., et al. “The importance of genes does not predict their mutation rate.” PLOS Biology, vol. 24, no. 11, November 12, 2026, pp. e3004089 Public Library of Science (PLoS), doi: 10.1371/journal.pbio.3004089. <https://doi.org/10.1371/journal.pbio.3004089>.

Cite this page:

Sharma, Priya. “Gene Importance Fails to Predict Mutation Rates in Major Microbes.” BioScience. BioScience ISSN 2521-5760, 13 November 2026. <https://www.bioscience.com.pk/en/subject/biology/the-importance-of-genes-does-not-predict-their-mutation-rate>. Sharma, P. (2026, November 13). “Gene Importance Fails to Predict Mutation Rates in Major Microbes.” BioScience. ISSN 2521-5760. Retrieved November 13, 2026 from https://www.bioscience.com.pk/en/subject/biology/the-importance-of-genes-does-not-predict-their-mutation-rate Sharma, Priya. “Gene Importance Fails to Predict Mutation Rates in Major Microbes.” BioScience. ISSN 2521-5760. https://www.bioscience.com.pk/en/subject/biology/the-importance-of-genes-does-not-predict-their-mutation-rate (accessed November 13, 2026).
End of the article