PySTACHIO: Python Single-molecule TrAcking stoiCHiometry Intensity and simulatiOn, a flexible, extensible, beginner-friendly and optimized program for analysis of single-molecule microscopy data

• Published in: Computational and structural biotechnology journal Vol. 19; pp. 4049 - 4058
• Main Authors: Shepherd, Jack W., Higgins, Ed J., Wollman, Adam J.M., Leake, Mark C.
• Format: Journal Article
• Language: English
• Published: Netherlands Elsevier B.V 01.01.2021; Research Network of Computational and Structural Biotechnology; Elsevier
• Subjects: Analysis software; biochemical compounds; biotechnology; cameras; computer software; Diffusion coefficients; dyes; fluorescence; Image analysis; light microscopy; Molecular stoichiometry; SMLM; stoichiometry; Super-resolution microscopy; Super-resolution microscopy; Analysis software; Image analysis; SMLM; Diffusion coefficients; Molecular stoichiometry
• ISSN: 2001-0370; 2001-0370
• DOI: 10.1016/j.csbj.2021.07.004

[Display omitted] •We present PySTACHIO, a refined version of our spot tracking algorithm.•We demonstrate highly improved performance over previous MATLAB versions.•PySTACHIO can accurately estimate stoichiometries and 2D diffusion coefficients.•Performance is comparable to state-of-the-art packages on challenge data.•PySTACHIO has both GUI and command line interfaces and can be hosted as a web app. As camera pixel arrays have grown larger and faster, and optical microscopy techniques ever more refined, there has been an explosion in the quantity of data acquired during routine light microscopy. At the single-molecule level, analysis involves multiple steps and can rapidly become computationally expensive, in some cases intractable on office workstations. Complex bespoke software can present high activation barriers to entry for new users. Here, we redevelop our quantitative single-molecule analysis routines into an optimized and extensible Python program, with GUI and command-line implementations to facilitate use on local machines and remote clusters, by beginners and advanced users alike. We demonstrate that its performance is on par with previous MATLAB implementations but runs an order of magnitude faster. We tested it against challenge data and demonstrate its performance is comparable to state-of-the-art analysis platforms. We show the code can extract fluorescence intensity values for single reporter dye molecules and, using these, estimate molecular stoichiometries and cellular copy numbers of fluorescently-labeled biomolecules. It can evaluate 2D diffusion coefficients for the characteristically short single-particle tracking data. To facilitate benchmarking we include data simulation routines to compare different analysis programs. Finally, we show that it works with 2-color data and enables colocalization analysis based on overlap integration, to infer interactions between differently labelled biomolecules. By making this freely available we aim to make complex light microscopy single-molecule analysis more democratized.