Europe PMC

This website requires cookies, and the limited processing of your personal data in order to function. By using the site you are agreeing to this as outlined in our privacy notice and cookie policy.

Abstract 


Background

Legg-Calvé-Perthes disease (LCPD) is a childhood hip disorder thought to be caused by disruption of blood supply to the developing femoral head. There is potential for imaging to help assess revascularization of the femoral head.

Purpose

To investigate whether quantitative susceptibility mapping (QSM) can detect neovascularization in the epiphyseal cartilage following ischemic injury to the developing femoral head.

Study type

Prospective.

Animal model

Right femoral head ischemia was surgically induced in 6-week-old male piglets. The animals were sacrificed 48 hours (n = 3) or 4 weeks (n = 7) following surgery, and the operated and contralateral control femoral heads were harvested for ex vivo MRI.

Field strength/sequence

Preclinical 9.4T MRI to acquire susceptibility-weighted 3D gradient echo (GRE) images with 0.1 mm isotropic spatial resolution.

Assessment

The 3D GRE images were used to manually segment the cartilage overlying the femoral head and were subsequently postprocessed using QSM. Vessel volume, cartilage volume, and vessel density were measured and compared between operated and control femoral heads at each timepoint. Maximum intensity projections of the QSM images were subjectively assessed to identity differences in cartilage canal appearance, location, and density.

Statistical tests

Paired t-tests with Bonferroni correction were used (P < 0.008 considered significant).

Results

Increased vascularity of the epiphyseal cartilage following ischemic injury was clearly identified using QSM. No changes were detected 48 hours after surgery. Vessel volume, cartilage volume, and vessel density were all increased in the operated vs. control femoral heads 4 weeks after surgery (P = 0.001, 0.002, and 0.001, respectively). Qualitatively, the increase in vessel density at 4 weeks was due to the formation of new vessels that were organized in a brush-like orientation in the epiphyseal cartilage, consistent with the histological appearance of neovascularization.

Data conclusion

QSM can detect neovascularization in the epiphyseal cartilage following ischemic injury to the femoral head.

Level of evidence

1 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2019;50:106-113.

Free full text 


Logo of nihpaLink to Publisher's site
J Magn Reson Imaging. Author manuscript; available in PMC 2020 Jul 1.
Published in final edited form as:
PMCID: PMC7249674
NIHMSID: NIHMS1580176
PMID: 30556613

Quantitative Susceptibility Mapping Detects Neovascularization of the Epiphyseal Cartilage After Ischemic Injury in a Piglet Model of Legg-Calvé-Perthes Disease

Casey P. Johnson, PhD,1,2,* Luning Wang, PhD,1,2 Ferenc Tóth, DVM, PhD,3 Olumide Aruwajoye, PhD,4 Brooke Kirkham, BA,1 Cathy S. Carlson, DVM, PhD,5 Harry K.W. Kim, MD, MSc,4,6 and Jutta M. Ellermann, MD, PhD1,2

Abstract

Background:

Legg-Calvé-Perthes disease (LCPD) is a childhood hip disorder thought to be caused by disruption of blood supply to the developing femoral head. There is potential for imaging to help assess revascularization of the femoral head.

Purpose:

To investigate whether quantitative susceptibility mapping (QSM) can detect neovascularization in the epiphyseal cartilage following ischemic injury to the developing femoral head.

Study Type:

Prospective.

Animal Model:

Right femoral head ischemia was surgically induced in 6-week-old male piglets. The animals were sacrificed 48 hours (n = 3) or 4 weeks (n = 7) following surgery, and the operated and contralateral control femoral heads were harvested for ex vivo MRI.

Field Strength/Sequence:

Preclinical 9.4T MRI to acquire susceptibility-weighted 3D gradient echo (GRE) images with 0.1 mm isotropic spatial resolution.

Assessment:

The 3D GRE images were used to manually segment the cartilage overlying the femoral head and were sub-sequently postprocessed using QSM. Vessel volume, cartilage volume, and vessel density were measured and compared between operated and control femoral heads at each timepoint. Maximum intensity projections of the QSM images were subjectively assessed to identity differences in cartilage canal appearance, location, and density.

Statistical Tests:

Paired t-tests with Bonferroni correction were used (P < 0.008 considered significant).

Results:

Increased vascularity of the epiphyseal cartilage following ischemic injury was clearly identified using QSM. No changes were detected 48 hours after surgery. Vessel volume, cartilage volume, and vessel density were all increased in the operated vs. control femoral heads 4 weeks after surgery (P = 0.001, 0.002, and 0.001, respectively). Qualitatively, the increase in vessel density at 4 weeks was due to the formation of new vessels that were organized in a brush-like orientation in the epiphyseal cartilage, consistent with the histological appearance of neovascularization.

Data Conclusion:

QSM can detect neovascularization in the epiphyseal cartilage following ischemic injury to the femoral head.

