Minimizing detection bias of somatic mutations in a highly heterozygous oak genome

Somatic mutations are particularly relevant for long-lived organisms. Sources of somatic mutations include imperfect DNA repair, replication errors, and exogenous damage such as ultraviolet radiation. A previous study estimated a surprisingly low number of somatic mutations in a 234-year-old individ...

Descripción completa

Detalles Bibliográficos
Autores: Xian, Wenfei, Carbonell-Bejerano, Pablo [0000-0002-7266-9665], Rabanal, Fernando A., Bezrukov, Ilja, Reymond, Philippe, Weigel, Detlef
Tipo de recurso: artículo
Estado:Versión publicada
Fecha de publicación:2025
País:España
Institución:Universidad de La Rioja (UR)
Repositorio:RIUR. Repositorio Institucional de la Universidad de La Rioja
OAI Identifier:oai:portal.dialnet.es:doc/68b33d7d2f359d6e8f724e4f
Acceso en línea:https://investigacion.unirioja.es/documentos/68b33d7d2f359d6e8f724e4f
Access Level:acceso abierto
Palabra clave:Quercus robur
HiFi assembly
Somatic mutations
Descripción
Sumario:Somatic mutations are particularly relevant for long-lived organisms. Sources of somatic mutations include imperfect DNA repair, replication errors, and exogenous damage such as ultraviolet radiation. A previous study estimated a surprisingly low number of somatic mutations in a 234-year-old individual of the pedunculate oak (Quercus robur), known as the Napoleon Oak. It has been suggested that the true number of somatic mutations was underestimated due to gaps in the reference genome and too conservative filtering of potential mutations. We therefore generated new high-fidelity long-read data for the Napoleon Oak (n = 12) to produce both a pseudo-haploid genome assembly and a partially phased diploid assembly. The high heterozygosity allowed for complete reconstruction of phased and gapless centromeres for 22 of the 24 chromosomes. On the other hand, the high heterozygosity posed challenges for short-read alignments. Use of only the pseudo-haploid assembly as a reference led to potential misalignments, while use of only the diploid assembly reduced variant detection sensitivity. Since most somatic mutations are layer-specific, the fraction of reads covering a specific somatic mutation is expected to be relatively low, even where all cells in a single layer contain a specific mutation. To address this challenge, we employed a read assignment strategy, selecting the appropriate reference sequence (pseudo-haploid or diploid) based on alignment score and mapping quality. Ultimately, we identified 198 high-confidence somatic mutations, compared with 17 somatic mutations identified before with the same set of short reads. Our approach thus increased the total estimated annual mutation rate by a factor of 5.