The effects of 12 weeks of static stretch training on the functional, mechanical, and architectural characteristics of the triceps surae muscle–tendon complex

Purpose We investigated the effects of 12 weeks of passive static stretching training (PST) on force-generating capacity, passive stiffness, muscle architecture of plantarflexor muscles. Methods Thirty healthy adults participated in the study. Fifteen participants (STR, 6 women, 9 men) underwent 12-week plantarflexor muscles PST [(5 × 45 s-on/15 s-off) × 2exercises] × 5times/week (duration: 2250 s/week), while 15 participants (CTRL, 6 women, 9 men) served as control (no PST). Range of motion (ROM), maximum passive resistive torque (PRTmax), triceps surae architecture [fascicle length, fascicle angle, and thickness], passive stiffness [muscle–tendon complex (MTC) and muscle stiffness], and plantarflexors maximun force-generating capacity variables (maximum voluntary contraction, maximum muscle activation, rate of torque development, electromechanical delay) were calculated Pre, at the 6th (Wk6), and the 12th week (Wk12) of the protocol in both groups. Results Compared to Pre, STR ROM increased (P < 0.05) at Wk6 (8%) and Wk12 (23%). PRTmax increased at Wk12 (30%, P < 0.05), while MTC stiffness decreased (16%, P < 0.05). Muscle stiffness decreased (P < 0.05) at Wk6 (11%) and Wk12 (16%). No changes in triceps surae architecture and plantarflexors maximum force-generating capacity variables were found in STR (P > 0.05). Percentage changes in ROM correlated with percentage changes in PRTmax (ρ = 0.62, P = 0.01) and MTC stiffness (ρ = − 0.78, P = 0.001). In CTRL, no changes (P > 0.05) occurred in any variables at any time point. Conclusion The expected long-term PST-induced changes in ROM were associated with modifications in the whole passive mechanical properties of the ankle joint, while maximum force-generating capacity characteristics were preserved. 12 weeks of PST do not seem a sufficient stimulus to induce triceps surae architectural changes.


