Whole-body vibration of mice induces articular cartilage degeneration with minimal changes in subchondral bone

Whole-body vibration of mice induces articular cartilage degeneration with minimal changes in subchondral bone

Accepted Manuscript Whole-body vibration of mice induces articular cartilage degeneration with minimal changes in subchondral bone Matthew R. McCann, ...

1MB Sizes 0 Downloads 112 Views

Accepted Manuscript Whole-body vibration of mice induces articular cartilage degeneration with minimal changes in subchondral bone Matthew R. McCann, Cynthia Yeung, Michael A. Pest, Anusha Ratneswaran, Steven I. Pollmann, David W. Holdsworth, Frank Beier, S. Jeffrey Dixon, Dr. Cheryle A. Séguin PII:

S1063-4584(16)30387-9

DOI:

10.1016/j.joca.2016.11.001

Reference:

YJOCA 3885

To appear in:

Osteoarthritis and Cartilage

Received Date: 2 March 2016 Revised Date:

29 July 2016

Accepted Date: 2 November 2016

Please cite this article as: McCann MR, Yeung C, Pest MA, Ratneswaran A, Pollmann SI, Holdsworth DW, Beier F, Dixon SJ, Séguin CA, Whole-body vibration of mice induces articular cartilage degeneration with minimal changes in subchondral bone, Osteoarthritis and Cartilage (2016), doi: 10.1016/j.joca.2016.11.001. This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.

ACCEPTED MANUSCRIPT

1 2

Title: Whole-body vibration of mice induces articular cartilage degeneration with minimal changes in subchondral bone

3 4 5 6 7

Running Title: WBV leads to cartilage degeneration

RI PT

a

SC

Department of Physiology and Pharmacology, Schulich School of Medicine & Dentistry, University of Western Ontario, London, Ontario, Canada. N6A 5C1 b Dentistry, Schulich School of Medicine & Dentistry, University of Western Ontario, London, Ontario, Canada. N6A 5C1 c Imaging Research Laboratories, Robarts Research Institute; Department of Medical Biophysics, and Department of Surgery, Schulich School of Medicine & Dentistry, University of Western Ontario, London, Ontario, Canada. N6A 5C1 d Bone and Joint Institute, University of Western Ontario, London, Ontario, Canada. N6A 5C1

M AN U

8 9 10 11 12 13 14 15 16

Authors: Matthew R. McCanna,d, Cynthia Yeunga, d, Michael A. Pesta, d, Anusha Ratneswarana, d, Steven I. Pollmannd, David W. Holdsworthc,d, Frank Beiera,d, S. Jeffrey Dixona,b,d, Cheryle A. Séguina,d

17

30 31 32

TE D

*Corresponding Author: Dr. Cheryle A. Séguin Department of Physiology and Pharmacology Schulich School of Medicine & Dentistry The University of Western Ontario London, Ontario, Canada, N6A 5C1 Email: [email protected]

EP

Keywords: Whole Body Vibration, Knee Joint, Cartilage Degeneration, Subchondral Bone

AC C

18 19 20 21 22 23 24 25 26 27 28 29

33 34 35

1

ACCEPTED MANUSCRIPT

ABSTRACT

37

Objective: Low-amplitude, high-frequency whole-body vibration (WBV) has been

38

adopted for the treatment of musculoskeletal diseases including osteoarthritis; however,

39

there is limited knowledge of the direct effects of vibration on joint tissues. Our recent

40

studies revealed striking damage to the knee joint following exposure of mice to WBV.

41

The current study examined the effects of WBV on specific compartments of the murine

42

tibiofemoral joint over 8 weeks, including microarchitecture of the tibia, to understand

43

the mechanisms associated with WBV-induced joint damage.

44

Design: Ten-week-old male CD-1 mice were exposed to WBV (45 Hz, 0.3 g peak

45

acceleration; 30 min/day, 5 days/week) for 4 weeks, 8 weeks, or 4 weeks WBV followed

46

by 4 weeks recovery. The knee joint was evaluated histologically for tissue damage.

47

Architecture of the subchondral bone plate, subchondral trabecular bone, primary and

48

secondary spongiosa of the tibia was assessed using micro-CT.

49

Results: Meniscal tears and focal articular cartilage damage were induced by WBV; the

50

extent of damage increased between 4 and 8-week exposures to WBV. WBV did not alter

51

the subchondral bone plate, or trabecular bone of the tibial spongiosa; however, a

52

transient increase was detected in the subchondral trabecular bone volume and density.

53

Conclusions: The lack of WBV-induced changes in the underlying subchondral bone

54

suggests that damage to the articular cartilage may be secondary to the meniscal injury

55

we detected. Our findings underscore the need for further studies to assess the safety of

56

WBV in the human population to avoid long-term joint damage.

AC C

EP

TE D

M AN U

SC

RI PT

36

57 58 59

2

ACCEPTED MANUSCRIPT

INTRODUCTION

60

Mechanical loading is critical for maintaining joint tissue homeostasis. Articular

62

cartilage health is paramount to overall joint function and is typically well maintained

63

under normal physiological loading [1]. Moderate physical activity has been shown to

64

sustain cartilage health in human [2] and animal studies [3]. In contrast, insufficient

65

mechanical loading can lead to articular cartilage degeneration; progressive thinning of

66

the articular cartilage is detected in patients with paralysis compared to healthy mobile

67

individuals [4]. Cartilage loss can be detected relatively quickly following changes in

68

mechanical loading; partial load bearing associated with walking with crutches resulted in

69

atrophy of knee joint articular cartilage in just 7 weeks [5]. Conversely, mechanical

70

overloading is perhaps the best-studied initiator of joint degeneration. Articular cartilage

71

is highly susceptible to osteoarthritis (OA) if subjected to acute [6] or chronic [7] impact

72

loading. Altered mechanical loading following injury to the knee joint, modelled in

73

rodents using surgically-induced joint destabilization (e.g. destablization of the medial

74

menicus, anterior cruciate ligament transection), leads to OA-like phenotypes [8].

75

Although the pathological processes mediating the progression of OA remain the focus of

76

ongoing investigation, OA is now generally considered a disease of the entire joint in

77

which articular cartilage loss, meniscus and ligament damage and the periarticular bone

78

changes contribute to the overall disease [9].

SC

M AN U

TE D

EP

AC C

79

RI PT

61

