Fast, flexible gene cluster family delineation with IGUA

Prokaryotic genomes harbor a variety of functional elements encoded as contiguous multi-gene clusters, with biosynthetic gene clusters (BGCs, genetic determinants of secondary metabolite biosynthesis) serving as a notable example. In a typical workflow, BGCs are clustered into Gene Cluster Families (GCFs), units that group BGCs encoding similar biosynthetic pathways together. However, existing methods cannot readily scale to massive datasets and cannot be used for GCF delineation tasks beyond BGC clustering. Here, we present IGUA (Iterative Gene clUster Analysis; https://github.com/zellerlab/IGUA), a scalable, flexible GCF delineation method for genomic segments with multi-gene architectures. On a BGC clustering task, IGUA is [≥]10x faster than the state-of-the-art (BiG-SCAPE/BiG-SLiCE), without sacrificing accuracy. To highlight its scalability, we use IGUA to cluster >2.8 million BGCs from {approx}1 million prokaryotic genomes in View larger version (18K): org.highwire.dtl.DTLVardef@1f27216org.highwire.dtl.DTLVardef@203ccborg.highwire.dtl.DTLVardef@7766c8org.highwire.dtl.DTLVardef@fcd4b5_HPS_FORMAT_FIGEXP M_FIG C_FIG