Introduction
Passive stretching is widely performed in sport and rehabilitation mainly to improve joint range of motion (ROM) and muscle performance. However, evidence exists that passive stretching, when performed acutely, may induce negative changes in muscle function, such as a reduced force-generating capacity (Power et al. 2004;Kay and Blazevich 2012;Longo et al. 2014;Behm et al. 2016;Trajano et al. 2017;Cè et al. 2020) and a depressed rate of torque development (RTD) (Simic et al. 2013;Trajano et al. 2019). These alterations are often accompanied by a reduction in the amplitude of the surface electromyographic (sEMG) signal from the contracting muscle after stretching (Behm et al. 2001(Behm et al. , 2016Cramer et al. 2005). These impairments may be ascribed to both neuromuscular, such as a reduced activation (Behm et al. 2016;Trajano et al. 2017), and mechanical mechanisms, such as alteration in the viscoelastic properties of the muscle-tendon complex (MTC) (Magnusson et al. 1996;Morse et al. 2008;Longo et al. 2014).
Whether or not a long-term (chronic) passive stretching training (PST) could affect muscle function still remains controversial. On one side, indeed, some studies reported an increase in maximum voluntary dynamic (Kokkonen et al. 2007;Nelson et al. 2012) and isometric (LaRoche et al. 2008) contraction, as well as a reduction in whole-joint (Kubo et al. 2002;Guissard and Duchateau 2004;Nakamura et al. 2012) and muscle stiffness (Blazevich et al. 2014;Nakamura et al. 2020) after chronic stretching. On the other side, other works reported no effects of PST on maximum muscle strength (Akagi and Takahashi 2014;Konrad and Tilp 2014;Blazevich et al. 2014;Sato et al. 2020), amplitude of sEMG detected during maximum voluntary force production (Blazevich et al. 2014) and whole-joint (Konrad and Tilp 2014;Blazevich et al. 2014) and muscle stiffness (Konrad and Tilp 2014). This apparent discrepancy could be justified by differences in methodological approach, PST duration, number and duration of weekly stretching sessions, and stretch intensity (Freitas et al. 2018).
After PST, architectural adaptations (fascicle length, fascicle angle, and muscle thickness) in the stretched muscle are still a matter of debate (Medeiros and Lima 2017;Nunes et al. 2020). Since changes in these variables would influence muscle contraction characteristics (Narici et al. 2016), it may be of great interest to investigate the potential PSTinduced architectural adaptations. Animal studies showed that chronically stretched muscles undergo an increase in fibre size/muscle mass and fibre length (Sola et al. 1973;Holly et al. 1980; Barnett et al. 1980), possibly due to mechanical signalling-induced increase in protein synthesis and addition of series sarcomeres (Goldspink et al. 1995(Goldspink et al. , 2002. Moreover, regional differences seem to occur in the chronically stretched muscle, since a different remodelling was observed between its middle and distal portions (Dix and Eisenberg 1990). However, the architectural responses to PST in humans are not clear (Nakamura et al. 2012;Akagi and Takahashi 2014;Konrad and Tilp 2014;Blazevich et al. 2014;Freitas and Mil-Homens 2015;Simpson et al. 2017;Sato et al. 2020;Beltrão et al. 2020) and more studies may be required to evaluate the impact of a long-term PST program on muscle architecture (Medeiros and Lima 2017).
Therefore, the aims of this study was to assess whether or not changes in ROM induced by long-term PST are accompanied by changes in muscle force-generating capacity, passive whole-joint and muscle stiffness, and architectural MTC characteristics. The possible correlations between changes in ROM and changes in force-generating capacity, passive stiffness, and architectural variables were also examined. Hypothesis was made that, with sufficient training volume, PST could reduce passive stiffness while increasing fascicle length and muscle thickness. Therefore, the reduction in passive stiffness would impair force transmission to the tendon insertion point, thus reducing the muscle force-generating capacity and RTD. Nevertheless, the hypothesised modifications in architectural characteristics would lead to an increase in force production. The net effect of these two balancing mechanisms would result into unchanged maximum force expression.

Participants
Based on pilot testing, the sample size was computed using statistical software (G-Power 3.1, Dusseldorf, Germany) expecting at least a moderate Cohen's d effect size (0.60) 1 3 in ROM changes. Considering α = 0.05 and a required power (1 − β) = 0.80, the desired sample size resulted in 24 participants. To ensure sufficient statistical power, 30 healthy volunteers (12 women) were enrolled in the study (mean ± standard deviation: age 22.7 ± 1.8 years, body mass 68.4 ± 9.4 kg, stature 1.74 ± 0.1 m), and were randomly assigned to a stretching training group (STR; N = 15, 6 women, 9 men; age 22.3 ± 0.8 years, body mass 68.5 ± 9.4 kg, stature 1.74 ± 0.08 m) and a control group (CTRL: N = 15, 6 women, 9 men; age 23.4 ± 0.8 years, body mass 67.4 ± 9.5 kg, stature 1.73 ± 0.08 m) with the same number of women in each group. Each participant received a full explanation of the aim of the study, the experimental procedures, and signed a written informed consent. The volunteers were recreationally active university students engaged in regular sports activities (2.0 ± 1.0 h/week). Exclusion criteria were: (i) the presence of musculoskeletal injury within the past 6 months; (ii) neurological deficit affecting their ability to stretch; (iii) being regularly involved in stretching training; and (iv) the presence of all those circumstances in which stretching training is contraindicated (e.g., joint and tissue laxity). Participants were asked to abstain from ergogenic beverages or similar in the 24 h preceding the test sessions and to report to the laboratory without any form of heavy intensity physical exercise in the previous 48 h. The study was approved by the local University Ethical Committee (CE 27/17) and had been performed in accordance with the latest principles of the Declaration of Helsinki.

Experimental design
A randomised pre-post parallel group design was adopted to study the effects of a 12-week PST programme on the functional, mechanical and architectural properties of the triceps surae MTC. The study lasted 15 weeks. In the first two weeks, participants underwent two familiarisation sessions interspersed by at least 48 h to get acquainted with all testing procedures. A third session was attended to collect baseline data (Pre). Subsequently, participants enrolled in the STR group were instructed about the exercises to be performed on the dominant limb (right for all participants) included in the training programme and familiarised with the exercise intensity. Thereafter, STR started the 12 weeks of PST, whereas participants in CTRL continued their habitual activities. All individuals were tested at the 6th week (Wk6) and at the 12th week (Wk12) of the protocol. In males, the last evaluation occurred within a week from the end of the PST in STR. To minimise the effects of the menstrual cycle on the assessments, female participants recorded their menstrual cycle in a personal diary at the beginning and throughout the study. Knowing the occurrence of the early follicular phase allowed women to be tested around the same menstrual days (3 ± 3 days from the early follicular phase). However, this choice implied that tests may have occurred with a ± 5 day-dispersion from the exact testing week at Wk6 evaluation .
All experiments were carried out in a room at constant temperature (22 ± 1 °C) and relative humidity (50 ± 5%). During tests, participants lay prone on a custom-made ergometer ( Fig. 1) (Longo et al. 2017), with a mobile metal platfrom for consistency with Fig. 1 connected to a previously calibrated load cell (mod. SM-2000 N, Interface, UK; operating linearly between 0 and 2000 N). The ankle of the dominant limb was firmly attached to the mobile metal plate by a Velcro ® strap (Velcro Industries Inc., Willemstad, Netherlands Antilles) to minimize heel displacement during assessments. Hip and shoulders were also firmly secured to the ergometer. The load cell was constantly kept in line with the axis of force. Force signal was driven to an A/D converter (mod. UM 150, Biopac, Biopac System Inc., Santa Barbara, CA, USA), sampled at 10,240 Hz, directed to an auxiliary input of the electromyography amplifier (mod. EMG-USB, OtBioelettronica, Turin, Italy) and stored on a personal computer. Torque was calculated by multiplying the force output by the distance between the apical aspect of the external malleolus and the force application point. A previously calibrated bi-axial angle transducer (mod. TSD 130A, Biopac System, CA, USA) was positioned on the external face of the fibula and on the calcaneum to monitor the changes in ankle ROM.
First, ultrasound images of the triceps surae muscle were obtained at rest with the ankle joint fixed at 0° (perpendicularity between the tibia and the longitudinal axis of the foot). Thereafter, maximum ROM together with passive force exerted during the manoeuvre was evaluated. Passive MTC mechanical properties were subsequently assessed by combining passive force and ultrasound techniques. Lastly, the foot was re-positioned at 0° to determine the peak isometric torque. Myoelectric activity was detected during passive movement and muscle contraction.

Ultrasound imaging
Since previous studies showed that the different triceps surae muscles can undergo different strain during stretch exercise (Hirata et al. 2016) and within the same muscle different portions can be unequally affected by the stretch manoeuvre (Andrade et al. 2020), we examined the possible PST-induced architectural adaptations in all triceps surae heads at different regions.
Muscle architecture was assessed in vivo at rest for gastrocnemius medialis (GM), gastrocnemius lateralis (GL) and soleus (SOL) muscles by B-mode ultrasound (LOGIQS7, 1 3 GE©, Fairfield, Connecticut, USA) with a 5-cm linear-array probe (mod. 9L, 3.1-10.0 MHz) in extended-field-of-view mode (LOGIQ view ). The participants lay prone on the examination bed with hip and knee joint extended and the ankle fixed at 0°. Transmission gel was applied to improve acoustic coupling. All muscles were inspected before extendedfield-of-view acquisition. The proximal and distal ends (i.e., myotendinous junction -MTJ) of each head of the triceps surae were identified by moving the ultrasound probe along the longitudinal axis of each muscle belly. Thereafter, each muscle was marked on the skin at 50% length. At this location, medial and lateral muscle boundaries were also identified for each muscle, and the 50% width was drawn on the skin. A single snapshot was collected at this site for each muscle to obtain muscle thickness (MT). Subsequently, after ultrasound inspection, the medial and lateral GM and GL muscle boundaries were drawn from their respective MTJ along muscle length. Then, the probe was moved from each MTJ towards the proximal end of the respective muscle belly. When necessary, transducer manipulation occurred so that the fascicles and both aponeuroses remained continuous and visible. This path was assumed to be the best fascicle plane that could be followed, and it was marked on the skin (Franchi et al. 2020). After the inspection, a continuous single view was taken for GM and GL heads starting from the MTJ and by moving slowly the probe along the drawn line ensuring that the extended-field-of-view image exceeded the 50% muscle length. Due to its short fascicles, SOL images were captured as single snapshots at the MTJ and mid-belly sites. The expert operator performing ultrasound scanning ensured minimal pressure was applied.
The images were analysed offline using an open-source computer program (ImageJ 1.44b, National Institutes of Health, USA). For each muscle, three clearly visible muscle fascicles were identified at 50% muscle length (MID) and close to the MTJ (DIST-i.e., about 2 cm along the deep aponeurosis from the MTJ) (Fig. 2a). Muscle fascicle length (L f ) was measured by drawing a line along the three fascicles between the deep and superficial aponeurosis. Any fascicle curvature was taken into account. On the same highlighted fascicles, their insertion angle into the deep aponeurosis was measured as fascicle angle (θ). The three measured fascicles and angles were averaged and used for the analysis. MT of each triceps surae head was measured as the distance between the superficial and deep aponeurosis at three points where the muscle belly was the widest and the two aponeuroses parallel (Fig. 2b) from the single snapshot taken at MID. These three measures were then averaged.

ROM and maximum passive resisting torque
As described in the "Experimental design" section, the ankle was securely fixed to the custom-made ergometer for ROM and maximum passive resisting torque (PRT max ) detection. After 10 conditioning passive ankle movements, the joint ROM was determined starting with the ankle at its resting position (~ − 20° of dorsiflexion) and manually dorsiflexed at slow speed to avoid reflex activations monitored by surface electromyography (sEMG) until the participant's point of discomfort was reached. At maximum ROM, the ankle joint was mechanically fixed for the recording of PRT max , which was considered as index of stretch tolerance .
The difference between the ankle at 0° and the end ROM was considered as joint ROM. Reliability was previously reported (intraclass correlation coefficient-ICC = 0.94; standard error of measurement as percentage-SEM% = 1.1) (Longo et al. 2014). PRT max was measured within a 250-ms epoch ).

Passive MTC and GM muscle stiffness
To allow measurements of passive force and GM MTJ position at different joint angles, the metal plate was manually fixed at 0°, 10°, 20° of ankle dorsiflexion, and at end ROM (Longo et al. 2014). Joint positioning was executed at slow speed to avoid reflex muscle activation (monitored by sEMG). PRT exerted by plantarflexor muscles was recorded at each angle as the average of force values during the first 5 s after ankle positioning. Such a short time allowed the operator to minimize the influence of the static position on the MTC viscoelastic properties. The passive torque-angle curve between 0° and 20° of dorsiflexion was fitted with the best polynomial regression model (Mizuno et al. 2013;Longo et al. 2014), and the slope of this curve at 20° of dorsiflexion (maximum common angle for all participants) represented passive MTC stiffness.
While passive MTC stiffness was assessed, GM MTJ displacement was also captured by means of B-mode ultrasound imaging using a 5-cm linear-array probe. A hypoechoic tape was applied on the skin close to the MTJ for monitoring possible shifting during measurement. GM MTJ was visualised in a continuous sagittal plane at 0°, 10°, 20° of dorsiflexion, at end ROM, and digitised offline relative to the hypoechoic marker. As previously described (Morse et al. 2008;Cè et al. 2015), possible artefacts due to skin shift during MTJ displacement assessment were corrected accordingly.
GM muscle stiffness was calculated by dividing the changes in PRT between 0° and 20° by the corresponding MTJ displacement. MTJ displacement between 0° and end ROM (ΔMTJ max ) represented maximum GM elongation.

Maximum voluntary contraction
After a standardized warm-up (10 × 2-s contractions of increasing intensity from 50% maximum voluntary contraction determined during familiarization up to maximum), plantarflexors maximum voluntary isometric contraction (MVC) was assessed at 0° (anatomical neutral position). Participants were instructed to contract "as fast and hard as possible" (Maffiuletti et al. 2016), to maintain the contraction for 3 s and then to relax. Two trials were performed separated by 3 min of rest. In case of a between-trial difference > 5%, a third trial was executed. The best trial was used as MVC.

Electromyographic assessment
sEMG signal was collected during MVC measurements in GM, GL and SOL. The myoelectric activity was detected by a linear array of eight electrodes (mod. KITAD008, OtBioelettronica, Turin, Italy; probe: 45 mm × 20 mm; electrode length: 2 mm; inter-electrode distance: 5 mm) fixed to the skin by dual-adhesive foams (mod. AD004, OtBioelettronica, Turin, Italy) and filled with conductive gel (Cogel, Comedical, Trento, Italy). The skin area under the sEMG electrodes was cleaned with ethyl alcohol, abraded gently with fine sandpaper and prepared with a conductive cream (Nuprep, Weaver and Co., Aurora, USA) to achieve an inter-electrode impedance below 2000 Ω. For each muscle, the sEMG array was placed over the muscle belly along the direction of the muscle fibres, in accordance with the European recommendations for surface EMG (Hermens et al. 1999). sEMG was acquired by a multichannel amplifier with a sampling rate of 10,240 Hz (mod. EMG-USB, OtBioelettronica, Turin, Italy; input impedance: > 90 MΩ; CMRR: > 96 dB), amplified (gain × 1000) and filtered (filter type: 4th order Butterworth filter; bandwidth: 10-500 Hz).
The analysis was performed by OtBiolab+ software (OtBioelettronica, Turin, Italy). The sEMG signal epochs were aligned with the same force signal epochs. The sEMG signal of each muscle was analysed in time domain within the same 1-s period detected in the middle of the MVC plateau. The sEMG root mean square (RMS) was calculated in consecutive 250-ms time windows and then averaged.

RTD
Rate of torque development was measured using the raw force trace while assessing MVC. A threshold of three standard deviations above the baseline signal for three consecutive points within a 100-ms interval of the resting condition immediately preceding the contraction was used as the onset of force signal. Due to large variability in the early time window, rate of torque development was calculated as the average change in torque per time interval from the torque onset to 200 ms (RTD 200 ) (Maffiuletti et al. 2016;Trajano et al. 2019). Reliability for RTD has been previously reported (ICC = 0.91, SEM% = 6.5) (Longo et al. 2014).

Electromechanical delay
The electromechanical delay (EMD) was calculated as the time between the onset of the band-passed, rectified GM sEMG signal and torque onset during MVC (Longo et al. 2017). The same criteria for torque onset identification (three standard deviations above the baseline signal preceding the contraction) were used for detecting sEMG onset. Reliability analysis for EMD was reported previously (ICC = 0.97, SEM% = 1.1) (Longo et al. 2017).

PST programme
The PST programme consisted of 12 weeks of training, 5 sessions per week (60 sessions in total). This choice was based on recent studies pointing out that investigations providing the participants with an adequate training volume (i.e., weekly frequency > 3 times/week, total duration > 8 weeks) were needed to clarify the effects of PST on the functional, mechanical and architectural characteristics of the MTC (Freitas et al. 2018;Nakamura et al. 2020;Andrade et al. 2020). Each session included two manoeuvres for plantarflexor muscles (Fig. 3), 45 s elongation and 15 s recovery in the starting position repeated for five times (Longo et al. 2014). The exercises were performed on the right limb and stretching duration was 450 s each session, giving a total stretching time of 2250 s per week. For the first exercise, the participants stood erect with arms supporting the body and the right foot in dorsiflexion against a board (Fig. 3a). For the second exercise, the participants lay supine starting with the right knee straightened and an elastic band around the right foot positioned at maximum dorsiflexion. Using the arms while holding the elastic band, the right hip was flexed until maximum point of discomfort (Fig. 3b). The exercises were chosen to represent a stretch modality commonly performed in recreational (e.g., gyms) and sports activities. As stated above, the CTRL group received no PST. To promote participants' compliance, daily classes were held at different day time (morning and afternoon) at the University Sports Centre gym. Each class was supervised by an expert operator, which monitored the attendance, the correct exercise execution, and the intensity exerted during the exercise. Participants were instructed to stretch around the maximally tolerable stretch within the pain limit (Blazevich et al. 2014;Cè et al. 2020). Individuals not attending at least the 80% of classes were excluded from the study, and new participants were recruited to substitute the drop out.

Statistical analysis
Statistical analysis was performed using a statistical software package (IBM SPSS Statistics v. 26, Armonk, NY, USA). The Shapiro-Wilk test was used to check the normal distribution of the sampling. To determine possible baseline differences between STR and CTRL before intervention, the unpaired Student's t test was applied for each variable of interest. A two-way mixed-model analysis of variance (ANOVA) [withingroup factor: time, 3 levels (Pre, Wk6, Wk12); between-groups factor: group, 2 levels (STR, CTRL)] was used to check for differences between groups over time. To assess betweengroups differences in the changes over time, an analysis of co-variance (ANCOVA) was applied, using values at Pre as covariate. Multiple comparisons were perfomed applying the Bonferroni's correction. For every variable of interest, the percentage variation was calculated for each participant at each time point. Thereafter, the mean percentage difference of the sample with 95% confidence interval (95% CI) was calculated. Since the Pre-Wk12 percentage changes in PRT max and GM muscle stiffness of the STR group were not normally distributed, the possible correlations between percentage changes in ROM and PRT max , MTC stiffness, and GM muscle stiffness were assessed by Spearman's Rho coefficient.

Participants' compliance
Attendance was about 90% (54/60 training sessions). Four participants dropped out because of personal reasons not linked to the study. They were immediately replaced to maintain the sample size.

Muscle architecture and thickness
Mean values of all parameters are presented in Table 2. No significant differences were found between groups in any variable at Pre (t-range = − 1.49-0.99, P-range = 0.14-0.88, g-range = 0.03-0.54, trivial-small). The ANOVA did not reveal any effect of time in any muscle or location (P-range = 0.12-0.98, small-medium) or interaction (P = 0.10-0.91, small-medium) ( Table 2).

Discussion
The main results of this investigation were that the expected PST-induced increase in ROM and PRT max was accompanied by a decrease in GM muscle and whole MTC stiffness, in line with the experimental hypothesis. The percentage changes in ROM at the end of the training protocol correlated with those in PRT max and MTC but not with GM muscle stiffness. Therefore, an adequate PST volume can induce whole-joint passive mechanical alterations, with GM muscle stiffness playing only a limited role in PST-induced ROM adaptations. Contrary to the experimental hypothesis, though, PST did not induce alterations in muscle architecture at any location, suggesting that this training protocol was not sufficient to induce structural adaptations detectable by ultrasound imaging technique. Interestingly, plantarflexor muscles force-generating capacity variables (MVC, sEMG RMS, RTD 200 , EMD) were not affected by PST, indicating that despite the reduction in stiffness, this training protocol did not compromise the mechanisms underpinning force transmission during MVC.

Preliminary considerations
As expected, PST induced an increase in ankle ROM by ~ 9% at Wk6 and ~ 24% at Wk12 compared to Pre. Concomitantly, PRT max (an index of stretch tolerance) increased at both time points despite being significantly higher only at Wk12 (~ 30%). A correlation was found between the percentage changes in ROM and PRT max at Wk12, suggesting that one of the factors explaining the increase in ankle ROM could be attributable to the increase in stretch tolerance, in line with previous findings (Guissard and Duchateau 2004;Blazevich et al. 2014;Kay et al. 2016;Nakamura et al. 2017). Although the exact origin of this phenomenon still needs to be fully clarified, the increase in stretch tolerance could be ascribed to a change in the afferent input from nociceptive nerve endings and mechanoreceptors, as well as to the increased participants' willingness to tolerate the greater loading at ROM endpoint (Law et al. 2009;Weppler and Magnusson 2010). week (Wk12) of study protocol. § § Significantly different from Pre with P < 0.01; *significantly different between groups with P < 0.05; **significantly different between groups with P < 0.01 ▸ 1 3

Neuromuscular activation and force-generating capacity variables after PST
Twelve weeks of PST did not alter MVC, sEMG RMS, RTD 200 , and EMD. Although the administration of one bout of passive stretching led to acute changes in these parameters (Longo et al. 2014(Longo et al. , 2017, the same was not observed after long-term stretching training. Interestingly, despite the PST-induced reduction in passive stiffness, the active neuromuscular and mechanical behaviour of the contracting muscle during MVC was not influenced by passive mechanical changes. Indeed, as hypothesised, a reduction in MTC and GM muscle stiffness would have implicated a delayed force transmission from the muscle to the tendon insertion point, a slower RTD, and, possibly, a reduction in MVC (Longo et al. 2017). However, the present results are not consistent with this potential occurrence. A possible explanation could be provided by two concurring mechanisms: unaltered motor drive towards the muscles involved in PST, and possible unchanged PST-induced tendon stiffness during an active contraction, which are consistent with similar sEMG RMS, RTD 200 and EMD during MVC after intervention compared to baseline. Similarly, previous studies found no PST-induced changes in maximum muscle strength (Akagi and Takahashi 2014;Konrad and Tilp 2014;Blazevich et al. 2014;Sato et al. 2020) and RTD (Guissard and Duchateau 2004;Blazevich et al. 2014), supporting our outcomes of no effects of passive mechanical alterations on active muscle mechanical characteristics. However, one investigation found a lengthening in EMD during MVC after 8 weeks of knee flexors PST (Minshull et al. 2014). Therefore, further studies are needed to clarify whether or not potential functional PST-induced adaptations are muscle-group-dependent.

Passive MTC and muscle mechanical properties after PST
Besides stretch tolerance, the PST-induced increase in ROM was accompanied by a decrease in MTC stiffness at Wk12 (~ 17%), an increase in GM muscle elongation (~ 26% at Wk6 and ~ 42% at Wk12), and a decrease in GM muscle stiffness (~ 11% at Wk6 and ~ 16% at Wk12). These results suggest that the increase in ROM could be associated with changes in passive mechanical properties of the structures and tissues surrounding the ankle joint, as well as GM muscle mechanical characteristics. Furthermore, a correlation was found only between the percentage changes in ROM and MTC stiffness after PST, suggesting that the increase in ROM could be explained by a reduction in whole-joint stiffness. Nonetheless, despite a similar time course in ROM and GM muscle stiffness alterations, the lack of correlation between percentage modifications in these two variables highlights that changes in the muscle component of stiffness cannot be the main factor explaining ROM changes over time.
The present results are in accordance with previous studies, in which a decrease in MTC (Kubo et al. 2002;Guissard and Duchateau 2004;Nakamura et al. 2012) and GM muscle stiffness (Blazevich et al. 2014;Nakamura et al. 2017Nakamura et al. , 2020 have been observed after a PST program, and a correlation between percentage changes in ROM and MTC stiffness has been found (Guissard and Duchateau 2004). In contrast, other studies did not find changes in triceps surae MTC (Konrad and Tilp 2014;Blazevich et al. 2014) and GM muscle stiffness (Konrad and Tilp 2014). The discrepancy can be attributable to different training modality, such as intensity, frequency, number of exercises, and overall duration (Freitas et al. 2018). Concerning GM muscle stiffness, another possible explanation could be that PST affected differently the stiffness of the individual triceps surae muscles and/or the regional stiffness within each muscle, as observed by shear-wave ultrasound imaging technique (Andrade et al. 2020). Nonetheless, it seems that combining several weeks of training with high weekly exercise frequency, as it has been done in the present study, can lead to MTC and GM muscle stiffness adaptations. Interestingly, between Wk6 and Wk12, we observed a non-significant decrease in muscle stiffness. Hence, it can be speculated that after Wk6, the changes in ROM could be more related to changes in other structures than muscle stiffness.
MTC stiffness could be related to the intrinsic stiffness of muscles, tendons, and connective tissues surrounding the whole MTC and joint (e.g., fascia, ligaments, joint capsule, bursa) and/or to neural mechanisms (Guissard and Duchateau 2004;Nakamura et al. 2012). Previous studies suggested that within the muscle component, a PST programme can change the compliance of the intra-and extracellular structures involved in the passive tension generation during a stretch manoeuvre, particularly: (i) resting filamentary tension, caused by the stretching of the stable cross-link between the actin and myosin filaments (Proske and Morgan 1999); (ii) stretching of the no-contractile proteins, such as titin and desmin proteins (Magid and Law 1985;Trombitás et al. 1998); and (iii) deformation of the connective tissues located within and surrounding the muscle (in particular, the perimysium) (Borg and Caulfield 1980;Rowe 1981). The present results in GM muscle stiffness suggest that mechanical modifications occurred at muscle level possibly involving one or more of the above-mentioned contributors. In contrast, the influence of long-term PST on the tendon component is still controversial. Indeed, as summarised in a recent review (Freitas et al. 2018), previous studies using different PST approaches did not find changes in tendon passive mechanical properties. Finally, despite not being measured in the present investigation, PST-induced neural adaptations explaining changes in MTC stiffness cannot be excluded (Guissard and Duchateau 2004).

Triceps Surae architecture after PST
Contrary to our hypothesis, the present results demonstrated a lack of architectural adaptations to 12 weeks of PST at all regional levels. The rationale of changes in muscle size and architecture (i.e., longitudinal muscle growth and/or fascicle elongation) was based on several animal studies indicating that joint immobilization in a lengthened position (i.e., stretch) for long periods led to muscle growth and increase in muscle length due to sarcomere addition (Goldspink 1977;Goldspink et al. 1995). Whether or not these adaptations can be obtained also in humans is still controversial. Indeed, Simpson et al. (2017) reported an increase in muscle thickness and GM and GL fascicle length at both mid-belly and MTJ locations after 6 weeks of overloaded static stretching training performed 5 times a week. Similarly, Andrade et al. (2020) observed an increase in mid-belly GM fascicle length after 12 weeks of PST performed 5 times a week, involving similar exercises and training protocol than as used in the present investigation. In contrast, albeit shorter duration (i.e., 3-6 weeks), other studies did not report such adaptations (Nakamura et al. 2012;Konrad and Tilp 2014;Blazevich et al. 2014;Sato et al. 2020). Compared to the study by Simpson et al. (2017), the use of unloaded stretch exercises in the present investigation can explain the differences between studies in the architectural adaptations. However, we do not have a clear explanation for the discrepancy between the present study and the one by Andrade et al. (2020). Hypothetically, this could be ascribed to the use of different ultrasound approaches for determining muscle architecture (i.e., extended field of view versus B-mode ultrasound, respectively). Nonetheless, the present lack of architectural changes extends the results of two recent reviews highlighting the absence of structural and architectural adaptations induced by unloaded static stretching training up to 8 weeks (Freitas et al. 2018;Nunes et al. 2020).
As final consideration, in light of previous studies (Simpson et al. 2017;Andrade et al. 2020), it can be speculated that combining overloaded stretching training with high training volume could lead to possible structural adaptations in the stretched muscles at rest. Future studies are needed to explore this hypothesis.

Study limitations
This study comes with some limitations. First, the studied population was fairly homogenous (healthy, physically active adults); therefore, generalisation to other populations is limited. Second, possible differences between sexes in the response to PST could not be investigated due to the small sample size; future studies are needed to explore this possible occurrence. Third, it has been recently pointed out that ROM changes can be also explained by nerve stiffness changes due to stretching intervention (Andrade et al. 2020). However, due to the lack of shear-wave elastography-equipped ultrasound, this variable was not measured. Last, the assessment of muscle activity during stretching sessions would have provided more information about the level of muscle activation during elongation. Nonetheless, all of our participants were instructed to reach their own perceived level of maximum stretch within the pain limit, which led to an increase in ROM over the training period.

Conclusion
The present findings demonstrated that with 12 weeks of PST, the expected increase in ROM and stretch tolerance was accompanied by passive stiffness reduction at both MTC and GM muscle level. Moreover, changes in muscle stiffness occurred already at an earlier stage (Wk6). However, changes in ROM were explained by changes in stretch tolerance and whole MTC stiffness, but not by changes in GM muscle stiffness, indicating that other tissues contributing to the entire joint stiffness could play a major role in determining joint ROM variations. Interestingly, no changes in muscle force-generating capacity variables and architectural features at rest were found. These results provide evidence that the stretching protocol used in the present study induced modifications in the passive mechanical properties of the ankle joint without compromising the plantarflexor muscles force-generating capacity. As an important implication, practitioners can therefore utilize stretching routines chronically without impairing sport or exercise maximum performance. Future studies are needed to assess whether or not PST of longer duration and/or with different protocols could induce resting architectural adaptations.

Declarations
The Authors would like to thank all the volunteers that participated in the study.
Funding Open access funding provided by Università degli Studi di Milano within the CRUI-CARE Agreement. The present study was funded by a dedicated grant (PSR Fondi Linea 2-Tipologia A-2016 CE) from the Department of Biomedical Sciences for Health, Università degli Studi di Milano.

Compliance with ethical standards
Conflict of interest No conflicts of interest, financial or otherwise, are declared by the author(s).
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons .org/licen ses/by/4.0/.