# The magnetization transfer ratio of the post-mortem canine intervertebral disk is positively correlated to Pfirrmann grading on high field 3.0T MRI: a pilot study

Antonin Martenne-Duplan, Corentin Tricou, Marlène Finck, Benjamin Cartiaux, Germain Arribarat, Giovanni Mogicato

► To cite this version:

Antonin Martenne-Duplan, Corentin Tricou, Marlène Finck, Benjamin Cartiaux, Germain Arribarat, et al.. The magnetization transfer ratio of the post-mortem canine intervertebral disk is positively correlated to Pfirrmann grading on high field 3.0T MRI: a pilot study. Frontiers in Veterinary Science, 2024, 11, pp.1335331. 10.3389/fvets.2024.1335331. hal-04696509

HAL Id: hal-04696509

https://ut3-toulouseinp.hal.science/hal-04696509v1

Submitted on 13 Sep 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for updates

OPEN ACCESS

EDITED BY
Mindaugas Tamosiunas,
University of Latvia, Latvia

REVIEWED BY
Hakyoung Yoon,
Jeonbuk National University, Republic of Korea
Renaldas Raisutis,
Kaunas University of Technology, Lithuania

*CORRESPONDENCE
Antonin Martenne-Duplan
✉ antonin.martenne@gmail.com

RECEIVED 08 November 2023
ACCEPTED 05 February 2024
PUBLISHED 14 February 2024

**CITATION**
Martenne-Duplan A, Tricou C, Finck M, Cartiaux B, Arribarat G and Magicato G (2024) The magnetization transfer ratio of the post-mortem canine intervertebral disk is positively correlated to Pfirrmann grading on high field 3.0T MRI: a pilot study.
*Front. Vet. Sci.* 11:1335331.
doi: 10.3389/fvets.2024.1335331

**COPYRIGHT**
© 2024 Martenne-Duplan, Tricou, Finck, Cartiaux, Arribarat and Magicato. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

# The magnetization transfer ratio of the post-mortem canine intervertebral disk is positively correlated to Pfirrmann grading on high field 3.0T MRI: a pilot study

Antonin Martenne-Duplan¹*, Corentin Tricou², Marlène Finck¹, Benjamin Cartiaux³, Germain Arribarat⁴ and Giovanni Magicato³

¹Centre Hospitalier Vétérinaire Massilia, Animedis, Marseille, France, ²Anatomy and Diagnostic Imaging Department of Ecole Nationale Vétérinaire de Toulouse (ENVT), Université de Toulouse, Toulouse, France, ³Toulouse Neuroimaging Center (ToNIC), University of Toulouse Paul Sabatier-INSERM-ENVT, Toulouse, France, ⁴Toulouse Neuroimaging Center (ToNIC), University of Toulouse Paul Sabatier-INSERM, Toulouse, France

**Objective:** Intervertebral disk (IVD) degeneration usually occurs earlier in chondrodystrophic dog breeds than in other breeds. Spinal cord compression secondary to IVD degeneration is the most common cause of myelopathy in these dogs. Standard magnetic resonance imaging (MRI) sequences permit the identification of IVD degeneration and its consequences on adjacent neurological structures. In human medicine, quantitative MRI sequences, such as magnetization transfer ratio (MTR) sequences, are developed and used to detect early IVD degeneration. This prospective randomized post-mortem comparative study aimed to evaluate the correlation between a qualitative Pfirrmann MRI grading and the MTR values of the IVD in chondrodystrophic dogs.

**Materials and methods:** Vertebral columns of eight canine cadavers were frozen and thawed prior to imaging with T2-weighted and MTR sequences using a 3.0 T high-field MRI. These sequences were reviewed by two observers. A Spearman correlation coefficient was calculated in order to compare the MTR values with the Pfirrmann grade. Pearson correlation coefficients were calculated to evaluate the inter-observer agreement of the delineation of the region of interest (ROI) around the NP and the MTR values. A Wilcoxon-Mann-Whitney test was used to conclude on the significance of the correlation between the MTR values and the Pfirrmann grades.

**Results:** There were 138 intervertebral disks analyzed: 29/138 (21.0%) IVD were grade I, 74/138 (53.6%) grade II, and 35/138 (25.4%) grade III. No grades IV and V were present in this study. Inter-observer agreement for delineation of IVD ROI was fair (r = 0.54) but inter-observer agreement of mean MTR value within the ROI was very good (r = 0.89). Mean MTR values were 16.459% (10.0305–21.0950%) for grade I, 18.888% (10.0750–27.2400%) for grade II, and 22.813% (12.5700–31.7600%) for grade III. The mean MTR value was significantly different between each Pfirrmann grade: between grades I and II (*p* < 0.005), grades II and III (*p* < 0.05), and grades I and III (*p* < 0.005). There was a significant moderate positive correlation between Pfirrmann grading and mean MTR values (r = 0.516).
**Conclusion:** The magnetization transfer ratio seems to be an objective method to detect early intervertebral disk degeneration via quantitative analysis.

