BioMed Research International
Volume 2019 (2019), Article ID 4250940, 12 pages
https://doi.org/10.1155/2019/4250940

Research Article

Histomorphometric Analysis of Callus Formation Stimulated by Axial Dynamisation in a Standardised Ovine Osteotomy Model

K. M. Reich,1,2 S. Tangl,1,2 P. Heimel,2,3 S. Lettner,1,2 A. Ignatius,⁴ L. E. Claes,⁴ J. Pfeil,⁵ A. Janousek,⁶ and H. Redl2,3

¹Karl Donath Laboratory for Hard Tissue and Biomaterial Research, Department of Oral Surgery, University Clinic of Dentistry, Medical University of Vienna, Vienna, Austria
²Austrian Cluster for Tissue Regeneration, Vienna, Austria
³Ludwig Boltzmann Institute for Experimental and Clinical Traumatology, AUVA Research Centre, Vienna, Austria
⁴Institute of Orthopedic Research and Biomechanics, Center of Musculoskeletal Research Ulm, University Hospital Ulm, Ulm, Baden-Württemberg, Germany
⁵Department of Orthopaedics and Orthopaedic Surgery, St. Josefs Hospital, Wiesbaden, Germany
⁶Lorenz Boehler Trauma Center, Vienna, Austria

Correspondence should be addressed to S. Tangl; stefan.tangl@meduniwien.ac.at

Received 10 September 2018; Revised 7 December 2018; Accepted 16 January 2019; Published 12 February 2019

Academic Editor: Mirella Falconi

Copyright © 2019 K. M. Reich et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The cyclic axial dynamisation of a stabilised fracture is intended to promote callus formation and bone healing. Most studies focused on biomechanical properties or the quantity of new bone formation. Far less is known about the quality of newly formed callus tissues, such as tissue distribution and arrangement within the callus. The aim of this current study was to investigate the effect of cyclic, axial dynamisation on the quantity and quality of callus in an established delayed fracture healing model. In 41 sheep transverse osteotomies with a gap size of 3 mm were stabilised with a unilateral external fixator. In 32 of these, fracture ends were axially stimulated with displacement amplitudes of 0.8 mm, 0.4 mm, 0.2 mm, or 0.0 mm, respectively, for six weeks. In the remaining 9 sheep of the control group, an additional external fixator was mounted to achieve almost total rigidity. Animal material originating from a past animal experiment was reanalysed in this study. Histological thin-ground sections were histomorphometrically analysed regarding the histological structure and composition of the defect region. A slight tendency towards an increase in size of total callus area, area of new bone (nB.Ar), and cartilage (Cg.Ar) was detected with increasing displacement amplitudes compared to the control group. At the anterior callus side nB.Ar and Cg.Ar were significantly larger than at the posterior side in all groups independent of treatment. Regarding the quality of callus, areas of very compact bone were predominant in the treatment groups whereas in the control group a slight shift to more porous bone was observed. No difference of callus compactness was observed between the anterior and the posterior side. The established method to assess the local compactness of callus areas is a useful tool to quantitatively determine the spatial distribution of new bone tissue within the callus. The application of this method in combination with biomechanical testing might reveal interesting relations between tissue distribution and bone strength that, with traditional histomorphometry, cannot be identified.

1. Introduction

In clinics, most biomechanical and patient-related factors effecting fracture healing are unmodifiable such as the type of (comminuted) fracture, severity of soft tissue injury, blood supply, and also patients’ age and pre-existing (inflammatory) diseases [1]. What can be adapted is the rigidity of fracture stabilisation. The use of external fixators for fracture stabilisation allows for the individual adjustment of the mechanical environment which is mainly composed of fracture end stability and mechanical stress particularly applied through weight bearing during the fracture healing process [2]. Rigid devices reduce interfragmentary movement (IFM) and thus mechanical stress on fracture ends, whereas flexible devices enable limited axial loading during weight bearing. With the application of an additional stimulation module the external fixator can stimulate the fracture ends in a controlled manner [3].

The principal idea to mechanically stimulate fractured bone in order to enhance healing is based on Wolff’s decisive conclusion on the Mechanotransduction of bone [4]. Due to the ability of bone to sense differences in mechanical strain and to translate this information into signals to stimulate bone remodeling, bone is capable to adapt its mechanical properties according to its present needs. Also, in fracture healing, the biomechanical environment is critical [5, 6]. The external mechanical stimulation of a stabilised fracture site is intended to promote bone healing by simulating loads that naturally occur in healthy bone underweight bearing, for instance. The force, frequency, timing, and duration of externally applied mechanical stimulation via so-termed micromovements can be varied.

In search of optimized stimulation parameters diverse delayed fracture healing models were introduced differing in the choice of animal model, bone region to be studied, osteotomy gap size, type of fixator, magnitude of applied force, timing, frequency, and direction of interfragmentary movement (compression versus distraction)[7, 8]. This variety of studied parameters and the partly diverging results make conclusions on the effect of external mechanical stimulation difficult [7, 9, 10].

It is generally accepted that moderate mechanical stimulation via such applied cyclic axial micromovements is potent to increase callus formation [11–15]. Unlike cyclic loading, static loading seemed to have no effect on bone healing [16]. Shear forces do not impede but seem to prolong the healing process [17]. Large gap size and large interfragmentary strains [9] especially when applied at the wrong time point during the healing process are detrimental to fracture healing [12, 18]. The early phase of fracture healing was shown to be particularly sensitive to mechanical stimulation [12, 19]: some found prolonged healing and an extended chondral phase when the external fixation allowed early IFM [12, 20]. Recently, the focus was put on the modulation of stiffness of fixators during the course of healing [21, 22]. Claes et al. found that primary stability in terms of rigid fixation for at least 3 weeks and following late dynamisation result in superior healing in rats compared to rats with constant flexible or constant rigid fixation [23]. In contrast, the concept of “inverse” or “reverse dynamisation” recommends early flexible fixation during the early healing stages and rigid fixation during the later phases of healing [19, 21, 22, 24]. In these studies, priority was mainly given to clinically relevant parameters such as biomechanical properties (bending stiffness and torsional stiffness) or to the quantity of callus formation in terms of area or volume of newly formed bone [14, 25]. Although the quantity of the callus alone is not a useful predictor for mechanical stability of bone healing [9, 23, 26], far less is known about its quality, that is, the “density” or “compactness” of the newly formed bone tissue and the distribution of tissues within the callus [27, 28].

The aim of this present study was to broaden the current knowledge of the quality and quantity of callus formation by analyzing a well-established delayed fracture healing model in sheep [25, 29, 30] with detailed semiautomatic histomorphometry. Sheep tibiae that have already been employed for two publications [29, 30] were reevaluated for this purpose. Emphasis was put on the arrangement of tissues, that is, the composition and alignment of newly formed bone and cartilage within the callus. Therefore, tibial osteotomies with a gap size of 3 mm were stabilised with unilateral external fixators and stimulated with cyclic axial compression for six weeks. Since best stimulating effects were reported for strains in the range of 0.2 mm to 1 mm [31] we focused in our study on four groups with displacement amplitudes of 0.0 mm, 0.2 mm, 0.4 mm, and 0.8 mm and compared them with a control group equipped with a more rigid two-plane osteotomy fixation device.