Whole-body vibration (WBV) training was originally implemented clinically as a

80

therapeutic approach to increase bone mass for the treatment of osteoporosis [10], and is

81

capable of increasing muscle strength and power [11] [12]. Underlying these biological

82

effects is the principle that WBV increases bone remodeling [13] and motor unit

3

ACCEPTED MANUSCRIPT

recruitment thresholds in muscle [14]. Contributing to the increasing popularity of WBV

84

training in the fitness industry, studies have reported beneficial effects of WBV on

85

neuromuscular performance in healthy individuals, including flexibility and strength [15].

86

However, recent studies suggest these effects may vary based on age, gender or health

87

status [16].

RI PT

83

Clinical trials have recently evaluated the therapeutic effect of WBV on joint pain

89

and functional performance in patients with chronic knee OA [17]. A number of studies

90

report beneficial effects including reduced pain intensity, increased muscle strength, and

91

reductions in the plasma concentrations of inflammatory markers in OA patients

92

following WBV training [18, 19]. Conversely, some trials reported no significant

93

improvements in pain intensity or other parameters among knee OA patients following

94

WBV [19]. Importantly, in these studies joint health was not directly assessed.

95

Quantification of the direct effects of WBV on knee joint tissues are limited; one study

96

reported the ability of WBV to prevent loss of tibial articular cartilage thickness

97

(determined by MRI) following prolonged immobilization in healthy male subjects [20].

98

A recent meta-analysis evaluating the use of WBV in knee OA management concluded

99

that there were no significant differences in functional performance or patient self-

100

reported pain following WBV, likely related to the diversity of the protocols implemented

101

[19]. Together, these findings suggest the need for for further studies evaluating the

102

safety and efficacy of WBV prior to its therapeutic implementation for the treatment of

103

OA.

AC C

EP

TE D

M AN U

SC

88

104

To address these issues, recent studies from our group evaluated the effect of

105

WBV on joint health using the mouse as a preclinical model. Our findings demonstrated

4

ACCEPTED MANUSCRIPT

that 4 weeks of repeated exposure of healthy young mice to WBV, using protocols

107

reported to be beneficial to bone, induced OA-like knee joint damage marked by

108

degeneration of both the meniscus and articular cartilage as well as activation of MMP-

109

mediated matrix degradation [21]. Building from these findings, the present study was

110

designed to determine: i) if knee joint degeneration was progressive with extended

111

exposure to WBV or was potentiated by cessation of vibration, and ii) if WBV promoted

112

changes in the subchondral bone that would constitute a potential mechanism for the OA-

113

like pathology detected in the knee.

M AN U

SC

RI PT

106

METHODS

114

Whole-Body Vibration. All procedures were approved by the Animal Use Committee at

116

The University of Western Ontario and this study was conducted in accordance the

117

ARRIVE guidelines [22]. Based on protocols of WBV used in our previous studies to

118

model those used clinically [21], 10-week-old male CD-1 mice were exposed to vertical

119

sinusoidal vibrations (frequency of 45 Hz, peak-to-peak amplitude of 74 µm, and 0.3 g

120

peak acceleration) for 30 min/day, 5 days/week for 4 weeks, 8 weeks, or 4 weeks WBV

121

followed by 4 weeks recovery, using a previously described vibration platform [23, 24].

122

Age-matched controls were housed in identical chambers on a sham (non-vibrated)

123

platform to replicate handling and environmental conditions. Following WBV, mice were

124

returned to conventional housing and monitored daily. Twenty-four hours after the final

125

exposure to WBV, mice were euthanized by a lethal dose of sodium pentobarbital. Spinal

126

tissues were isolated to characterize effects of WBV on the intervertebral disc and

127

vertebral bone, detailed in [25].

AC C

EP

TE D

115

128

5

ACCEPTED MANUSCRIPT

Micro-computed Tomography (micro-CT). Right hindlimbs (mid-femur to ankle) were

130

isolated and fixed in 4% paraformaldehyde (PFA) for 24 h and embedded in 1% agarose

131

in PBS in 50 mL conical tubes. Tissues were scanned using a laboratory micro-CT

132

scanner (eXplore Locus, GE Healthcare Biosciences) with an internal calibrating

133

phantom composed of air, water, and synthetic bone-mimicking epoxy (SB3 2.8 x 3.4 x 8

134

mm; Gammex RMI). The scanning protocol consisted of 900 projection images acquired

135

at a 0.4º angle increment, obtained over 165 minutes of gantry rotation. The x-ray tube

136

potential was 80 kVp with a tube current of 450 µA, an exposure time of 4,500 msec, and

137

2 frames per view angle averaged. Images were reconstructed into 3D volumes with 20

138

µm isotropic voxel spacing and linearly rescaled into Hounsfield units using the internal

139

calibration standards, as previously described [26]. Using three-dimensional visualization

140

and analysis software (MicroView, GE Healthcare Biosciences) full volumes were

141

cropped to contain only the tibia and fibula, and reoriented to the same axes. Regions of

142

interest (ROI) were manually outlined to define the subchondral bone plate, subchondral

143

trabecular bone, primary and secondary spongiosa of the tibia, where each compartment

144

was marked within a series of 2D planes and splined together to generate 3D ROIs. Bone

145

morphometry was quantified using the Bone Analysis tool in MicroView, as previously

146

reported [27].

SC

M AN U

TE D

EP

AC C

147

RI PT

129

148

Histological Assessment and Scoring. Following micro-CT, right hind limbs were

149

decalcified for 10 days in 5% EDTA in PBS (pH 7.0), processed, embedded in paraffin

150

wax, and 5 µm coronal sections were obtained. Serial sections were taken every 50 µm to

151

encompass all weight-bearing areas of the femoro-tibial joint. Using established protocols

6

ACCEPTED MANUSCRIPT

[23, 28], sections were stained with 0.1% Safranin O/0.02% fast green to detect

153

proteoglycans and glycosaminoglycans. Sections were imaged on a Leica DM1000

154

microscope with Leica Application Suite (Leica Microsystems). The four quadrants of

155

the knee joint were scored by three independent blinded observers for degeneration using

156

the murine Osteoarthritis Research Society International (OARSI) recommended

157

histopathologic scale [29].

RI PT

152

SC

158

Statistical Analyses. Data are presented from experiments conducted with 5-6 mice per

160