KEYWORDS

dog, MRI, nucleus pulposus, degeneration, MTR

# 1 Introduction

The intervertebral disk (IVD) is composed of an incompressible nucleus pulposus (NP) surrounded by an annulous fibrosus (AF), and by the vertebral endplates on its cranial and caudal aspects (1). The role of the IVD, with the ligamentous structures, is to limit vertebral column motion in any plane and to resist compressive forces along the long axis of the vertebral column that are created by the weight of the head and the abdominal muscles (2).

IVD degeneration is a common disorder in dogs (3). It occurs by two main pathophysiologic mechanisms, chondroid and fibroid degenerations, also named Hansen type 1 and Hansen type 2, respectively, (2). IVD degeneration is a common finding in clinically neurological normal dogs, especially in chondrodystrophic breeds with frequent and early chondroid degeneration of the NP. In veterinary medicine, magnetic resonance imaging (MRI) is the gold standard for evaluating IVD degeneration and spinal cord consequences due to its excellent tissular contrast in the case of intervertebral disk herniation (4). A 5-category Pfirrmann system is used in human medicine to classify lumbar IVD degeneration. This scoring takes into account the intensity of the signal and the height of the IVD. This scoring has been previously used in veterinary medicine and has shown a great correlation to histopathological degeneration (4). However, this grading is a qualitative method with limitations especially concerning inter-operator agreement.

In human medicine, MRI quantitative IVD evaluation has been developed using quantitative MRI sequences such as magnetization transfer ratio (MTR). MTR sequence provides a measurement of the rate of magnetization exchange between the free hydrogen protons in fluid and those that are bound to macromolecules in tissue and is highly predictive of histopathological changes (5). The MT is assessed by applying an off-resonance RF pulse at a different frequency than the Larmor frequency of free-water protons. This RF only saturates the protons bound to macromolecules, such as collagen. Collagen has been found to produce a significant magnetization transfer effect (6). Collagen is the main macromolecule of IVD AF. Collagen is also present in the physiological NP in lower proportion. This macromolecule is produced by chondrocyte-like cells during chondroid IVD degeneration of the NP. A previous study using a canine model showed that early degenerative or traumatic IVD changes are detectable with both T2WI and MTR MR imaging and are well correlated to histopathological findings (7). This same study showed that MTR increases with the severity of T2WI IVD signal abnormalities and IVD degeneration degree on histopathology. Another human study showed a moderate and significant linear correlation between MTR and Pfirrmann grades suggesting that collagen relative density increases with the degeneration of the NP. It showed that MTR MRI may serve as a non-invasive diagnostic tool for IVD degeneration (5). However, to our knowledge, IVD

degeneration grade evaluation using MTR has never been attempted in veterinary research.

So the aim of this study was to evaluate the correlation between the Pfirrmann intervertebral disk degeneration grading on T2WI images with the quantitative magnetization transfer ratio values of post-mortem intervertebral disks on a 3.0 T high-field MRI. It also aimed to evaluate the inter-observer agreement of the delineation of the region of interest (ROI) around the NP and the MTR values.

Our first hypothesis was that quantitative MTR would be significantly positively correlated to Pfirrmann grading. The second hypothesis was that MTR values have a high inter-observer agreement.

# 2 Materials and methods

This study was a prospective randomized post-mortem comparative study. The experimental protocol and euthanasia conditions were approved by the Ethical Committee 'Sciences et Santé Animale 115 - Ecole Nationale Vétérinaire de Toulouse with the authorization number APAFIS# 21559-2019071917392588v3. The dogs were raised in the Avogadro Laboratory (Fontenilles, France) and were euthanized for pedagogic purposes unrelated to this study at the Ecole Nationale Vétérinaire de Toulouse (ENVT) in the Anatomy Department. Cervical and thoracic vertebral columns were harvested from the cadavers and frozen at -20°C for approximately 5 months. Four days before the MRI acquisition, vertebral columns were placed in a cold room at 4°C to achieve slow defrosting. Once completely thawed, the vertebral columns were stored at 4°C. On the day of the MRI, vertebral columns were stored in a cooler and transported from the Pavillon Baudot of the Center Hospitalier Universitaire de Purpan (Toulouse, France) to perform the acquisition.