We hypothesised that mechanical stimulation with different displacement amplitudes applied daily over a period of 6 weeks has a stimulatory effect not only on the size of callus formation but also on callus distribution and compactness. Based on several studies describing a side-dependent difference in the amount of medial and lateral bone formation within the callus [20, 27, 32], we examined if such regional differences also exist between the anterior and posterior side of the callus. Since mechanical properties of the callus are dependent not only on the quantity but also on the quality of newly formed bone [9, 26, 28, 33], including its porosity, the distribution of areas of low, medium, and high bone density within the callus was analysed.

2. Material and Methods

This study reuses sheep tibiae that have already been employed for two publications focusing on the biomechanical and radiological examination of osteogenesis in the cortical bone [29, 30]. Here we addressed new aspects of the same animal model that have not been studied so far and raised novel scientific questions concerning the histological structure and composition of the defect region. While the first two publications concentrated on the question if the method is clinically feasible or not [29, 30], this current study aimed to describe the healing processes at the histological level. New and additional techniques were used for the acquisition of the present data.

Forty-one skeletally mature female merino sheep (average age: 24 months, weighting 50-60 kg) underwent a transverse osteotomy of the left tibia. An osteotomy gap size of 3 mm was produced to simulate a delayed fracture healing model [25]. In all animals osteotomy gaps were laterally stabilised with a custom-made unilateral external fixator comprising a telescoping shaft for controlled axial interfragmentary movement when being unlocked. Sheep were randomly divided into 5 groups: in one group no mechanical stimulation was applied (0.0 mm group). In 3 groups the healing zone of the osteotomy fragments was axially stimulated by means of an attached stimulation microprocessor for 20 min/day with a frequency of 1 Hz and displacement amplitudes of 0.2 mm, 0.4 mm, and 0.8 mm, respectively. Animals of the control group received a unilateral external fixator combined with an additional monotube external fixator (anterolaterally) to increase interfragmentary rigidity and thus disable movement. Table 1 summarises the treatment groups (including the number of evaluated and excluded animals) and the respective stimulation regimens via the motor-driven external fixator.

Table 1: Treatment groups and stimulation regimen.

Group
Evaluated Animals (n)
Excluded Animals (n)
Fixation Device
Mechanical Stimulation

0.0 mm group
6
0
unilateral external fixator
no stimulation

0.2 mm group
6
1
unilateral external fixator
1200 cycles/day with 1 Hz, amplitude 0.2 mm

0.4 mm group
7
2
unilateral external fixator
1200 cycles/day with 1 Hz, amplitude 0.4 mm

0.8 mm group
8
2
unilateral external fixator
1200 cycles/day with 1 Hz, amplitude 0.8 mm

control group
4
5
unilateral external fixator AND unilateral monotube external fixator
no stimulation

Total
31
10

The study protocol was approved by the National Commission for Animal Experiments in South Africa (No. 0947, Biocon Research Institute, Pretoria, South Africa). All animal experiments were performed in Pretoria (South Africa) and were in accordance with the institutional guidelines for care and use of animals. Subsequent histological processing was performed in the Karl Donath Laboratory for Hard Tissue and Biomaterial Research, School of Dentistry, Medical University of Vienna (Austria).

2.1. Surgery

Detailed surgical protocols, characterisation of the construction, and working mechanisms of the external fixators have been described in greater detail by Wolf et al. [29] and Krischak et al. [30]. In short, a custom-made unilateral external fixator (F. Haas, Künzelsau, Germany) was attached to the exposed lateral tibial diaphysis of all 41 sheep using four 5 mm half-screws (Synthes, Umkrich, Germany). The periosteum was resected 2 mm proximal and distal of the osteotomy line to avoid irritation during sawing. Standardised transverse osteotomy was performed under general anaesthesia using an oscillatory saw and a saw guide, being constantly irrigated with saline solution. Bone fragments were distracted to a gap size of 3 mm and stabilised with the unilateral external fixator. Single knot technique was used to close the superficial layer and skin.

In 9 of the 41 sheep an additional unilateral monotube external fixator (Synthes, Umkrich, Germany) was attached to the anterolateral side of the same tibia with additional four 5 mm half-screws to achieve increased rigidity. These sheep being equipped with a more rigid two-plane fixation configuration are referred to as “control group” throughout the text.

During surgery a single dose of ampicillin (1 g) was given for antibiotic prophylaxis. Pin care was performed daily in all animals. Sheep were allowed full weight bearing immediately after surgery.

2.2. Stimulation via External Cyclic Compression

In the 32 sheep being exclusively equipped with the custom-made unilateral external fixator, a custom stimulation module was mounted on the telescoping rod of the fixator which was electromechanically driven and controlled by a microprocessor. Animals were randomly divided into 4 treatment groups with different displacement amplitudes of cyclic compression (0.0 mm, 0.2 mm, 0.4 mm, and 0.8 mm). The 0.0 mm group received no mechanical stimulation. The other groups, that is, 0.2 mm group, 0.4 mm group, and 0.8mm group, were stimulated 20 minutes per day at a frequency of 1 Hz. Stimulation was started 12 days post-op. To improve legibility, the 0.0 mm, 0.2 mm, 0.4 mm, and 0.8 mm groups are referred to as “treatment groups” throughout the text.

The rest of the day the telescoping mechanism of the external fixator was locked with screws in order to impede uncontrolled axial interfragmentary movement. Mechanical properties of the external fixators are listed in Table 2. After a healing period of 6 weeks animals were killed with an overdose of pentobarbitone. However, due to deep pin infections and fractures through the screw holes, ten sheep were not able to achieve full weight bearing and thus had to be excluded from the study.

Table 2: Mechanical properties of external fixator configurations [30].

One-plane fixation  
(Stimulation groups)
Two-plane fixation  
(Control group)

unilateral external fixator
unilateral external fixator AND unilateral monotube ext. fixator

Axial stiffness
183 N/mm
388 N/mm

Torsional stiffness
2.5 Nm/deg
2.5 Nm/deg

Bending stiffness in mediolateral direction
2Nm/deg  
(167 N/mm )
3 Nm/deg  
(250 N/mm)

Bending stiffness in anterior-posterior direction
31 N/mm
170 N/mm

Estimated value. Comparative data for bending stiffness in mediolateral direction given in N/mm is missing in Krischak et al. (2002)

2.3. Histology

After sacrificing the animals the external fixators were carefully removed from the dissected tibia. For the preparation of longitudinal undecalcified thin-ground sections bone blocks of the osteotomised tibiae were removed, fixed in neutral-buffered formalin, and embedded in methyl methacrylate (PMMA, Heraeus Kulzer GmbH, Wehrheim, Germany). Longitudinal slices were cut with a band saw, mounted on plastic slides, and ground to a thickness of approximately 100 μm. The polished thin ground sections were finally stained with Paragon.

2.4. Histomorphometry