treatment group. Parameters from mice exposed to WBV were compared to those from

161

age-matched non-vibrated sham controls using a parametric, two-tailed, unpaired t-test

162

with a Welch’s correction (4 week time point), or one-way ANOVA with a Tukey post-

163

hoc test (8 week time point). P < 0.05 was considered significant.

M AN U

159

165

TE D

164

RESULTS

To assess the effects of daily WBV on the knee joint and determine if protocols of

167

WBV used clinically [30] induce long-term joint damage, 10-week-old male CD-1 mice

168

were exposed to sinusoidal WBV (30 min/day, 5 days/week, at 45 Hz, 0.3 g peak

169

acceleration) for 4 weeks, 8 weeks, or 4 weeks of WBV followed by 4 weeks without

170

WBV (recovery) (Figure 1). Consistent with our recent study [21], analysis of Safranin

171

O/fast green stained histological sections demonstrated that 4 weeks of WBV induced

172

knee joint damage in 3 of 6 mice (OARSI score >2), localized to the medial joint

173

compartment in 2 mice and the lateral femoral condyle in the third mouse. In 1 of these

174

mice, articular damage was detected on both the femoral and tibial condyles. The OARSI

AC C

EP

166

7

ACCEPTED MANUSCRIPT

semi-quantitative histologic scoring system assigns a grade between 0-6 to each joint

176

quadrant, where a grade of 2 is characterized by loss of articular cartilage with vertical

177

clefts extending beyond the superficial layer and some loss of surface lamina. WBV-

178

induced damage included evidence of meniscal tears, focal articular cartilage defects with

179

erosion down to the tide mark, loss of glycosaminoglycans in the superficial layer, and

180

osteophyte formation on the medial tibial plateau and femoral condyle (Figure 2A). No

181

joint damage (OARSI >2) was detected in any age-matched sham controls. Despite these

182

findings, no significant differences were detected in the average maximum OARSI scores

183

between mice exposed to 4 weeks of WBV and sham controls.

M AN U

SC

RI PT

175

To determine if WBV-induced joint damage was progressive with continued daily

185

WBV or potentiated by inclusion of a recovery period, mice were exposed to either 8

186

weeks of WBV or 4 weeks of WBV followed by 4 weeks without (recovery). Following

187

8 weeks of WBV, 5 of 6 mice demonstrated articular cartilage damage (OARSI score

188

>2), localized to the medial joint compartment in 2 mice, the lateral joint compartment in

189

1 mouse, and both the medial and lateral joint compartments in 1 mouse. In 3 of these

190

mice, articular damage was detected on both the femoral and tibial condyles. Following 8

191

weeks WBV, joint damage was characterized by severe loss of meniscal tissue, focal

192

articular

193

glycosaminoglycans in the superficial layer of articular cartilage adjacent to focal defects,

194

and osteophyte formation on the medial tibial plateau and femoral condyle (Figure 2B).

195

Focal articular cartilage damage (OARSI score >2) was also detected in 2 of 6 age-

196

matched sham controls, localized to the medial tibial plateau. Exposure of mice to 8

AC C

EP

TE D

184

cartilage

lesions

with

erosion

8

to

the

subchondral

bone,

loss

of

ACCEPTED MANUSCRIPT

197

weeks of WBV induced a significant increase in the average maximum OARSI score in

198

the lateral tibial plateau compared to age-matched sham controls. Articular cartilage damage (OARSI score >2) was detected in only 1 of 5 mice

200

subjected to 4 weeks of WBV followed by 4 weeks of recovery (WBV/REC). Of note

201

however, in this animal the meniscal and articular cartilage damage detected was similar

202

to that of mice exposed to 4 weeks WBV (Figure 2B), suggesting that inclusion of a

203

recovery period failed to mitigate the damage induced by WBV in this animal.

SC

RI PT

199

We next examined if WBV promoted changes to the subchondral bone

205

microarchitecture that would correlate with WBV-induced soft tissue damage. Using

206

high-resolution micro-CT, we first assessed the subchondral bone plate subjacent to the

207

articular cartilage of the tibial plateau since this region showed the most pronounced

208

WBV-induced articular cartilage degeneration (Figure 3A). Subchondral bone plate

209

thickness was measured within a region of interest measuring 500 µm in mediolateral

210

length and 1,150 µm ventrodorsal length, centred on the medial tibial plateau, based on

211

previously reported protocols [31]. No significant differences in bone plate thickness

212

(mm), bone mineral content (BMC, mg), or bone mineral density (BMD, mg/cm3) were

213

detected in the subchondral bone of mice exposed to 4 weeks of WBV (Figure 3B), 8

214

weeks of WBV, or 4 weeks of WBV followed by 4 weeks recovery (Figure 3C).

TE D

EP

AC C

215

M AN U

204

Bone microarchitecture was next assessed in the epiphyseal subchondral

216

trabecular bone in the medial compartment of the tibia. For each specimen, a region of

217

interest was manually outlined in serial planes to capture trabecular bone within this area;

218

the lowest point of the extensor sulcus was extended as a boundary, and the subchondral

9

ACCEPTED MANUSCRIPT

bone plate, growth plate and cortical bone were excluded (Figure 4A). Exposure of mice

220

to 4 weeks of WBV induced a significant increase in the subchondral trabecular bone

221

volume fraction (BVF) and trabecular BMD (mg HA/cm3) compared to non-vibrated

222

sham controls (Figure 4B). In contrast, no significant differences were detected in BVF

223

or BMD (mg HA/cm3) of the subchondral trabecular bone in mice exposed to 8 weeks of

224

WBV or 4 weeks of WBV followed by 4 weeks recovery compared to age-matched sham

225

controls (Figure 4C).

SC

RI PT

219

Previous studies suggest the primary and secondary spongiosa within long bones

227

as sites that demonstrate robust adaptive bone remodelling in response to dynamic

228

mechanical loading [32]. We therefore assessed bone microarchitecture within both the

229

primary and secondary spongiosa of the tibia, and quantified bone properties within each

230

separately to delineate region-specific changes. A region of interest was defined to

231

include the primary spongiosa, manually outlined in serial planes to capture the trabecular

232

bone within 50 µm distal to the growth plate and exclude the cortical bone (Figure 5A).

233

No significant differences were detected in trabecular BVF or trabecular BMD in the

234

