A workflow for rapid unbiased quantification of fibrillar feature alignment in biological images

• Published in: Frontiers in computer science (Lausanne) Vol. 3
• Main Authors: Marcotti, Stefania, de Freitas, Deandra Belo, Troughton, Lee D, Kenny, Fiona N, Shaw, Tanya J, Stramer, Brian M, Oakes, Patrick W
• Format: Journal Article
• Language: English
• Published: Switzerland Frontiers Media S.A 01.10.2021
• Subjects: alignment; anisotropy; cytoskeleton; extracellular matrix (ECM); fast Fourier transform (FFT); fibers; extracellular matrix (ECM); alignment; Fast Fourier Transform (FFT); cytoskeleton; fibres; anisotropy
• ISSN: 2624-9898; 2624-9898
• DOI: 10.3389/fcomp.2021.745831

Measuring the organisation of the cellular cytoskeleton and the surrounding extracellular matrix (ECM) is currently of wide interest as changes in both local and global alignment can highlight alterations in cellular functions and material properties of the extracellular environment. Different approaches have been developed to quantify these structures, typically based on fibre segmentation or on matrix representation and transformation of the image, each with its own advantages and disadvantages. Here we present , a workflow to quantify the alignment of fibrillar features in microscopy images exploiting 2D Fast Fourier Transforms (FFT). Using pre-existing datasets of cell and ECM images, we demonstrate our approach and compare and contrast this workflow with two other well-known ImageJ algorithms to quantify image feature alignment. These comparisons reveal that has a number of advantages due to its grid-based FFT approach. 1) Flexibility in defining the window and neighbourhood sizes allows for performing a parameter search to determine an optimal length scale to carry out alignment metrics. This approach can thus easily accommodate different image resolutions and biological systems. 2) The length scale of decay in alignment can be extracted by comparing neighbourhood sizes, revealing the overall distance that features remain anisotropic. 3) The approach is ambivalent to the signal source, thus making it applicable for a wide range of imaging modalities and is dependent on fewer input parameters than segmentation methods. 4) Finally, compared to segmentation methods, this algorithm is computationally inexpensive, as high-resolution images can be evaluated in less than a second on a standard desktop computer. This makes it feasible to screen numerous experimental perturbations or examine large images over long length scales. Implementation is made available in both MATLAB and Python for wider accessibility, with example datasets for single images and batch processing. Additionally, we include an approach to automatically search parameters for optimum window and neighbourhood sizes, as well as to measure the decay in alignment over progressively increasing length scales.