Specimens were photographed with a digital camera (DXM1200; Nikon Cooperation, Tokyo, Japan) mounted on a Microphot FXA microscope (Nikon, Cooperation) at a resolution of 229 pixel/mm. Single images were assembled to create large overview pictures (Lucia G 4.71, Laboratory Imaging Ltd., Praha, Czech Republic) which were further processed in Adobe Photoshop 7.0® (Adobe, San Jose, CA, USA).

The region of interest was defined as the area of the callus 1.5 cm proximal and 1.5 cm distal from the centre of the osteotomy gap, comprising all tissues within the newly formed callus area (new bone, cartilage, and soft/fibrous tissues) but excluding the marrow space surrounding the callus as well as the original proximal and distal autochthonous cortical bone. The lateral margins of the ROI were defined by the maximal extension of the callus.

Based on an algorithm established with Definiens Image Analysis Software® (Definiens AG, Munich, Germany) newly formed bone, cartilage, soft/fibrous tissue, and background were automatically segmented and classified within the region of interest. Falsely classified areas were corrected manually under microscopic control.

Histomorphometric evaluation focused on four aspects:

(1)Total callus area (Cl.Ar in mm²) comprising newly formed bone, cartilage, and fibrous tissue encapsulated within the callus present within the given region of interest.(2)Total area of newly formed bone tissue (nB.Ar in mm²) and the total area of cartilage (Cg.Ar in mm²) were measured within the total callus area. Besides, new bone volume fraction (nBV/TV in %) and cartilage volume fraction (CgV/TV in %) were calculated within the total callus area, that is, tissue volume (TV).(3)The region of interest was longitudinally divided into the anterior and the posterior callus area. Newly formed bone area and cartilage area were measured within the anterior (An.) and the posterior (Pt.) callus area, respectively, [An.nB.Ar, Pt.nB.Ar, An.Cg.Ar, Pt.Cg.Ar (in mm²)]. Moreover, new bone volume fraction (nBV/TV in %) and cartilage volume fraction (CgV/TV in %) were calculated within these regions.(4)Since the mechanical properties of the callus are dependent not only on the quantity but also on the quality of newly formed bone including its porosity [9, 26], the distribution of new bone within the callus was measured. Therefore, three classes of new Bone Density (nB.Dn) were defined by sorting the bone into three brackets based upon their local “bone area density” (not to be confused with bone mineral density as measured by densitometry). Areas where less than 33% were occupied with new bone were defined as regions with low bone density, representing bone with high porosity. Areas with 33% to 66% new bone were defined as regions of medium bone density areas and areas exceeding the threshold of 66% were classified as high bone density regions, representing very compact bone. Morphometrically speaking, the local bone density of each pixel was measured as the bone area fraction within a distance of 0.8 mm around the pixel. Size and distribution of low, medium, and high nB.Dn areas were analysed (Figures 2 and 3).

2.5. Statistics

Statistical calculations were conducted using R version 3.4.0 [34]. Mixed models [35] for dependent variables (nB.Ar and Cg.Ar) including a random intercept for animal ID, fixed effects for treatment and location (anterior and posterior), and an interaction term were fitted. For percentages of density areas, an additional level of random effects was included to account for all densities simultaneously. Kenward-Roger approximation for F-tests [36] was used to test for the hypotheses of no effect of several fixed effects.

3. Results

In total five sheep of the treatment groups and five sheep of the control group had to be excluded from the study as they did not achieve full weight bearing due to deep pin infections and/or fractures through the screw holes. Consequently, callus formation was analysed in six sheep of each, the 0.0 mm group and the 0.2 mm group, seven sheep of the 0.4 mm group, eight sheep of the 0.8 mm group, and four sheep of the control group. Because of these unexpected dropouts statistical power of our analysis was very low. Detailed histomorphometric results for the remaining sheep are presented in box plots (Figure 1 and Tables 3 and 4).

Table 3: nB.Ar and Cg.Ar were assessed separately for the anterior and the posterior callus side. Anterior new bone formation was up to twice as high as posterior. Also cartilage area was significantly higher at the anterior callus side. Data are presented as mean values (mean), standard deviations (SD), median, minimum (Min), and maximum (Max) in mm.

Location
Treatment Group
Mean
SD
Median
Min
Max

nB.Ar
anterior
0.0 mm
88.1
21.2
92.6
47.0
107.8

0.2 mm
70.9
32.3
74.4
18.6
109.8

0.4 mm
111.7
38.4
99.5
65.4
179.3

0.8 mm
103.8
36.4
111.9
42.3
158.5

Control
70.7
44.0
65.6
27.0
124.7

posterior
0.0 mm
48.8
27.3
44.4
22.3
98.1

0.2 mm
63.1
25.4
56.6
34.9
102.1

0.4 mm
58.8
28.1
56.4
17.4
99.1

0.8 mm
56.8
25.5
54.2
27.1
96.0

Control
34.4
9.1
35.4
23.0
43.6

Total
0.0 mm
136.9
16.0
146.4
116.0
147.9

0.2 mm
134.0
20.4
131.3
109.2
156.4

0.4 mm
170.5
53.1
164.5
106.1
260.8

0.8 mm
160.6
55.8
158.2
85.5
230.7

Control
105.1
43.6
98.6
67.1
156.2

Cg.Ar
anterior
0.0 mm
5.2
4.4
4.3
1.0
10.4

0.2 mm
5.9
4.6
7.1
0.0
11.0

0.4 mm
8.7
9.3
7.6
0.3
25.1

0.8 mm
5.6
3.0
5.4
1.9
11.7

Control
0.3
0.3
0.2
0.0
0.7

posterior
0.0 mm
1.3
2.4
0.2
0.0
6.0

0.2 mm
3.5
1.9
3.9
1.0
5.6

0.4 mm
2.4
3.3
0.4
0.0
8.5

0.8 mm
2.4
2.1
1.8
0.4
6.6

Control
0.0
0.0
0.0
0.0
0.0

Total
0.0 mm
6.6
4.5
7.4
1.0
11.9

0.2 mm
9.4
5.7
10.4
2.0
16.7

0.4 mm
11.1
10.1
9.3
0.3
30.2

0.8 mm
7.9
4.4
8.1
2.4
15.2

Control
0.3
0.3
0.2
0.0
0.7

Table 4: nBV/TV and CgV/TV were assessed separately for the anterior and the posterior callus side. Although the absolute area of callus tissues is significantly different between the anterior and the posterior side, the relative amount of new bone within the anterior and the posterior callus area (nBV/TV) is very similar. Data are presented as mean values (mean), standard deviations (SD), median, minimum (Min), and maximum (Max) in %.

Location
Treatment Group
Mean
SD
Median
Min
Max

nB.Ar
anterior
0.0 mm
67.6
5.0
67.8
59.4
74.6

0.2 mm
64.0
11.4
62.9
50.0
82.7

0.4 mm
66.4
9.9
67.1
50.8
78.8

0.8 mm
60.2
11.2
61.1
37.3
74.9

Control
59.7
2.7
59.1
57.6
63.3

posterior
0.0 mm
64.4
7.1
63.8
53.2
72.3

0.2 mm
62.5
7.4
62.3
54.0
75.1

