Growing evidence suggests an association of lumbar paraspinal muscle morphology with low back pain (LBP) and lumbar pathologies. Unilateral spinal disorders provide unique models to study this association, with implications for diagnosis, prognosis, and management. Statistical shape analysis is a technique that can identify signature shape variations related to phenotypes but has never been employed in studying paraspinal muscle morphology. We present the first investigation using this technique to reveal disease-related paraspinal muscle asymmetry, using MRIs of patients with a single posterolateral disc herniation at the L5-S1 spinal level and unilateral leg pain. Statistical shape analysis was conducted to reveal disease- and phenotype-related morphological variations in the multifidus and erector spinae muscles at the level of herniation and the one below. With the analysis, shape variations associated with disc herniation were identified in the multifidus on the painful side at the level below the pathology while no pathology-related asymmetry in cross-sectional area (CSA) and fatty infiltration was found in either muscle. The results demonstrate higher sensitivity and spatial specificity for the technique than typical CSA and fatty infiltration measures. Statistical shape analysis holds promise in studying paraspinal muscle morphology to improve our understanding of LBP and various lumbar pathologies.