Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2022 Oct 4;32(4):10736.
doi: 10.4081/ejtm.2022.10736.

New indexes for myofibril linearity in muscle image analysis

Ettore Rocchi  1 , Sara Peluso  2 , Stefano Amatori  3 , Davide Sisti  4
Affiliations

New indexes for myofibril linearity in muscle image analysis

Ettore Rocchi et al. Eur J Transl Myol. .

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.

PubMed Disclaimer

Conflict of interest 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.

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. proposed an index that can formally be described as follows:

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 (lsi) 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 lht), 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):

The second proposal is based on circular statistics, taking a cue from the same αi proposed by Cisterna et al.; circular statistics has been used earlier to study muscle cell alignment. Ideally, for each angle αi, 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):

where sin αi and cas αi are respectively the mean of the sine and cosine components of the n αi angles.

(b) Much more important would be the length of the mean vector r:

the mean vector length r can range between 0 (representing perfect isotropy – or a circular uniform distribution – i.e., the maximum possible misalignment) and 1 (representing perfect anisotropy, i.e., the maximum possible alignment). This approach, as opposed to the one proposed by Cisterna et al., 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) 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:

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). Moreover, a concentration parameter k 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):

Nevertheless, for n . 15 or r < 0.45, this estimation would be corrected as follows:

In this way, the concentration parameter k 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), 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.,

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. - PubMed
    1. 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. - PubMed
    1. 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.
    1. Batschelet E. Circular statistics in biology. Academic Press. 111 Fifth Ave., New York, NY 10003; 1981.
    1. 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 in Statistics-Simulation and Computation. 1981; 10(5), 493-502.