## 2.1 Image acquisition

MR Imaging was performed using a high field 3.0 Tesla magnet (Achieva, Philips, Amsterdam, The Netherlands). For imaging, the vertebral column was placed with the cervical part facing the MR gantry and the ventral aspect of the vertebral bodies facing the patient table to mimic a “sternal recumbency”-like position. “Dorsal recumbency”-like position was not achievable due to instability of the vertebral column on the table during acquisition. A 16-channel flexible elbow coil (Philips) was placed around the vertebral column. The field of view included the cervical and thoracic spine. T2-weighted images were obtained in the sagittal plane with the detailed MR acquisition parameters displayed in Table 1. Then a magnetization transfer (MT) sequence was performed in the sagittal plane using the same placement of slices as the T2-weighted sequence with detailed acquisition parameters also displayed in Table 1. The
<table><caption>TABLE 1 Image acquisition parameters for the T2-weighted images and for the magnetization transfer ratio (MTR) sequence.</caption><thead><tr><th>Sequence</th><th>T2WI sagittal</th><th>MTR sagittal (M0)</th><th>MTR sagittal (MSAT)</th></tr></thead><tbody><tr><td>Repetition time TR (ms)</td><td>3,000</td><td>72</td><td>72</td></tr><tr><td>Echo time TE (ms)</td><td>90</td><td>4,5</td><td>4,5</td></tr><tr><td>Angle flip (°)</td><td>18°</td><td>18°</td><td>18°</td></tr><tr><td>Matrix</td><td>400*400*16</td><td>400*400*16</td><td>400*400*16</td></tr><tr><td>Voxel size (mm)</td><td>0,7*0,7</td><td>0,7*0,7</td><td>0,7*0,7</td></tr><tr><td>Slice thickness (mm)</td><td>2</td><td>2</td><td>2</td></tr><tr><td>Interslice gap (mm)</td><td>0</td><td>0</td><td>0</td></tr><tr><td>Number of slices</td><td>16</td><td>16</td><td>16</td></tr><tr><td>Echo train (slice)</td><td>1</td><td>1</td><td>1</td></tr><tr><td>Acquisition number</td><td>3</td><td>3</td><td>3</td></tr><tr><td>Acquisition time (min)</td><td>3</td><td>27</td><td>27</td></tr><tr><td>Band width (Hz)</td><td>1,200</td><td>1,200</td><td>1,200</td></tr><tr><td>Pre-saturation impulsion RF</td><td>-</td><td>-</td><td>1,000 Hz for 2.5 ms</td></tr></tbody></table>

MTR data were obtained using a sagittal gradient echo sequence with dual acquisition and collected with or without the application of a magnetization transfer pre-pulses. One acquisition was performed with an off-resonance pulse applied at 1000 Hz below the free water proton resonance frequency for 2.5 ms (MSAT) and the second one without the off-resonance pulse (M0). Then, a third set of images consisting of the Magnetization Transfer Ratio images was calculated on a pixel-by-pixel algorithm using the formula: MTR = (M0-MSAT)/M0. The signal value of MTR can be expressed as the ratio value (0-1) or as the percentage of saturation (0-100%) calculated with the formula: MTR (%) = [(M0-MSAT)/M0]x100. In the result section below, the authors decided to express the results as the percentage of saturation.

## 2.2 Qualitative analysis - Pfirrmann grading

T2-weighted images were reviewed by a first-year ECVDI veterinary radiology resident (AMD) and a ECVDI-certified veterinary radiologist (MF) using a Digital Imaging and Communication in Medicine (DICOM)-reader system (3.0, Horos Project, Geneva, Switzerland) in a high-resolution screen workstation (iMac, Apple, USA). To grade each IVD space, the image was adjusted to display the highest height and length of the NP hyperintensity. The image used for grading could be different in the same vertebral column depending on the IVD of interest. The observers graded the IVD by consensus using the 5-category Pfirrmann system to grade the degree of IVD degeneration. Due to the post-mortem nature of this study, the comparison between NP and cerebrospinal fluid signals was not evaluable for the grading. The IVD was considered grade I if the

disc was homogeneous with bright hyperintense signal and normal disc height; grade II if the disc was inhomogeneous but with preserved hyperintense signal, clear differentiation between the NP and the AF and normal height; grade III if the disk was inhomogeneous with intermittent gray signal intensity, unclear differentiation between the NP and the AF and normal to slightly decreased height; grade IV if the disk was inhomogeneous with hypointense dark gray signal intensity, absence of distinction between the NP and the AF and a slightly to moderately decreased size and grade V if the disk was inhomogeneous with hypointense black signal intensity, absence of distinction between the NP and the AF and a collapsed disk space. Both observers were familiar with the Pfirrmann grading procedure and had trained together prior to this study on various human and dog T2WI images to homogenize their grading.

