Gene annotations for wild alpine reindeer (Rangifer tarandus tarandus)

Here we provide the gene annotations for wild alpine reindeer (also called mountain reindeer; Rangifer tarandus tarandus). In addition, we provide the gene annotations for an updated genome assembly of Svalbard reindeer (Rangifer tarandus platyrhynchus). We provide these for both convenience and because some of the functional annotations of genes/proteins are removed when we prepare these for uploading to ENA. We also provide the FASTA files for the assemblies we have made.

We annotated the genome assemblies using a pre-release version of the EBP-Nor genome annotation pipeline (https://github.com/ebp-nor/GenomeAnnotation). First, AGAT (https://zenodo.org/record/7255559) agat_sp_keep_longest_isoform.pl and agat_sp_extract_sequences.pl were used on the GRCh38 genome assembly and annotation to generate one protein (the longest isoform) per gene. Miniprot (Li, 2023) was used to align the proteins to the curated assemblies. UniProtKB/Swiss-Prot (Consortium et al., 2023) release 2024_04 in addition to the Vertebrata part of OrthoDB v11 (Kuznetsov et al., 2022) were also aligned separately to the assemblies. Red (Girgis, 2015) was run via redmask (https://github.com/nextgenusfs/redmask) on the assemblies to mask repetitive areas. GALBA (Brůna et al., 2023; Buchfink et al., 2015; Hoff and Stanke, 2018; Li, 2023; Stanke et al., 2006) was run with the GRCh38 proteins using the miniprot mode on the masked assemblies. The funannotate-runEVM.py script from Funannotate was used to run EvidenceModeler (Haas et al., 2008) on the alignments of GRCh38 proteins, UniProtKB/Swiss-Prot proteins, Vertebrata proteins and the predicted genes from GALBA. The resulting predicted proteins were compared to the protein repeats that Funannotate distributes using DIAMOND blastp, and the predicted genes were filtered based on this comparison using AGAT. The filtered proteins were compared to the UniProtKB/Swiss-Prot release 2024_04 using DIAMOND (Buchfink et al., 2015) blastp to find gene names, and InterProScan was used to discover functional domains. AGATs agat_sp_manage_functional_annotation.pl was used to attach the gene names and functional annotations to the predicted genes.

List of files provided here and their description:

mRanTar1.2.hap1.fa.gz  - genome assembly of Svalbard reindeer (hap1)

mRanTar1.2.hap1.gff.gz - genome annotation of Svalbard reindeer (hap1)

mRanTar1.2.hap1.proteins.fa.gz - predicted proteins of Svalbard reindeer (hap1)

mRanTar1.2.hap2.fa.gz - genome assembly of Svalbard reindeer (hap2)

mRanTar1.2.hap2.gff.gz - genome annotation of Svalbard reindeer (hap2)

mRanTar1.2.hap2.proteins.fa.gz - predicted proteins of Svalbard reindeer (hap2)

mRanTar2.1.hap1.fa.gz - genome assembly of wild alpine reindeer (hap1)

mRanTar2.1.hap1.gff.gz - genome annotation of wild alpine reindeer (hap1)

mRanTar2.1.hap1.proteins.fa.gz - predicted proteins of wild alpine reindeer (hap1)

mRanTar2.1.hap2.fa.gz - genome assembly of wild alpine reindeer (hap2)

mRanTar2.1.hap2.gff.gz - genome annotation of wild alpine reindeer (hap2)

mRanTar2.1.hap2.proteins.fa.gz - predicted proteins of wild alpine reindeer (hap2)