Systematic assessment of the growth plates of the wrist in young gymnasts: development and validation of the Amsterdam MRI assessment of the Physis (AMPHYS) protocol.

OBJECTIVES: To develop and validate a protocol for MRI assessment of the distal radial and ulnar periphyseal area in gymnasts and non-gymnasts.
METHODS: Twenty-four gymnasts with wrist pain, 18 asymptomatic gymnasts and 24 non-gymnastic controls (33 girls) underwent MRI of the wrist on a 3T scanner. Sequences included coronal proton density-weighted images with and without fat saturation, and three-dimensional water-selective cartilage scan and T2 Dixon series. Skeletal age was determined using hand radiographs. Three experienced musculoskeletal radiologists established a checklist of possible (peri)physeal abnormalities based on literature and clinical experience. Five other musculoskeletal radiologists and residents evaluated 30 MRI scans (10 from each group) using this checklist and reliability was determined using the intraclass correlation coefficient (ICC) and Fleiss' kappa. A final evaluation protocol was established containing only items with fair to excellent reliability.
RESULTS: Twenty-seven items were assessed for reliability. Intra-rater and inter-rater agreement was good to excellent (respective ICCs 0.60-0.91 and 0.60-0.78) for four epiphyseal bone marrow oedema-related items, physeal signal intensity, metaphyseal junction and depth of metaphyseal intrusions. For physeal thickness, thickness compared with proximal physis of first metacarpal, metaphyseal intrusions, physeal connection of intrusions and metaphyseal bone marrow signal intensity, intra-rater agreement was fair to excellent (ICC/kappa 0.55-0.85) and inter-rater agreement was fair (ICC/kappa 0.41-0.59). Twelve items were included in the final protocol.
CONCLUSION: The Amsterdam MRI assessment of the Physis protocol facilitates patient-friendly and reliable assessment of the (peri)physeal area in the radius and ulna.

The concise Amsterdam MRI assessment of the Physis (AMPHYS) protocol contains 12 items derived from the literature and clinical experience that can be used to assess characteristics of the periphyseal area in the distal radius and ulna.

All 12 items of the AMPHYS evaluation protocol showed good content validity and fair to excellent intra-rater and inter-rater agreement.

The accompanying 15 min MRI protocol includes coronal three-dimensional frequency-selective, fat-suppressed gradient-echo imaging, a T2-weighted two-point Dixon sequence and proton density-weighted sequences with and without fat suppression.

The AMPHYS protocol contributes to uniform assessment and reporting of suspected physeal stress injury on MRI, with minimal patient and clinician burden.

Use of the AMPHYS protocol in research and clinical practice provides more insight in morphological differences between stress-induced and maturity-related changes in the periphyseal area.

Introduction

In young athletes, physeal injury can occur as traumatic fractures or as stress injuries caused by repetitive microtrauma. The latter are commonly located at physes of the distal radius (gymnast wrist), the distal humerus (Little League elbow) and the proximal humerus (Little League shoulder). In the physis, or growth plate, cartilage cells are generated, proliferate, hypertrophy, and eventually calcify into bone, vascularised by metaphyseal and epiphyseal vessels. Injury to the multilayered physis or its vascularisation can cause (partial) physeal growth arrest, sometimes resulting in permanent growth disturbances and damage to surrounding structures.

Wrist pain and overuse wrist injury occur frequently in young gymnasts, and early diagnosis of physeal stress injury allows timely intervention and recovery to prevent long-term sequelae. Radiographic signs include distal radial physeal irregularity and widening. These characteristics are incorporated in a radiographic grading system, irregularity representing grade 1 injury and widening severe (grade 3) injury. MRI can, in addition, depict non-osseous tissues like physeal cartilage, the highly vascular primary metaphyseal spongiosa and traumatic signs like bone marrow oedema (BMO). Cartilage-sensitive and fluid-sensitive MRI sequences such as three-dimensional (3D) frequency-selective, fat-suppressed gradient-echo images and fat-saturated T2-weighted images have been recommended for detailed imaging of physeal cartilage and stress-induced BMO. Dixon chemical shift MRI can be used to achieve uniform fat suppression in the hand and wrist.