## 2.3 Quantitative analysis - magnetization transfer ratio

For the comparative study, ROI was drawn separately by a first-year ECVDI veterinary radiology resident (AMD) and a final-year veterinary student (CT) using a DICOM reading system (MRICron, 2007, USA) for each IVD. The ROI was drawn over the hyperintense signal of the NP in the T2WI images (See figure in Supplementary materials file). Each ROI was saved and then transferred to another software (Sisyphe, France) using the “copy/paste” function on the corresponding image of the MTR sequence (8, 9). For each ROI, the following data were recorded: ROI volume, minimal MTR value (MTRmin), mean MTR value (MTRmean), median MTR value (MTRmed), and maximal MTR value (MTRmax). For the comparative study, for each ROI, the average value of both operators’ MTRmean was used (mMTR). Then, for all the IVD in the same Pfirrmann grade, an average of all the IVD values was calculated (mean MTR).

## 2.4 Comparative study - statistical analysis

Results were assessed by using commercially available softwares: Excel 2018 (16.16.27 version, Microsoft Corporation, USA), R 2004-2016 (3.6.3 version, R Core Team, New Zealand) and Prism 10 (10.1.1 version, GraphPad Software, USA). Pearson’s correlation coefficients were used to assess the inter-operator agreement for the ROIs volume and the MTR mean of the nucleus pulposus. The agreement strength was classified as very poor (<0.19), poor (0.2-0.39), fair (0.4-0.59), good (0.6-0.79), and very good (>0.8). A Spearman correlation coefficient was also calculated in order to compare the MTR values with the Pfirrmann grade. The significance threshold was set to value of *p*<0.05 for all statistical sets.

In addition, box plots of MTR values as a function of Pfirrmann grade were produced using the Excel software (16.16.27 version, Microsoft Corporation, USA). The average value of the different MTR for each grade is represented by a black cross within the box.

Finally, to conclude on the significance of the correlation between the MTR values and the Pfirrmann grades, a Wilcoxon-Mann-Whitney test was used.
## 3 Results

Eight vertebral columns were included in this study. All dogs were beagle, considered as a chondrodystrophic breed. Dogs had a median age of 2 years and 7 months (range 23–36 months) and a median weight of 14.4 kg (range 12.8–16.0 kg) at the time of euthanasia. Five were neutered males and three were spayed females. Dogs presented no abnormalities on clinical and neurological examinations. The number of IVD visible on the MRI sequences was variable for each dog. A total of 138 IVD were analyzed in this study (range: 15–19 IVD per dog).

### 3.1 Pfirrmann grading and MTR sequences

After T2WI assessment, 29/138 IVD were grade I (21.01%), 74/138 grade II (53.62%), and 35/138 were grade III (25.36%). No grades IV and V were present in this study.

Four types of images were obtained from the magnetization transfer sequences: image without the off-resonance RF with the signal M0, image with the off-resonance saturation RF with the signal MSAT, image calculated with the MTR pixel-by-pixel algorithm (MTR = (M0-MSAT)/M0) and image of the MTR color map. These images are displayed in Figure 1.

A greater signal saturation in the MTR sequences was observed with increased Pfirrmann grade: the mean MTR value for Pfirrmann

grade I IVD degeneration was 16.459% (10.0305–21.0950%), the mean MTR value for Pfirrmann grade II IVD degeneration was 18.888% (10.0750–27.2400%), and the mean MTR value for Pfirrmann grade III IVD degeneration was 22.813% (12.5700–31.7600%). The overlaid MTR color maps (Figure 2) offer a clear visualization of the NP. In grade I (Figure 2A white arrays), the NP appears as a homogeneous low-saturated signal (red), and around the NP the well-defined annulus fibrosus of mildly increased saturation (yellow). In grade II (Figure 2A asterisk), the findings were variable with heterogeneity of the NP and/or increased saturation of the NP center, and the annulus fibrosus is less well demarcated than in grade I. In grade III (Figure 2B), the NP appears moderately heterogeneous with increased saturation and the annulus fibrosus is barely visible.

### 3.2 Statistical analysis

The Pearson's correlation coefficients for the inter-operator agreement for the measurement of ROIs volume and the MTR mean of the nucleus pulposus were fair (0.54) and very good (0.89), respectively (Figures 3, 4).

