Representation and detection of inversions in pangenome graphs
- 1Univ Rennes, Inria, CNRS , IRISA—UMR 6074, F-35000, Rennes ,
- 2GenPhySE, Université de Toulouse, INRAE, ENVT , 31326, Castanet Tolosan ,
- 3IGEPP, INRAE, Institut Agro, University of Rennes , F-35653, Le Rheu ,
Abstract
Pangenome graphs are increasingly used in genetic diversity analyses because they reduce reference bias in read mapping and enhance variant discovery and genotyping from SNPs to structural variants. In pangenome graphs, variants appear as bubbles, which can be detected by dedicated bubble calling tools. Although these tools report essential information on the variant bubbles, such as their position and allele walks in the graph, they do not annotate the type of the detected variants. While simple SNPs, insertions, and deletions are easily distinguishable by allele size, large balanced variants like inversions are harder to differentiate among the large number of unannotated bubbles and remain underexplored in pangenome graph benchmarks and analyses. In this work, we focused on inversions, which have been drawing renewed attention in evolutionary genomics studies in the past years, and aimed to assess how this type of variant is handled by state of the art pangenome graph pipelines. We identified two distinct topological motifs for inversion bubbles: one path-explicit and one alignment-rescued, and developed a tool to annotate them from bubble-caller outputs. We constructed pangenome graphs with both simulated data and real data using four state of the art pipelines, and assessed the impact of inversion size, genome divergence, and variant density on inversion representation and accuracy. Our results reveal substantial differences between pipelines in simulated graphs, with some inversions either misrepresented or lost. In addition, recovery rates are strikingly low in real human datasets, highlighting major challenges in analyzing inversions through pangenomic approaches.
Related
- Amira: detection of single and multi-copy antimicrobial resistance genes directly from long reads using gene-space de Bruijn graphs
- Chromosomal inversions modulate 3D chromatin interactions and regulatory architecture in Anopheles gambiae
- From plant pangenome assemblies to operational graph references
- MineGraph facilitates plastid and mitochondrial graph pangenome analysis in plants
- The Emirati pangenome featuring a diploid chromosome-level reference
BibTeX
@article{romain2026representation,
title = {{Representation and detection of inversions in pangenome graphs}},
author = {Sandra Romain and Siegfried Dubois and Fabrice Legeai and Claire Lemaitre},
journal = {NAR Genomics and Bioinformatics},
year = {2026},
month = oct,
doi = {10.1093/nargab/lqag125},
url = {https://academic.oup.com/nargab/article/doi/10.1093/nargab/lqag125/8916833},
note = {CC BY 4.0},
}
RIS
TY - JOUR
TI - Representation and detection of inversions in pangenome graphs
AU - Sandra Romain
AU - Siegfried Dubois
AU - Fabrice Legeai
AU - Claire Lemaitre
PY - 2026
DA - 2026/10/15
JO - NAR Genomics and Bioinformatics
DO - 10.1093/nargab/lqag125
UR - https://academic.oup.com/nargab/article/doi/10.1093/nargab/lqag125/8916833
AB - Pangenome graphs are increasingly used in genetic diversity analyses because they reduce reference bias in read mapping and enhance variant discovery and genotyping from SNPs to structural variants. In pangenome graphs, variants appear as bubbles, which can be detected by dedicated bubble calling tools. Although these tools report essential information on the variant bubbles, such as their position and allele walks in the graph, they do not annotate the type of the detected variants. While simple SNPs, insertions, and deletions are easily distinguishable by allele size, large balanced variants like inversions are harder to differentiate among the large number of unannotated bubbles and remain underexplored in pangenome graph benchmarks and analyses. In this work, we focused on inversions, which have been drawing renewed attention in evolutionary genomics studies in the past years, and aimed to assess how this type of variant is handled by state of the art pangenome graph pipelines. We identified two distinct topological motifs for inversion bubbles: one path-explicit and one alignment-rescued, and developed a tool to annotate them from bubble-caller outputs. We constructed pangenome graphs with both simulated data and real data using four state of the art pipelines, and assessed the impact of inversion size, genome divergence, and variant density on inversion representation and accuracy. Our results reveal substantial differences between pipelines in simulated graphs, with some inversions either misrepresented or lost. In addition, recovery rates are strikingly low in real human datasets, highlighting major challenges in analyzing inversions through pangenomic approaches.
SN - 2631-9268
ER -