Semi-Automated Pipeline Counts Muscle Fiber Nuclei and Satellite Cells at Scale

Skeletal muscle is a remarkably plastic tissue. Its individual fibers—each one a single multinucleated muscle cell—change in size, nuclear content, and even molecular identity in response to exercise, disease, disuse, and aging. For decades, researchers have tracked these changes by staining biopsy samples with fluorescent antibodies and laboriously counting structures by hand under the microscope. That work, while foundational, is slow, subjective, and difficult to scale to the hundreds or even thousands of fibers contained in a single human muscle cross-section. Now, a team at the University of Oregon has developed and validated a semi-automated image analysis pipeline that brings high-throughput quantification of muscle histology within reach of any laboratory equipped with a standard widefield fluorescence microscope.

The new workflow, described in Physiological Reports, measures five key endpoints simultaneously: the number of myonuclei per muscle cell, the number of Pax7-positive satellite cells per cell, the cross-sectional area of each fiber, the perimeter of each fiber, and the fiber’s type, classifying it as slow-twitch type I or fast-twitch type II. According to the authors, this is the first published pipeline to combine all five measurements in a single workflow compatible with widely available widefield microscopy, rather than requiring specialized imaging platforms. Each of these endpoints carries real biological weight. Satellite cells are the resident stem cells of skeletal muscle, and their expansion and integration into fibers as new myonuclei underpins theories such as the myonuclear domain hypothesis, which holds that each nucleus governs a finite volume of cytoplasm. Emerging evidence even suggests the domain may track fiber perimeter more closely than cross-sectional area, making accurate perimeter measurement newly important.

The pipeline begins before any computation, with careful tissue handling and staining. Human biopsies were taken from the vastus lateralis, the lateral thigh muscle that is a standard site in exercise physiology studies. Tissue was embedded in optimal cutting temperature compound, flash-frozen in liquid-nitrogen-cooled isopentane, and sectioned at seven micrometers on a cryostat. Immunohistochemistry then labeled three targets: laminin, a protein that outlines the sarcolemma surrounding each fiber; type I myosin heavy chain, which marks slow fibers; and Pax7, the canonical marker of quiescent and activated satellite cells. Slides were imaged at twenty-fold magnification on a Leica widefield microscope across four channels—DAPI for all nuclei, green fluorescence for Pax7, yellow for myosin heavy chain I, and far-red for laminin—producing the raw material for analysis.

A crucial early stage is manual image cleaning in Fiji/ImageJ. Rather than treating every pixel as trustworthy data, the protocol calls for brightness and contrast optimization on each channel, followed by manual removal of artifacts such as tissue folds, blood vessels, blurred or overlapping laminin signal, and out-of-focus nuclei. The authors argue this step, though the most time-consuming part of the workflow, is essential: it strips out the noise that would otherwise generate false positives downstream and brings the analyzed regions into line with what a careful human analyst would measure. Cleaned images are then exported with a strict naming convention that encodes subject, timepoint, leg, and set, enabling fully automated batch processing.

The computational core pairs two open-source tools. CellProfiler, the modular image-analysis platform maintained at the Broad Institute, orchestrates segmentation and measurement, while Cellpose, a machine-learning model developed at HHMI’s Janelia Research Campus, performs deep-learning-based detection of cellular boundaries. Fiber segmentation uses Cellpose’s cyto2 model on the laminin channel to trace sarcolemmal borders. Nuclear segmentation is trickier, because muscle cross-sections contain nuclei from many cell types, and the default Cellpose nuclei model is highly sensitive. The pipeline handles this with an iterative subtraction strategy: a Gaussian filter smooths the DAPI signal, an initial Cellpose run identifies candidate nuclei, background intensity statistics are computed, the seventeenth percentile of pixel intensity is subtracted to suppress non-specific fluorescence, and a refined second run produces the final nuclear masks.

Notably, deep learning proved less suitable for the rarest target. Pax7-positive satellite cells are sparse—on average only about one per ten fibers in the validation dataset—and their signal is faint against the sarcoplasmic background. For these, the pipeline falls back on CellProfiler’s classical IdentifyPrimaryObjects module, using Otsu thresholding with size, circularity, and intensity filters tuned to the amplified Pax7 signal. The authors note the modular design means this component could be swapped to detect other sparse cell populations; they have already adapted the same pipeline to identify macrophages and atrophied fibers in an elderly clinical cohort.