0.4 mm
65.8
8.8
65.5
57.5
77.8

0.8 mm
61.9
12.9
64.9
32.4
75.3

Control
55.0
8.4
57.3
42.9
62.5

Total
0.0 mm
66.4
4.2
65.1
62.7
74.0

0.2 mm
61.6
7.1
62.0
51.2
72.6

0.4 mm
66.0
8.9
63.3
52.8
78.2

0.8 mm
60.6
12.0
62.3
34.6
75.0

Control
58.6
3.6
59.4
53.9
61.7

Cg.Ar
anterior
0.0 mm
3.7
2.6
3.6
0.8
6.5

0.2 mm
3.9
2.7
4.8
0.0
6.3

0.4 mm
4.3
4.3
3.3
0.3
10.9

0.8 mm
3.4
1.8
3.2
1.3
6.3

Control
0.2
0.3
0.2
0.0
0.5

posterior
0.0 mm
1.1
1.6
0.3
0.0
3.8

0.2 mm
3.7
2.3
3.1
1.1
6.9

0.4 mm
2.0
2.3
1.0
0.0
5.2

0.8 mm
2.8
2.6
2.2
0.4
7.8

Control
0.0
0.0
0.0
0.0
0.0

Total
0.0 mm
3.0
1.9
3.3
0.6
5.3

0.2 mm
4.0
2.1
4.6
1.2
6.5

0.4 mm
3.8
3.2
2.9
0.2
9.2

0.8 mm
3.2
1.9
3.0
1.0
5.5

Control
0.2
0.2
0.1
0.0
0.3

(a)

(b)

(c)

Figure 1: (a) Total Callus Area (Cl.Ar). Total callus area (comprising newly formed bone, cartilage, and fibrous tissue within the ROI) increased with the increase of displacement amplitudes in the treatment groups. Largest Cl.Ar was found in the 0.4 mm group and the 0.8 mm group, smallest callus area in the control group. (b) New Bone Area (nB.Ar) within the total callus. Formation of new bone within the total callus area was lowest in the control group. Highest nB.Ar was found in the 0.4 mm group and the 0.8 mm group. (c) Cartilage Area (Cg.Ar) within the total callus. Presence of cartilage in the total callus area was very low in all groups. Variation within the groups was very high since cartilage was not present in all specimens. In the treatment groups Cg.Ar was slightly increased compared to the control group.

Figure 2: Size of Different Bone Density Areas (nB.Dn). Mineralised callus areas were classified as low nB.Dn (areas less than 33% new bone), medium nB.Dn (areas with 33 to 66% new bone), and high nB.Dn (areas exceeding 66% new bone). In this box plot diagram, the average size of low, medium, and high bone density areas is depicted for the 0.0 mm, 0.2 mm, 0.4 mm, and 0.8 mm group as well as for the control group.

Figure 3: Areas of Different Bone Densities (nB.Dn) within the Callus. Histological images were semiautomatically segmented and areas of different bone densities were colour encoded. Mineralised callus areas with less than 33% new bone were defined as regions with low nB.Dn (red), areas with 33 to 66% new bone as regions of medium nB.Dn (turquoise), and areas exceeding 66% were classified as high nB.Dn (light blue). For better visualisation, histological images are superimposed by the respective false colour images. The distribution of areas with different bone densities over the total callus area was very similar between the groups and showed a consistent inherent pattern. Besides, treatment groups (a-d) show a very similar proportional distribution: a relatively small area of low nB.Dn (5-10%), approximately 30% medium nB.Dn, and 50-65% high nB.Dn. In the control group (e), a shift to less dense bone was observed: the upwards trend observed in the treatment groups levelled off at the medium density level in the control group (approximately 45% medium nB.Dn). Areas with high nB.Dn were less common (approximately 42%) as compared to the treatment groups. (a) 0.0 mm group, (b) 0.2 mm group, (c) 0.4 mm group, (d) 0.8 mm group, and (e) control group. Anterior side is on the left and posterior side is on the right.

3.1. Total Callus Area (Cl.Ar)

Total callus area increased with higher displacement amplitudes (Figure 1(a)). The largest Cl.Ar was found in the 0.8 mm group (268.5 mm², SD 89.8) and 0.4 mm group (262.6 mm², SD 82.9), the smallest callus area in the control group (182.0 mm², SD 83.9). The difference between the groups was statistically not significant (p= 0.342).

3.2. nB.Ar in the Total Callus Area

Formation of new bone within the total callus area was lowest in the control group showing a nB.Ar of 105.1 mm² (SD 43.6). In the treatment groups nB.Ar was slightly increased in the 0.0 mm group (136.9 mm², SD 16.0) and the 0.2 mm group (134.0 mm², SD 20.4) when compared with the control group. Highest new bone formation was found in the 0.4 mm group and the 0.8 mm group with nB.Ar of 170.5 mm² (SD 53.1) and 160.6 mm² (SD 55.8), respectively. Data are presented in Figure 1(b). The difference between the groups failed to reach statistical significance (p=0.198).

3.3. Cg.Ar in the Total Callus Area

Presence of cartilage in the total callus area was very low in all groups (Figure 1(c)). Lowest Cg.Ar was measured in the control group (0.3 mm², SD 0.3). In the treatment groups Cg.Ar was slightly higher ranging from 6.6 mm² (SD 4.5) in the 0.0 mm group to 11.1 mm² (SD 10.1) in the 0.4 mm group. Variation within the groups was very high since cartilage was not present in all specimens. The difference between the groups was statistically not significant (p=0.206).

3.4. nB.Ar in the Anterior and Posterior Callus Area

When dividing the callus axially in an anterior and a posterior region of interest a statistically significant difference in nB.Ar between the regions was detected (p=0.001). Detailed results are presented in Table 3.

At the anterior side of the callus new bone formation was significantly higher compared to the posterior side in all, the control and the treatment groups. Again there was a slight tendency towards an increase of nB.Ar in the groups with higher stimulation amplitudes in both regions.

The lowest values for An.nB.Ar and Pt.nB.Ar were found in the control group (An.nB.Ar 70.7 mm², SD 44.0 and Pt.nB.Ar 34.4 mm², SD 9.1). Highest An.nB.Ar was measured in the 0.8 mm group (103.8 mm², SD 36.4) and the 0.4 mm group (111.7 mm², SD 38.4). In contrast, Pt.nB.Ar was highest in the 0.2 mm group (63.1 mm², SD 25.4). The interaction between region and treatment group was not statistically significant (p=0.403).

3.5. Cg.Ar in the Anterior and Posterior Callus Area

Cg.Ar was significantly higher in the anterior callus side compared to the posterior callus side (p=0.029). Anterior Cg.Ar was found to be smallest in the control group with 0.3 mm² (SD 0.3). Highest An.Cg.Ar was found in the 0.4 mm group (8.7 mm², SD 9.3). Posteriorly, no cartilage was present in any animal of the control group. In the 0.2 mm group, Pt.Cg.Ar was highest with 3.5 mm² (SD1.9). In the other groups Pt.Cg.Ar was very close ranging from 1.3 mm² (SD 2.4) in the 0.0 mm group to 2.4 mm² (SD 3.3) and 2.4 mm² (SD 2.1) in the 0.4 mm and 0.8 mm group, respectively. See Table 3. The interaction between region and treatment group was not statistically significant (p=0.536).