LEGG-CALVÉ-PERTHES DISEASE (LCPD) is a serious childhood hip disorder that can result in femoral head deformity and early-onset osteoarthritis.1 LCPD is caused by the interruption of blood supply to the developing femoral head, which can lead to extensive ischemic injury to the secondary ossification center (SOC) and the overlying subarticular epiphyseal cartilage (ie, growth cartilage).2 Over time, the affected femoral head spontaneously revascularizes, which initiates reparative and resorptive processes that weaken the femoral head and make it susceptible to deformation and collapse.2 The timing and rate of revascularization are important factors in determining the progression of femoral head healing and to predict clinical outcome. However, there is a lack of imaging technologies to noninvasively and serially assess revascularization patterns in both the SOC and the epiphyseal cartilage following ischemic injury.3

The vascular supply to the developing femoral head is provided by blood vessels passing through cartilage canals in the epiphyseal cartilage.3,4 Unlike articular cartilage, which is an avascular tissue, epiphyseal cartilage is only present during childhood and is extensively vascularized via a network of cartilage canal vessels.5 Epiphyseal cartilage serves as a growth cartilage responsible for the volumetric growth of the SOC. As the femoral head grows, the epiphyseal cartilage is gradually replaced by the advancing SOC until the calcified layer of the cartilage is formed and only the articular cartilage remains.5 In LCPD, there is a growth arrest of the SOC due to the cessation of blood flow and necrosis of the deep layer of the epiphyseal cartilage.1 It has been shown in a piglet model of LCPD, using micro-computed tomography (CT) angiography (μCTA), that revascularization of the femoral head following ischemic injury is initiated by invasion of new vessels (neovascularization) from the femoral neck through the epiphyseal cartilage.6 These vessels express high levels of hypoxia-inducible factor-1 and vascular endothelial growth factor.6 The pattern of neovascularization observed in the piglet model is consistent with the findings from a serial perfusion study in patients with LCPD using contrast-enhanced magnetic resonance imaging (MRI), which demonstrated that revascularization in the SOC is first observed at the periphery of the femoral head.3 Interestingly, younger children with LCPD tend to have better clinical outcomes than older children,7,8 which may be related to their greater epiphyseal cartilage volume and, thus, greater capacity to repair and remodel the femoral head following ischemic injury.6 Clinically, it would be useful to image the vascularity of the epiphyseal cartilage in LCPD to better understand the recovery of the epiphyseal cartilage from ischemic injury, to identify restoration of epiphyseal growth, and to assess the efficacy of new angiogenic treatments. Currently, contrast-enhanced MRI is the only available in vivo imaging technique to visualize vascularity in the developing femoral head.9,10 However, recent findings of gadolinium deposition in the brain has renewed concerns regarding its potential long-term effects in the pediatric patient population.11

Recently, an MRI technique based on susceptibility-weighted imaging (SWI) was introduced12 and further developed1319 to image vascularity in the epiphyseal cartilage without the administration of exogenous contrast agents. SWI has been applied to visualize epiphyseal cartilage vascularity in the distal femur ex vivo1217 and, more recently, in vivo in children,18,19 demonstrating its potential for clinical use. SWI is sensitive to magnetic susceptibility differences between the cartilage canals and the surrounding cartilage matrix.12 SWI images can be further processed using quantitative susceptibility mapping (QSM) to improve visualization of the cartilage canals and enable quantification.14,17 QSM of cartilage canals has been extensively validated histologically and provides detailed visualization of the 3D vascular architecture of epiphyseal cartilage in the distal femur.13,16,17

The purpose of this study was to investigate whether QSM can detect neovascularization following ischemic injury to the developing femoral head in a piglet model of LCPD. We hypothesized that neovascularization in the epiphyseal cartilage can be detected in the ischemic femoral head during the early repair phase.

Materials and Methods

Animals

Our prospective study was conducted from November 2015 to February 2018 in accordance with approved Institutional Animal Care and Use Committee protocols. Eleven 6-week-old male piglets were purchased from one of two commercial providers (Change of Pace in Aubrey, TX; or Manthei Hog Farm in Elk River, MN) and underwent surgery to induce whole right femoral head ischemia. Male piglets were chosen because LCPD is five times more common in boys than girls,20 and the size of the femoral head of a 6-week-old piglet is similar to that of a 4- to 6-year-old child, which is a typical age of LCPD onset.21 The animal model is well established and involves placing a tight ligature about the femoral neck and transecting the ligamentum teres to interrupt all blood flow to the developing femoral head.22,23 The surgeries were performed by either a board-certified pediatric orthopedic surgeon (H.K.W.K.) or a board-certified large animal orthopedic surgeon (F.T.) with respectively 20 years and 1 year of experience using the model. The piglets were then euthanized either 48 hours (n = 3) or 4 weeks (n = 8) after surgery, and the operated (right) and contralateral control (left) femoral heads were harvested and frozen at −20°C for later imaging. The 48-hour timepoint corresponds to acute ischemic damage to the femoral head, whereas the 4-week timepoint corresponds to the early repair stage including revascularization of the epiphyseal cartilage.6 One of the surgeries at the 4-week postsurgical timepoint failed to induce ischemic injury to the femoral head, as confirmed by histological assessment and was excluded. Thus, a total of n = 7 femoral head pairs were evaluated at the 4-week timepoint.

