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.

Recommendations for extending the GFF3 specification for improved interoperability of genomic data

[A consensus on the standardization of the next generation sequencing process for the diagnosis of genetic diseases (3) - Data analysis].

Recommendations for the FAIRification of genomic track metadata

Gfastats: conversion, evaluation and manipulation of genome sequences using assembly graphs

Efficient indexing and querying of annotations in a pangenome graph

Standards recommendations for the Earth BioGenome Project

TextFormats: Simplifying the definition and parsing of text formats in bioinformatics

FAIR Header Reference genome: a TRUSTworthy standard

DATA RESOURCES AND ANALYSES FAIR Header Reference genome: A TRUSTworthy standard

Best practices for genetic and genomic data archiving

Genomics and Biological Big Data: Facing Current and Future Challenges around Data and Software Sharing and Reproducibility

OME-NGFF: a next-generation file format for expanding bioimaging data-access strategies

Glbase: a Framework for Combining, Analyzing and Displaying Heterogeneous Genomic and High-Throughput Sequencing Data

The GIAB genomic stratifications resource for human reference genomes

GVAF: generalized, flexible filtering software for annotated variant files

OmniGenBench: Automating Large-scale in-silico Benchmarking for Genomic Foundation Models

Updates to the Alliance of Genome Resources Central Infrastructure

Functional equivalence of genome sequencing analysis pipelines enables harmonized variant calling across human genetics projects

Assessing and assuring interoperability of a genomics file format

The changing face of genome assemblies: Guidance on achieving high‐quality reference genomes

Variant interpretation using population databases: lessons from gnomAD