Abstract
Cataloguing the distribution of genes within natural bacterial populations is essential for understanding evolutionary processes and the genetic basis of adaptation. Advances in whole genome sequencing technologies have led to a vast expansion in the amount of bacterial genomes deposited in public databases. There is a pressing need for software solutions which are able to cluster, catalogue and characterise genes, or other features, in increasingly large genomic datasets. Results: Here we present a pangenomics toolbox, PIRATE (Pangenome Iterative Refinement and Threshold Evaluation), which identifies and classifies orthologous gene families in bacterial pangenomes over a wide range of sequence similarity thresholds. PIRATE builds upon recent scalable software developments to allow for the rapid interrogation of thousands of isolates. PIRATE clusters genes (or other annotated features) over a wide range of amino acid or nucleotide identity thresholds and uses the clustering information to rapidly identify paralogous gene families and putative fission/fusion events. Furthermore, PIRATE orders the pangenome using a directed graph, provides a measure of allelic variation, and estimates sequence divergence for each gene family. Conclusions: We demonstrate that PIRATE scales linearly with both number of samples and computation resources, allowing for analysis of large genomic datasets, and compares favorably to other popular tools. PIRATE provides a robust framework for analysing bacterial pangenomes, from largely clonal to panmictic species.
Original language | English |
---|---|
Article number | giz119 |
Journal | GigaScience |
Volume | 8 |
Issue number | 10 |
DOIs | |
Publication status | Published - 1 Oct 2019 |
Bibliographical note
Funding Information:This work has been supported by BBSRC/NERC grant BB/M026388/1 awarded to E.F and MRC grant MR/L015080/1 awarded to S.S.
Publisher Copyright:
© 2019 The Author(s) 2019. Published by Oxford University Press.
Keywords
- bioinformatics
- microbial genomics
- next-generation sequencing
- pangenomics