DOI: 10.17044/scilifelab.28606814.v1
title: Gene annotation of <i>Blastobotrys mokoenaii</i>, <i>Blastobotrys illinoisensis</i>, and <i>Blastobotrys malaysiensis</i>

description: This dataset contains the gene annotation data for three species of <i>Blastobotrys</i> yeats: <i>B. mokoenaii</i>, <i>B. illinoisensis</i>, and <i>B. malaysiensis</i>.The genome assemblies for <i>B.</i><i> mokoenaii </i>(NRRL Y-27120)<i> </i>and <i>B. malaysiensis </i>(NRRL Y-6417) were publicly available on the National Center for Biotechnology Information (NCBI) under accessions GCA_003705765.3 and GCA_030558815.1, respectively.The genome assembly for <i>B. illinoisensis</i> (NRRL YB-1343) was generated by SciLifeLab's National Genomics Infrastructure (NGI) using PacBio long-read data and deposited in the European Nucleotide Archive (ENA) under accession GCA_965113335.1.<b>File description</b>bmokoenaii_annotation.gff<br>This file contains the gene models predicted for <i>B. mokoenaii </i>(GCA_003705765.3).billinoisensis_annotation.gff<br>This file contains the gene models predicted for <i>B. illinoisensis </i>(GCA_003705765.3).bmalaysiensis_annotation.gff<br>This file contains the gene models predicted for <i>B. malaysiensis </i>(GCA_030558815.1).<b>Gene annotation methods</b>Repeat MaskingPrior to annotation, a repeat library was built for each species using RepeatModeler2 v2.0.2 and the genomes were soft-masked using RepeatMasker v4.1.5.<br>$ RepeatModeler -database ${DB} -engine ncbi -pa 16<br>$ RepeatMasker -dir . -gff -u -no_is -xsmall -e ncbi -lib ${LIBRARY} -pa 16 genome.fastaStructural Annotation<br>Structural annotation was performed on the soft-masked genomes using Braker3 v3.0.3 incorporating external evidence in the form of all fungal proteins from OrthoDB v11 (available at https://bioinf.uni-greifswald.de/bioinf/partitioned_odb11).<br>$ braker.pl --genome="$genome" \--prot_seq=${protein} --workingdir=${PWD} \<br>--gff3 --threads=16 --verbosity=3 \<br>--nocleanup --species=${i}<br>Functional Annotation<br>The predicted genes were functionally annotated using the National Bioiformatics Infrastructure Sweden (NBIS) <code>functional_annotation</code> nextflow pipeline<code> </code>v2.0.0 (https://github.com/NBISweden/pipelines-nextflow). Briefly, this pipeline performs similarity searches between the annotated proteins and the UniProtKB/Swiss-Prot database (downloaded on 2023-12) using the Basic Local Alignment Search Tool (BLAST). Then it uses InterProScan to query the proteins against InterPro v59-91 databases, and merges results using AGAT v1.2.0.tRNAs and rRNAs<br>Transfer RNA (tRNA) and ribosomal RNA (rRNA) genes were annotated using tRNAscan-SE v2.0.12 and barrnap<code> </code>v0.9, respectively. Other ncRNAs, such as SRP RNA, RNase P RNA, spliceosomal ncRNAs etc. have not been predicted. Finnally, the functionally annotated protein-coding genes, tRNAs, and rRNAs were combined into a single GFF file using AGAT v1.2.0.$ tRNAscan-SE -E --gff ${output}_trnas.gff --thread 16 ${genome}.fasta<br>$ barrnap --kingdom euk --threads 6 ${genome}.fasta &gt; ${output}_rrna.gffAnnotation integrationFinnally, the functionally annotated protein-coding genes, tRNAs, and rRNAs were combined into a single GFF file using AGAT v1.2.0.$ agat_sp_complement_annotations.pl --ref ${protein_coding} --add ${trna} --add ${rrna} --out full_annotation.gff