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.

Figeno: multi-region genomic figures with long-read support

Analysis and Visualization of ChIP-Seq and RNA-Seq Sequence Alignments Using Ngs.plot.

Cluster efficient pangenome graph construction with nf-core/pangenome

genomation: a toolkit to summarize, annotate and visualize genomic intervals

D-GENIES: dot plot large genomes in an interactive, efficient and simple way

BRGenomics for analyzing high-resolution genomics data in R

Efficient indexing and querying of annotations in a pangenome graph

gggenomes: effective and versatile visualizations for comparative genomics

HAPPE: A Tool for Population Haplotype Analysis and Visualization in Editable Excel Tables

genomepy: genes and genomes at your fingertips

Panacus: fast and exact pangenome growth and core size estimation

GenoVi, an open-source automated circular genome visualizer for bacteria and archaea

PyamilySeq: A Python Tool for Interpretable Gene (Re)Clustering and Pangenomic Inference Across Species and Genera

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

Giraffe: a tool for comprehensive processing and visualization of multiple long-read sequencing data

Scalepopgen: bioinformatics workflow resources implemented in Nextflow for comprehensive population genomic analyses

Building pangenome graphs

pyHiM: a new open-source, multi-platform software package for spatial genomics based on multiplexed DNA-FISH imaging

GREAP: a comprehensive enrichment analysis software for human genomic regions

funkea: Functional Enrichment Analysis in Python

Unbiased pangenome graphs