MRI

Femoral head specimens were individually thawed at 4°C and imaged using a preclinical 9.4T MRI system (Agilent Technologies; Santa Clara, CA) equipped with a Varian console and millipede radiofrequency coil (Varian NMR Systems; Palo Alto, CA). To avoid air-tissue interface artifacts at the cartilage surface of the femoral head, the specimens were immersed in Fomblin (Specialty Fluids; Castaic, CA). Imaging included a high-spatial-resolution 3D gradient echo (GRE) sequence with magnetic susceptibility weighting to provide detailed visualization of the epiphyseal cartilage vascularity. Typical sequence parameters were: field-of-view = 38.4 × 38.4 × 38.4 mm3; sampling matrix = 384 × 384 × 384; acquired spatial resolution = 0.1 × 0.1 × 0.1 mm3; repetition time / echo time (TR/TE) = 40/15 msec; flip angle = 15°; bandwidth (BW) = 37 Hz/px; and scan time = 100 minutes. Fat saturation was not applied. These scan parameters were chosen to match those used in prior studies.12,14,17 Imaging was performed by MRI scientists (C.P.J. and L.W.) with 12 and 10 years of experience, respectively.

Histology

Immediately after imaging, the specimens were bisected with a bone saw in a plane including the insertion point of the ligamentum teres and the apex of the greater trochanter and fixed in 10% neutral buffered formalin. The specimens were then decalcified in 10% EDTA, processed in paraffin, sectioned at 5 μm in one location oriented along the bisection plane, and stained with hematoxylin and eosin. Selected sections were additionally stained with toluidine blue and safranin O. Histological analysis was performed by a board-certified orthopedic veterinary pathologist (C.S.C.) with 30 years of experience to qualitatively evaluate the morphology of sites from both the operated and control femoral heads. The pathologist was blinded to the surgery to validate that ischemic injury was induced in the operated femoral heads. Additionally, regions of reparative changes were noted, including evidence of neovascularization in the epiphyseal cartilage.

Data Analysis

Each 3D susceptibility-weighted GRE MRI image was postprocessed using QSM to better visualize and quantify the vascularity of the epiphyseal cartilage.14,17 The workflow is summarized in Fig. 1. QSM postprocessing was performed in MatLab (R2013b; Math-Works; Natick, MA) using a script and methods provided by Bilgic et al,24 which included Laplacian phase unwrapping, sophisticated harmonic artifact reduction on phase (SHARP) filtering, and calculation of magnetic susceptibility values. First, the 3D GRE magnitude image was used to manually segment the cartilage overlying the femoral head with ITK-SNAP software,25 yielding a binary cartilage mask. Second, the cartilage mask was applied to the 3D GRE phase image, and the phase was unwrapped using a Laplacian filter.26 Third, SHARP filtering was applied to remove the background field.27 The SHARP filtering used an adaptive kernel size (9 × 9 × 9 voxels to 3 × 3 × 3 voxels, corresponding to filter lengths of ~0.9–0.3 mm) and the cartilage mask was eroded by two pixels (0.2 mm) to reduce artifacts at the edge of the cartilage.24 A truncation factor of 0.25 was selected based on visual quality. Finally, QSM values were calculated using a closed-form L2-regularized solution, which is one of four algorithms provided by Bilgic et al.24

An external file that holds a picture, illustration, etc.
Object name is nihms-1580176-f0001.jpg

Overview of QSM postprocessing steps to visualize and quantify cartilage canals of the femoral head. (a) Original 3D GRE magnitude image (single slice, cropped to a field-of-view of 23.4 × 23.4 mm2 for display). (b) Cartilage mask generated after segmentation of the 3D GRE magnitude image using ITK-SNAP. This mask was used to measure cartilage volume. (c) Original 3D GRE phase image. (d) Laplacian unwrapped phase of the segmented cartilage. (e) Field map after SHARP filtering. (f) QSM map generated after magnetic susceptibility calculation using an L2-regularized algorithm. (g) Cartilage canal vessel mask generated after applying a threshold value of 0.025 ppm to the QSM map. This mask was used to measure vessel volume and vessel density. (h) Vessel density map generated by calculating the local density of vessels in a sphere with a 2.1 mm diameter.