The Spearman correlation showed a significant positive moderate correlation (rho = 0.516; *p* < 0.001) between the Pfirrmann grades and the MTR values of the IVD (Figure 5).

The statistical analysis showed that there were significant differences between all the Pfirrmann grades and their correlated

<figure><img src="image_3.png" /></figure>

FIGURE 1

Different sets of images of the MTR sequence: Image without the off-resonance RF with the signal M0 (A). Image with the off-resonance saturation RF with the signal MSAT (B). Image calculated with the MTR pixel-by-pixel algorithm (MTR = (M0-MSAT)/M0) (C). Image of the MTR color map (D).
<figure><img src="image_4.png" /><figcaption>**FIGURE 2**
Overlaid MTR Color maps ROI of intervertebral disks. (A) – white arrows: grade I, the NP appears as a homogeneous low-saturated signal (red), and around the NP the well-defined annulus fibrosus of mildly increased saturation (yellow). (A) – asterisk: grade II, heterogeneity of the NP and the annulus fibrosus is less well demarcated than in grade I. (B): grade III, the NP appears moderately heterogeneous with increased saturation and the annulus fibrosus is barely visible.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**FIGURE 3**
Scatter diagram depicting the delineation of IVD ROI between observers. The Pearson correlation coefficient (0.54) indicates a fair inter-observer difference.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**FIGURE 4**
Scatter diagram illustrating the inter-observer MTR mean of the nucleus pulposus. The Pearson correlation coefficient (0.89) indicates a very good inter-observer agreement.</figcaption></figure>

MTR values: between grades I and II ($p < 0.005$), grades II and III ($p < 0.05$), and grades I and III ($p < 0.005$). These results are depicted in box plots in Figure 6.

## 4 Discussion

The standard assessment for IVD degeneration is based on T2-weighted images of the NP and also on the Pfirrmann grading depending on the radiologist's preference. The Pfirrmann grading only relies on qualitative observations such as the signal intensity and the shape of the NP in T2-weighted images (4, 10). This is prone to intrapersonal and interpersonal biases and as the Pfirrmann grading is a numerological grading, some IVD degeneration may share features of two adjacent grades. In contrast, the MTR images are based on quantitative measurement with specificity of certain macromolecules content, providing a more objective assessment than the previously mentioned methods (11). Therefore, this study evaluated correlations between the IVD NP MTR values and Pfirrmann grade in dogs.

In human studies, quantitative MRI techniques have been employed to identify IVD degeneration by linking MRI signals with the biochemical components of the NP and AF. Superior to semi-quantitative morphological MRI methods like traditional T2-weighted imaging, quantitative MRI techniques, particularly MTR imaging, exhibited heightened accuracy (12–14). A recent study with 67 human patients affected by IVD degeneration at the cervicothoracic junction revealed a positive correlation between MTR values of the IVD NP and Pfirrmann grades (15). In our study, we observed similar results with comparable MTR values.

IVD degeneration stands out as one of the most prevalent neurologic disorders in dogs (3). Traditional diagnostic measures for IVD degeneration encompass clinical signs, neurological manifestations, and imaging observations. While various diagnostic imaging techniques exist, MRI holds the status of the preferred method for delineating the stage of disc degeneration, the extent of herniation, and concurrent spinal cord parenchymal lesions in dogs.
(16–18). Conventional MRI protocols in canines typically involve two-dimensional (2D) fast spin-echo (FSE) sequences replicated across multiple planes, and more recently, three-dimensional (3D) FSE (19). However, quantitative MRI techniques, especially Magnetization Transfer Ratio (MTR) imaging, remain underutilized. A recent study has shown modifications of the spinal cord MTR values secondary to IVD herniation with significantly decreased MTR values in chronic cases compared to acute cases (20). Another experiment using a canine model demonstrated the detectability of early traumatic or degenerative changes of the IVD through MTR, confirmed by histologic and biochemical analysis (7). Nevertheless, as of now, no study has assessed the application of MTR imaging in correlation with Pfirrmann grading for IVD degeneration.

In our study, the MTR values within the ROI showed a strong inter-operator agreement despite a low inter-operator agreement for the ROI volume. IVD of grade I appear in MTR images as a homogeneous low-saturated signal with a well-defined border of mildly increased saturation. The higher saturation of the ventral and dorsal aspect of the NP might be secondary to volume averaging on T2-weighted images leading to inclusion within the ROI of the inner part of the AF or of the transitional zone with a higher collagen content than the NP resulting in higher saturation. In grade III, the NP appears moderately heterogeneous with increased saturation. The technique used to draw the ROI in our study has already been described in previous studies (21, 22) and the strong inter-observer agreement of the MTR values within the ROI in our study showed that the method of ROI drawing may be a reliable method.

