vafSOMA: Detecting cell-layer-fixed somatic mutations in plants
- 1Technical University of Munich
Abstract
Somatic mutations are a source of heritable genetic variation in plants and have long been exploited by breeders through selection and vegetative propagation of bud sports. They often arise in specific meristem layers and can become fixed as descendant cell lineages expand during development. Bulk sequencing of plant organs mixes DNA from these layers, so the expected variant allele frequency of a layer-fixed mutation depends on that layer's tissue contribution, which complicates their detection. Here we present vafSOMA, a statistical model that uses tissue composition to detect fixed somatic mutations, assign them to layers and obtain counts for mutation-rate estimation. It accommodates bulk and layer-enriched samples, including those that retain DNA from other layers after enrichment. In simulations, vafSOMA performed close to the best variant allele frequency cutoff chosen in hindsight, without requiring a user-selected cutoff. Its diagnostics flagged layers for which the data did not support separate mutation counts or rates. Analysing replicate organs substantially improved detection of mutations fixed in the branch meristem and yielded estimates close to the simulated stem-cell mutation rate. Application to apricot sequencing data showed that most reference mutations fixed in the dominant layer were recovered from fruit mixtures at high depth, whereas epidermal mutations remained unresolvable in fruit and leaves. These findings suggest that published bulk somatic sequencing studies mainly capture mutations fixed in the dominant tissue layer and may therefore underestimate the combined input of fixed mutations by roughly two- to threefold when layer-specific rates are similar. vafSOMA is available at https://github.com/jlab-code/vafSOMA.
Related
- Age-Related Patterns of Selected Protein-Altering Somatic Mutations and Overall Survival in Colorectal Cancer: A Retrospective Comparative Cohort Study
- Somatic mutations in neurodegeneration: from detection to disease mechanisms
- Clonal selection of somatic mutations shapes NK cell memory
- Parent-of-origin phasing of somatic mutations shows equal mutation burden between parental genomes in human cancers
- Somatic mutations and single-nucleus transcriptomics reveal uniquely human properties of neuronal aging
BibTeX
@misc{zhou2026vafsoma,
title = {{vafSOMA: Detecting cell-layer-fixed somatic mutations in plants}},
author = {M. Zhou and F. Johannes},
howpublished = {bioRxiv},
year = {2026},
month = nov,
doi = {10.64898/2026.11.02.763622},
url = {https://www.biorxiv.org/content/10.64898/2026.11.02.763622v1},
note = {CC BY},
}
RIS
TY - UNPB
TI - vafSOMA: Detecting cell-layer-fixed somatic mutations in plants
AU - M. Zhou
AU - F. Johannes
PY - 2026
DA - 2026/11/02
PB - bioRxiv
DO - 10.64898/2026.11.02.763622
UR - https://www.biorxiv.org/content/10.64898/2026.11.02.763622v1
AB - Somatic mutations are a source of heritable genetic variation in plants and have long been exploited by breeders through selection and vegetative propagation of bud sports. They often arise in specific meristem layers and can become fixed as descendant cell lineages expand during development. Bulk sequencing of plant organs mixes DNA from these layers, so the expected variant allele frequency of a layer-fixed mutation depends on that layer's tissue contribution, which complicates their detection. Here we present vafSOMA, a statistical model that uses tissue composition to detect fixed somatic mutations, assign them to layers and obtain counts for mutation-rate estimation. It accommodates bulk and layer-enriched samples, including those that retain DNA from other layers after enrichment. In simulations, vafSOMA performed close to the best variant allele frequency cutoff chosen in hindsight, without requiring a user-selected cutoff. Its diagnostics flagged layers for which the data did not support separate mutation counts or rates. Analysing replicate organs substantially improved detection of mutations fixed in the branch meristem and yielded estimates close to the simulated stem-cell mutation rate. Application to apricot sequencing data showed that most reference mutations fixed in the dominant layer were recovered from fruit mixtures at high depth, whereas epidermal mutations remained unresolvable in fruit and leaves. These findings suggest that published bulk somatic sequencing studies mainly capture mutations fixed in the dominant tissue layer and may therefore underestimate the combined input of fixed mutations by roughly two- to threefold when layer-specific rates are similar. vafSOMA is available at https://github.com/jlab-code/vafSOMA.
ER -