Cartilage volume, vessel volume, and vessel density were then quantitatively measured. Cartilage volume was measured as the total number of segmented cartilage voxels prior to QSM postprocessing (ie, using the cartilage mask). Vessel volume was measured as the total number of voxels remaining after thresholding the postprocessed QSM images; all voxels with a QSM value ≥0.025 ppm were assumed to belong to cartilage canals. This threshold value was visually found to robustly separate the cartilage canal voxels (which have high QSM values on the order of 0.040 ppm17) from the surrounding cartilage matrix. Vessel density was then calculated by dividing the vessel volume by the total cartilage volume in the QSM postprocessed image.

3D maps of vessel density were generated for each specimen by counting the number of cartilage canal voxels in the postprocessed QSM image (ie, those voxels with a QSM value ≥0.025 ppm) within a sphere of diameter 21 voxels (ie, volume = 4.9 mm3). Each voxel of the vessel density map was assigned the ratio of the number of cartilage canal voxels to the total number of voxels in the sphere (ie, 4945 voxels) centered at the same location. In cases where the sphere extended outside of the cartilage volume, the density calculation only included those voxels within the cartilage volume.

Lastly, the QSM images were subjectively assessed by two independent observers (C.P.J., an MRI scientist with 12 years of experience; and J.M.E., a board-certified musculoskeletal radiologist and MRI scientist with 30 years of experience) to identify changes in cartilage canal architecture. To visualize the cartilage canals, maximum intensity projections (MIPs) of the 3D QSM images were generated at multiple projection angles. The general appearance, location, and density of the cartilage canals were noted.

Statistical Analysis

Measurements of vessel volume, cartilage volume, and vessel density were compared between the operated and control femoral heads at the 48-hour and 4-week timepoints using two-tailed paired t-tests conducted using R statistical software (v. 3.3.1; R Foundation for Statistical Computing; Vienna, Austria). Given a target Type I error rate of α = 0.05 and six hypotheses (three measures × two timepoints), P < 0.008 was considered statistically significant after conservative Bonferroni correction for multiple comparisons. Secondary post-hoc two-tailed Student’s t-tests were also performed to compare measurements between the 48-hour and 4-week timepoints, with P < 0.05 (uncorrected) considered hypothesis generating.

Results

At the 48-hour timepoint, no differences in vessel volume (P = 0.479), cartilage volume (P = 0.503), or vessel density (P = 0.147) were detected between the operated and control femoral heads (Table 1 and Fig. 2). The paired specimens had very similar vessel volume (mean difference [Δ] = 1.0 ± 2.0 mm3), cartilage volume (Δ = −46 ± 99 mm3), and vessel density (Δ = 0.12 ± 0.09%). GRE and QSM images and vessel density maps from a representative pair of control and operated femoral heads at the 48-hour timepoint are shown in Fig. 3. 3D rotating MIPs of the QSM images for each femoral head are shown in Supplemental Movie 1. In both the control and operated femoral heads, the cartilage canals were oriented linearly and were evenly distributed across the extent of the epiphyseal cartilage. These observations were consistent across all three pairs of femoral heads at the 48-hour timepoint (see Supplemental Fig. S1).

An external file that holds a picture, illustration, etc.
Object name is nihms-1580176-f0002.jpg

Quantitative analysis of vessel volume, cartilage volume, and vessel density of paired control and operated femoral heads. Plotted are the means and standard deviations for the three measurements at both the 48-hour (gray bars) and 4-week (black bar) timepoints. Results are shown for the control femoral heads, operated femoral heads, and the paired difference (Δ) between the operated and control femoral heads. No differences were detected between the operated and control femoral heads 48 hours following surgical induction of complete femoral head ischemia. In contrast, 4 weeks following surgery, all three measures were significantly increased in the operated vs. control femoral heads (*P < 0.008 for the paired differences).

An external file that holds a picture, illustration, etc.
Object name is nihms-1580176-f0003.jpg

Images of a representative pair of control and operated femoral heads 48 hours after surgical induction of ischemia. Shown are: (i) single slices of the 3D GRE magnitude images (cropped to a field-of-view of 23.9 × 17.8 mm2 for display), on which the cartilage canals appear dark against a bright cartilage background; (ii) thin-slab MIPs (2.0 mm) of the 3D QSM images at the same location as the GRE images, on which the cartilage canals appear bright (ie, higher magnetic susceptibility values compared to the cartilage background); (iii) full-volume MIPs of the 3D QSM images, which show the 3D architecture of the cartilage canals; and (iv) full-volume MIPs of 3D vessel density maps derived from the QSM images, which show the regional vessel density across the femoral head. The two femoral heads have a similar appearance on all images, including similar size of the secondary ossification center, thickness of the cartilage, and vessel architecture and density, which demonstrates that the QSM method is not sensitive to the acute effects of ischemia. Also see Supplemental Movie 1 and Supplemental Figure S1.

TABLE 1.