primary spongiosa following 4 weeks of WBV (Figure 5B), 8 weeks of WBV or 4 weeks

235

of WBV followed by 4 weeks recovery (Figure 5C). We then moved distally and

236

measured the secondary spongiosa of the tibia, manually outlined in serial planes to

237

capture the trabecular bone within an area extending 1mm from the most distal point of

238

the primary spongiosa, and exclude the cortical bone (Figure 6A). No significant

239

differences were detected in trabecular BVF or trabecular BMD in the secondary

240

spongiosa following 4 weeks of WBV (Figure 6B), 8 weeks of WBV or 4 weeks of

241

WBV followed by 4 weeks recovery (Figure 6C).

AC C

EP

TE D

M AN U

226

10

ACCEPTED MANUSCRIPT

DISCUSSION

242

The reported beneficial effects of WBV on bone mass and muscle strength have

244

led to increased use of WBV devices both clinically and in the sports and fitness industry

245

to treat pain, increase joint mobility, and enhance neuromuscular performance. However,

246

recent conflicting reports have highlighted inconsistencies in the anabolic effects of WBV

247

on bone, in both animal models [21, 33] and human clinical trials [30, 34]. In this report,

248

we investigated the mechanisms associated with our previous reports of WBV-induced

249

damage to the murine femoro-tibial joint (knee) following exposure of mice to protocols

250

of WBV that model those used clinically [21]. Based on histological characterization, our

251

findings suggest progressive degeneration of the meniscus and articular cartilage

252

following exposure to WBV. Between 4 and 8 weeks, progressive degeneration was

253

marked by severe loss of meniscus tissue, an increase in the both the number of mice

254

showing articular cartilage defects and the extent of the articular cartilage damage,

255

resulting in a significant increase in the average maximum OARSI score. Based on

256

detailed micro-CT based examination of the tibia, we show that WBV-induced joint

257

damage is not associated with either thickening of the subchondral bone plate or changes

258

in the trabecular bone of the primary or secondary spongiosa.

SC

M AN U

TE D

EP

AC C

259

RI PT

243

Although loss of articular cartilage is the hallmark pathological feature in OA,

260

disease progression involves biomechanical and biochemical alterations in the

261

composition, structure and functional properties of the articular cartilage and periarticular

262

bone [35]. Within the bone compartment, OA-associated changes include increases in

263

subchondral cortical bone thickness, early bone remodeling in the subchondral trabeculae

264

resulting in decreased trabecular bone mass, formation of osteophytes at the joint margins

11

ACCEPTED MANUSCRIPT

and expansion of the zone of calcified cartilage leading to thinning of the articular

266

cartilage [36, 37]. Studies from our group and others have demonstrated that rodent

267

models of OA likewise demonstrate hallmark bone changes, including thickening of the

268

subchondral bone plate and decreased BMD in the subchondral trabecular bone,

269

particularly in regions adjacent to affected articular cartilage [38-40]. Given the interplay

270

between these tissues, the current study used high-resolution micro-CT to determine if

271

focal defects in articular cartilage induced by WBV were associated with changes to the

272

periarticular bone. Using manually defined regions of interest, we did not detect

273

significant differences in bone microarchitecture in either the subchondral bone plate or

274

primary and secondary spongiosa of the tibia. These findings are consistent with

275

quantification of bone microarchitecture in the vertebral bone of these mice, which

276

likewise demonstrated no change following exposure to WBV [25]. Our findings did

277

however demonstrate a transient increase in subchondral trabecular bone mass after 4

278

weeks of WBV, although no differences were detected between mice exposed to 8 weeks

279

WBV compared to controls. Interestingly, we noted that following 4 weeks of WBV,

280

mice with the highest OARSI scores also showed the highest bone plate thickness, bone

281

mineral content and bone mineral density. This association was not however seen in mice

282

following 8 weeks of WBV and may therefore reflect a transient adaptive response to

283

vibration loading.

SC

M AN U

TE D

EP

AC C

284

RI PT

265

Sclerostin is a marker of osteocyte differentiation that serves as a negative

285

regulator of bone formation by inhibiting the osteo-anabolic Wnt signalling pathway. In

286

humans, loss of function mutations of the human SOST gene (encoding sclerostin) lead to

287

sclerosteosis [41, 42], and deletion of Sost in mice results in a high bone mass phenotype

12

ACCEPTED MANUSCRIPT

[43]. Despite increased subchondral bone, Sost-deficient mice do not demonstrate

289

changes to the articular cartilage in the knee associated with OA [44]. Interestingly, using

290

a femur osteotomy model of fracture healing in aged mice, a recent study demonstrated

291

that WBV increased the expression of Sost in osteocytes, associated with impaired

292

fracture healing [45]. These findings may be in keeping with recent studies by our group

293

[21] and others [46] highlighting inconsistencies in the previously reported anabolic

294

effects of WBV on bone, as well as the confounding effects of age and hormone levels.

295

Furthermore, these findings highlight the potential role of Wnt signaling as a candidate

296

regulator of WBV mechanotransduction that should be investigated in joint tissues such

297

as meniscus and articular cartilage.

M AN U

SC

RI PT

288

Our previous study reported meniscal damage in mice following 2 weeks of

299

WBV, changes that preceded the detection of overt articular cartilage damage [21].

300

Together with the current study, these findings suggest a correlation between meniscus

301

damage and the induction of articular cartilage degeneration following WBV. As such,

302

we speculate that WBV-induced articular cartilage damage may occur following joint

303

destabilization, rather than as a result of changes in the mechanical environment resulting

304

from increased subchondral bone characteristic of OA. Meniscal damage is often

305

associated with pain and joint dysfunction [47], and loss of meniscal function leads to

306

rapid and progressive osteoarthritis [48]. Damage or partial loss of the meniscus

307

influences the displacement and deformation of the meniscus under load [49], resulting in

308

increased focal-point stresses on the articular cartilage surface of the tibial plateau [50]

309

that we propose lead to the formation of the focal cartilage defects we detected in mice

310

exposed to WBV.

AC C

EP

TE D

298

13

ACCEPTED MANUSCRIPT

An interesting aspect of this and our previous study is the fact that while severe

312

knee joint damage was induced by WBV (including severe damage to the meniscus, focal

313

articular cartilage defects extending to the subchondral bone, and osteophyte formation

314

on both the femoral and tibial condyles), damage was only detected in a subset of mice at

