buscolite 23.10.26

buscolite: busco analysis for gene predictions

BUSCOlite: simplified BUSCO analysis for gene prediction

BUSCOlite can run the miniprot/Augustus mediated genome predictions as well as the pyhmmer HMM predictions using the BUSCO v9 or v10 databases. It also provides a python API to run busco analysis from within python, ie to be used inside the eukaryotic gene prediction pipeline Funannotate.

This tool is not meant to be a replacment of BUSCO, for most general use cases you should continue to use BUSCOv5

BUSCO models/lineages can be downloaded from the BUSCO site: v5, v4. BUSCOlite does not provide an internal method to do this, as it is trivial to download the lineage you need from your organism(s) by following these links.

There are limited dependencies with BUSCOlite:

• augustus (note: many versions on conda have non-functional PPX/--proteinprofile mode)
• miniprot
• pyhmmer
• pyfastx
• natsort

Why?

Funannotate uses BUSCO to find core conserved marker genes that it uses as a basis to train several ab-initio gene predictors. When BUSCO v2 came out it was python3 only and at that time funannotate was still python2, so I modified the BUSCOv2 source code to be compatible with python2 so it could be run within funannotate. Now BUSCOv5 is the current release, that has numerous bells and whistles that funannotate does not need (no knock against bells and whistles) but the real problem is that due to the large number of dependencies associated with these extra tools is that I cannot build a conda image that includes funannotate and BUSCOv5. So I re-wrote BUSCOv2 here so that it has limited dependencies and will make it easier to incorporate as a dependency of funannotate. A side note is that the metaeuk method that BUSCOv5 now uses as default does not produce complete gene models, in fact the protein sequences it outputs have lowercase sequences that are actually not found in your genome at all. So for training ab-initio predictors, the metaeuk method is not useful -- however, it is faster to get your simple stats on "how complete is my genome assembly".

To install release versions use the pip package manager, like so:

python -m pip install buscolite

To install the most updated code in master you can run:

python -m pip install git+https://github.com/nextgenusfs/buscolite.git