Quantitative Measurements of Epiphyseal Cartilage Vascularity in Paired Control and Operated Femoral Heads 48 Hours and 4 Weeks Following Surgical Induction of Complete Femoral Head Ischemia

48 Hours4 Weeks
ControlOperatedPaired Diff.ControlOperatedPaired Diff.
Vessel Volume (mm3)15.8 ± 4.516.8 ± 5.91.0 ± 2.013.6 ± 3.846.9 ± 17.033.3 ± 15.5 *
Cartilage Volume (mm3)1296 ± 1141250 ± 80−46 ± 991689 ± 2662271 ± 323582 ± 289 *
Vessel Density (%)1.21 ± 0.311.33 ± 0.400.12 ± 0.090.80 ± 0.162.02 ± 0.551.22 ± 0.55 *

Values are reported as mean ± standard deviation.

*Statistically significant (P < 0.008).

In contrast, significant differences between the operated and control femoral heads were detected at the 4-week timepoint. Quantitatively, the operated femoral heads had significantly greater vessel volume (P = 0.001; Δ = 33.3 ± 15.5 mm3), cartilage volume (P = 0.002; Δ = 582 ± 289 mm3), and vessel density (P = 0.001; Δ = 1.22 ± 0.55%) than their paired controls (Table 1 and Fig. 2). Compared to the collective operated and control findings at the 48-hour timepoint, the operated femoral heads at the 4-week timepoint had greater vessel volume (P = 0.001; Δ = 30.6 ± 17.6 mm3), cartilage volume (P < 0.001; Δ = 998 ± 336 mm3), and vessel density (P = 0.014; Δ = 0.75 ± 0.64%). By comparison, the control femoral heads at the 4-week timepoint had similar vessel volume to the collective findings at 48 hours (P = 0.273; Δ = −2.7 ± 6.1 mm3) but greater cartilage volume (P = 0.004; Δ = 416 ± 282 mm3), resulting in reduced vessel density (P = 0.006; Δ = −0.47 ± 0.36%). Thus, whereas the control femoral heads had reduced vessel density over time, the operated femoral heads had increased vessel density.

Qualitatively, increased vessel density in the operated femoral heads at the 4-week timepoint was primarily localized to the epiphyseal cartilage near the metaphyseal physis. GRE and QSM images and vessel density maps from a representative pair of operated and control femoral heads at the 4-week timepoint are shown in Fig. 4. Comparative rotating MIPs of their 3D QSM images are shown in Supplemental Movie 2. Clusters of brush-like vessels were seen near the metaphyseal physis in the operated femoral head, whereas the control femoral head had the same linear architecture observed at the 48-hour timepoint. Additionally, the operated femoral head had thicker cartilage and smaller SOC due to the cessation of blood flow, which produces a growth arrest of the SOC. The increased density of vessels near the metaphyseal physis was further evident in the vessel density map. This pattern of increased vessel density was consistent across all pairs of femoral heads at 4-week timepoint (see Supplemental Fig. S2). The histological appearance of cartilage canals for a pair of control and operated femoral heads are shown in Fig. 5a. While the cartilage canals of the control femoral head appeared singly in cross-section on histology, the cartilage canals of the operated femoral head were composed of a large central vessel surrounded by many smaller vessels, which is characteristic of neovascularization. For comparison, targeted MIPs of 3D QSM images from control and operated femoral heads are shown in Fig. 5b, which have a cartilage canal architecture consistent with the histological findings.

An external file that holds a picture, illustration, etc.
Object name is nihms-1580176-f0004.jpg

Images of a representative pair of control and operated femoral head 4 weeks after surgical induction of ischemia. Images are the same as described in Fig. 2. The field-of-view of the images was cropped to 32.2 × 24.0 mm2 for display. In contrast to the results at 48 hours, the control and operated femoral heads have significant differences at the 4-week timepoint. In addition to thicker cartilage and a smaller secondary ossification center in the operated femoral head (due to lack of development), increased vessel density is evident in the epiphyseal cartilage near the metaphyseal physis (arrows). On the 3D QSM full-volume MIP, this dense vasculature can be seen to have a brush-like architecture that is different from that seen in the control femoral head or at the 48-hour timepoint (cf. Fig. 2). The vessel density map shows regional density near the metaphyseal physis exceeding 15%, which is much greater than that seen in the control femoral head or the results at 48 hours (cf. Fig. 2). Also see Supplemental Movie 2 and Supplemental Figure S2.

An external file that holds a picture, illustration, etc.
Object name is nihms-1580176-f0005.jpg