315

each time point examined. As such, statistical analysis of OARSI scores demonstrated

316

significant WBV-induced damage only following 8 weeks of WBV in a single joint

317

compartment. The absence of a significant increase in OARSI scores following 4 weeks

318

WBV, or in multiple joint compartments at either time point may be associated with the

319

limited sample size used in this study. Of note, the individual cages in which mice are

320

housed during WBV allow for free range of movement. We have observed that during

321

vibration, mice alternate through a variety of postures including rearing on hind limbs,

322

free movement on all 4 limbs, grooming and sleeping. It is tempting to speculate that

323

increased time spent in specific postures may be directly associated with increased joint

324

damage since the magnitude of vibration experienced at skeletal sites likely varies due to

325

muscle contraction, soft tissue damping, or resonance effects [51]. We have recently

326

developed non-invasive techniques to quantify vibration amplitude in the skeleton of live

327

mice using radiographic images of surgically implanted metal marker beads [24]. Future

328

studies will utilize this approach to quantify the vibration amplitude and deformations at

329

specific sites and in various postures. This will allow us to directly relate the biological

330

response of joint tissues to WBV to its biomechanical effects.

AC C

EP

TE D

M AN U

SC

RI PT

311

331

While this and previous studies by our group [21, 25] point to alarming potential

332

for damage induced by whole-body vibration on joint tissues, the preclinical model used

333

is associated with inherent limitations. First, these studies were conducted using a single

14

ACCEPTED MANUSCRIPT

protocol of WBV that does not effectively capture the wide range of parameters currently

335

being used. The variables that control the intensity of vibration are the frequency and

336

amplitude (peak-to-peak displacement), which together determine peak acceleration.

337

Clinical protocols of WBV vary widely in frequency (2.5-90 Hz) and acceleration (0.3-

338

3.6 g), and are implemented for periods ranging from one week to one year. Second, we

339

examined the effects of WBV using young, skeletally mature male CD-1 mice. Previous

340

studies have demonstrated that the ability of WBV to increase bone mineral density and

341

bone formation is dependent on age and genetic background in mouse models [52, 53].

342

As such, physiological factors including age, sex, and genetic background may likewise

343

affect the response of joint tissues to WBV. In fact, a recent study from our group

344

reported a lack of WBV-induced degeneration in either the knee or intervertebral discs of

345

C57Bl/6 mice exposed to the same protocols of WBV [54]. Taken together, these studies

346

underscore the need for further rigorous research to assess the relative contribution of

347

vibration parameters and physiological factors in modulating the effects of WBV on joint

348

health.

TE D

M AN U

SC

RI PT

334

In conclusion, WBV results in progressive degeneration of the meniscus and

350

articular cartilage, damage that may not be repaired following removal of the stimulus.

351

Furthermore, we demonstrate that WBV promotes a transient increase in subchondral

352

trabecular bone; however, no change was detected in morphometry of the subchondral

353

bone plate or primary and secondary spongiosa of the tibia.

AC C

EP

349

354 355

ACKNOWLEDGEMENTS

15

ACCEPTED MANUSCRIPT

356

We thank Dr. Kim Beaucage for her contributions, in particular for her assistance with

357

micro-CT training and protocols.

358

AUTHOR CONTRIBUTIONS

360

All authors were involved in drafting the article or revising it critically for content, and

361

all authors approved the final version. Drs. McCann and Séguin had access to all of the

362

data in the study and take responsibility for the integrity of the data and the accuracy of

363

the data analysis.

364

Study conception and design: M.R.M., M.A.P., A.R., D.W.H., F.B., S.J.D., and C.A.S.

365

Acquisition of data: M.R.M., C.Y., M.A.P., A.R.

366

Analysis and interpretation of data: M.R.M., C.Y., M.A.P., A.R., S.I.P., D.W.H., F. B.,

367

S.J.D., and C.A.S.

SC

M AN U

TE D

368

EP

Role of the funding source: This work was funded by the Natural Sciences and Engineering Research Council of Canada [371546] and Canadian Institutes of Health Research [132377, 312615, 86574, 133575]. Infrastructure was provided in part by the Canada Foundation for Innovation-Leaders Opportunity Fund [25086 to C.A.S]. M.R.M. & M.A.P were supported by CIHR Doctoral Awards. M.R.M., M.A.P. & A.R. were supported in part by the Joint Motion Program - A CIHR Training Program in Musculoskeletal Health Research and Leadership. A.R. was supported by a Doctoral Award from The Arthritis Society. F.B. is the Canada Research Chair in Musculoskeletal Research. C.A.S. is supported by a CIHR New Investigator Award and Early Researcher Award from the Ontario Ministry of Research and Innovation.

AC C

369 370 371 372 373 374 375 376 377 378 379 380 381 382 383

RI PT

359

Competing interest statement: The authors declare that they have no competing interests.

16

ACCEPTED MANUSCRIPT

REFERENCES

384 1.

RI PT

2.