3.6. nBV/TV and CgV/TV

Although the absolute size of the callus area was significantly different between the anterior and the posterior side, the volume fraction of newly formed bone within these given regions was very similar (p=0.421). Detailed histomorphometric results for total nBV/TV and total CgV/TV as well as for the anterior and posterior volume fractions are presented in Table 4.

The parameter nBV/TV, however, only reflects the overall density of the newly formed bone within the callus area. In order to gather information about the actual distribution of the new bone within the callus area, the new parameter new Bone Density (nB.Dn) was introduced as follows.

3.7. Size of Different Bone Density Areas (nB.Dn)

Differences in the size of different bone density areas are presented in Figure 2. Their distribution within the ROI is illustrated in Figures 3(a)–3(e).

Areas of low bone density (low nB.Dn) were defined as callus areas of high porosity, that is, callus areas comprising less than 33% newly formed bone. Low bone density areas were very rare in all groups. In the treatment groups low nB.Dn areas ranged from 10.0 mm² (SD 5.5) in the 0.0 mm group to 29.2 mm² (SD 32.1) in the 0.8 mm group. In the control group, low bone density areas had a size of 27.4 mm² (SD 23.1).

Areas of medium bone density (medium nB.Dn), that is, callus areas comprising 33-66% new bone, were more frequent than low bone density areas in all groups. In the 0.0 mm group, the 0.2 mm group and the 0.4 mm group medium bone density areas were very similar with 61.4 mm² (SD 24.9), 65.1 mm² (SD 19.0), and 71.4 mm² (SD 45.1), respectively. Largest medium nB.Dn areas were found in the 0.8 mm group (100.5 mm², SD 62.9). Medium nB.Dn was the commonest density in the control group averaging 81.3 mm² (SD 37.8).

Areas of high bone density (high nB.Dn), that is, very compact callus areas comprising more than 66% newly formed bone, were predominant in all treatment groups. In the 0.0 mm group, the 0.2 mm group, and the 0.8 mm group high nB.Dn areas were very similar, averaging 129.2 mm² (SD 13.6), 122.1 mm² (SD 27.3), and 130.8 mm² (SD 60.4), respectively. Largest high bone density areas were measured in the 0.4 mm group (164.3 mm², SD 59.0). In the control group high nB.Dn areas were smallest (73.1 mm², SD 27.6).

When looking at the ratio of low, medium, and high density areas within the callus, treatment groups showed a very similar proportional distribution: a relatively small area of low nB.Dn (5-10%), approximately 30% of medium nB.Dn, and 50-65% of high nB.Dn. However, in the control group a shift to less dense bone was observed: the upward tendency observed in the treatment groups levelled off at the medium density level in the control group (approximately 45% are medium nB.Dn). Areas with highly dense bone were less common (approximately 42%) as compared to the treatment groups.

Due to the small sample size the difference failed to reach statistical significance between the groups (p=0.556). No difference between the anterior and the posterior side of the callus concerning the proportion of low, medium, and high bone density areas was detected (p=0.371).

3.8. Regional Distribution of Different Bone Density Areas

The distribution of areas with different bone density over the total callus area was very similar between the groups and showed a consistent inherent pattern (Figures 3(a)–3(e)).

Areas of highly dense bone accounted for the majority of the total callus area and were predominantly located adjacent to the periosteal and endosteal surfaces of the tibial diaphysis and in the periphery of the periosteal callus. This finding was regardless of whether the gap was already bridged periosteally, endosteally, or intercortically. In specimens with an osseous bridged gap, highly dense bone was also found in the intercortical site and the endosteal compartment adjacent to the former gap.

Areas of medium dense bone most frequently occurred in specimens with a less developed osseous bridge in the intercortical region, that is, between the proximal and distal cortices. Independent of the grade of intercortical bridging, medium dense bone was found to be the predominant density bridging the anterior and the posterior endosteal surface within the marrow space. In addition, medium bone density areas were present at the outer margin of the periosteal callus at the level of osteotomy gap, and rarely between the cortical fragments. If cartilage was present (mainly in the periosteal gap), low and medium dense bone areas were surrounding cartilage remnants. Bone was observed to be denser with increasing distance to cartilage.

In all groups some individuals showed mid-sized to large areas of fibrous tissues within the callus.

4. Discussion

In this present study, a well-established delayed fracture healing model in the sheep [25] was used to analyse the effect of cyclic axial stimulation on the size of callus tissues and the spatial distribution of newly formed bone within the callus in terms of three density classes. A slight tendency towards an increase in size of total Cl.Ar, nB.Ar, and Cg.Ar was detected with rising displacement amplitudes. In addition, a side-specific difference of nB.Ar and Cg.Ar was identified between the anterior and the posterior callus side in all groups independent of treatment. Regarding the quality of callus, areas of very compact bone were predominant in the treatment groups whereas in the control group a slight shift to more porous bone was observed. Statistical power was very low due to the high complication rate in the control group.

4.1. Quantity of Callus Tissues

Our findings obtained by detailed semiautomatic histomorphometry are in line with other studies performed in sheep reporting increased new bone area in dynamically stabilised osteotomy gaps compared to more rigidly-stabilised gaps measured by various methods [9, 11, 15, 31]. A recent μCT study of Tufekci et al. [19] demonstrated that 1 mm compressive IFMs, that were applied within the first 3 weeks, result in significantly increased bone volume compared to osteotomy gaps that were completely isolated from functional loading and IFM.

Considering the previous studies and the present results, displacement amplitudes in the range of 0.5-1 mm appear to have a potent stimulatory effect on bone healing in terms of increased callus formation. However, the timeframe of stimulation seems to be critical since bone is assumed to respond differently to mechanical strains depending on the stage of the bone healing process [12, 19, 23]. Recommendations concerning the optimal point in time for mechanical stimulation/dynamisation are conflicting [19, 21, 23, 24, 37]: Tufekci et al. demonstrated that the early stimulation within the first 3 weeks (starting 5 days after surgery) resulted in increased callus bone volume and torsional strength after 9 weeks in comparison with ongoing stimulation with decreasing displacement amplitudes. This so-termed “inverse stimulation” started with displacement amplitudes of 1 mm; stimulation was then continuously reduced in 0.25 mm increments until it was stopped after six weeks [19]. The authors demonstrated that, during the later stage of bone healing, mechanical stimulation is not needed. By contrast, it rather tends to compromise healing compared with early stimulation [19].

This is in clear contrast to Claes et al. who demonstrated superior healing with primary rigid (for at least 3 weeks) and subsequent flexible stabilisation of osteotomised femora of rats [23]. Compared to the constant flexible fixation group, callus area was decreased and bending rigidity was increased after 5 weeks. Late dynamisation was hypothesised to accelerate bone remodeling and consolidation of the callus [23]. In the current study, cyclic axial compressive dynamisation was also started comparatively late (12 days after surgery) and continued until 6 weeks postoperatively which is in between the dynamisation protocols of Tufekci et al. [19] and Claes et al. [23]. After 14 days, Epari et al. reported that periosteal woven bone formation is just about to start while soft callus composed of fibrous tissues is still present [20]. According to this, the present study reflects the effect of dynamisation on the beginning of hard callus formation.