A custom Python script built on OpenCV, pandas, and NumPy then fuses the separate segmentation masks into a single biological dataset. Each fiber mask is decomposed by its unique grayscale label, and every nucleus or satellite cell is assigned to the fiber with which it shares the greatest bitwise pixel overlap—a ‘greatest overlap’ rule that resolves objects straddling two fibers. Fiber area and perimeter are computed in physical units using a pixel-to-micron scaling factor, and fiber type is assigned automatically by measuring mean type I myosin heavy chain fluorescence within each fiber mask after rolling-ball background subtraction. The intensity threshold separating type I from type II fibers, set at a pixel value of twenty, was derived empirically from an analysis of 3,722 cells across twelve human biopsies. Output is a unified Excel workbook listing fiber ID, area, perimeter, nuclear count, satellite cell count, and typing for every cell, alongside diagnostic overlay images that allow researchers to audit the automated calls visually.

Validation against traditional manual analysis is where the pipeline earns its credentials. The team compared both methods head-to-head on 310 matched cells from a human biopsy. For myonuclei per cell, manual counting averaged 2.471 nuclei versus 2.474 for the pipeline, a difference of just -0.003 nuclei, with no statistically significant difference by paired t-test. Satellite cell counts likewise agreed, averaging 0.110 manually versus 0.119 automatically. Fiber typing was the strongest performer: the pipeline classified 49 type I and 261 type II fibers against the manual count of 48 and 262, with a correlation of r = 0.939 between methods. Cross-sectional area and perimeter correlated even more highly, at r = 0.959 and r = 0.935 respectively, though with systematic offsets that the authors explain candidly.

Those offsets trace to boundary conventions rather than measurement failure. Manual outlines were drawn with a ten-pixel brush slightly inside the laminin border, while the automated segmentation places its one-pixel outline on or marginally outside it. Because area and perimeter are computed from within the drawn boundary, this small difference inflates automated values—the pipeline averaged 3,532 versus 3,005 square micrometers for area, and 264 versus 249 micrometers for perimeter. Similar discrepancies between manual and automated fiber outlines have been reported before, and the authors stress that users simply need to apply consistent methodology. Other limitations include the mandatory image-cleaning step, strict file-naming requirements, and occasional detection failures on very sparse or dim fields of view, for which the authors suggest alternative Cellpose detection modes or the newer Cellpose3.

The practical payoff is throughput. Published power calculations indicate that comparisons of satellite cells per fiber require at least fifty type I and seventy-five type II fibers, and many labs have settled on a minimum of 150 cells per biopsy as a working threshold. A single human biopsy can yield well over a thousand fibers, and the new pipeline can process all of them, dramatically increasing sample representation while reducing both analysis time and operator dependence. Every cell receives a unique identifier, every segmentation decision is saved as a reviewable overlay, and the modular CellProfiler architecture invites adaptation to new antibodies, tissues, and imaging conditions. As interest grows in satellite cell dynamics across the human lifespan—from athletic performance to sarcopenia and rehabilitation—the authors argue this workflow offers muscle biologists a scalable, transparent, and reproducible foundation for turning microscope images into meaningful numbers.

Subject of Research: A semi-automated microscopy pipeline for quantifying satellite cells, myonuclei, and muscle fiber morphology in human skeletal muscle cross-sections

Article Title: A semi‐automated pipeline for quantitation of Pax7+, myonuclei, and cross‐sectional area by fiber type

Article References: Megowan, H. G., Luu, M., Shuaib, A., Augienello, K. B., Fries, A. C., Searcy, J., & Dreyer, H. C. (2026). A semi‐automated pipeline for quantitation of Pax7+, myonuclei, and cross‐sectional area by fiber type. Physiological Reports, 14(18), Article e71097. https://doi.org/10.14814/phy2.71097

Image Credits: AI Generated

DOI: 10.14814/phy2.71097

Keywords: skeletal muscle, satellite cells, Pax7, myonuclei, cross-sectional area, fiber type, image analysis, CellProfiler, Cellpose, immunohistochemistry, muscle biology, widefield microscopy