Histological appearance of cartilage canals in control and operated femoral heads 4 weeks following surgical induction of complete femoral head ischemia. (a) Histological sections stained with safranin-O show distinct differences in the appearance of cartilage canals in a control vs. operated femoral head. In the images, the cartilage is stained red (articular cartilage is on the right side), the trabecular bone of the secondary ossification center is stained blue/green, and the clear spaces in the cartilage represent vessels. While the cartilage canals are singular and have a circular appearance in the control femoral head (arrowheads), they are multiplied (with a large central vessel surrounded by many smaller vessels) in the operated femoral head (arrows). (b) The distinct cartilage canal morphology seen histologically is consistent with the 3D appearance of the cartilage canals using QSM, from which targeted MIPs of cartilage canals are shown on a similar scale as the histology images. Although cartilage canals in control femoral heads have a normal linear appearance (arrowheads), cartilage canals in operated femoral heads have dense brush-like architecture (arrows).

Discussion

Our results demonstrate that QSM can detect vascularity in the epiphyseal cartilage of the developing femoral head without the use of an exogeneous contrast agent in a piglet model of LCPD. Furthermore, QSM can detect increased density of vessels in the epiphyseal cartilage following ischemic injury to the femoral head, which can be attributed to neovascularization during the early repair phase. Histological assessment confirmed the finding.

Our findings at the 4-week timepoint support the hypothesis that QSM can detect neovascularization in the epiphyseal cartilage. Comparison of the control and operated femoral heads revealed significantly greater vessel density in the operated femoral heads despite having significantly greater cartilage volume due to interruption of endochondral ossification.6 The greater vessel density was primarily localized to the epiphyseal cartilage near the metaphyseal physis, which is where vessels supplying the SOC are known to originate.6,10 Importantly, the vessel density of the operated femoral heads was also greater at 4 weeks vs. 48 hours postsurgery. This supports the previously reported finding in the piglet model that the epiphyseal cartilage expresses increased levels of angiogenic factors that promote new vessel ingrowth into the cartilage after ischemic injury.6 The dense, brush-like appearance of the cartilage canals in the operated femoral heads at the 4-week timepoint were fundamentally different than the linear appearance of cartilage canals for all other examined femoral heads. This brush-like orientation of vessels was confirmed histologically to be indicative of neovascularization. A similar appearance of neovascularization was recently observed in a goat model of osteochondritis dissecans, another developmental disease that affects epiphyseal vascularity.17

Our results at the 4-week timepoint are consistent with the μCTA findings of Kim et al6 using the same piglet model and a similar experimental design. Kim et al also found that both vessel volume and cartilage volume were significantly greater in the operated vs. control femoral heads. Kim et al also concluded, based on differences between the measurements in operated femoral heads at 48 hours and 4 weeks postsurgery, that the increase in vessel volume (despite a concomitant increase in cartilage volume) is a result of neovascularization. However, there were some discrepancies between our results and those of Kim et al at the 4-week timepoint. While cartilage volume measurements were similar using μCTA and QSM, vessel volume measurements using μCTA were about a factor of 10 lower than those using QSM. This discrepancy may be due to several factors: 1) μCTA images the vessel lumen of the cartilage canals, whereas QSM images the cartilage canal structure enclosing the lumen and may overestimate the size of the cartilage canals due to magnetic field disturbances near the cartilage canals; 2) QSM images both patent and occluded vessels, whereas μCTA images only patent vessels (ie, those that are accessible via injection of iodinated contrast material); and 3) our QSM images appear to show greater vessel detail than the μCTA images, which suggests that the QSM images had greater spatial resolution and/or sensitivity. Another discrepancy is that Kim et al found that the operated and control femoral heads had similar vessel density, while we found the operated femoral heads had significantly greater vessel density. This is primarily due to detection of a greater difference in vessel volume between the operated and control femoral heads using QSM. As noted above, this may due to QSM detecting finer vessel detail and/or the presence of both patent and occluded vessels. In either case, our finding of increased vessel density is consistent with the growth or invasion of new vessels.

The lack of a difference between the operated and control femoral heads at the 48-hour timepoint is a result of the sensitivity of QSM to the presence of cartilage canal vessels independent of blood flow. In contrast, the μCTA study by Kim et al did not detect vessels in the epiphyseal cartilage 48 hours postsurgery,6 since iodine contrast material injected into the extraosseous vessels feeding the femoral head was unable to reach the epiphyseal cartilage vessels due to the surgically induced vascular occlusion. Thus, while QSM has the advantage of being a noninvasive technique that may be suitable for in vivo imaging of children, its potential to detect vascular occlusion (eg, by measuring vessel oxygenation levels28,29) has not yet been explored in the epiphyseal cartilage. However, it has been previously observed that cartilage canal profiles can become less distinct and even vanish over an extended period in response to ischemic injury to epiphyseal cartilage,17 thus degenerative as well as proliferative changes in the cartilage canals may eventually be detectable.