The diverging results between these studies concerning the timing of dynamisation might be attributable to differences in the studied animal model, the gap size, the type of mechanical stimulation/dynamisation, the force and displacement amplitude of IFMs via loading or of external dynamisation, and the treatment duration, respectively. A standardised comparison of these dynamisation protocols is needed.

4.2. Quality of Callus Tissues

Since the quantity of the callus alone is not a useful predictor for mechanical stability of bone healing [9, 23, 26, 28], a second focus of this present study was set on a parameter of callus quality: the “density” or “compactness” of newly formed bone tissue or, in other words, the porosity of the calcified callus. In fact, mechanical axial dynamisation might have a slight effect on the density of newly formed callus: in the treatment groups, areas of high bone density were predominant whereas in the control group a slight shift to less dense bone was observed. The distribution of areas with different bone densities did not show an anterior- /posterior-specific pattern. Further studies with a larger study sample are needed to verify the observed tendency of this study. In addition, it has to be elucidated how the local distribution of areas with different compactness is correlated with bone strength and thus the resistance to refracture.

To our knowledge only a few studies were dealing with callus density. In these cases, density was defined as bone mineral density (BMD) of the callus or as bone area or bone volume per tissue area (BV/TV) which does not allow us to draw any conclusions about the distribution of newly formed bone within the callus and its actual porosity [20, 27, 33, 38, 39]. In this present study, we established a method to quantitatively describe the density or porosity of callus areas and the spatial distribution of new bone tissue within the callus. As there is no valid definition in the literature, we split densities into thirds of < 33% (low), 33-66% (medium), and > 66% bone tissue per callus area (high bone density). Despite the fact that the classification is not based on existing hard data, it provides interesting quantitative information about the distribution of bone tissue within the calcified callus that was qualitatively described in previous studies [20].

When interpreting the density or “compactness” of the callus, the stage of fracture healing and also the region (periosteal, endosteal, or intracortical) has to be considered. In this current study, moderate resorptive activity was predominantly present at the periosteal surface of the callus in all groups indicating that the ignition of external callus remodeling was rather at the same level in all groups. Accordingly, the lower bone density and also the lower nBV/TV in the control group do not seem to be the result of advanced resorption due to remodeling. From a theoretical point of view, they might rather be interpreted as a reduced need of bone mass in this area. In fact, the control group was equipped with two external fixators which provided higher stability. Most of the mechanical strains, which, under physiological conditions, are sustained by the bone itself, were now carried by the fixators. In consequence, in the control group, bone was subjected to reduced stresses when compared to the treatment groups. In other words, bone might have been understimulated [19] or stress shielded [40, 41]. Consequently, it might have simply been unnecessary for bone to invest in the formation of a larger and denser callus because stability was already provided. New bone might have only been formed to such an extent that stability is warranted in conjunction with the fixator. The mechanical phenomenon of “understimulation” or “stress shielding” of bone which occurs in areas where stability is already provided by an implant or fixator is mainly described for joint prosthesis in the hip or the knee [40, 42]. However, “stress shielding” was recently also described as being a potential risk of stiff osteotomy or fracture fixation allowing insufficient axial compression via external fixators [43–46]. Studies dealing with the effect of stress shielding mainly report bone resorption of bone areas adjacent to prosthesis or implants that are no longer subjected to mechanical strains due to changes in load transfer after implant integration [42]. Since the local strain distribution at an osteotomy site is assumed to provide mechanobiological signals that lead to a specific pattern of tissue distribution [47] it is conceivable that stress shielding also might have a direct effect on the bone formation process in osteotomy and fracture healing resulting in a reduced callus area and callus “density.”

In order to draw conclusions about the potential role of stress shielding in the context of external fixation of osteotomies, long-term studies are required analysing the course of bone formation from its beginning up to that point when an equilibrium between bone formation and bone resorption is achieved. The bone defect model established by Tufekci et al. (2018) which isolates an osteotomy gap from any kind of IFM might be a good tool to investigate the effect of “understimulation” on bone formation and callus density over the course of time.

4.3. Anterior/Posterior Difference in Callus Formation

The third focus of this present study was put on the regional distribution of newly formed tissues between the anterior and the posterior side of the callus. Several studies reported a side-dependent difference in callus formation between the medial and the lateral callus side when the osteotomies were stabilised with rigid or flexible fixators mounted on the medial aspect of the tibia [12, 20]. In this present study, we examined if differences exist also in the anterior-posterior plane when the fixator is mounted on the lateral side of the tibia. In fact, anterior nB.Ar was almost twice as high as posterior nB.Ar in all groups. The reasons for this obvious bias towards one side are not completely understood. When looking for potential explanations of this side-difference we came across the stiffness properties of the unilateral external fixator. In fact, 4-point bending stiffness of the unilateral fixator was significantly reduced in the anterior-posterior direction (31 N/mm) compared to the mediolateral direction (167 N/mm) in this present study as previously described by Wolf et al. [29]. This suggests the assumption that, in the treatment groups, interfragmentary movement in the anterior-posterior direction was possible during weight bearing at least to some extent. However, the same pattern of bone distribution was observed in the control group in which stiffness in the anterior-posterior direction was increased to 170 N/mm by the application of the second external fixator.

As reported by others, the orientation of an external unilateral fixator is potent to cause nonuniform interfragmentary movements and additional shear movements within the osteotomy gap [12, 15, 17, 48] which in turn might influence bone healing in one way. Beyond the mounting plane of the external fixators and the associated changes in shear forces, muscle strains, and gait patterns, also the greater soft tissue coverage at the anterior side and one-sided damage of blood vessels might contribute to differences in cell activation and thus in the temporal and spatial callus formation [12, 20, 49–51]. Further examinations are needed to determine if the observed anterior-posterior difference also exists in tibias being stabilised with medially mounted external fixators.

4.4. Limitations

Unfortunately, five sheep of the four treatment groups and five sheep of the control group had to be excluded from the study due to severe complications (Table 1). Pin track infections and fractures through the screw holes proximally and distally of the osteotomy gap occurred in more than 55% of the animals of the control group which led to a very small group size and weak power of our statistics.

A potential reason for the high drop-out rate in the control group might be the localisation and the increased number of screws that was necessary to fix the second external fixator [30]. The four additional screws undoubtedly resulted in an additional local tissue trauma and thus in an increased infection risk at the pin insertion sites. Also the anterolateral localisation of those screws might have been unfavourable due to its vicinity to muscle tissue. Weight bearing during walking leads to relative movements between the soft tissue and the screws, thus provoking muscle and soft tissue irritation, which is reported to be the most common reason for pin infections. [52] Similar complication rates were reported for studies in sheep and rabbits due to pin infections, pin loosening, or fractures [9, 53]. Also, in the clinical setting, pin site infections and fracture susceptibility of the screw holes are common complications [46, 54] but the better hygienic conditions facilitate prevention and healing in human patients.

