Abstract
Background:
Few studies have explored relationships between low back pain (LBP) and physical characteristics (physique, muscle strength, tightness, and flexibility) in female rhythmic gymnasts (RGs).
Hypothesis:
Due to their extreme flexibility, modifiable physical factors for LBP in RGs are not related to muscle tightness.
Study Design:
Cross-sectional study.
Level of Evidence:
Level 4.
Methods:
College female RGs (n = 95) were categorized into LBP and non-LBP groups based on questionnaires and orthopaedic surgeon interview. Physical assessments included whole-body bone mineral content and density, trunk lean body mass, and scoliosis presence by dual-energy X-ray absorptiometry scans, presence of abnormal findings and spine alignment by magnetic resonance imaging scans, hip muscle strength testing (flexion, extension, and abduction), range of motion (ROM), and flexibility testing. LBP and non-LBP groups were compared, and multivariate regression analysis performed.
Results:
RGs with LBP exhibited significantly longer practice time, lower hip flexion muscle strength on the nondominant side versus dominant side, lower active straight leg raise on the nondominant side versus dominant side, and lower hip external rotation (ER) ROM on the nondominant side versus dominant side. A history of LBP was also associated with current LBP. Multivariate logistic regression analysis revealed that a history of LBP (odds ratio [OR], 6.33; 95% CI, 1.56-25.62), differences in hip flexion strength (OR, 0.96; 95% CI, 1.00-1.08), and differences in hip ER ROM (OR, 1.12; 95% CI, 1.03-1.19) were factors associated with LBP.
Conclusion:
LBP was experienced by 27.3% RGs and associated not only with a history of LBP and a deficit in hip flexor strength on the nondominant side but also a deficit in hip ER ROM on the nondominant side.
Clinical Relevance:
Focusing on hip flexor strength and hip ER ROM on the nondominant side may be of value for LBP rehabilitation and prevention in RGs.
Rhythmic gymnastics is characterized by combining classical ballet movements with strength and fitness derived from artistic gymnastics and modern dance. Rhythmic gymnastics has a shorter history than artistic gymnastics, having been adopted as an official event for individual competitions at the 1984 Los Angeles Games and for team competitions at the 1996 Atlanta Games. Rhythmic gymnasts (RGs) perform programs with music that include intricate body maneuvers coordinated with manipulation of handheld apparatus, such as ropes, hoops, balls, clubs, and ribbons. These programs feature leaps, turns, balances, and acrobatic maneuvers, demanding exceptionally high flexibility and muscle strength. 30 Gymnastics, including rhythmic gymnastics, is characterized in specialized sports from an early period, with training typically commencing before adulthood.4,16,22,32 To attain elite status, RG athletes undergo gymnastics-specific interventions, such as extreme stretching, from an early age to acquire advanced flexibility and balance skills. 8
RG athletes have been reported to be at risk of low back pain (LBP) due repetitive hyperextension of the lumbar spine during competition. 29 Hutchinson 12 reported that 86% of elite-level RGs experience LBP during practice and competition, with LBP being the sole recorded injury necessitating time away from the sport. Cupisti et al 6 found that injuries occurred in 70% of RG athletes, with back injuries constituting 22.2%, after ankle and foot injuries (38.9%). Cugusi et al 5 summarized risk factors for LBP in female artistic gymnasts, identifying increased age, weight, and gymnastics practice time as factors associated with LBP. Although RGs experiencing more overuse injuries than artistic gymnasts and therefore would benefit from reducing RGs’ practice time, it is unclear whether this is practical. 28 Therefore, identifying modifiable factors associated with LBP in RGs is essential. Mizoguchi et al 24 reported that decreased flexibility of hip and shoulder flexors on the dominant side was associated with current LBP of high school volleyball players in case control setting. In a prospective study, Enoki et al 7 also found that lower passive hip flexion angle was an associated factor for chronic LBP in pole vaulters and decathletes. However, it remains unclear whether these factors apply to RGs with their heightened flexibility.
This study aimed to identify modifiable physical characteristics associated with LBP in RGs, particularly the range of motion (ROM) and muscle strength of the lumbar region and adjacent joints to the lumbar, which are clinically closely associated with LBP. Given their flexibility, we hypothesized that factors other than muscle tightness are associated with the characteristics of LBP in RG athletes.
Methods
Participants
This study included 95 collegiate female RGs from 3 universities. All participants had started rhythmic gymnastics before entering college and were competitive RGs who competed in national to international competitions. Detailed explanation of the experimental procedures was provided to the participants before obtaining written informed consent. The study received approval from the ethics committee of Waseda University (2022-094).
Assessment Procedures
The assessment comprised a questionnaire survey, dual-energy X-ray absorptiometry (DXA) scan, magnetic resonance imaging (MRI) scan, hip muscle strength tests, muscle tightness and flexibility tests, and spinal alignment evaluation, all conducted on the same day. Participants first completed the questionnaire, followed by measurements of height and weight (taken clothed and without shoes) using a stadiometer (YG200DN, Yagami Co) to the nearest 0.1 cm and a digital scale (BC622, TANITA Co) to the nearest 0.1 kg, respectively. The questionnaire survey covered age at beginning competition, age at menarche, training hours per day, days per week, and current and previous history of LBP. All participants reported no lumbar pain during testing.
DXA Measurement
DXA (Horizon, Hologic Inc) was used to measure whole-body bone mineral content and density as well as the regional lean body mass of the trunk. The participants were placed in the supine position so that the midline of the body coincided with the centerline of the examination table. The supine posture involved the following: the shoulder joint at 45°; the forearm pronated; the hip and knee fully extended; and the feet fixed in way that the toes were in contact and did not move during scanning. The trunk region was defined as extending from the first cervical vertebra to the femoral neck, excluding the arm region from the axillary fossa to the fingertips. 35 Scoliosis presence was determined from the DXA image by an orthopaedic surgeon with >40 years of experience, considering the low radiation dose following the method of Jamaludin et al. 15
MRI Scan and Image Acquisition
Spinal alignment and lumbar pathologies, such as lumbar spondylolysis, lumbosacral facet syndrome, and herniated discs, were assessed using short tau inversion recovery (STIR) and 3-dimensional (3-D) LAVA imaging on a 3T MRI scanner (SIGNA Premier 3.0T; GE Healthcare). STIR scanning parameters were as follows: repetition time/echo time (TR/TE), 1800/minimum (60) ms; slice thickness, 2.0 mm; field of view (FOV), 24.0 cm; matrix, 320 × 224; inversion time, 150 ms. The lumbar spine, including the lumbosacral vertebrae, was scanned using 3-D LAVA based on the previous study with the following parameters 2 : TR/TE, 4.5/2.1 ms; slice thickness, 0.80 mm; FOV, 24.0 cm; matrix, 384 × 38. Participants were supine during scanning, with their ankles in 90° dorsiflexion and knees extended using a fixture.
Hip Muscle Strength Testing
Maximal isometric hip muscle strength was measured for hip flexion, extension, and abduction using handheld dynamometry (HHD; Mobie, SAKAI Medicine). Two examiners, one a physical therapist with >5 years of experience operating the HHD and the other ensuring proper posture, conducted the test. Hip flexion strength was measured in a half-sitting position (adjusted for participant height), with the hip on the measurement side in flexion and the opposite hip and knee in extension (Figure 1a). Hip extension strength was performed in the prone position, with the hip on the measurement side in extension and the contralateral hip and knee in flexion (Figure 1b). Hip abduction strength was conducted in the lateral decubitus, with the hip and knee on the measurement side in extension and the contralateral hip and knee joints in mild flexion (Figure 1c). Participants were instructed to exert maximum force for 3 s, with the HHD device placed proximal to the knee joint with maximum force without leg and pelvis rotation. Maximal hip muscle torque was calculated as hip muscle strength by multiplying force by lever arm and normalized using body weight (Nm/kg): force (kg) × lever arm (m)/body mass (kg). Lever arm length was measured from the superior aspect of the greater trochanter to the lateral epicondyle on the DXA image.