The QSM technique may play an important role in the clinical management of patients with LCPD, as it potentially can help determine initiation of neovascularization and restoration of the epiphyseal cartilage function to sustain femoral head growth. The technique may also be useful for determining the efficacy of proangiogenic treatments to accelerate epiphyseal cartilage neovascularization and subsequent epiphyseal healing.6 However, further work is needed to translate this method for in vivo imaging at clinical magnetic field strengths. The use of long scan times and a preclinical 9.4T MRI system enabled detailed depiction of the vascularity of the epiphyseal cartilage, but image spatial resolution and quality will need to be reduced to accommodate in vivo imaging. This is potentially feasible, as QSM has recently been applied to examine the vascularity of the distal femoral condyle in children.18,19 Additionally, the distinct appearance of neovascularization compared with normal vasculature may aid in its detection at lower spatial resolution.

Scientific limitations of our study include its cross-sectional design, a somewhat limited sample size, and the use of an ex vivo animal model. Future work will focus on translation of these methods for in vivo studies to longitudinally assess changes in the vascularity of the epiphyseal cartilage during normal femoral head development and in response to ischemic injury.

In conclusion, QSM can detect neovascularization in the epiphyseal cartilage following ischemic injury to the femoral head. This noncontrast-enhanced method may prove useful for understanding the pathophysiology of LCPD and related childhood hip disorders, clinical staging of these diseases, and evaluating the efficacy of proangiogenic treatments.

Supplementary Material

supp fig 1

movie 1

movie 2

supp fig 2

Acknowledgments

Funding for this study was provided by the National Institute for Arthritis and Musculoskeletal and Skin Diseases (K01AR070894), the National Institute for Biomedical Imaging and Bioengineering (P41EB015894), the W. M. Keck Foundation, and Texas Scottish Rite Hospital for Children.

We thank Paula Overn for assistance with the preparation of the histological sections.

Footnotes

Additional supporting information may be found in the online version of this article.

References