Leong DJ, Hardin JA, Cobelli NJ, Sun HB. Mechanotransduction and cartilage integrity. Ann N Y Acad Sci 2011; 1240: 32-37. Manninen P, Riihimaki H, Heliovaara M, Suomalainen O. Physical exercise and risk of severe knee osteoarthritis requiring arthroplasty. Rheumatology (Oxford) 2001; 40: 432-437. Musumeci G, Castrogiovanni P, Trovato FM, Imbesi R, Giunta S, Szychlinska MA, et al. Physical activity ameliorates cartilage degeneration in a rat model of aging: A study on lubricin expression. Scand J Med Sci Sports 2014. Vanwanseele B, Eckstein F, Knecht H, Stussi E, Spaepen A. Knee cartilage of spinal cord-injured patients displays progressive thinning in the absence of normal joint loading and movement. Arthritis Rheum 2002; 46: 2073-2078. Hinterwimmer S, Krammer M, Krotz M, Glaser C, Baumgart R, Reiser M, et al. Cartilage atrophy in the knees of patients after seven weeks of partial load bearing. Arthritis Rheum 2004; 50: 2516-2520. Borrelli J, Jr., Silva MJ, Zaegel MA, Franz C, Sandell LJ. Single high-energy impact load causes posttraumatic OA in young rabbits via a decrease in cellular metabolism. J Orthop Res 2009; 27: 347-352. Ko FC, Dragomir C, Plumb DA, Goldring SR, Wright TM, Goldring MB, et al. In vivo cyclic compression causes cartilage degeneration and subchondral bone changes in mouse tibiae. Arthritis Rheum 2013; 65: 1569-1578. Fang H, Beier F. Mouse models of osteoarthritis: modelling risk factors and assessing outcomes. Nat Rev Rheumatol 2014; 10: 413-421. Loeser RF, Goldring SR, Scanzello CR, Goldring MB. Osteoarthritis: a disease of the joint as an organ. Arthritis Rheum 2012; 64: 1697-1707. Wysocki A, Butler M, Shamliyan T, Kane RL. Whole-body vibration therapy for osteoporosis: state of the science. Ann Intern Med 2011; 155: 680-686, W206-613. Wilcock IM, Whatman C, Harris N, Keogh JW. Vibration training: could it enhance the strength, power, or speed of athletes? J Strength Cond Res 2009; 23: 593-603. Lam FM, Lau RW, Chung RC, Pang MY. The effect of whole body vibration on balance, mobility and falls in older adults: a systematic review and meta-analysis. Maturitas 2012; 72: 206-213. Totosy de Zepetnek JO, Giangregorio LM, Craven BC. Whole-body vibration as potential intervention for people with low bone mineral density and osteoporosis: a review. J Rehabil Res Dev 2009; 46: 529-542. Pollock RD, Woledge RC, Martin FC, Newham DJ. Effects of whole body vibration on motor unit recruitment and threshold. J Appl Physiol (1985) 2012; 112: 388-395. Karatrantou K, Gerodimos V, Dipla K, Zafeiridis A. Whole-body vibration training improves flexibility, strength profile of knee flexors, and hamstrings-to-quadriceps strength ratio in females. J Sci Med Sport 2013; 16: 477-481. Segal NA, Glass NA, Shakoor N, Wallace R. Vibration platform training in women at risk for symptomatic knee osteoarthritis. PM R 2013; 5: 201-209; quiz 209. Salmon JR, Roper JA, Tillman MD. Does acute whole-body vibration training improve the physical performance of people with knee osteoarthritis? J Strength Cond Res 2012; 26: 2983-2989. Park YG, Kwon BS, Park JW, Cha DY, Nam KY, Sim KB, et al. Therapeutic effect of whole body vibration on chronic knee osteoarthritis. Ann Rehabil Med 2013; 37: 505-515.

3.

4.

8. 9. 10. 11. 12.

13.

14.

TE D

7.

EP

6.

M AN U

SC

5.

AC C

385 386 387 388 389 390 391 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431

15.

16. 17.

18.

17

ACCEPTED MANUSCRIPT

23.

24.

25.

26.

27.

28.

29.

RI PT

SC

22.

M AN U

21.

TE D

20.

Li X, Wang XQ, Chen BL, Huang LY, Liu Y. Whole-Body Vibration Exercise for Knee Osteoarthritis: A Systematic Review and Meta-Analysis. Evid Based Complement Alternat Med 2015; 2015: 758147. Liphardt AM, Mundermann A, Koo S, Backer N, Andriacchi TP, Zange J, et al. Vibration training intervention to maintain cartilage thickness and serum concentrations of cartilage oligometric matrix protein (COMP) during immobilization. Osteoarthritis Cartilage 2009; 17: 1598-1603. McCann MR, Patel P, Pest MA, Ratneswaran A, Lalli G, Beaucage KL, et al. Repeated exposure to high-frequency low-amplitude vibration induces degeneration of murine intervertebral discs and knee joints. Arthritis Rheumatol 2015; 67: 21642175. Kilkenny C, Browne WJ, Cuthill IC, Emerson M, Altman DG. Improving bioscience research reporting: the ARRIVE guidelines for reporting animal research. PLoS Biol 2010; 8: e1000412. McCann MR, Patel P, Beaucage KL, Xiao Y, Bacher C, Siqueira WL, et al. Acute vibration induces transient expression of anabolic genes in the murine intervertebral disc. Arthritis Rheum 2013; 65: 1853-1864. Hu Z, Welch I, Yuan X, Pollmann SI, Nikolov HN, Holdsworth DW. Quantification of mouse in vivo whole-body vibration amplitude from motion-blur using x-ray imaging. Phys Med Biol 2015; 60: 6423-6439. McCann MR. Veras M, Young C, Lalli G, Patel P, Leitch KM, Holdsworth DW, Beier F, Dixon SJ, Sé guin CA. Whole-body vibration of mice induces progressive degeneration of intervertebral discs associated with increased expression of IL-1β and multiple matrix degrading enzymes. Osteoarthritis Cartilage; Manuscript Number: OAC6969R1. Du LY, Umoh J, Nikolov HN, Pollmann SI, Lee TY, Holdsworth DW. A quality assurance phantom for the performance evaluation of volumetric micro-CT systems. Phys Med Biol 2007; 52: 7087-7108. Ulici V, Hoenselaar KD, Agoston H, McErlain DD, Umoh J, Chakrabarti S, et al. The role of Akt1 in terminal stages of endochondral bone formation: angiogenesis and ossification. Bone 2009; 45: 1133-1145. Pest MA, Russell BA, Zhang YW, Jeong JW, Beier F. Disturbed cartilage and joint homeostasis resulting from a loss of mitogen-inducible gene 6 in a mouse model of joint dysfunction. Arthritis Rheumatol 2014; 66: 2816-2827. Glasson SS, Chambers MG, Van Den Berg WB, Little CB. The OARSI histopathology initiative - recommendations for histological assessments of osteoarthritis in the mouse. Osteoarthritis Cartilage 2010; 18 Suppl 3: S17-23. Slatkovska L, Alibhai SM, Beyene J, Hu H, Demaras A, Cheung AM. Effect of 12 months of whole-body vibration therapy on bone density and structure in postmenopausal women: a randomized trial. Ann Intern Med 2011; 155: 668-679, W205. Botter SM, van Osch GJ, Clockaerts S, Waarsing JH, Weinans H, van Leeuwen JP. Osteoarthritis induction leads to early and temporal subchondral plate porosity in the tibial plateau of mice: an in vivo microfocal computed tomography study. Arthritis Rheum 2011; 63: 2690-2699. Sugiyama T, Price JS, Lanyon LE. Functional adaptation to mechanical loading in both cortical and cancellous bone is controlled locally and is confined to the loaded bones. Bone 2010; 46: 314-321. Wenger KH, Freeman JD, Fulzele S, Immel DM, Powell BD, Molitor P, et al. Effect of whole-body vibration on bone properties in aging mice. Bone 2010; 47: 746-755.