Hansen's study published approximately 70 years ago described the pathological features of IVD degeneration in dogs and classified it into two distinct forms each typically occurring in different types of dogs (2). Despite general support of Hansen's original description, more recent histologic analysis suggested that there is less difference

between degeneration in chondrodystrophic and non-chondrodystrophic dogs than was previously assumed. In fact, advanced AF degeneration in non-chondrodystrophic dogs has been shown to be associated with the replacement of notochord-like cells by chondrocyte-like cells similar to IVD degeneration in chondrodystrophic dogs (23, 24). In our study only the MTR values of the MP were analyzed with a ROI drawn over the hyperintense signal of the NP in the T2WI images. So it would be interesting to study the MTR values of the AF. Moreover only Beagle dogs, considered as a chondrodystrophic breed, were imaged. The collagen content of the IVD has been shown to vary with the breed and more specifically with chondrodystrophic status (25). A further study is necessary to know if our results might be transposable to non-chondrodystrophic dogs as well.

<figure><img src="image_7.png" /><figcaption>FIGURE 5
Distribution of MTR values according to Pfirrmann grades. The Spearman correlation revealed a statistically significant positive moderate correlation (rho = 0.516; *p* < 0.001) between Pfirrmann grades and MTR values of the IVD.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>FIGURE 6
Variation of the MTR values in function of the Pfirrmann grading represented in plot boxes. The black cross represents the average of MTR within one Pfirrmann grade.</figcaption></figure>
The Pfirrmann grading system focuses on IVD characteristic changes such as signal intensity, disk structure, differentiation between AF and NP, and IVD height but not on changes in the tissues surrounding the IVD. These changes in the surrounding tissues should be included to obtain a complete picture of the status of IVD. In fact, to accurately evaluate the IVD degeneration in dogs, the Pfirrmann grading system should be used in conjunction with the presence or absence of disk herniation and abnormalities of the other structures composing the IVD space (17). To obtain this information, transverse and midsagittal MRI images are needed. In our study, this modified Pfirrmann wasn't usable because of the absence of CSF due to the transaction of the central nervous system at the level of the foramen magnum.

Several limitations are acknowledged with this study design. The first main limitation of this study is the frozen post-mortem *ex vivo* protocol. In human literature, the freezing process appears to be a good method of tissue conservation when the sample is rapidly imaged after defrosting. This technique showed minimal alteration of the image quality in MRI (26). In other species, it has been shown that frozen-thawed samples have a significantly increased MTR compared to lived samples, however the alteration of the MTR signal was homogeneous in the entire scan field of view (27). In our study, all the vertebral columns were homogeneously frozen and thawed. For this reason, difference of MTR between different Pfirrmann grades in our samples was considered interpretable in frozen-thawed cadavers. One limitation of this study is the restriction to the cervical and thoracic areas. In fact, a human study showed that thoracic IVD NP had a higher collagen content compared to the lumbar IVD NP, resulting in higher MT saturation in the thoracic region than in the lumbar region (5, 28). Another limitation of this study is the absence of advanced IVD degeneration of grades IV and V using the Pfirrmann grading. Further studies of the entire *in vivo* vertebral column with severely degenerated IVD are necessary to confirm our results in live animals.

Finally, one of the main clinical limitations of this study is the acquisition time of the MTR sequence of 27 min. This acquisition time seems non-compatible with daily clinical use and may only be feasible in clinically stable animals in whom significant prolongation of anesthesia time would not result in severely increased anesthetic risk. Furthermore, only sagittal images were acquired in our study but transverse images may provide new diagnostic information and lengthen even more the anesthetic time.

In conclusion, we statistically demonstrated a significant moderate positive correlation between Pfirrmann grading and the MTR values, showing that there is an increased magnetization transfer within the NP with IVD degeneration. The MTR MRI seems to be an objective method to detect early IVD degeneration via quantitative analysis. Nonetheless, additional investigations are required to validate our findings in canine patients from both chondrodystrophic and non-chondrodystrophic breeds.

## Data availability statement

The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.

## Ethics statement

The animal study was approved by Sciences et Santé Animale 115 - Ecole Nationale Vétérinaire de Toulouse with the authorization number APAFIS# 21559-2019071917392588v3. The study was conducted in accordance with the local legislation and institutional requirements.

## Author contributions

