Apis mellifera graph genome

<strong>AmelGraph 1.1.0</strong> We aligned 5 different publicly available assemblies with the cactus pangenome workflow (v2.0.5): <strong>Use </strong> <strong>Species </strong> <strong>Genome ID </strong> <strong>Accession </strong> <strong>Graph ID </strong> <strong>Size (Mb) </strong> Reference <em>A. mellifera</em> (DH4) Amel_HAv3.1 GCF_003254395.2 DH4 225.2 Derivate <em>A. m. mellifera </em> INRA_AMelMel_1.0 GCA_003314205.1 mellifera 227.0 Derivate <em>A. m. carnica </em> ASM1384124v2 GCA_013841245.2 carnica 226.0 Derivate <em>A. m. caucasica </em> ASM1384120v1 GCA_013841205.1 caucasica 224.8 Derivate <em>A. m. ligustica </em> ASM1932182v1 GCA_019321825.1 ligustica 231.1 We used the cactus¹ pangenome workflow to generate a 5-ways pangenome alignment, masking in blocks of 10Kb, and considering all the sequences in the different data. Due to incompatibilities with the pangenome workflow generation of the indexes, we regenerated the giraffe indexes from the output GFA/VCF file. We also downloaded annotations available for three genomes (DH4, <em>A. m. carnica</em> and <em>A. m. caucasica</em>), modifying the naming of the contigs to include the subspecies name (i.e. 'LG1' for <em>A. m. caucasica</em> modified to 'caucasica.LG1'). We then ran the vg autoindex function: <pre><code>vg autoindex -w giraffe -o pangenome -t 8 -T ./TMP -x carnica.gff -x caucasica.gff -R XG -x DH4.gff</code></pre> We compared this 'full' cactus graph to one built using only sequence on the linkage groups, and to one generated using the PGGB workflow. We evaluated the different graphs using PGGE. This full cactus graph returned the highst aligned identity, sequence matches, and unique alignments, while having the lowest number of multiple mapping and missing alignments. <strong>Parameter </strong> <strong>CACTUS (LG) </strong> <strong>CACTUS (FULL) </strong> <strong>PGGB (norm) </strong> HAv3.1 size 221,626,419 225,250,884 225,250,884 Graph length (bp) 243,077,200 246,675,144 315,672,860 Extra sequence 21,450,781 21,424,260 90,421,976 # nodes 17,120,383 11,664,080 17,952,597 # edges 21,385,634 15,990,593 21,702,958 <strong>References</strong> 1. Armstrong, J. <em>et al</em>. Progressive Cactus is a multiple-genome aligner for the thousand-genome era. <em>Nature</em> <strong>587</strong>, (2020).