GENCODEHub

GENCODEHub is a Bioconductor AnnotationHub package for GENCODE human gene annotation. It provides GENCODE release 50 for GRCh38.p14 and Ensembl release 116 as separate GRanges resources and as a structural TxDb database.

The annotation files are hosted in Zenodo record 21784647. AnnotationHub downloads only the resource selected by the user and keeps a local cached copy for later use.

1. Install the package

if (!requireNamespace("BiocManager", quietly = TRUE)) {
    install.packages("BiocManager")
}
BiocManager::install("GENCODEHub")

2. List the available resources

library(AnnotationHub)
library(GENCODEHub)

hub <- AnnotationHub()
gencodeV50Resources(hub)

The collection contains the complete GTF annotation, separate resources for gene, transcript, exon, CDS, UTR, start_codon, stop_codon, and Selenocysteine features, and a TxDb database.

3. Retrieve the resource needed for an analysis

genes <- gencodeV50("gene", hub = hub)
transcripts <- gencodeV50("transcript", hub = hub)
txdb <- gencodeV50("TxDb", hub = hub)

Use the gene resource for gene-body overlaps or gene-level transcription start sites. Use the transcript resource when transcript-specific coordinates are required. Use the TxDb for relationships among genes, transcripts, exons and coding sequences.

4. Check coordinates and identifiers

genes
head(genes$gene_id)
S4Vectors::mcols(genes)

The GRanges resources retain the versioned GENCODE identifiers and imported GTF attributes. Coordinates are one-based, closed intervals on GRCh38.p14. GRCh37 or hg19 data must be converted to GRCh38 before direct annotation.

Source and citation

The resources are derived from GENCODE Human Release 50 using the official gencode.v50.annotation.gtf.gz file. Cite the applicable GENCODE publication, the Zenodo record and GENCODEHub when using the resources. Use citation("GENCODEHub") for the package citation.