AM-D: Conceptualization, Investigation, Methodology, Writing – original draft. CT: Conceptualization, Formal analysis, Validation, Writing – review & editing. MF: Conceptualization, Investigation, Methodology, Supervision, Validation, Writing – review & editing. BC: Conceptualization, Investigation, Supervision, Validation, Writing – review & editing. GA: Conceptualization, Methodology, Supervision, Validation, Writing – review & editing. GM: Conceptualization, Funding acquisition, Investigation, Methodology, Supervision, Validation, Writing – review & editing.

## Funding

The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.

## Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

## Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

## Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2024.1335331/full#supplementary-material
## References

1. Evans H, Christensen G. Joints and ligaments In: L Alexander, editor. Miller's *anatomy of the dog. 2nd ed* ed. Philadelphia: WB Saunders Co (1979). 225-68.

2. Hansen HJ. A pathologic-anatomical study on disc degeneration in dog, with special reference to the so-called enchondrosis intervertebralis. *Acta Orthop Scand Suppl.* (1952) 23:1-130. doi: 10.3109/ort.1952.23.suppl-11.01

3. Fluehmann G, Doherr MG, Jaggy A. Canine neurological diseases in a referral hospital population between 1989 and 2000 in Switzerland. *J Small Anim Pract.* (2006) 47:582-7. doi: 10.1111/j.1748-5827.2006.00106.x

4. Bergknut N, Auriemma E, Wijsman S, Voorhout G, Hagman R, Lagerstedt AS, et al. Evaluation of intervertebral disk degeneration in chondrodystrophic and nonchondrodystrophic dogs by use of Pfirrmann grading of images obtained with low-field magnetic resonance imaging. *Am J Vet Res.* (2011) 72:893-8. doi: 10.2460/ajvr.72.7.893

5. Wang C, Witschey W, Goldberg A, Elliott M, Borthakur A, Reddy R. Magnetization transfer ratio mapping of intervertebral disc degeneration. *Magn Reson Med.* (2010) 64:1520-8. doi: 10.1002/mrm.22533

6. Antoniou J, Mwale F, Demers CN, Beaudoin G, Goswami T, Aebi M, et al. Quantitative magnetic resonance imaging of enzymatically induced degradation of the nucleus pulposus of intervertebral discs. *Spine*. (2006) 31:1547-54. doi: 10.1097/01.brs.0000221995.77177.9d

7. Chen C, Jia Z, Han Z, Gu T, Li W, Li H, et al. Quantitative T2 relaxation time and magnetic transfer ratio predict endplate biochemical content of intervertebral disc degeneration in a canine model. *BMC Musculoskelet Disord.* (2015) 16:157. doi: 10.1186/s12891-015-0610-6

8. Tensaouti F, Lotterie JA, Clarisse P, Berry I. Sisyphe-neuroimaging software toolbox. *Eur Soc Magn Reson Med Bio.* (2008) 2008:2-4.

9. Tensaouti F, Lotterie JA, Clarisse P, Berry I. Sisyphe-software for neuroimaging In: K Saetzler, editor. *First international workshops on image processing theory, tools and applications (IPTA08)*. New York: Springer (2008)

10. Pfirrmann CW, Metzdorf A, Zanetti M, Hodler J, Boos N. Magnetic resonance classification of lumbar intervertebral disc degeneration. *Spine*. (2001) 26:1873-8. doi: 10.1097/00007632-200109010-00011

11. Gray ML, Burstein D, Lesperance LM, Gehrke L. Magnetization transfer in cartilage and its constituent macromolecules. *Magn Reson Med.* (1995) 34:319-25. doi: 10.1002/mrm.1910340307

12. Kim DK, Ceckler TL, Hascall VC, Calabro A, Balaban RS. Analysis of water-macromolecule proton magnetization transfer in articular cartilage. *Magn Reson Med.* (1993) 29:211-5. doi: 10.1002/mrm.1910290209

13. Baum T, Joseph GB, Karampinos DC, Jungmann PM, Link TM, Bauer JS. Cartilage and meniscal T2 relaxation time as non-invasive biomarker for knee osteoarthritis and cartilage repair procedures. *Osteoarthr Cartil.* (2013) 21:1474-84. doi: 10.1016/j.joca.2013.07.012

14. Xin H, Zhang C, Wang D, Shi Z, Gu T, Wang C, et al. Tissue-engineered allograft intervertebral disc transplantation for the treatment of degenerative disc disease: experimental study in a beagle model. *Tissue Eng Part A.* (2013) 19:143-51. doi: 10.1089/ten.tea.2012.0255