Considering these complications and the high exclusion rate of animals it would have been very desirable to conduct additional animal experiments especially in order to increase the low number of control animals. As mentioned, the current study reevaluated sheep tibiae from an animal experiment performed about twenty years ago [29, 30]. We found that under these conditions it would not have been possible to perform a repetition of the surgery with the required degree of standardisation a posteriori. Therefore we decided not to supplement this current study with any new animal experimentation consciously accepting the low number of control animals.

Beyond the small group sizes, the main limitation of this study is that callus formation was analysed only after six weeks which makes it impossible to draw conclusions about potential differences between the groups in the early healing period. It remains open if the observed slight differences in callus densities are related to the timing of dynamisation. To analyse the progress of callus formation multiple time points should be considered.

5. Conclusion

In brief, axial dynamisation tended to slightly increase the volume of the newly formed callus. Besides, mechanical stimulation might have a slight effect on the density or compactness of the calcified callus. It could be demonstrated that the established method to assess the density or porosity of callus areas is a useful tool to quantitatively determine the spatial distribution of new bone tissue within the callus. The application of this method in combination with biomechanical testing might reveal interesting relations between tissue distribution of healing bone and bone strength that, with the traditional histomorphometric parameters BV/TV or B.Ar, cannot be identified.

In addition to the medial/lateral differences in new bone volume described in the literature, this study demonstrated that there exists also a difference between the anterior and posterior side. If this side difference also persists if the mounting plane of the external fixators is changed is issue of further analysis.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Disclosure

No external funding, apart from the support of the authors’ institution, was obtained for this study.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this article.

References

L. Claes, S. Recknagel, and A. Ignatius, “Fracture healing under healthy and inflammatory conditions,” Nature Reviews Rheumatology, vol. 8, no. 3, pp. 133–143, 2012.

S. M. Perren, “Physical and biological aspects of fracture healing with special reference to internal fixation,” Clinical Orthopaedics and Related Research, pp. 175–196, 1979.

A. E. Goodship, “Mechanical stimulus to bone,” Annals of the Rheumatic Diseases, vol. 51, no. 1, pp. 4–1, 1992.

J. Wolff, Das Gesetz der Transformation der Knochen, A. Hirschwald, Ed., Berlin, 1892.

E. Santolini, R. West, and P. V. Giannoudis, “Risk factors for long bone fracture non-union: a stratification approach based on the level of the existing scientific evidence,” Injury, vol. 46, Supplement 8, pp. S8–S19, 2015.

P. V. Giannoudis, T. A. Einhorn, and D. Marsh, “Fracture healing: the diamond concept,” Injury, vol. 38, no. 4, pp. S3–S6, 2007.

M. Jagodzinski and C. Krettek, “Effect of mechanical stability on fracture healing--an update,” Injury, vol. 38, supplement 1, pp. S3–S10, 2007.

J. C. Reichert, S. Saifzadeh, M. E. Wullschleger et al., “The challenge of establishing preclinical models for segmental bone defect research,” Biomaterials, vol. 30, no. 12, pp. 2149–2163, 2009.

L. Claes, P. Augat, G. Suger, and H.-J. Wilke, “Influence of size and stability of the osteotomy gap on the success of fracture healing,” Journal of Orthopaedic Research, vol. 15, no. 4, pp. 577–584, 1997.

J. Reifenrath, N. Angrisani, M. Lalk, and S. Besdo, “Replacement, refinement, and reduction: Necessity of standardization and computational models for long bone fracture repair in animals,” Journal of Biomedical Materials Research Part A, vol. 102, no. 8, pp. 2884–2900, 2014.

T. Yamaji, K. Ando, S. Wolf, P. Augat, and L. Claes, “The effect of micromovement on callus formation,” Journal of Orthopaedic Science, vol. 6, no. 6, pp. 571–575, 2001.

P. Klein, H. Schell, F. Streitparth et al., “The initial phase of fracture healing is specifically sensitive to mechanical conditions,” Journal of Orthopaedic Research, vol. 21, no. 4, pp. 662–669, 2003.

J. Kenwright and A. E. Goodship, “Controlled mechanical stimulation in the treatment of tibial fractures,” Clinical Orthopaedics and Related Research, no. 241, pp. 36–47, 1989.

A. E. Goodship, J. L. Cunningham, and J. Kenwright, “Strain rate and timing of stimulation in mechanical modulation of fracture healing,” Clinical Orthopaedics and Related Research, no. 355, pp. S105–S115, 1998.

L. Claes, H.-J. Wilke, P. Augat, S. Rübenacker, and K. Margevicius, “Effect of dynamization on gap healing of diaphyseal fractures under external fixation,” Clinical Biomechanics, vol. 10, no. 5, pp. 227–234, 1995.

L. E. Lanyon and C. T. Rubin, “Static vs dynamic loads as an influence on bone remodelling,” Journal of Biomechanics, vol. 17, no. 12, pp. 897–905, 1984.

H. Schell, D. R. Epari, J. P. Kassi, H. Bragulla, H. J. Bail, and G. N. Duda, “The course of bone healing is influenced by the initial shear fixation stability,” Journal of Orthopaedic Research, vol. 23, no. 5, pp. 1022–1028, 2005.

M. J. Gardner, M. C. H. van der Meulen, D. Demetrakopoulos, T. M. Wright, E. R. Myers, and M. P. Bostrom, “In vivo cyclic axial compression affects bone healing in the mouse tibia,” Journal of Orthopaedic Research, vol. 24, no. 8, pp. 1679–1686, 2006.

P. Tufekci, A. Tavakoli, C. Dlaska et al., “Early mechanical stimulation only permits timely bone healing in sheep,” Journal of Orthopaedic Research, vol. 36, no. 6, pp. 1790–1796, 2018.

D. R. Epari, H. Schell, H. J. Bail, and G. N. Duda, “Instability prolongs the chondral phase during bone healing in sheep,” Bone, vol. 38, no. 6, pp. 864–870, 2006.

V. Glatt, M. Miller, A. Ivkovic et al., “Improved healing of large segmental defects in the rat femur by reverse dynamization in the presence of bone morphogenetic protein-2,” The Journal of Bone & Joint Surgery—American Volume, vol. 94, no. 22, pp. 2063–2073, 2012.

D. R. Epari, T. Wehner, A. Ignatius, M. A. Schuetz, and L. E. Claes, “A case for optimising fracture healing through inverse dynamization,” Medical Hypotheses, vol. 81, no. 2, pp. 225–227, 2013.

L. Claes, R. Blakytny, J. Besse, C. Bausewein, A. Ignatius, and B. Willie, “Late dynamization by reduced fixation stiffness enhances fracture healing in a rat femoral osteotomy model,” Journal of Orthopaedic Trauma, vol. 25, no. 3, pp. 169–174, 2011.

N. Bartnikowski, L. E. Claes, L. Koval et al., “Modulation of fixation stiffness from flexible to stiff in a rat model of bone healing,” Acta Orthopaedica, vol. 88, no. 2, pp. 217–222, 2017.