EP

19.

AC C

432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481

30.

31.

32.

33.

18

ACCEPTED MANUSCRIPT

40.

41.

42.

43.

44.

45.

46.

RI PT

39.

SC

37. 38.

M AN U

36.

TE D

35.

Liphardt AM, Schipilow J, Hanley DA, Boyd SK. Bone quality in osteopenic postmenopausal women is not improved after 12 months of whole-body vibration training. Osteoporos Int 2015; 26: 911-920. Goldring SR. Alterations in periarticular bone and cross talk between subchondral bone and articular cartilage in osteoarthritis. Ther Adv Musculoskelet Dis 2012; 4: 249-258. Patel V, Issever AS, Burghardt A, Laib A, Ries M, Majumdar S. MicroCT evaluation of normal and osteoarthritic bone structure in human knee specimens. J Orthop Res 2003; 21: 6-13. Goldring MB, Goldring SR. Osteoarthritis. J Cell Physiol 2007; 213: 626-634. Holland JC, Brennan O, Kennedy OD, Mahony NJ, Rackard S, O'Brien FJ, et al. Examination of osteoarthritis and subchondral bone alterations within the stifle joint of an ovariectomised ovine model. J Anat 2013; 222: 588-597. Loeser RF, Olex AL, McNulty MA, Carlson CS, Callahan M, Ferguson C, et al. Disease progression and phasic changes in gene expression in a mouse model of osteoarthritis. PLoS One 2013; 8: e54633. McErlain DD, Appleton CT, Litchfield RB, Pitelka V, Henry JL, Bernier SM, et al. Study of subchondral bone adaptations in a rodent surgical model of OA using in vivo micro-computed tomography. Osteoarthritis Cartilage 2008; 16: 458-469. Brunkow ME, Gardner JC, Van Ness J, Paeper BW, Kovacevich BR, Proll S, et al. Bone dysplasia sclerosteosis results from loss of the SOST gene product, a novel cystine knot-containing protein. Am J Hum Genet 2001; 68: 577-589. Balemans W, Ebeling M, Patel N, Van Hul E, Olson P, Dioszegi M, et al. Increased bone density in sclerosteosis is due to the deficiency of a novel secreted protein (SOST). Hum Mol Genet 2001; 10: 537-543. Li X, Ominsky MS, Niu QT, Sun N, Daugherty B, D'Agostin D, et al. Targeted deletion of the sclerostin gene in mice results in increased bone formation and bone strength. J Bone Miner Res 2008; 23: 860-869. Roudier M, Li X, Niu QT, Pacheco E, Pretorius JK, Graham K, et al. Sclerostin is expressed in articular cartilage but loss or inhibition does not affect cartilage remodeling during aging or following mechanical injury. Arthritis Rheum 2013; 65: 721-731. Wehrle E, Liedert A, Heilmann A, Wehner T, Bindl R, Fischer L, et al. The impact of low-magnitude high-frequency vibration on fracture healing is profoundly influenced by the oestrogen status in mice. Dis Model Mech 2015; 8: 93-104. Kiel DP, Hannan MT, Barton BA, Bouxsein ML, Sisson E, Lang T, et al. Low-Magnitude Mechanical Stimulation to Improve Bone Density in Persons of Advanced Age: A Randomized, Placebo-Controlled Trial. J Bone Miner Res 2015; 30: 1319-1328. McNulty AL, Guilak F. Mechanobiology of the meniscus. J Biomech 2015; 48: 14691478. Fairbank TJ. Knee joint changes after meniscectomy. J Bone Joint Surg Br 1948; 30B: 664-670. Yao J, Snibbe J, Maloney M, Lerner AL. Stresses and strains in the medial meniscus of an ACL deficient knee under anterior loading: a finite element analysis with imagebased experimental validation. J Biomech Eng 2006; 128: 135-141. Kurosawa H, Fukubayashi T, Nakajima H. Load-bearing mode of the knee joint: physical behavior of the knee joint with or without menisci. Clin Orthop Relat Res 1980: 283-290.

EP

34.

AC C

482 483 484 485 486 487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529

47.

48.

49.

50.

19

ACCEPTED MANUSCRIPT

52.

53.

54.

Zhao L, Dodge T, Nemani A, Yokota H. Resonance in the mouse tibia as a predictor of frequencies and locations of loading-induced bone formation. Biomech Model Mechanobiol 2013. Xie L, Jacobson JM, Choi ES, Busa B, Donahue LR, Miller LM, et al. Low-level mechanical vibrations can influence bone resorption and bone formation in the growing skeleton. Bone 2006; 39: 1059-1066. Xie L, Rubin C, Judex S. Enhancement of the adolescent murine musculoskeletal system using low-level mechanical vibrations. J Appl Physiol (1985) 2008; 104: 1056-1062. Kerr G.J. McCann MR, Branch JK., Ratneswaran A, Pest MA., Beier F, Holdsworth DW, Dixon SJ, Seguin CA. C57BL/6 mice are resistant to whole-body vibration-induced joint degeneration. In Revison at Osteoarthritis Cartilage (Manuscript Number OAC7134).

RI PT

51.

SC

530 531 532 533 534 535 536 537 538 539 540 541 542

AC C

EP

TE D

M AN U

543

20

ACCEPTED MANUSCRIPT

FIGURE LEGENDS

545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 572

Figure 1: Schematic representation of experimental design. Whole-body vibration was delivered by an electromagnetic shaker regulated by an open-loop controller set to parameters of 45 Hz, 0.3 g peak acceleration, with a peak-to-peak amplitude of 74 µm. Ten-week-old, wild-type male mice (n=5-6 mice/group) were exposed to 30 mins of WBV/day, 5 days/week for 4 weeks, or 8 weeks, or 4 weeks of WBV followed by 4 weeks recovery (WBV/REC). 24 h following the last exposure to WBV, mice were sacrificed and knee joint tissues were assessed compared to age and gender-matched nonvibrated controls.

SC

M AN U

TE D

EP

574 575 576 577 578 579 580 581 582 583 584 585 586