1. Kim HK. Pathophysiology and new strategies for the treatment of Legg-Calve-Perthes disease. J Bone Joint Surg Am 2012;94:659–669. [Abstract] [Google Scholar]
2. Kim HK, Herring JA. Pathophysiology, classifications, and natural history of Perthes disease. Orthop Clin North Am 2011;42:285–295. [Abstract] [Google Scholar]
3. Kim HK, Burgess J, Thoveson A, Gudmundsson P, Dempsey M, Jo CH. Assessment of femoral head revascularization in Legg-Calve-Perthes disease using serial perfusion MRI. J Bone Joint Surg Am 2016;98: 1897–1904. [Abstract] [Google Scholar]
4. Chung SM. The arterial supply of the developing proximal end of the human femur. J Bone Joint Surg Am 1976;58:961–970. [Abstract] [Google Scholar]
5. Blumer MJ, Longato S, Fritsch H. Structure, formation and role of cartilage canals in the developing bone. Ann Anat 2008;190:305–315. [Abstract] [Google Scholar]
6. Kim HK, Bian H, Aya-ay J, Garces A, Morgan EF, Gilbert SR. Hypoxia and HIF-1alpha expression in the epiphyseal cartilage following ischemic injury to the immature femoral head. Bone 2009;45:280–288. [Abstract] [Google Scholar]
7. Stulberg SD, Cooperman DR, Wallensten R. The natural history of Legg-Calve-Perthes disease. J Bone Joint Surg Am 1981;63:1095–1108. [Abstract] [Google Scholar]
8. Herring JA, Kim HT, Browne R. Legg-Calve-Perthes disease. Part II: Prospective multicenter study of the effect of treatment on outcome. J Bone Joint Surg Am 2004;86–A:2121–2134. [Abstract] [Google Scholar]
9. Barnewolt CE, Shapiro F, Jaramillo D. Normal gadolinium-enhanced MR images of the developing appendicular skeleton: Part I. Cartilaginous epiphysis and physis. AJR Am J Roentgenol 1997;169:183–189. [Abstract] [Google Scholar]
10. Jaramillo D, Villegas-Medina OL, Doty DK, et al. Age-related vascular changes in the epiphysis, physis, and metaphysis: normal findings on gadolinium-enhanced MRI of piglets. AJR Am J Roentgenol 2004;182: 353–360. [Abstract] [Google Scholar]
11. Soares BP, Lequin MH, Huisman T. Safety of contrast material use in children. Magn Reson Imaging Clin N Am 2017;25:779–785. [Abstract] [Google Scholar]
12. Nissi MJ, Toth F, Zhang J, et al. Susceptibility weighted imaging of cartilage canals in porcine epiphyseal growth cartilage ex vivo and in vivo. Magn Reson Med 2014;71:2197–2205. [Europe PMC free article] [Abstract] [Google Scholar]
13. Toth F, Nissi MJ, Zhang J, et al. Histological confirmation and biological significance of cartilage canals demonstrated using high field MRI in swine at predilection sites of osteochondrosis. J Orthop Res 2013;31: 2006–2012. [Europe PMC free article] [Abstract] [Google Scholar]
14. Nissi MJ, Toth F, Wang L, Carlson CS, Ellermann JM. Improved visualization of cartilage canals using quantitative susceptibility mapping. PLoS One 2015;10:e0132167. [Europe PMC free article] [Abstract] [Google Scholar]
15. Toth F, Nissi MJ, Ellermann JM, et al. Novel application of magnetic resonance imaging demonstrates characteristic differences in vasculature at predilection sites of osteochondritis dissecans. Am J Sports Med 2015; 43:2522–2527. [Europe PMC free article] [Abstract] [Google Scholar]
16. Martel G, Kiss S, Gilbert G, et al. Differences in the vascular tree of the femoral trochlear growth cartilage at osteochondrosis-susceptible sites in foals revealed by SWI 3T MRI. J Orthop Res 2016;34:1539–1546. [Abstract] [Google Scholar]
17. Wang L, Nissi MJ, Toth F, et al. Quantitative susceptibility mapping detects abnormalities in cartilage canals in a goat model of preclinical osteochondritis dissecans. Magn Reson Med 2017;77:1276–1283. [Europe PMC free article] [Abstract] [Google Scholar]
18. Dymerska B, Bohndorf K, Schennach P, Rauscher A, Trattnig S, Robinson SD. In vivo phase imaging of human epiphyseal cartilage at 7 T. Magn Reson Med 2018;79:2149–2155. [Abstract] [Google Scholar]
19. Kolb A, Robinson S, Stelzeneder D, et al. Vessel architecture in human knee cartilage in children: an in vivo susceptibility-weighted imaging study at 7 T. Eur Radiol 2018;28:3384–3392. [Europe PMC free article] [Abstract] [Google Scholar]
20. Perry DC, Hall AJ. The epidemiology and etiology of Perthes disease. Orthop Clin North Am 2011;42:279–283. [Abstract] [Google Scholar]
21. Perry DC, Bruce CE, Pope D, Dangerfield P, Platt MJ, Hall AJ. Legg-Calve-Perthes disease in the UK: geographic and temporal trends in incidence reflecting differences in degree of deprivation in childhood. Arthritis Rheum 2012;64:1673–1679. [Abstract] [Google Scholar]
22. Kim HK, Su PH, Qiu YS. Histopathologic changes in growth-plate cartilage following ischemic necrosis of the capital femoral epiphysis. An experimental investigation in immature pigs. J Bone Joint Surg Am 2001;83–A:688–697. [Abstract] [Google Scholar]
23. Kim HK, Su PH. Development of flattening and apparent fragmentation following ischemic necrosis of the capital femoral epiphysis in a piglet model. J Bone Joint Surg Am 2002;84–A:1329–1334. [Abstract] [Google Scholar]
24. Bilgic B, Fan AP, Polimeni JR, et al. Fast quantitative susceptibility mapping with L1-regularization and automatic parameter selection. Magn Reson Med 2014;72:1444–1459. [Europe PMC free article] [Abstract] [Google Scholar]
25. Yushkevich PA, Piven J, Hazlett HC, et al. User-guided 3D active contour segmentation of anatomical structures: significantly improved efficiency and reliability. Neuroimage 2006;31:1116–1128. [Abstract] [Google Scholar]
26. Li W, Wu B, Avram AV, Liu C. Magnetic susceptibility anisotropy of human brain in vivo and its molecular underpinnings. Neuroimage 2012; 59:2088–2097. [Europe PMC free article] [Abstract] [Google Scholar]
27. Schweser F, Deistung A, Lehr BW, Reichenbach JR. Quantitative imaging of intrinsic magnetic tissue properties using MRI signal phase: an approach to in vivo brain iron metabolism? Neuroimage 2011;54: 2789–2807. [Abstract] [Google Scholar]
28. Rauscher A, Sedlacik J, Barth M, Haacke EM, Reichenbach JR. Noninvasive assessment of vascular architecture and function during modulated blood oxygenation using susceptibility weighted magnetic resonance imaging. Magn Reson Med 2005;54:87–95. [Abstract] [Google Scholar]
29. Shen Y, Kou Z, Kreipke CW, Petrov T, Hu J, Haacke EM. In vivo measurement of tissue damage, oxygen saturation changes and blood flow changes after experimental traumatic brain injury in rats using susceptibility weighted imaging. Magn Reson Imaging 2007;25:219–227. [Abstract] [Google Scholar]

Citations & impact 


Impact metrics

Jump to Citations

Citations of article over time

Article citations


Go to all (8) article citations

Data 


Data behind the article

This data has been text mined from the article, or deposited into data resources.

Similar Articles 


To arrive at the top five similar articles we use a word-weighted algorithm to compare words from the Title and Abstract of each citation.


Funding 


Funders who supported this work.

NIAMS NIH HHS (1)

NIBIB NIH HHS (1)

NIH HHS (1)

W. M. Keck Foundation