ο€  Myofibril linearity in muscle image analysis Eur J Transl Myol 32 (4): 10736, 2022 doi: 10.4081/ejtm.2022.10736 - 1 - New indexes for myofibril linearity in muscle image analysis Ettore Rocchi (1), Sara Peluso (1), Stefano Amatori (2), Davide Sisti (2) (1) Independent researcher, Urbino, Italy; (2) Department of Biomolecular Sciences, University of Urbino Carlo Bo, Urbino, Italy. This article is distributed under the terms of the Creative Commons Attribution Noncommercial License (CC BY-NC 4.0) which permits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited. Abstract The endeavor to evaluate the linearity of myofibrillar structures and their potential deviation from a straight line is a fascinating problem in muscle tissue image analysis. In this Letter, we suggest two different strategies for solving the same challenge. The first strategy is based on an alignment index, which could be derived by comparing the sum of the lengths of the individual sarcomeres with the distance between the "head" of the first and the "tail" of the last sarcomere. The second strategy relies on circular statistics, which takes a cue from an already suggested method. Our proposed methods are alternatives: the former has the advantage of simplicity; the latter is certainly more elegant and gives greater substance to statistical analysis, but in contrast, it also has greater computational complexity. Key Words: Circular statistics; myofibrillar structures; sarcomeres. Eur J Transl Myol 32 (4): 10736, 2022 doi: 10.4081/ejtm.2022.10736 An intriguing problem in analysing the image of muscle tissue is represented by the attempt to evaluate the linearity of myofibrillar structures and their possible deviation from a straight line. By indicating with Ξ±i the angles between two successive sarcomeres along the same myofibril, Cisterna et al.1 proposed an index that can formally be described as follows: 𝛼𝛼� = βˆ‘ |𝛼𝛼𝑖𝑖|𝑛𝑛 𝑖𝑖=1 𝑛𝑛 where n is the number of angles (i.e., the number of sarcomeres minus 1); the absolute value was introduced to avoid that random fluctuations around a hypothetical central axis could lead to an average equal to zero (i.e., the same result we would obtain in case of perfect alignment). In this paper, we propose two alternative approaches to the same problem. An alignment index could be obtained by comparing the sum of the lengths of the individual sarcomeres (𝑙𝑙𝑙𝑙𝑖𝑖) with the distance between the "head" of the first sarcomere (indicated by h) and the "tail" of the last sarcomere (indicated by t) (we will call this distance π‘™π‘™β„Žπ‘‘π‘‘), so that if the sarcomeres are perfectly aligned, the two quantities coincide, and their ratio is therefore 1; on the contrary, the more the behaviour is "zigzag", the more this ratio increases (> 1): βˆ‘ 𝑙𝑙𝑙𝑙𝑖𝑖𝑛𝑛 𝑖𝑖=1 π‘™π‘™β„Žπ‘‘π‘‘ The second proposal is based on circular statistics, taking a cue from the same 𝛼𝛼𝑖𝑖 proposed by Cisterna et al.;1 circular statistics has been used earlier to study muscle cell alignment.2 Ideally, for each angle 𝛼𝛼𝑖𝑖, it is possible to construct a unit vector having the base at the centre of a goniometric circumference and the vertex on the circumference from which to calculate the resulting vector, which, divided by n (that is the number of angles or the number of sarcomeres minus 1), gives the resulting mean vector. From here, two parameters can be obtained (a) The direction of the mean vector π›Όπ›ΌοΏ½βˆ— (if it is equal to 0, it means that the displacements to the right balance those to the left; or, as an extreme case, that the sarcomeres are perfectly aligned): π›Όπ›ΌοΏ½βˆ— = ⎩ βŽͺ ⎨ βŽͺ ⎧ π‘‘π‘‘π‘‘π‘‘π‘›π‘›βˆ’1 𝑙𝑙𝑠𝑠𝑛𝑛 𝛼𝛼𝑖𝑖 𝑐𝑐𝑐𝑐𝑙𝑙 𝛼𝛼𝑖𝑖 , 𝑐𝑐𝑐𝑐𝑙𝑙 𝛼𝛼𝑖𝑖 > 0 πœ‹πœ‹ + π‘‘π‘‘π‘‘π‘‘π‘›π‘›βˆ’1 𝑙𝑙𝑠𝑠𝑛𝑛 𝛼𝛼𝑖𝑖 𝑐𝑐𝑐𝑐𝑙𝑙 𝛼𝛼𝑖𝑖 , 𝑙𝑙𝑠𝑠𝑛𝑛 𝛼𝛼𝑖𝑖 < 0 where 𝑙𝑙𝑠𝑠𝑛𝑛 𝛼𝛼𝑖𝑖 and 𝑐𝑐𝑐𝑐𝑙𝑙 𝛼𝛼𝑖𝑖 are respectively the mean of the sine and cosine components of the n 𝛼𝛼𝑖𝑖 angles. (b) Much more important would be the length of the mean vector r: π‘Ÿπ‘Ÿ = 1 𝑛𝑛 οΏ½(βˆ‘ 𝑙𝑙𝑠𝑠𝑛𝑛 𝛼𝛼𝑖𝑖𝑛𝑛 𝑖𝑖=1 )2 + (βˆ‘ 𝑐𝑐𝑐𝑐𝑙𝑙 𝛼𝛼𝑖𝑖𝑛𝑛 𝑖𝑖=1 )2; the mean vector length r can range between 0 (representing perfect isotropy – or a circular uniform distribution – i.e., the maximum possible misalignment) Myofibril linearity in muscle image analysis Eur J Transl Myol 32 (4): 10736, 2022 doi: 10.4081/ejtm.2022.10736 - 2 - and 1 (representing perfect anisotropy, i.e., the maximum possible alignment). This approach, as opposed to the one proposed by Cisterna et al.,1 is based on the hypothesis of circularity (and non-linearity of the angles). It is also possible to perform a statistical test using the mean vector length r as test statistics: Rayleigh test (1919)3 is the best known in circular statistics; its null hypothesis is the uniform circular distribution of the angles, and the alternative hypothesis is a generic anisotropy. Likewise, the V-test is also based on the mean vector length r, which uses the projection V of the mean vector over an a priori direction (in our case, we can think of it as the direction of the muscle fibre) as test statistics: V=r cosΟ†, where πœ‘πœ‘ is the angle between the mean vector and the a priori direction. Both the Rayleigh test and the V-test have appropriate tables of critical values (see Batschelet 1981).4 Moreover, a concentration parameter π‘˜π‘˜ can be estimated [using the Maximum Likelihood approach (ML)] starting from the mean vector length r, under the hypothesis of a Von Mises distribution (an analogue of Normal distribution for angular data); the following approximated formula was proposed by Best and Fisher (1981)5: π‘˜π‘˜οΏ½π‘€π‘€π‘€π‘€ = ⎩ βŽͺ ⎨ βŽͺ ⎧ 2π‘Ÿπ‘Ÿ + π‘Ÿπ‘Ÿ3 + 5 6 π‘Ÿπ‘Ÿ5 , π‘Ÿπ‘Ÿ < 0.53 βˆ’0.4 + 1.39π‘Ÿπ‘Ÿ + 0.43 1 βˆ’ π‘Ÿπ‘Ÿ , 0.53 ≀ π‘Ÿπ‘Ÿ < 0.85 1 π‘Ÿπ‘Ÿ3 βˆ’ 4π‘Ÿπ‘Ÿ2 + 3π‘Ÿπ‘Ÿ , π‘Ÿπ‘Ÿ β‰₯ 0.85 Nevertheless, for n ≀ 15 or r < 0.45, this estimation would be corrected as follows: π‘˜π‘˜οΏ½ = ⎩ βŽͺ ⎨ βŽͺ βŽ§π‘šπ‘šπ‘‘π‘‘π‘šπ‘š οΏ½π‘˜π‘˜οΏ½π‘€π‘€π‘€π‘€ βˆ’ 2 π‘˜π‘˜οΏ½π‘€π‘€π‘€π‘€ , 0οΏ½ , π‘˜π‘˜οΏ½π‘€π‘€π‘€π‘€ < 2 (𝑛𝑛 βˆ’ 1)3π‘˜π‘˜οΏ½π‘€π‘€π‘€π‘€ 𝑛𝑛3 + 𝑛𝑛 , π‘˜π‘˜οΏ½π‘€π‘€π‘€π‘€ β‰₯ 2 In this way, the concentration parameter π‘˜π‘˜οΏ½ should also be used as a linearity index, useful for our purposes: the higher the k value, the higher the alignment. In conclusion, our first method represents an alternative to the one previously proposed by Cisterna et al. (2021),1 equivalent in terms of potential but with the practical advantage of not having to measure angles, but only lengths, which makes it much more "convenient”. The second method we propose, on the other hand, uses angles but uses circular analysis techniques instead of linear analysis methods, which makes it more elegant and gives greater substance to statistical analysis, but in contrast, it also has greater computational complexity. Our methods have potential use in several sarcomere- related conditions by providing a quantitative definition of myofibril linearity in skeletal muscle.6,7 List of acronyms 𝑙𝑙𝑙𝑙𝑖𝑖 - lengths of the individual sarcomeres ML - Maximum Likelihood approach Contributions of Authors All authors have read and approved the final edited typescript. Acknowledgments None Funding The authors received no specific funding for this work. Conflict of Interest The authors declare no financial, personal, or other conflicts of interest. Ethical Publication Statement We confirm that we have read the Journal’s position on issues involved in ethical publication and affirm that this report is consistent with those guidelines. Corresponding Author Stefano Amatori, Department of Biomolecular Sciences, University of Urbino Carlo Bo, Piazza Rinascimento 7, 61029 Urbino, Italy ORCID ID: 0000-0001-7497-755X E-mail: stefano.amatori1@uniurb.it E-mails and ORCID iD of co-authors Ettore Rocchi: ettoreroc@gmail.com ORCID id: 0000-0002-7612-2819 Sara Peluso: sarapeluso99@gmail.com ORCID id: 0000-0003-1924-1093 Davide Sisti: davide.sisti@uniurb.it ORCID id: 0000-0002-7925-7495 References 1. Cisterna B, Malatesta M, Zancanaro C, Boschi F. A computational approach to quantitatively define sarcomere dimensions and arrangement in skeletal muscle. Comput Methods Programs Biomed. 2021 Nov; 211:106437. doi: 10.1016/j.cmpb.2021. 106437. Epub 2021 Sep 24. 2. Coletti D, Teodori L, Albertini MC, Rocchi M, PristerΓ  A, Fini M, Molinaro M, Adamo S. Static magnetic fields enhance skeletal muscle differentiation in vitro by improving myoblast alignment. Cytometry A. 2007 Oct;71(10):846-56. doi: 10.1002/cyto.a.20447. 3. Rayleigh L. XXXI. On the problem of random vibrations, and of random flights in one, two, or three dimensions. The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science. 1919; 37(220), 321-347. 4. Batschelet E. Circular statistics in biology. Academic Press. 111 Fifth Ave., New York, NY 10003; 1981. 5. Best DJ & Fisher NI. The BIAS of the maximum likelihood estimators of the von Mises-Fisher concentration parameters: the BIAS of the maximum likelihood estimators. Communications mailto:stefano.amatori1@uniurb.it mailto:ettoreroc@gmail.com mailto:sarapeluso99@gmail.com mailto:davide.sisti@uniurb.it Myofibril linearity in muscle image analysis Eur J Transl Myol 32 (4): 10736, 2022 doi: 10.4081/ejtm.2022.10736 - 3 - in Statistics-Simulation and Computation. 1981; 10(5), 493-502. 6. Edmunds KJ, GΓ­slason MK, Arnadottir ID, Marcante A, Piccione F, Gargiulo P. Quantitative Computed Tomography and Image Analysis for Advanced Muscle Assessment. Eur J Transl Myol. 2016 Jun 22;26(2):6015. doi: 10.4081/ejtm.2016.6015. eCollection 2016 Jun 13. 7. Coletti C, Acosta GF, Keslacy S, Coletti D. Exercise-mediated reinnervation of skeletal muscle in elderly people: An update. Eur J Transl Myol. 2022 Feb 28;32(1). doi: 10.4081/ejtm.2022.10416. Disclaimer All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher. Submission: August 21, 2022 Revision received: August 26, 2022 Accepted for publication: August 26, 2022