A. E. Goodship and J. Kenwright, “The influence of induced micromovement upon the healing of experimental tibial fractures,” The Journal of Bone & Joint Surgery, vol. 67, no. 4, pp. 650–655, 1985.

P. Augat, J. Merk, H. K. Genant, and L. Claes, “Quantitative assessment of experimental fracture repair by peripheral computed tomography,” Calcified Tissue International, vol. 60, no. 2, pp. 194–199, 1997.

M. D. Markel, M. A. Wikenheiser, and E. Y. Chao, “A study of fracture callus material properties: Relationship to the torsional strength of bone,” Journal of Orthopaedic Research, vol. 8, no. 6, pp. 843–850, 1990.

M. Casanova, A. Schindeler, D. Little, R. Müller, and P. Schneider, “Quantitative phenotyping of bone fracture repair: a review,” BoneKEy Reports, vol. 3, 2014.

D.-I. S. Wolf, A. Janousek, J. Pfeil et al., “The effects of external mechanical stimulation on the healing of diaphyseal osteotomies fixed by flexible external fixation,” Clinical Biomechanics, vol. 13, no. 4-5, pp. 359–364, 1998.

G. D. Krischak, A. Janousek, S. Wolf, P. Augat, L. Kinzl, and L. E. Claes, “Effects of one-plane and two-plane external fixation on sheep osteotomy healing and complications,” Clinical Biomechanics, vol. 17, no. 6, pp. 470–476, 2002.

L. E. Claes, C. A. Heigele, C. Neidlinger-Wilke et al., “Effects of mechanical factors on the fracture healing process,” Clinical Orthopaedics and Related Research, no. 355, pp. S132–S147, 1998.

F. Witt, A. Petersen, R. Seidel, A. Vetter, R. Weinkamer, and G. N. Duda, “Combined in vivo/in silico study of mechanobiological mechanisms during endochondral ossification in bone healing,” Annals of Biomedical Engineering, vol. 39, no. 10, pp. 2531–2541, 2011.

F. C. den Boer, J. A. Bramer, P. Patka et al., “Quantification of fracture healing with three-dimensional computed tomography,” Archives of Orthopaedic and Trauma Surgery, vol. 117, no. 6-7, pp. 345–350, 1998.

R Core Team, “R: A Language and Environment for Statistical Computing., R Foundation for Statistical Comput- ing, Vienna, Austria,” 2017, http://www.R-project.org.

D. Bates, M. Mächler, B. M. Bolker, and S. C. Walker, “Fitting linear mixed-effects models using lme4,” Journal of Statistical Software, vol. 67, no. 1, pp. 1–48, 2015.

U. Halekoh and S. Højsgaard, “A kenward-roger Ap- proximation and parametric bootstrap methods for tests in linear mixed models – The R package pbkrtest,” Journal of Statistical Software, vol. 59, no. 9, pp. 1–30, 2014.

L. Claes, R. Blakytny, M. Gockelmann, M. Schoen, A. Ignatius, and B. Willie, “Early dynamization by reduced fixation stiffness does not improve fracture healing in a rat femoral osteotomy model,” Journal of Orthopaedic Research, vol. 27, no. 1, pp. 22–27, 2009.

E. F. Morgan, Z. D. Mason, K. B. Chien et al., “Micro-computed tomography assessment of fracture healing: relationships among callus structure, composition, and mechanical function,” Bone, vol. 44, no. 2, pp. 335–344, 2009.

B. Preininger, S. Checa, F. L. Molnar, P. Fratzl, G. N. Duda, and K. Raum, “Spatial-temporal mapping of bone structural and elastic properties in a sheep model following osteotomy,” Ultrasound in Medicine & Biology, vol. 37, no. 3, pp. 474–483, 2011.

R. Huiskes, H. Weinans, and B. Van Rietbergen, “The relationship between stress shielding and bone resorption around total hip stems and the effects of flexible materials,” Clinical Orthopaedics and Related Research, vol. 274, pp. 124–134, 1992.

D. M. O'Doherty, S. P. Butler, and A. E. Goodship, “Stress protection due to external fixation,” Journal of Biomechanics, vol. 28, no. 5, pp. 575–586, 1995.

D. R. Sumner and J. O. Galante, “Determinants of stress shielding: design versus materials versus interface,” Clinical Orthopaedics and Related Research, vol. 274, pp. 202–212, 1992.

P. Låtman, O. S. Nilsson, O. Brosjö, and L. Strömberg, “Stress shielding by rigid fixation studied in osteotomized rabbit tibiae,” Acta Orthopaedica, vol. 60, no. 6, pp. 718–722, 1989.

H. T. Aro, P. J. Kelly, D. G. Lewallen, and E. Y. Chao, “The effects of physiologic dynamic compression on bone healing under external fixation,” Clinical Orthopaedics and Related Research, vol. 256, pp. 260–273, 1990.

D. Paley, B. Fleming, M. Catagni, T. Kristiansen, and M. Pope, “Mechanical evaluation of external fixators used in limb lengthening,” Clinical Orthopaedics and Related Research, no. 250, pp. 50–57, 1990.

A. T. Fragomen and S. R. Rozbruch, “The mechanics of external fixation,” HSS Journal, vol. 3, no. 1, pp. 13–29, 2007.

P. Augat, U. Simon, A. Liedert, and L. Claes, “Mechanics and mechano-biology of fracture healing in normal and osteoporotic bone,” Osteoporosis International, vol. 16, supplement 2, pp. S36–S43, 2005.

G. N. Duda, K. Eckert-Hübner, R. Sokiranski, A. Kreutner, R. Miller, and L. Claes, “Analysis of inter-fragmentary movement as a function of musculoskeletal loading conditions in sheep,” Journal of Biomechanics, vol. 31, no. 3, pp. 201–210, 1997.

B. McKibbin, “The biology of fracture healing in long bones,” The Journal of Bone & Joint Surgery (British Volume), vol. 60, no. 2, pp. 150–162, 1978.

A. L. Wallace, E. R. Draper, R. K. Strachan, I. D. Mccarthy, and S. P. Hughes, “The vascular response to fracture micromovement,” Clinical Orthopaedics and Related Research, vol. 301, pp. 281–290, 1994.

F. Zhao, Z. Zhou, Y. Yan et al., “Effect of fixation on neovascularization during bone healing,” Medical Engineering & Physics, vol. 36, no. 11, pp. 1436–1442, 2014.

F. Behrens and W. Johnson, “Unilateral external fixation. Methods to increase and reduce frame stiffness,” Clinical Orthopaedics and Related Research, no. 241, pp. 48–56, 1989.

T. Matsushita and T. Kurokawa, “Comparison of cyclic compression, cyclic distraction and rigid fixation: Bone healing in rabbits,” Acta Orthopaedica, vol. 69, no. 1, pp. 95–98, 1998.

J.-J. Wu, H. S. Shyr, E. Y. S. Chao, and P. J. Kelly, “Comparison of osteotomy healing under external fixation devices with different stiffness characteristics,” The Journal of Bone & Joint Surgery, vol. 66, no. 8, pp. 1258–1264, 1984.