15. Zhang C, Lin Y, Han Z, Gao L, Guo R, Shi Q, et al. Feasibility of T2 mapping and magnetic transfer ratio for diagnosis of intervertebral disc degeneration at the cervicothoracic junction: a pilot study. *Biomed Res Int.* (2019) 2019:1-9. doi: 10.1155/2019/6396073

16. Mayhew PD, McLean RC, Ziemer LS, Culp WTN, Russell KN, Shofer FS, et al. Risk factors for recurrence of clinical signs associated with thoracolumbar intervertebral disk herniation in dogs: 229 cases (1994-2000). *J Am Vet Med Assoc.* (2004) 225:1231-6. doi: 10.2460/javma.2004.225.1231

17. Besalti O, Pekcan Z, Sirin YS, Erbas G. Magnetic resonance imaging findings in dogs with thoracolumbar intervertebral disk disease: 69 cases (1997-2005). *J Am Vet Med Assoc.* (2006) 228:902-8. doi: 10.2460/javma.228.6.902

18. Longo S, Gomes SA, Briola C, Duffy K, Targett M, Jeffery ND, et al. Association of magnetic resonance assessed disc degeneration and late clinical recurrence in dogs treated surgically for thoracolumbar intervertebral disc extrusions. *J Vet Intern Med.* (2021) 35:378-87. doi: 10.1111/jvim.15989

19. Lee S, Hwang J, Ko J, Jeong J, Eom K, Kim J. Comparison between T2-weighted two-dimensional and three-dimensional fast spin-echo MRI sequences for characterizing thoracolumbar intervertebral disc disease in small-breed dogs. *Vet Radiol Ultrasound.* (2022) 63:216-23. doi: 10.1111/vru.13049

20. Shinn RL, Pancotto TE, Stadler KL, Werre SR, Rossmeisl JH. Magnetization transfer and diffusion tensor imaging in dogs with intervertebral disk herniation. *J Vet Intern Med.* (2020) 34:2536-44. doi: 10.1111/jvim.15899

21. Hoppe S, Quirbach S, Mamisch TC, Krause FG, Werlen S, Benneker LM. Axial T2 mapping in intervertebral discs: a new technique for assessment of intervertebral disc degeneration. *Eur Radiol.* (2012) 22:2013-9. doi: 10.1007/s00330-012-2448-8

22. Ying J, Han Z, Zeng Y, Du Y, Pei S, Su L, et al. Evaluation of intervertebral disc regeneration with injection of mesenchymal stem cells encapsulated in PEGDA-microcryogel delivery system using quantitative T2 mapping: a study in canines. *Am J Transl Res.* (2019) 11:2028-41.

23. Bergknut N, Smolders LA, Grinwis GCM, Hagman R, Lagerstedt AS, Hazewinkel HAW, et al. Intervertebral disc degeneration in the dog. Part 1: anatomy and physiology of the intervertebral disc and characteristics of intervertebral disc degeneration. *Vet J.* (2013) 195:282-91. doi: 10.1016/j.tvjl.2012.10.024

24. Hansen T, Smolders LA, Tryfonidou MA, Meij BP, Vernooij JCM, Bergknut N, et al. The myth of fibroid degeneration in the canine intervertebral disc: a histopathological comparison of intervertebral disc degeneration in Chondrodystrophic and Nonchondrodystrophic dogs. *Vet Pathol.* (2017) 54:945-52. doi: 10.1177/0300985817726834

25. Ghosh P, Taylor TK, Braund KG, Larsen LH. The collagenous and non-collagenous protein of the canine intervertebral disc and their variation with age, spinal level and breed. *Gerontology.* (1976) 22:124-34. doi: 10.1159/000212129

26. Zarb F, McNulty J, Gatt A, Formosa C, Chockalingam N, Evanoff MG, et al. Comparison of in vivo vs. frozen vs. Thiel cadaver specimens in visualisation of anatomical structures of the ankle on proton density magnetic resonance imaging (MRI) through a visual grading analysis (VGA) study. *Radiography.* (2017) 23:117-24. doi: 10.1016/j.radi.2016.12.004

27. Chang EY, Bae WC, Statum S, Du J, Chung CB. Effects of repetitive freeze-thawing cycles on T2 and T2 of the Achilles tendon. *Eur J Radiol.* (2014) 83:349-53. doi: 10.1016/j.ejrad.2013.10.014

28. Eyre DR, Muir H. Quantitative analysis of types I and II collagens in human intervertebral discs at various ages. *Biochim Biophys Acta.* (1977) 492:29-42. doi: 10.1016/0005-2795(77)90211-2