MRI has been proposed as the imaging method of choice for evaluation and therapeutic decision-making in patients with physeal stress injury. However, some injury manifestations may resemble normal growth and physeal development on MRI, and BMO is often present in wrists of asymptomatic children. This renders uniform assessment and diagnosis of physeal stress injury challenging, especially at an early stage. In addition, the number of potential injury signs described on radiographs, and translated to MRI, is extensive and may require lengthy MRI examinationsand evaluations. A concise and reliable procedure for MRI-based assessment of the physis can aid in identifying physeal injury in a patient-friendly manner. This study aimed to develop and validate a standardised protocol for MRI evaluation of the distal radial and ulnar periphyseal area to improve uniform assessment of the wrist growth plates.

Methods

Study design

This study was performed according to the Declaration of Helsinki and approved by our institution’s review board (reference no 2014_382#B2015303). It consisted of multiple phases: prospective MRI collection, literature-based and expert-based protocol development, and reliability testing by different experts. MRI acquisition took place at the Academic Medical Center, Amsterdam, between June 2015 and November 2017. Participants and their parent(s) or guardian(s) provided written informed consent to participate.

Study population

The study cohort consisted of 24 gymnasts with wrist pain, 18 asymptomatic gymnasts and 24 non-gymnasts, aged 12–18 years. Symptomatic gymnasts were referred by their (sports) physician. Sex-matched and skeletal age-matched asymptomatic gymnasts and non-gymnasts were recruited through gymnastics clubs, the bring-a-friend strategy and notices within our institution. Gymnasts had performed their sport for at least 1 year and up to 6 months or less prior to participation. The symptomatic group consisted of gymnasts with clinically suspected physeal injury and wrist pain in the 6 months prior to inclusion. The non-gymnast group consisted of healthy children without wrist pain. Exclusion criteria were the following: history of fracture, wrist surgery or infection, growth disorders, systemic or oncological disease involving the musculoskeletal system (eg, juvenile idiopathic arthritis) and closed growth plate on hand radiograph. Participants filled out a questionnaire on demographic information, sports participation and wrist pain.

Imaging

Standard posterior-anterior radiographs with focus–detector distance of 1.30 m were obtained of one (symptomatic) hand and wrist. Skeletal age was determined using a computerised method (BoneXpert, v1.0; Visiana, Holte, Denmark; www.BoneXpert.com) validated in a healthy Dutch population. MRI of the (symptomatic) wrist was performed in a feet-first, supine position with both arms resting alongside the body on a 3.0 Tesla MRI scanner (Ingenia; Philips Healthcare, Best, The Netherlands) using a dedicated eight-channel, receive-only wrist coil.

The MRI protocol included coronal turbo spin-echo (TSE) proton density (PD)-weighted sequences with and without fat suppression (spectral attenuated inversion recovery) as part of our institution’s standard clinical protocol. Two additional sequences were used: a coronal T2-weighted two-point Dixon sequence and a coronal 3D water-selective cartilage scan (3D WATSc) (online supplementary appendix 1). The field of view was centred on the distal radial and ulnar physes and included the proximal physis of the first metacarpal bone (MC-1) on all images (online supplementary appendix 1).

Development of Amsterdam MRI assessment of the Physis (AMPHYS) protocol

Following guidelines for achieving good content validity, potentially relevant growth plate characteristics and physeal (stress) injury signs were collected from the literature (figure 1). The expert group in the development phase consisted of three experienced musculoskeletal radiologists (MCdJ, EHGO, MM) from different institutions and two physicians experienced in research on musculoskeletal imaging (LSK, RBJK). The senior author added items to the literature-derived list, based on the experience in evaluating physeal stress injuries on standard clinical MRI. A standardised scoring form with these items was created, allowing radiologists to indicate a finding’s presence in any image series, its extent or severity, its location in the radius and/or ulna, and which sequences provided optimal assessability.

Figure 1

Flow chart showing the development process of the Amsterdam MRI assessment of the Physis evaluation protocol.

After the introduction of the sequences comprising the MRI protocol, the radiologists individually evaluated a random sample of blinded MRI scans from each of the participant groups. Image evaluation was performed on a PC workstation with high-resolution monitor using IMPAX software version 6.6.1.4024 (AGFA HealthCare N.V., Mortsel, Belgium). One radiologist (MM) assessed scans of 60 participants, the others (MCdJ, EHGO) of at least 15 participants. The expert group evaluated the reading directly afterwards, discussing disagreements in scoring until consensus was reached. Items considered not relevant or not assessable on MRI were removed. An initial MRI scoring checklist was formed, including instructions and images illustrating abnormalities.

Validation of AMPHYS protocol

The expert group in the validation phase was not involved in protocol development and consisted of three musculoskeletal radiologists and two musculoskeletal radiology residents (KFvD, FFS, MPT, SJ, RH) from three institutions, representing daily radiological practice. To assure consensus on interpretation of scoring instructions, the group discussed the newly developed MRI scoring checklist and examples, and assessed multiple trial cases prior to image evaluation. Using the same workstations and software as the development group, all raters evaluated 30 blinded MRI scans, that is, 10 from each participant group (figure 1). Four weeks after the first session, one rater (RH) re-evaluated all 30 MRI scans to determine intra-rater agreement.

Statistical analysis

Inter-rater agreement of scoring items was determined by calculating intraclass correlation coefficients (ICCs) for absolute agreement between raters using a two-way random analysis of variance (ANOVA) model (case 2, ICC(2,1)) for ordinal variables, and the unweighted Fleiss’ kappa for binary variables. The set of 30 double assessments by one rater was used to calculate the intra-rater ICC for absolute agreement for each item, using a two-way random ANOVA model (case 2, ICC(2,1)) for ordinal variables, and the unweighted Fleiss’ kappa for dichotomous variables. Agreement was expressed by the ICC or kappa and its CI, and levels of agreement measured by ICC were defined according to Cicchetti (<0.4, poor; 0.4 to 0.59, fair; 0.6 to 0.74, good; ≥0.75, excellent), and for kappa, according to Landis and Koch (<0, poor; 0 to 0.2, slight; 0.21 to 0.4, fair; 0.41 to 0.6, moderate; 0.61 to 0.8, substantial; 0.81 to 1.0, almost perfect). Only items with fair to excellent inter-rater and intra-rater agreement were considered to show sufficient agreement and variability within the population and included in the final protocol. Based on an expected ICC ≥0.8 and preferred 95% CI of 0.6 to 1.0, the preferred sample size was eight images per group to be rated by each rater. Data analyses were performed using SPSS V24.0.

Results

Participant characteristics

The cohort consisted of 33 boys and 33 girls. MRI scans for inter-rater and intra-rater agreement were of five girls and five boys in the symptomatic gymnast group, six girls and four boys in the asymptomatic gymnast group and seven girls and three boys in the non-gymnast group. Respective mean calendar and skeletal ages were 14.7 (SD 1.3) and 13.0 (SD 0.7) years in symptomatic gymnasts, 14.1 (SD 1.0) and 12.2 (SD 1.0) years in asymptomatic gymnasts and 13.5 (SD 1.1) and 13.0 (SD 1.8) years in non-gymnasts.

Content validity

The development phase rendered 16 relevant items from the literature and 14 items from clinical experience. After the first reading and consensus meeting, three items were excluded (physeal depressions, metaphyseal cysts and striations). The remaining 27 items were considered relevant and dichotomous or ordinal variables assessable on MRI. Items were divided into three categories: epiphysis (ie, various characteristics of BMO), physis (ie, thickness, signal intensity, disruptions and epiphyseal border) and metaphysis (ie, physeal border, intrusions, BMO, widening, periosteal bone formation and sclerosis) (table 1). Preferred sequences were defined per item and visibility of BMO on the cartilage-specific 3D WATSc sequence was identified separately to potentially identify very severe BMO. Figure 2 shows the normal multilayered aspect of the physis, and figure 3 shows examples of epiphyseal BMO, physeal thickening, metaphyseal intrusions and metaphyseal BMO.

Table 1

Scoring items and grades for physeal characteristics in the distal radius and ulna

Scoring item	Grades
Epiphysis
Bone marrow oedema
Extent*	No oedema	<50% of epiphyseal volume					>50% of epiphyseal volume
Location10 32	No oedema	Oedema adjacent to physis					Oedema not adjacent to physis
Signal intensity*	No oedema	1	2	3	4	5
Visibility on 3D WATSc*	No oedema	Oedema not visible					Oedema visible
Physis
Thickness9 20	Normal	Increased
Location of thickness20	No increased thickness	Increased on radial side					Increased on ulnar side
Thickness compared with proximal physis of MC-1*	Not increased	Twice		Three times			Four times
Signal intensity8	1	2	3	4	5
Disruptions32 33	No disruptions	1	2	>2 disruptions
Width of disruptions*	No disruptions	<2 mm		≥2 mm
Physeal border on epiphyseal side25	Undulating	Irregular
Depth of irregularities*	No irregularities	<Thickness of physis					>Thickness of physis
Metaphysis
Physeal border on metaphyseal side25	Undulating	Slightly irregular					Distinctly irregular
Metaphyseal intrusions32 34	Absent	Present
Signal intensity intrusions*	No intrusions	Less than physis		Same as physis			Higher than physis
Connection of intrusion with physis*	No intrusions	Connected with physis					Not connected with physis
Depth of intrusions*	No intrusions	<2 mm					>2 mm
Bone marrow oedema
Presence10 32	Present	Absent
Depth*	No oedema	Area <2 cm from physis					Area ≥2 cm from physis
Location10 32	No oedema	Oedema adjacent to physis					Oedema not adjacent to physis
Signal intensity*	No oedema	1	2	3	4	5
Visibility on 3D WATSc*	No oedema	Oedema not visible					Oedema visible
Homogeneity*	No oedema	Homogeneous oedema					Inhomogeneous oedema
Shape25	Normal	Unilateral widening					Bilateral widening
Location of widening*	No widening	Radial					Ulnar
Periosteal bone formation10	Absent	Present
Sclerosis2 9	Absent	Present

*Item added based on clinical experience with evaluation of physeal stress injuries.

3D WATSc, three-dimensional water-selective cartilage scan; MC-1, first metacarpal bone.

Figure 2

Left: three-dimensional water-selective cartilage scan image showing the trilaminar appearance of the physis. Right: schematic overview of the three layers of the physis. A: cartilaginous part, B: zone of provisional calcification, C: primary spongiosa of metaphysis.

Figure 3

(A) T2 Dixon image showing a focal patch of epiphyseal bone marrow oedema in the radius, indicated by a white arrowhead. (B) Three-dimensional water-selective cartilage scan (3D WATSc) image showing diffuse thickening of the distal radial physis (located to the right of the white asterisk) compared with the proximal physis of the first metacarpal bone (located to the right of the black asterisk). (C) 3D WATSc image showing intrusions of physeal cartilage into the radial metaphysis, marked by white arrowheads. (D) T2 Dixon image showing extensive bone marrow oedema of the radial metaphysis, marked with a black asterisk.

Inter-rater agreement

In the radius, inter-rater agreement was excellent for BMO extent and signal intensity (ICC, 0.78 and 0.76, respectively) (table 2). For BMO location and visibility in the epiphysis on the 3D WATSc sequence, inter-rater agreement was good (ICC 0.60 and 0.70, respectively), as well as for physeal signal intensity (ICC 0.62), metaphyseal border of physis (ICC 0.60) and depth of intrusions into the metaphysis (ICC 0.69). Agreement was moderate for metaphyseal intrusions (kappa 0.59) and physeal thickness (kappa 0.47), and fair for physeal thickness compared with the proximal MC-1 physis (ICC 0.58), physeal connection of metaphyseal intrusions (ICC 0.41) presence (kappa 0.32) and signal intensity of metaphyseal BMO (ICC 0.51), physeal border on epiphyseal side (kappa 0.35) and periosteal bone formation (kappa 0.26).

Table 2

Intra-rater and inter-rater agreement for items relating to the radius and ulna

Categories and items	Radius				Ulna
	Intra-rater agreement		Inter-rater agreement		Intra-rater agreement		Inter-rater agreement
	ICC	95% CI	ICC	95% CI	ICC	95% CI	ICC	95% CI
Epiphysis
Bone marrow oedema
Extent	0.86	0.72 to 0.94	0.78	0.66 to 0.88	0.74	0.53 to 0.87	0.30	0.15 to 0.50
Location	0.60	0.32 to 0.79	0.60	0.43 to 0.75	0.69	0.44 to 0.84	0.26	0.11 to 0.45
Signal intensity	0.85	0.69 to 0.93	0.76	0.63 to 0.86	0.59	0.30 to 0.78	0.27	0.12 to 0.47
Visibility on 3D WATSc	0.90	0.79 to 0.95	0.70	0.55 to 0.83	0.28	0 to 0.57	0.34	0.18 to 0.54
Physis
Thickness*	0.80	0.62 to 0.90	0.47	0.35 to 0.59	0.44	0.09 to 0.69	0.22	0.10 to 0.34
Location of thickness	0	0 to 0.45	0.07	0 to 0.24	0	NA	0	0 to 0.14
Thickness compared with proximal physis of MC-1	0.83	0.67 to 0.91	0.58	0.42 to 0.73	0.52	0.19 to 0.74	0.21	0.07 to 0.40
Signal intensity	0.81	0.65 to 0.90	0.62	0.47 to 0.76	0.62	0.35 to 0.80	0.39	0.23 to 0.58
Disruptions	0.36	0 to 0.64	0.23	0.08 to 0.43	0.23	0 to 0.52	0.25	0.10 to 0.45
Width of disruptions	0.39	0.02 to 0.66	0.26	0.11 to 0.45	0.42	0.09 to 0.67	0.26	0.11 to 0.46
Physeal border on epiphyseal side*	0.67	0.40 to 0.83	0.35	0.23 to 0.47	0.67	0.42 to 0.83	0.20	0.08 to 0.32
Depth of irregularities	0	0 to 0.36	0	0 to 0.15	0	0 to 0.36		0 to 0.13
Metaphysis
Physeal border on metaphyseal side	0.84	0.69 to 0.92	0.60	0.43 to 0.75	0.58	0.28 to 0.78	0.43	0.27 to 0.61
Metaphyseal intrusions*	0.85	0.71 to 0.93	0.59	0.47 to 0.71	0.71	0.46 to 0.86	0.57	0.44 to 0.69
Signal intensity intrusions	0.75	0.53 to 0.87	0	0 to 0.15	0.73	0.46 to 0.87	0	0 to 0.15
Connection of intrusion with physis	0.85	0.71 to 0.93	0.41	0.24 to 0.60	0.56	0.25 to 0.77	0.43	0.26 to 0.62
Depth of intrusions	0.91	0.82 to 0.96	0.69	0.54 to 0.82	0.67	0.36 to 0.84	0.54	0.38 to 0.71
Bone marrow oedema
Presence*	0	0 to 0.36	0.32	0.21 to 0.44	0	NA	0.33	0.21 to 0.45
Depth	0	0 to 0.36	0.36	0.20 to 0.56	0	NA	0.28	0.12 to 0.49
Location	0	NA	0.34	0.18 to 0.53	0	NA	0.34	0.18 to 0.54
Signal intensity	0.55	0.22 to 0.76	0.51	0.35 to 0.69	0.46	0.10 to 0.71	0.43	0.26 to 0.61
Visibility on 3D WATSc	0.43	0.11 to 0.68	0.31	0.16 to 0.51	0.23	0 to 0.52	0.34	0.18 to 0.53
Homogeneity	0	0 to 0.04	0.34	0.14 to 0.58	0	0 to 0.02	0.15	0.03 to 0.34
Shape	0	0 to 0.36	0.33	0.17 to 0.52	0	0 to 0.34	0.12	0.01 to 0.28
Location of widening	0	NA	0.10	0 to 0.28	0	0 to 0.35	0	0 to 0.14
Periosteal bone formation*	0	NA	0.26	0.14 to 0.39	0	NA	0	0 to 0.09
Sclerosis*	0.23	0 to 0.54	0.08	0 to 0.22	0	NA	0	0 to 0.04

*Fleiss’ kappa is reported for dichotomous variables.

3D WATSc, three-dimensional water-selective cartilage scan; ICC, intraclass correlation coefficient; MC-1, first metacarpal bone; NA, not available (due to zero variance).

In the ulna, inter-rater agreement was moderate for metaphyseal intrusions (kappa 0.57). Agreement was fair for metaphyseal border of the physis (ICC 0.43), connection of metaphyseal intrusions with physis (ICC 0.43), intrusion depth (ICC 0.54), presence (kappa 0.33) and signal intensity of metaphyseal BMO (ICC 0.43) and physeal thickness (kappa 0.22). Agreement was poor or slight for the other items.

Intra-rater agreement

For the radius, intra-rater agreement was excellent or substantial for extent, signal intensity and visibility of epiphyseal BMO (ICC 0.86, 0.85 and 0.90, respectively), as well as for physeal thickness (kappa 0.80), thickness compared with the proximal MC-1 physis (ICC 0.83), signal intensity (ICC 0.81), metaphyseal border (ICC 0.84), metaphyseal intrusions (kappa, 0.85) and their physeal connection (ICC 0.85) and depth (ICC 0.91) (table 2). Agreement was good for epiphyseal BMO location (ICC 0.60) and epiphyseal border of physis (ICC 0.67), and fair for metaphyseal BMO signal intensity (ICC 0.55) and visibility on 3D WATSc (ICC 0.43).

Intra-rater agreement for ulnar items was good or substantial for extent (ICC 0.74) and location (ICC 0.69) of epiphyseal BMO, physeal signal intensity (ICC 0.62), and metaphyseal intrusion presence (kappa 0.71) and depth (ICC 0.67). Agreement was moderate or fair for epiphyseal BMO signal intensity (ICC 0.59), metaphyseal BMO (ICC 0.46), physeal thickness (kappa 0.44), thickness compared with the proximal MC-1 physis (ICC 0.52), metaphyseal border (ICC 0.58) and intrusion connection with the physis (ICC 0.56). The remaining items showed poor or slight agreement.

The final AMPHYS protocol consisted of 12 items for the radius and 5 for the ulna (table 3, online supplementary appendix 2) and a scoring form (online supplementary appendix 3).

Table 3

Amsterdam MRI assessment of the Physis protocol for (peri)physeal characteristics of the radius and ulna

Scoring item	MRI sequence	Radius	Ulna	Grade 0	Grade 1	Grade 2	Grade 3	Grade 4	Grade 5
Epiphysis
Bone marrow oedema
1.1. Extent	T2 Dixon, PD TSE SPAIR	+		No oedema	<50% of epiphyseal volume	>50% of epiphyseal volume
1.2. Location	T2 Dixon, PD TSE SPAIR	+		No oedema	Oedema adjacent to physis	Oedema not adjacent to physis
1.3. Signal intensity	T2 Dixon, PD TSE SPAIR	+		No oedema	1 (minimal signal intensity)	2	3	4	5
1.4. Visibility on 3D WATSc	3D WATSc	+		No oedema	Oedema not visible		Oedema visible
Physis
2. Thickness	3D WATSc, T2 Dixon, PD TSE			Normal	Increased
2.1. Compared with MC-1 proximal physis	3D WATSc, T2 Dixon, PD TSE	+		Not increased	Twice as thick as MC-1	Three times as thick as MC-1	Four times as thick as MC-1
3. Signal intensity	3D WATSc, T2 Dixon, PD TSE	+		Normal signal intensity	1 (minimal signal intensity)	2	3	4
Metaphysis
4. Physeal border on metaphyseal side	3D WATSc	+	+	Undulating	Slightly irregular	Distinctly irregular
5. Metaphyseal intrusions	3D WATSc	+	+	Absent	Present
5.1. Connection of intrusion with physis	3D WATSc	+	+	No intrusions	Connected with physis	Not connected with physis
5.2. Depth of intrusions	3D WATSc	+	+	No intrusions	<2 mm	>2 mm
Bone marrow oedema
5.3. Signal intensity	T2 Dixon, PD TSE, 3D WATSc, PD TSE SPAIR	+	+	No oedema	1 (minimal signal intensity)	2	3	4	5

3D WATSc, three-dimensional water-selective cartilage scan; MC-1, first metacarpal bone; PD, proton density; SPAIR, spectral attenuated inversion recovery; TSE, turbo spin-echo.

+

+

+

+

+

+

+

+

+

+

+

+

+

+

+

+

Discussion

The AMPHYS protocol contains 12 elements with good content validity to assess characteristics of the physis, epiphysis and metaphysis of the radius and ulna, with fair to excellent inter-rater and intra-rater agreement.

Bone marrow oedema

While BMO can indicate injury, it can be present in wrist bones of 40%–49% of asymptomatic children, especially during rapid skeletal maturation. Local areas of focal periphyseal oedema (FOPE) have been described as early sign of physiological physeal fusion. In young athletes, asymptomatic BMO may result from a physiological stress response to exercise. We therefore expected that periphyseal BMO would be present in all three groups, but its characteristics might be used to differentiate between physiological and injury-related BMO.

While five BMO-related items in the epiphysis and metaphysis showed fair to excellent intra-rater and inter-rater agreement, reliability was best for multiple characteristics of epiphyseal BMO (figure 3A). Periphyseal BMO was present in all groups, but its signal intensity and epiphyseal extent and location showed fair to excellent agreement for the ulna and the radius. Gradient-echo sequences commonly used for (physeal) cartilage imaging, such as 3D WATSc, are well-known to be insensitive for BMO, compared with fat-suppressed, T2-weighted fast spin-echo sequences. However, our results show that gradient-echo images may be useful in identifying severe cases of periphyseal BMO.

Physeal changes

Increased physeal thickness is a frequently described sign of stress injury on conventional radiographs and MRI. During normal growth, thickness remains constant almost until maturity, when the physis thins and eventually disappears. Physeal widening is thought to be caused by disrupted enchondral calcification, leading to accumulation of hypertrophied cartilage cells failing to ossify. This pathological process can occur focally with unilateral widening or affect the entire physis (figure 3B). We propose a method to assess distal radial and ulnar physeal thickness by comparison with the proximal MC-1 physis. Since this bone suffers less repetitive axial loading than the radius and ulna in gymnasts, we used this as within-person ‘reference physis’ that can easily be included in the field of view.

Other studies have also described bony bridges indicating physeal closure, sometimes prematurely, after injury. 3D volumetric reconstruction of the physis and bony bars was proposed to aid in treatment decision-making. We found poor reliability for physeal disruptions on MRI, suggesting these characteristics are less suitable for physeal assessment. Interpretation possibly depends on disruption size: small physeal disruptions that were seen in all three participant groups were not graded as such by all observers. These may be signs of normal physeal maturation, or susceptibility artefacts associated with gradient-echo MRI, likely caused by calcium depositions during local growth plate closure.

Physeal haziness or decreased radiolucency, another literature-derived characteristic of physeal stress injury, was excluded because it is mainly assessed on radiographs. Its MRI equivalent may be increased physeal signal intensity, reported after indirect physeal trauma. Although most of this study’s MRI sequences depict the physis as a relatively high signal structure, physeal signal intensity showed fair reliability and can be indicative of injury.

Metaphyseal changes

The physeal appearance on MRI is trilaminar: a hyperintense cartilaginous layer, a hypointense zone of provisional calcification and a hyperintense region of primary spongiosa or metaphyseal vascularisation (figure 2). The healthy physis can appear undulated. Irregularity of physeal borders is described on radiographs and MRI of gymnasts with stress injury of the distal radial physis. Our results show better inter-rater agreement for irregularity of the metaphyseal border compared with the epiphyseal side. Discontinuations of the metaphyseal provisional calcification zone are known injury signs. Most likely, the physeal–metaphyseal junction is more frequently affected by stress injury because of its role in chondrocyte calcification during longitudinal growth (figure 2).

Similarly, intrusions or physeal cartilage ‘tongues’ with signal intensities similar to the physis can extend into the metaphysis after physeal injury (figure 3C). Reliability was moderate to almost perfect for the presence of metaphyseal intrusions in both the radius and the ulna. Physeal connection and depth of these intrusions showed fair to good inter-rater agreement, while intra-rater agreement was excellent. Thus, interpretation of a high signal area’s physeal connection varies largely between observers and may in some cases include focal patches of high signal intensity also interpretable as metaphyseal BMO or FOPE. Intrusion extent or depth may therefore be more reliable injury signs.

Other metaphyseal changes associated with physeal stress injury include widening, cystic changes, sclerosis and striations on radiographs, and low signal and periosteal bone formation on MRI. These were excluded from the AMPHYS protocol because of poor reliability, likely caused by low prevalence in all groups.

Strengths and limitations

The concise AMPHYS protocol provides reliable physeal assessment while minimising scan time, and accompanying patient burden. We aimed to achieve good validity by involving a varied group of five observers for reliability assessment, and a heterogeneous sample of symptomatic and asymptomatic gymnasts and non-gymnasts to ensure sufficient population variety. However, some variability in abnormalities between groups may be caused by artefacts interpreted as (peri)physeal changes, like susceptibility artefacts on 3D WATSc and chemical shift artefacts on spin-echo PD images. In addition, the MC-1 physis is less reliable as reference for thickness when it has nearly fused. Nevertheless, assessment of radial physeal thickness by itself also showed moderate to substantial inter-rater and intra-rater agreement. Finally, sample size was based on the aim to assess the protocol’s reliability. For assessment of its diagnostic accuracy and score interpretation, a separate study in a larger sample is necessary.

Clinical implications and future directions

The patient-friendly and radiologist-friendly AMPHYS protocol is directly available for standardised and quick assessment of periphyseal changes in children with suspected physeal injury. MRI of the wrist on 3T is currently the standard of care in many clinical settings, and the prevalent sequences on which the evaluation protocol is based are supplied by most vendors and can likely even be modified for 1.5T scanners. In addition, patient burden is minimal, with absence of ionising radiation and a scan time of less than 15 min. The protocol can be used for initial injury assessment, treatment and recovery evaluation, and monitoring of (peri)physeal changes in children at risk of physeal injury. Uniform reporting of physeal stress changes on MRI can contribute to patient care and further research on related topics such as prognosis of injury recovery and potential complications.

This study provides data from gymnasts as the patient group most frequently affected with physeal stress injury of the wrist. Comparison with asymptomatic gymnasts is recommended because of changes that can be present due to physiological responses to exercise. Outcome scores need to be validated on a larger scale to provide diagnostic accuracy, grading interpretation and cut-off values for presence and severity of physeal stress injury. We will proceed to evaluate this grading system in a larger cohort study.

Conclusion

The AMPHYS protocol is a concise collection of radiographic and MRI-based characteristics of the periphyseal area of the radius and ulna that can be reliably assessed on 3T MRI after merely 15 minutes of scan time. Its 12 items include epiphyseal and metaphyseal BMO, physeal thickness and signal intensity, and metaphyseal intrusions and irregularities, with fair to excellent inter-rater and intra-rater agreement.