Figure 2: Histological appearance of mouse medial knee joints after exposure to WBV. A. Representative coronal sections of the medial compartment of the knee joint, stained with Safranin O/fast green from mice exposed to WBV for 4 weeks and nonvibrated sham controls. Images are oriented with the medial meniscus on the left, the femoral condyle on top, and the tibial plateau on the bottom. Microfissures were detected within the medial meniscus in 3 of 6 mice exposed to 4 weeks of WBV (arrows). No meniscal damage was detected in control mice not exposed to WBV. Focal defects in the articular cartilage were detected in 2 of the 5 mice subjected to WBV (arrow heads). B. Histological sections from mice exposed to 8 weeks of WBV or 4 weeks of WBV followed by 4 weeks of recovery and non-vibrated age-matched sham controls. Asterisks indicate developing osteophytes at the medial margin of the joints. No articular cartilage degeneration was detected in the controls. 8 weeks of WBV resulted in severe cartilage erosion and meniscal damage in 4 of 6 mice. 4 weeks of WBV followed by 4 weeks of recovery demonstrated cartilage degeneration but to a lesser degree than that detected following 8 weeks of WBV. Whole joints were evaluated using the OARSI scoring system to quantify the degree of joint degeneration, and maximum OARSI scores are presented corresponding to the medial femoral condyle (MFC), medial tibial plateau (MTP), lateral femoral condyle (LFC) and lateral tibial plateau (LTP). Bars indicate mean ± 95% confidence interval. (n = 5-6 mice per group).

Figure 3: Analysis of tibial bone plate microarchitecture following WBV A. Representative coronal micro-computed tomography (micro-CT) images of the proximal tibia in mice exposed to 4 weeks of WBV and sham controls. Subchondral bone plate thickness was measured within a manually defined 3D region of interest (black lines) measuring 500 µm in mediolateral length and 1,150 µm ventrodorsal length, centred on the medial tibial plateau. B. Morphometric analysis of the subchondral bone plate in mice exposed to 4 weeks of WBV shows no significant difference plate thickness, bone mineral content (BMC), or bone mineral density (BMD, mg/cm3) when compared to nonvibrated sham controls. C. Similarly, 8 weeks of WBV or 4 weeks of WBV followed by 4 weeks of recovery (WBV/REC) resulted in no change in subchondral bone plate thickness, bone mineral content, bone mineral density when compared age-matched nonvibrated sham controls. Bars indicate mean ± 95% confidence interval. (n = 5-6 mice per group).

AC C

573

RI PT

544

21

ACCEPTED MANUSCRIPT

587

RI PT

Figure 4: Analysis of subchondral trabecular bone morphometry in mice exposed to WBV. A. Representative coronal micro-computed tomography (micro-CT) images of the proximal tibia in mice exposed to 4 weeks of WBV and sham controls. Bone morphometry was assessed within a manually defined 3D region of interest (black lines) in the medial subchondral trabecular bone to capture trabecular bone within this area, the lowest point of the extensor sulcus was extended as a boundary, and the subchondral bone plate, growth plate and cortical bone were excluded. B. Mice subjected to 4 weeks of WBV demonstrate a significant increase in bone mineral content (BMC), bone mineral density (BMD, mg/cm3), bone volume fraction (BVF), and trabecular thickness (Th) compared to non-vibrated sham controls. C. Mice, subjected to 8 weeks of WBV or 4 weeks of WBV followed by 4 weeks of recovery (WBV/REC), demonstrate no significant difference in any of the reported parameters compared to age-matched sham controls. Bars indicate mean ± 95% confidence interval. (n = 5-6 mice per group).

SC

588 589 590 591 592 593 594 595 596 597 598 599 600

615 616 617 618 619 620 621 622 623 624 625 626 627

TE D

EP

614

Figure 5: Analysis of bone microarchitecture in the tibial primary spongiosa following WBV. A. Representative coronal micro-computed tomography (micro-CT) images of the proximal tibia in mice exposed to 4 weeks of WBV and sham controls. Bone morphometry was assessed within a manually defined 3D region of interest (white lines) to capture the trabecular bone within 50 µm distal to the growth plate and exclude the cortical bone. B. Mice subjected to 4 weeks of WBV demonstrate no significant change in trabecular bone volume fraction (BVF), bone mineral density (BMD, mg/cm3), or bone mineral content (BMC) when compared to age-matched non-vibrated sham controls. C. Similarly, mice exposed to 8 weeks of WBV or 4 weeks of WBV followed by 4 weeks of recovery (WBV/REC) demonstrated no significant difference in any of the reported parameters compared to age-matched non-vibrated sham controls. Bars indicate mean ± 95% confidence interval. (n = 5-6 mice per group).

Figure 6: Analysis of bone microarchitecture in the tibial secondary spongiosa following WBV. A. Representative coronal micro-computed tomography (micro-CT) images of the proximal tibia in mice exposed to 4 weeks of WBV and sham controls. Bone morphometry was assessed within a manually defined 3D region of interest (white lines) to capture the trabecular bone within an area 50 µm distal from the growth plate extending 970 µm distally and exclude the cortical bone. B. Mice subjected to 4 weeks of WBV demonstrate no significant change in trabecular bone volume fraction (BVF), bone mineral density (BMD, mg/cm3), or bone mineral content (BMC) compared to agematched non-vibrated sham controls. C. Similarly, mice exposed to 8 weeks of WBV or 4 weeks of WBV followed by 4 weeks of recovery (WBV/REC) demonstrated no significant difference in any of the reported parameters compared to age-matched nonvibrated sham controls. Bars indicate mean ± 95% confidence interval. (n = 5-6 mice per group).

AC C

602 603 604 605 606 607 608 609 610 611 612 613

M AN U

601

22

AC C

EP

TE D

M AN U

SC

RI PT

ACCEPTED MANUSCRIPT

AC C

EP

TE D

M AN U

SC

RI PT

ACCEPTED MANUSCRIPT

AC C

EP

TE D

M AN U

SC

RI PT

ACCEPTED MANUSCRIPT

AC C

EP

TE D

M AN U

SC

RI PT

ACCEPTED MANUSCRIPT

AC C

EP

TE D

M AN U

SC

RI PT

ACCEPTED MANUSCRIPT

AC C

EP

TE D

M AN U

SC

RI PT

ACCEPTED MANUSCRIPT