Abstract:The current reference genome is the backbone of diverse and rich annotations. Simple text formats, like VCF or BED, have been widely adopted and helped the critical exchange of genomic information. There is a dire need for tools and formats enabling pangenomic annotation to facilitate such enrichment of pangenomic references. The Graph Alignment Format (GAF) is a text format, tab-delimited like BED/VCF files, which was proposed to represent alignments. GAF could also be used to store paths representing annotations in a pangenome graph, but there are no tools to index and query them efficiently. Here, we present extensions to vg and HTSlib that provide efficient sorting, indexing, and querying for GAF files. With this approach, annotations overlapping a subgraph can be extracted quickly. Paths are sorted based on the IDs of traversed nodes, compressed with BGZIP, and indexed with HTSlib/tabix via our extensions for the GAF format. Compared to the binary GAM format, GAF files are easier to edit or inspect because they are plain text, and we show that they are twice as fast to sort and half as large on disk. In addition, we updated vg annotate, which takes BED or GFF3 annotation files relative to linear sequences and projects them into the pangenome. It can now produce GAF files representing these annotations' paths through the pangenome. We showcase these new tools on several applications. We projected annotations for all Human Pangenome Reference Consortium Year 1 haplotypes, including genes, segmental duplications, tandem repeats and repeats annotations, into the Minigraph-Cactus pangenome (GRCh38-based v1.1). We also projected known variants from the GWAS Catalog and expression QTLs from the GTEx project into the pangenome. Finally, we reanalyzed ATAC-seq data from ENCODE to demonstrate what a coverage track could look like in a pangenome graph. These rich annotations can be quickly queried with vg and visualized using existing tools like the Sequence Tube Map or Bandage.

Assessing genome conservation on pangenome graphs with PanSel

Wheat Panache - a pangenome graph database representing presence/absence variation across 16 bread wheat genomes

Wheat Panache: A pangenome graph database representing presence-absence variation across sixteen bread wheat genomes

Building pangenome graphs

Panacus: fast and exact pangenome growth and core size estimation

Unbiased pangenome graphs

PanGraph: scalable bacterial pan-genome graph construction

PanKmer: k-mer-based and reference-free pangenome analysis

Revisiting pangenome openness with -mers

PPanG: a precision pangenome browser enabling nucleotide-level analysis of genomic variations in individual genomes and their graph-based pangenome

PanGeneHome : A Web Interface to Analyze Microbial Pangenomes

Panx: Pan-Genome Analysis and Exploration

Exploring environmental intra‐species diversity through non‐redundant pangenome assemblies

Exploring gene content with pangene graphs

Compressive Pangenomics Using Mutation-Annotated Networks

Cluster efficient pangenome graph construction with nf-core/pangenome

PanPA: generation and alignment of panproteome graphs

Pan-evolutionary and regulatory genome architecture delineated by an integrated macro- and microsynteny approach

Efficient indexing and querying of annotations in a pangenome graph