Testing of physical characteristics. (a) Hip flexion muscle strength. (b) Hip extension muscle strength. (c) Hip abduction muscle strength. (d) Modified active SLR. (e) Thomas test. (f) Hip IR. (g) Hip ER. (h) Shoulder mobility test. (i) FFD. ER, external rotation; FFD, finger-floor distance; IR, internal rotation; SLR, straight leg raise.
Muscle Tightness and Flexibility Testing
Modified active and passive straight leg raise (SLR), hip extension, and hip external rotation (ER) and internal rotation (IR) were evaluated using a standard goniometer. Considering that RGs reach their maximum ROM in normal SLR, modified active/passive SLR was measured to assess hamstring tightness in RGs with excessive flexibility, with participants flexing the hip on the measurement side to maximum and actively/passively extending the knee from that position (Figure 1d). Hip extension ROM was measured by Thomas 21 test in which the participant positioned at the edge of the bed with the hip in extension on the measurement side (Figure 1e). Specifically, the participant supine on the testing table, actively and maximally flexing 1 hip joint, bringing the knee toward the chest to flatten the lumbar spine and stabilize the pelvis, which indicates the hip extension ROM. Hip IR and ER were performed with the hip in 90° of flexion (Figure 1f,g).
The shoulder mobility test, which includes a complex function of thoracic extension and shoulder flexion, required participants to rotate their upper limbs from front to back while holding a bar with both hands (Figure 1h). 23 Participants aimed to narrow the width between their hands holding the bar (100 cm in length and 2.5 cm in diameter), which had an attached scale for measuring the between-hands distance. The finger-floor distance (FFD) task involved participants standing on a platform, measuring the distance from the platform surface to the fingertips. FFD indicates that the higher the value shows the higher the flexibility. Measurements for shoulder mobility and FFD were collected in 0.1 cm increments.
Spinal Alignment Evaluation
The thoracic kyphosis angle, measured in a standing position by an examiner using an inclinometer (Bubble Inclinometer; Baseline Inc) over the spinous processes of the first and second thoracic spines and over the twelfth thoracic and first lumbar spines, was recorded once for each angle. The overall thoracic kyphosis angle was calculated by summing all angles. 27 The intratester and intertester reliability of this measurement was established (intraclass correlation coefficient, 0.95; 95% CI, 0.91-0.97). 20 Lumbar lordosis and sacral slope angles were obtained from the midsagittal plane in MRI scans between the superior edge of L1 and the inferior edge of L5 and between the superior edge of the sacrum and the horizontal line, respectively. 36
Assessment of LBP and MRI Findings
Inclusion criteria for LBP included selecting “yes” in the questionnaire survey to either of the questions “Have you often felt pain in your lower back during the past month?” or “Has the degree of pain interfered with your performance during competition (practice or games)?”, or responding “yes” in a direct interview conducted by an orthopaedic surgeon with >10 years of experience to questions regarding experiencing lumbar region pain. In addition, the presence or absence of abnormal findings in L1-L5 was recorded by the same orthopaedic surgeon from MRI scans, including STIR and 3-D LAVA mode. Determination of MRI scans was deidentified by a third party and the surgeon was blinded so as not to allow matching with other data.
Statistical Analysis
A statistical power analysis was conducted for sample size estimation. The study required a minimum of 18 athletes with LBP to achieve 90% power, an alpha of 0.05, based on a reference value of 20% prevalence of LBP as shown in previous studies.5,10 In addition, since this study is a cross-sectional design, we set the power higher than that of a prospective cohort study. 36
Descriptive statistics (means and standard deviations) were used to analyze the participants’ age, height, weight, arm span, arm span-height ratio, athletic careers, age at menarche, practice time, hip muscle strength, ROM, flexibility, and spinal alignment. Left and right hip muscle strength and ROM data were categorized by defining the dominant side (ie, the superior side for holding a 1-leg support posture). After confirming data normality, t tests were used to compare variables between the LBP and non-LBP groups. Effect size was calculated, with criteria for interpretation set at >0.2 (small), >0.6 (moderate), >1.2 (large), and >2.0 (very large). 11 The chi-square test was used to examine the relationships between current LBP and a history of LBP as well as scoliosis prevalence. Multivariate logistic regression, using variables with P < 0.05 from the 2-group comparisons, was employed to identify predictive variables for LBP. Only significant variables (P < 0.05) were retained in the model for each dataset. Models were developed using forward stepwise (conditional) regression. The significance was set at P < 0.05 in all cases, and data were analyzed using SPSS for Windows Version 28.0 (SPSS Inc).
Results
Among the 95 participants, (27.3%) were classified in the LBP group. The MRI findings determined by the orthopaedic surgeons are shown in Table 1. A comparison of basic characteristics between the LBP and non-LBP groups is presented in Table 2. The LBP group exhibited a significantly higher arm span-height ratio (1.03 ± 0.18 vs. 1.01 ± 0.02; P = 0.013) and longer practice time (6.7 ± 2.7 vs 5.7 ± 2.3 h/day; P = 0.05). In addition, a significant association was observed between the prevalence of LBP and a history of LBP (χ2 = 13.8; P < 0.001). No significant differences were noted in age, height, weight, bone mineral content, bone mineral density, trunk lean mass, athletic careers, age at menarche, or scoliosis between the 2 groups.
MRI findings in RGs with LBP
LBP, low back pain; MRI, magnetic resonance imaging; RGs, rhythmic gymnasts.
Characteristics of the LBP and non-LBP groups
LBP, low back pain; BMC, bone mineral content; BMD, bone mineral density. Bold P values indicate significance levels below 5%.
Table 3 presents the differences in hip muscle strength between the LBP and non-LBP group. The LBP group exhibited lower hip flexion muscle strength on the nondominant side (97.2 ± 18.1 vs. 105.7 ± 17.3 Nm/kg; P = 0.02), and the difference between the nondominant and dominant sides was greater compared with the non-LBP group (−9.1 ± 10.3 vs. −3.0 ± 16.3 Nm/kg; P = 0.04).
Differences in hip muscle strength between the LBP and non-LBP groups
Data are mean ± SD. LBP, low back pain. Differences were calculated by subtracting the nondominant side value from the dominant side value.
Differences in ROM, flexibility, and spinal alignment between the LBP and non-LBP groups are shown in Table 4. The LBP group had an increased difference between nondominant and dominant sides on modified active straight leg raise (m-ASLR) compared with the non-LBP group (−7.4° ± 11.2° vs −3.8° ± 7.7°; P = 0.04), a decreased hip ER ROM on the nondominant side compared with the dominant side (47.2° ± 5.3° vs 53.3° ± 7.7°; P < 0.001), and an increased difference between nondominant and dominant sides on hip ER ROM (−5.6° ± 6.1° vs −1.2° ± 6.6°; P = 0.002) than non-LBP group.
Differences in ROM, flexibility, and spine alignment between the LBP and non-LBP groups
Data are means ± SD. Bold P values indicate significance levels below 5%. Differences were calculated by subtracting the nondominant side value from the dominant side value. FFD, finger-floor distance; LBP, low back pain; m-ASLR, modified active straight leg raise; m-PSLR, modified passive straight leg raise; ROM, range of motion.
In multivariate logistic regression analysis (Table 5), 3 variables, namely LBP history (odds ratio [OR], 6.33; 95% CI, 1.56-25.62; P = 0.011), differences in hip flexion strength (OR, 1.04; 95% CI, 1.00-1.08; P = 0.040), and differences in hip ER ROM (OR, 1.12, 95% CI, 1.03-1.19; P = 0.013), emerged as the most significant factors associated with LBP in RGs.
Multivariate logistic regression analysis of variables associated with LBP
Model chi-square test (P < 0.001). Overall percentage correctly predicted, 78.3%. Pseudo-R2 (McFadden) = 0.406. Differences were calculated by subtracting the nondominant side value from the dominant side value. Bold P values indicate significance levels below 5%. ER, external rotation; LBP, low back pain; m-ASLR, modified active straight leg raise; OR, odds ratio.
Discussion
Joint hypermobility and lumbar spine hyperextension, which are specific physical characteristics, are recognized as related factors affecting LBP in RGs. However, limited findings on joint mobility, flexibility, muscle tightness, and muscle strength in RGs hinder identifying modifiable LBP risk factors. The present study aimed to determine the modifiable physical characteristics of RGs with LBP, with a focus on side-to-side differences. We hypothesized that characteristics other than flexibility are associated with LBP. Our results revealed associations with a history of LBP and side-to-side differences in both hip flexion strength and hip ER ROM.
In our study, LBP prevalence in RGs was 27.3%, with 57.7% of these participants exhibiting abnormal findings on MRI scans. The rate is slightly higher than the approximately 20% reported in previous studies,5,10 indicating a high LBP rate among RGs. Abnormal MRI scan findings were observed mostly in the lower lumbar spine, including lumbar spondylolysis at L5, intervertebral arthritis at L3/4, L4/5, and L5/S1, and lumbar disc herniation at L5/S1. Previous studies have shown that 90% of patients with LBP are diagnosed with nonspecific LBP, where the cause cannot be clinically identified. 18 In RGs, LBP is often associated with mechanical stress to the lumbar region. 12,30,33 In addition, a previous history of LBP was associated strongly with current LBP. Although the timing of lumbar spondylolysis found on MRI scans remains unclear, managing lumbar injuries during the growth period may be crucial for preventing LBP in collegiate RGs. Previous studies indicated that lumbar bone stress injury and spondylolysis are more likely to occur during skeletal immaturity,19,36 and a history of lumbar spondylolysis is an associated factor for LBP. 14
Frutuoso et al 9 reported that the nondominant leg predominately acts as the propulsion leg during most jumps in rhythmic gymnastics, indicating functional asymmetry. The physical characteristics of the LBP group had an increased difference between nondominant and dominant sides on m-ASLR, a decreased hip ER ROM on the nondominant side compared with the dominant side, and an increased difference between nondominant and dominant sides on hip ER ROM than the non-LBP group. Wu et al 37 also observed higher bone mass density at the proximal femur on the take-off side (ie, nondominant side) compared with the landing side (ie, dominant side) in female RGs. Given the higher load on the nondominant side, it is conceivable that the dysfunction of the nondominant side may be related to the LBP in RGs.
In the present study, the hip flexion muscle test involved extending the contralateral hip to place the hip on the measurement side in the final range of flexion. A previous study showed a significant contribution of the iliopsoas muscle with a greater hip flexion angle. 17 Therefore, RGs with LBP may have reduced iliopsoas muscle strength on the nondominant side. The iliopsoas muscle serves as a stabilizer linking the trunk, pelvis, and legs. 1 Bachrach 3 also suggested that weakness in the psoas muscles—part of the hip flexors—could lead to laxity and mechanical instability in the lumbosacral region, 13 stimulating nociceptors in the surrounding soft tissues and inducing pain. Modifying the function of the nondominant iliopsoas muscle is preferable for maintaining trunk stability in preventing and treating LBP in RGs. Notably, the peak torque of hip flexion in female RGs has been shown to be lower on the nondominant side compared with the dominant side, 9 as observed in the non-LBP group in the present study. Given the demonstrated association between hip flexion muscle strength and LBP, 26 understanding the degree of asymmetry in hip flexor strength is needed for LBP prevention and treatment in RGs.
The difference in hip ER ROM between nondominant and dominant side emerged as a factor associated with LBP in RGs. In addition, since the nondominant hip ER ROM of the LBP group was lower than that of the non-LBP group, this result would indicate a relative limitation of hip ER ROM on the nondominant side. Sadeghisani et al 31 identified a greater lumbopelvic rotation ROM during passive hip ER in patients with LBP regularly involved in rotation-related sports. In addition, Scholtes et al 34 suggested that patients with LBP who participate in rotation-related sports may exhibit greater and earlier lumbopelvic movement during active hip lateral rotation compared with those without LBP. Thus, limited hip ER ROM could potentially lead to compensatory lumbosacral movement. Therefore, given the limited hip ER ROM in in hip flexion, it would be recommended to assess and improve the tightness of the nondominant lateral hip external rotators, such as the gluteus maximus and deep external rotators, to address LBP in RG.
In the LBP group, the difference between the nondominant and dominant sides was more pronounced than that in the non-LBP group, indicating a relative lack of m-ASLR on the nondominant side. However, m-ASLR was not identified as an associated factor for LBP in the multivariate analysis. In addition, no difference in modified passive straight leg raise was observed between the LBP and non-LBP groups, confirming greater muscle flexibility in the posterior femoral muscles of RGs. Moltubakk et al 25 demonstrated that RGs with high hamstring flexibility had an advantage in flexion torque production in the knee-extended position compared with non-LBP group. In other words, m-ASLR may be a suitable assessment of hamstring flexibility in highly flexible RGs, indicating functional ROM reflecting muscle function in the terminal ROM of the joint angle. RGs are required to maintain a straight posture with straightened lower extremities for an attractive appearance during competition. Based on these observations, we speculate that the LBP group developed compensatory movements in the lumbar region due to the lack of m-ASLR.
This study had several limitations. First, given the cross-sectional nature of the study, we could not determine whether effects were due to LBP-associated disorders or characteristics of RGs who develop LBP. Second, LBP may not entirely align with abnormal findings in the lumbar region identified on MRI scans. In addition, as MRI coverage was limited to the first to the fifth lumbar spine, MRI-based findings for the thoracic spine, considered a potentially susceptible region in highly flexible RGs, could not be determined. This study provides valuable knowledge for RGs with high flexibility.
Conclusion
RGs had a high rate of LBP (27.3%), and LBP history was associated as a nonmodifiable parameter for LBP. In addition, deficiencies not only in nondominant hip flexor muscle strength but also in nondominant hip ER ROM were associated with LBP as a modifiable physical function. These results may suggest the usefulness of incorporating flexor strength exercises into training routines as well as implementing daily assessments and stretching for nondominant hip ER ROM for even highly flexible RGs to effectively manage LBP.
