Abstract
Skeletal muscle function is inferred from the spatial arrangement of muscle fiber architecture, which corresponds to myofiber molecular and metabolic features. Myofiber features are often determined using immunofluorescence on a local sampling, typically obtained from a median region. This median region is assumed to represent the entire muscle. However, it remains largely unknown to what extent this local sampling represents the entire muscle. We present a pipeline to study the architecture of muscle fiber features over the entire muscle, including sectioning, staining, imaging to image quantification and data-driven analysis with Myofiber type were identified by the expression of myosin heavy chain (MyHC) isoforms, representing contraction properties. We reconstructed muscle architecture from consecutive cross-sections stained for laminin and MyHC isoforms. Examining the entire muscle using consecutive cross-sections is extremely laborious, we provide consideration to reduce the dataset without loosing spatial information. Data-driven analysis with over 150,000 myofibers showed spatial variations in myofiber geometric features, myofiber type, and the distribution of neuromuscular junctions over the entire muscle. We present a workflow to study histological changes over the entire muscle using high-throughput imaging, image quantification, and data-driven analysis. Our results suggest that asymmetric spatial distribution of these features over the entire muscle could impact muscle function. Therefore, instead of a single sampling from a median region, representative regions covering the entire muscle should be investigated in future studies.