Skip to main content
  • Research article
  • Open access
  • Published:

Computer simulation of syringomyelia in dogs

Abstract

Background

Syringomyelia is a pathological condition in which fluid-filled cavities (syringes) form and expand in the spinal cord. Syringomyelia is often linked with obstruction of the craniocervical junction and a Chiari malformation, which is similar in both humans and animals. Some brachycephalic toy breed dogs such as Cavalier King Charles Spaniels (CKCS) are particularly predisposed. The exact mechanism of the formation of syringomyelia is undetermined and consequently with the lack of clinical explanation, engineers and mathematicians have resorted to computer models to identify possible physical mechanisms that can lead to syringes. We developed a computer model of the spinal cavity of a CKCS suffering from a large syrinx. The model was excited at the cranial end to simulate the movement of the cerebrospinal fluid (CSF) and the spinal cord due to the shift of blood volume in the cranium related to the cardiac cycle. To simulate the normal condition, the movement was prescribed to the CSF. To simulate the pathological condition, the movement of CSF was blocked.

Results

For normal conditions the pressure in the SAS was approximately 400 Pa and the same applied to all stress components in the spinal cord. The stress was uniformly distributed along the length of the spinal cord. When the blockage between the cranial and spinal CSF spaces forced the cord to move with the cardiac cycle, shear and axial normal stresses in the cord increased significantly. The sites where the elevated stress was most pronounced coincided with the axial locations where the syringes typically form, but they were at the perimeter rather than in the central portion of the cord. This elevated stress originated from the bending of the cord at the locations where its curvature was high.

Conclusions

The results suggest that it is possible that repetitive stressing of the spinal cord caused by its exaggerated movement could be a cause for the formation of initial syringes. Further consideration of factors such as cord tethering and the difference in mechanical properties of white and grey matter is needed to fully explore this possibility.

Background

Syringomyelia is a severe progressive pathological condition characterised by the presence of abnormal fluid-filled cavities (singular syrinx and plural syringes) within the spinal cord [1,2,3,4,5,6,7,8]. Clinical signs of syringomyelia include pain, loss of sensation, and motor dysfunction [2]. Syringomyelia affects both humans and animals (mainly dogs and cats) [9,10,11]. In humans, syringomyelia has a low prevalence of about one in 10,000 [3, 4]. In contrast, in some toy breed of dogs such as Cavalier King Charles Spaniels (CKCS), King Charles spaniels, Griffon Bruxellois, Maltese, Chihuahuas, and Pomeranians the prevalence is high. In the CKCS up to 70% of the total population is affected by the condition although in many dogs it may be asymptomatic [12]; in the Griffon Bruxellois, 50% of the American population was affected with again many dogs being asymptomatic [11].

Syringomyelia is associated with the obstruction of CSF channels and the subarachnoid space (SAS), particularly if the obstruction is at the foramen magnum [4]. The most common cause of foramen magnum obstruction in humans is Chiari malformation and in animals the analogous condition, Chiari-like malformation, to which brachycephalic toy breed dogs such as Cavalier King Charles Spaniels (CKCS) are particularly predisposed [13]. A common feature of Chiari-like malformation in several toy breeds is occipital bone hypoplasia and reduced caudal fossa volume associated with a compensatory increase in height of the cranial fossa [14]. Traits that increase the risk of syringomyelia in CKCS include brachycephaly with rostral displacement of the atlas and axis, acute angulation between the sphenoid and basioccipital bone, reduced occipital crest and increased cervical flexure and odontoid (dens) angulation [14,15,16,17]. Chiari-like malformation is analogous to Chiari type 1 and 0 malformation in humans which is also associated with syringomyelia [18].

The pathogenesis of syringomyelia is still unknown and this ultimately compromises the effectiveness of clinical (surgical) management of syringomyelia [2, 4, 8, 19]. It is widely believed that variation to the normal flow pathways of the cerebrospinal fluid (CSF) plays the main role in the formation and expansion of syringes [20]. CSF constantly moves in a pulsatile manner between the cranial and spinal SAS via the foramen magnum to compensate for transient shifts in the cranial blood volume associated with the arterial pulse [21, 22]. In Chiari malformation, the obstruction of the foramen magnum compromises the normal communication between the cranial and spinal CSF compartments [3, 23]. Similarly, stenoses in the SAS due to traumatic injury obstruct the flow of CSF [4]. In-vivo and in-vitro experiments [24,25,26,27] indicate that such blockages can result in pressure dissociation in the SAS and in excessive motion of the spinal cord [28]. Typically, syringes form caudally from a blockage. In dogs most affected with syringomyelia, the emergence and progression of cavities has a breed-related pattern. Thus, for example in Cavalier King Charles Spaniels the initial cavities usually first emerge in the cranial-cervical and thoracic-spinal cord segments (C1-C4 and T1-T4, respectively) [29]. The twelfth thoracic-second lumbar (T12-L2) segment is also a common site where syringes may form [29]. By comparison, in some breeds such as Griffon Bruxellois the cavities are often widest in the cervicothoracic junction and cranial thoracic region (unpublished data). In humans the analogous condition Chiari-I malformation is more likely to be associated with a cervical syringe [30].

Early theories on the origins of syringomyelia speculated that syringes form because bulk CSF flow is redirected to the central canal of the spinal cord which is connected with the ventricular system of the brain [1, 25]. However, these theories were abandoned since they could not account for non-caninular syringes (not continuous with the central canal) or syringes emerging in the absence of a patent central canal [20]. Current theories focus mainly on the accumulation of fluid into the spinal cord parenchyma [31,32,33,34,35,36,37] and on excessive stressing of the spinal cord tissue [38,39,40]. Most of the current theories are derived from in-vitro and computer models that can simulate propagation of CSF pulse waves in the spine under various scenarios [31,32,33,34,35, 38,39,40,41,42,43,44,45,46,47,48,49,50]. These models can capture the interaction between elastic structures (spinal cord and dura) and the CSF as well as predict spatial and temporal variation of parameters such as CSF pressure, CSF velocity, and stress in the spinal cord. Cirovic et al. [41, 49], and Bertram et al. [39, 51] used mathematical/computer models to identify four pulse modes in the spinal cord. Three of the four modes are involved in the movement of the CSF and have speeds of the order of several metres per second. More recently Støverud et al. predicted pulse speeds between 2 and 5 m/s using an anatomically realistic computer model of poroelastic spinal cord [52]. The findings from the computer modelling studies are consistent with the pulse speed data obtained from in-vivo human and animal experiments [24, 53, 54]. Elliot [40] further examined the pulse modes and concluded that stenoses in the SAS lead to elevated stress in the cord which is released with the emergence of a syrinx (branding this “a homeostatic mechanism for syringomyelia”). Bertram developed a fluid-structure interaction computer model of the human spinal cavity that included the cord with a syrinx, the SAS fluid compartment, and the elastic dura [38, 39, 45]. The model geometry was idealised to be axisymmetric neglecting the spinal curvature. He concluded that arachnoiditis can lead to excessive stressing of the cord, and that the sloshing of the fluid in an existing syrinx can induce stresses at the caudal end of the syrinx. Recently, Bertram and Heil [32, 35] extended Bertram’s original model to examine the case of a stenosis overlaying a large syrinx for a poro-elastic spinal cord. They concluded that this situation results in a “diastolic valve” that may pump small amounts of fluid from the SAS to the syrinx with each cardiac cycle. While these results are promising, they mainly apply to the enlargement of an existing syrinx or a “presyrinx” (an initial state of increased fluid in the spinal cord [55]) rather than with its emergence. None of the models convincingly explain why syringes emerge at specific anatomical locations even in the absence of arachnoiditis or an initial small “pre-syrinx”.

We hypothesize that preferential locations for the emergence of syringes may be related to the curvature of the spine. This cannot be examined with existing computer models of syringomyelia since they all approximate the spinal cord as a straight cylindrical structure. In order to investigate this possibility, we constructed a three-dimensional computer model of the canine spine. The CKCS spine was chosen for the study due to a high prevalence of syringomyelia in this breed and because of a highly regular pattern in the development of cavities in the spinal cord.

Methods

The dynamics of the canine spinal cord was investigated using a finite element model of a canine spinal cavity. The finite element method is an approximate method for solving engineering problems involving complex geometries and material properties that cannot be solved analytically. In this method, the domain of interest is “discretized” by breaking it into a large number of small geometric components called “elements”. For each element in the model, simplifying approximations are applied to the equations governing the mechanical behaviour of the system examined. The resulting set of equations involving all the elements is solved using computers to yield quantities such as displacement, speed, strain and stress in each element.

The model

The computational model was developed from a CKCS with a Chiari-like malformation and fully developed syrinx (C2-L3; Fig. 1a). The geometry of the spinal canal was extracted from T2-weighted transverse anatomical magnetic resonance imaging (MRI) scans. The imaging was obtained as part of the necessary diagnostic work and consent was obtained from the owner. Due to the size of the subject, two series of transverse scans were compiled to cover the full length of the spine (C1 – S1) at 5 mm intervals with a 10% interslice gap. During the MRI scan the dog was placed in dorsal recumbency with its pelvic limbs extended caudally on the scanning couch, with the spine being as straight as possible to avoid any external distortions of the images. The following T2-weighted imaging parameters were applied: Time to Repetition (TR) 3200, Echo Time (TE) 102, Field of View (FOV) 150, matrix resolution 384 × 320, and averages 5, using combination of spine table integrated coil and torso array coil for maximum signal collection.

Fig. 1
figure 1

a T2-weighted mid-sagittal MRI of a CKCS suffering from a large syrinx (syrinx ends indicated with arrows). b The geometry reconstruction process. 68 transverse plane scans were obtained (1). Each anatomical layer was segmented to create a mask (2). The mask was filtered out (3) and then binarised (4) to create a pure black and white image

In order to reduce the complexity of the model and create a high quality mesh, the components of the spinal canal were assumed to be circular and concentric. Previous studies have demonstrated that the CSF pressures are relatively unaffected by the concentricity and shape of the anatomic layers [56, 57]. The geometry of the spinal canal (C2 - S1) was manually segmented using the Simpleware software package (Simpleware Ltd., Exeter, UK) generating two-dimensional (2D) transverse masks for the following components of the spinal canal; the epidural space (EDS), the subarachnoid space, the spinal cord and the syrinx cavity (Fig. 1b). The 2D masks were exported into ImageJ, an open source image processing platform [58], to calculate the cross-sectional areas and the corresponding radii of each component. The centroids of the spinal cord masks were assumed to be the centres for all components. The reconstructed geometric data is summarised in Table 1 and in the Additional file 1.xlsx. A volumetric model was generated in ANSYS (ANSYS Inc., Pittsburgh USA) by connecting the transverse sections with splines. It was found that it is sufficient to use only 13 sections to construct a smooth three-dimensional geometry that conforms well with the reconstructed 68 transverse sections (see Table 1). The spine was assumed to be symmetric about the mid-sagittal plane and medio-lateral motion was neglected. The cranial end of the model was extended by 15 mm (transparent section in Fig. 2). This was done to remove spurious stress concentrations which would form in the cranial end of spine if the model was loaded directly at the C2 section of the model. Hexagonal “brick” elements of 0.5 mm overall size were used to discretize the model. The element size was chosen so that each layer has at least two elements across the thickness. The finite element model comprised of 903,232 elements in total. Figure 2 detail A shows all the layers in the final model, with the finite element mesh. The dura mater was modelled as a layer of shell elements with 1 mm thickness between the SAS and the EDS.

Table 1 Geometric parameters of the model
Fig. 2
figure 2

The finite element model of the CKCS suffering from a large syrinx. The anatomical layers are indicated in Detail A. The transparent cranial segment was added to the original geometry and the inputs were applied at its cranial end

The mechanical properties for the materials of the layers comprising the model were taken from the published literature [20]; the cord, dura, and the epidural fat were represented as linear elastic solids, and the CSF was represented as an inviscid fluid with the compressibility equal to that of water (Table 2). Treating CSF as inviscid in the first approximation is justified since the pulsatile flow of the CSF in the spine is dominated by inertia rather than viscous effects [56]. From a practical point of view, inviscid approximation was useful since it allowed for the simultaneous solution of solid and fluid equations using a Lagrangian approach. The model was constrained from moving at the outer edge of the EDS to account for the constraint imposed by the vertebrae. Constraint to movement was also imposed to the caudal ends of the spinal cord and the SAS.

Table 2 Material properties

Simulations

An input with amplitude and frequency representative of CSF movement generated by the arterial pulse in the cranium were used to excite the model. This was done by imposing a sinusoidal volumetric flow with a period of 0.4 s (150 beats per minute as the upper limit of heart rate for CKCS [59]) and an amplitude of 0.35 cm3/s. This resulted in a peak CSF speed of about 0.02 m/s which is within a physiologically plausible range [28, 53, 56, 60,61,62]. The volumetric input was applied at the cranial end of the extruded section of the model. Three extreme hypothetical situations were considered by varying the boundary conditions at the cranial end of the extruded section and by excluding or including the syrinx:

  1. 1)

    Cord without syrinx; cranial end of the cord constrained from movement and the CSF allowed to flow freely at the cranial end (representing the normal condition).

  2. 2)

    Cord without syrinx; cranial end of the cord allowed to move freely and the CSF blocked at the cranial end (representing Chiari with the syrinx still not formed).

  3. 3)

    Cord with syrinx; cranial end of the cord allowed to move freely and the CSF blocked at the cranial end (representing Chiari with the syrinx formed).

Prior to applying the above listed scenarios, a preliminary simulation was run in which the model was excited with a short transient pulse of 0.01 s duration and 0.35 cm3/s volumetric flow amplitude. This was done to test if the model can predict the pulse modes reported in the published theoretical and computational studies. All simulations were performed using the LS-DYNA finite element package for dynamic simulations (LS-DYNA Version 971, Livermore Software Technology Corp., Livermore, CA).

Results

Figure 3 shows the distribution of the axial (cranial-to-caudal) component of speed for a short pulse excitation. Only the syrinx, cord, and SAS layers are shown since the movement in the EDS was very small. At 5 ms after the excitation, there was a single disturbance at the cranial end of the model. At 25 ms that disturbance split into three waves that moved at different speeds and covered different distances. For the fastest wave modes, all three layers moved in the same (caudal) direction. For the second fastest mode, the fluids in the syrinx and SAS moved in the opposite direction from the spinal cord. For the slowest mode, the syrinx and the spinal cord moved in the opposite direction from the CSF in the SAS. From the distance covered and the time elapsed, the speeds of the three modes can be estimated to be approximately 1.5 m/s, 7 m/s, and 11 m/s. The speed of movement of the tissues in the anatomical layers of the spine generated by the pulse modes is of the order of 10− 2 m/s (see Fig. 3).

Fig. 3
figure 3

Axial (cranial-to-caudal) speed of tissue movement in the spinal column due to a short transient cranial excitation. a At t = 0.005 s from the onset of excitation an initial disturbance with the amplitude of 0.01 m/s (1 cm/s) formed at the cranial end. b At t = 0.025 s the initial disturbance split into three waves which travelled at speeds of 1.5 m/s (wave 1), 7 m/s (wave 2), and 11 m/s (wave 3). The waves generated tissue movement in the range of ±0.01 m/s

The “cardiac” excitation of the model generated an oscillatory motion of the SAS and/or spinal cord at the cranial end of the model. For the “normal condition” the peak speed of the CSF at the cranial end of the SAS was approximately 0.02 m/s. The amplitude of the speed in the SAS diminished in a linear way from the cranial to the caudal end. For the “Chiari” conditions, the peak speed of the cord at the cranial end was just below 0.01 m/s and the same linear tendency of amplitude decrease towards the caudal end was observed. In all three cases considered, the excitation resulted in harmonic pressure fluctuations in the SAS with a peak amplitude of approximately 400 Pa (about 3 mmHg). The pressure amplitude was approximately constant along the length of the spine; the pressure gradient between the cranial and the caudal ends did not exceed 0.5 mmHg over the cardiac cycle.

Figure 4 shows the distribution of stress in the cord along the axis of the spine for the three inputs. The displayed results were obtained in the following way: for each element comprising the discretized spinal cord, a stress-versus-time curve was obtained and the peak value of stress was determined. Next, the cord was divided into segments of one millimetre width in the cranial-to-caudal direction and the maximal value of peak stress recorded in the elements within each segment was determined. Thus, each point in a curve represents the maximal value of stress within a 1 mm thick slice of the cord over the whole simulated event. Please note that the normal stresses in the cord were entirely negative (compressive) and that the values shown represent the absolute maximal values. The radial stress perpendicular to the sagittal plane is shown in Fig. 4a. The stress was fairly uniform along the length of the cord for all three cases considered and approximately equal to the pressure in the SAS (400 Pa). The same was found for the radial stress in the sagittal plane (not shown in the figure). The radial stress was thus a reflection of the pressure in the SAS. For the normal case, the axial normal stress (cranial-to-caudal direction) showed little variation along the length of the spine and had the same amplitude as the radial stress. In contrast, for the Chiari case without a syrinx, the axial normal stress had two prominent peaks of just under 800 Pa, which were located at the level of C2 and C3 vertebrae and cervico-thoracic junction. Elsewhere, the stress values were about 200 Pa higher than for the normal condition. With the syrinx introduced, the stress levels remained elevated, the peaks were less distinct and there was a further stress increase in the lower thoracic portion of the spinal cord. This was also observed for the shear stress, as can be seen in Fig. 5c. For the normal condition, the shear stress was very low, an order of magnitude lower than the normal stress. For Chiari, the shear stress was considerably elevated with two peaks of approximately 200 Pa at the level of C2 and C3 vertebrae and cervico-thoracic junction. To get a more complete understanding as to the location of maximal stresses, we must look at the three-dimensional distribution. Figure 5 shows the distribution of the peak shear stress in the cranial portion of the spinal cord for the case of Chiari without a syrinx. The results displayed were obtained by assigning to each element the peak value of the shear stress recorded in that element during the simulated event. It can be seen that regions with the highest level of stress were on the ventral side of the cranial-cervical cord and at the dorsal side of the caudal-cervical/cranial thoracic spine. The same type of result was obtained for the normal axial stress. This suggests that regions of localised high stress are due to the bending deformation of the spinal cord.

Fig. 4
figure 4

The distribution of maximal stress in the axial (cranial-to-caudal) direction along the spinal cord. The axial coordinate is measured from C2; the shape of the spinal cord is shown at the top for reference. a Radial normal stress. b Axial normal stress. c Shear stress

Fig. 5
figure 5

The distribution of peak shear stress in the cervical and cranial thoracic part of the spinal cord (Chiari without a syrinx case). The regions in red are the locations of the highest shear stress

Discussion

The results for the short transient excitation show that the model can successfully predict the wave modes of the CSF pulse propagation reported in the open literature. The speeds of the three modes estimated from the results, as well as the associated movement of the anatomical layers are consistent with the theoretical models presented by Cirovic [41, 49]. They are also within the physiologically plausible range determined from in-vivo experiments on humans and animals [53, 54]. As the previous theoretical and computational studies neglected the curvature of the spine, it can be concluded that the curvature has a negligible effect on the propagation of the CSF pulse. The spinal curvature also had little effect on the pressure generated in the SAS; for all three cases considered, the pressure amplitude in the SAS was approximately 400 Pa, and it was uniform along the length of the spine. The predicted small spatial pressure gradient is consistent with findings of other published studies [63]. In other words, the spinal cavity behaved as a lumped compartment whose compliance was chiefly determined by the elasticity of the dural sac. This is because the time scales of the cardiac input are larger than the time needed for the CSF pulse to traverse the entire spinal cord. For all cases considered all the stress components in the cord were of the same order of magnitude as the pressure in the SAS (several hundreds of Pa). For the normal case where the cranial end of the cord was fixed, all the components of the normal stress were compressive and effectively equal to the pressure in the SAS. In other words, the pressure in the SAS generated an equal “hydrostatic pressure” (average normal stress) in the spinal cord. When the CSF was blocked and the cord was allowed to move, the radial normal stress did not change, but the axial (cranio-caudal) stress increased. We attribute this result to the fact that the axial movement of the cord resulted in bending in the sagittal plane. Bending deformation generates normal axial stresses that are largest at the perimeter of the section, and this is consistent with the predictions of the model. It should be noted that the maximal displacement of the cord in the caudal direction was small; about 1 mm in total. The axial normal stress was particularly prominent at the locations where the spinal cord had large curvature and the bending effect was pronounced. The same case was for the shear stress which is also generated by bending. Furthermore, these locations are also the locations where the cavities usually first start forming in the CKCS affected with syringomyelia. When a large syrinx was present, stress concentration moved to the lower thoracic region which is the next most likely site of syrinx formation [29]. This is an interesting result, and to our knowledge this is the first case where computer simulation predicts stress concentration at specific regions of the spinal cord in absence of stenosis and arachnoiditis. Although the stress magnitude at the most affected regions increased by as much as 100% from the normal values, the overall stress levels were generally quite low, and should not be able to cause an immediate damage to the spinal cord tissue. However, it is still feasible that exaggerated stressing repeated with each cardiac pulse does cause pathological changes to the spinal cord tissue over a long period of time. Actually, this is consistent with the fact that syringes take months or years to form [2, 8]. While the axial location of the increased cord stressing is coincident with the sites where the syringes tend to form, the model predicted that the stress was concentrated at the perimeter of the spinal cord rather than at the central portion. This is not in agreement with the fact that the cavities are usually located more centrally. Therefore, as it stands now, the results of this study cannot fully support the assumption that stressing of the cord near the maximum curvature locations may initiate the formation of syringes via direct mechanical damage. However, it should be considered that the model presented in this study does not include numerous factors that could modify the current result. Such factors include non-concentric realistic geometry of the anatomical layers, different material properties of the white and grey matter of the spinal cord, and the effect of denticulate ligaments and tethering in general. It is also possible that the syrinx initiation mechanism is less direct; it has been shown that the damage of microstructure at the surface of the cord may disrupt the blood-spinal cord barrier and affect fluid content of the spinal cord [37, 64, 65], thus leading to fluid accumulation.

Similar to other computer models, the model presented in this paper is based on a number of idealisations. Furthermore, the material properties of the spinal compartments are still not well characterised. We opted for choosing the material properties of the dura and the spinal cord that are most commonly used in the current (human) models of the spinal cavity [32, 35, 38, 39, 45]. Prior to making this choice, the model had been tested with the full range of spinal cord and dura material properties reported in the literature as summarised by Elliot et al. [20]. The tests showed that the dura properties listed in Table 2 lead to the most realistic results in terms of SAS pressure magnitude. When a lower stiffness was assigned to the spinal cord tissue, this resulted in proportionately lower stresses in the cord. However, the general trends remained the same; there was a significant increase in the axial and shear stresses at the upper cervical and thoracic spinal cord for the Chiari condition. This suggests that the conclusions for the current model are largely independent of the material constants, as long as they are within a plausible range. On the other hand, this also means that all the quantitative values of stress reported in this paper should not be taken literally, but rather in reference to the set of material properties used in the current model.

Virtually all computer models of syringomyelia end at the cranial end of the cervical spine or at best extend into the cranio-cervical junction. In this model, an ad hoc extruded cranial segment was added to avoid spurious stress concentrations at the cranial end of the cord. Effectively, this problem arises from the fact that the brain and the cranium are not included in the model. As Chiari deformation results in the brain being “crammed” in the cranium, it is possible that the lower portion of the brain and perhaps cranial section of the cervical spinal cord are under permanent stressing. Therefore, ultimately, a more complete model of syringomyelia should include both the cranial and spinal compartments.

Conclusions

A computer model was used to investigate the effect of excess movement of the spinal cord in dogs suffering from Chiari-like malformations with the goal of identifying possible factors that could promote the emergence of fluid-filled cavities in the spinal cord. It was found that small movements of the spinal cord associated with a cardiac pulse can result in a significant increase in shear and axial normal stresses and that this effect is most pronounced in the cervical and upper thoracic regions where the spinal curvature is high. The regions of high stress as predicted by the model are at the locations where syringes first form, but they are positioned near the perimeter, rather than at the central portion of the cord which is where the cavities usually are. We believe that the latter may be due to simplifying assumptions inherent to the model, and that this mechanism should be further investigated.

References

  1. Williams B. On the pathogenesis of syringomyelia: a review. J R Soc Med. 1980;73(11):798–806.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Levine DN. The pathogenesis of syringomyelia associated with lesions at the foramen magnum: a critical review of existing theories and proposal of a new hypothesis. J Neurol Sci. 2004;220(1–2):3–21.

    Article  PubMed  Google Scholar 

  3. Heiss JD, Patronas N, DeVroom HL, Shawker T, Ennis R, Kammerer W, Eidsath A, Talbot T, Morris J, Eskioglu E, et al. Elucidating the pathophysiology of syringomyelia. J Neurosurg. 1999;91(4):553–62.

    Article  CAS  PubMed  Google Scholar 

  4. Brodbelt AR, Stoodley MA. Post-traumatic syringomyelia: a review. J Clin Neurosci. 2003;10(4):401–8.

    Article  CAS  PubMed  Google Scholar 

  5. Stoodley MA, Jones NR, Yang L, Brown CJ. Mechanisms underlying the formation and enlargement of noncommunicating syringomyelia: experimental studies. Neurosurg Focus. 2000;8(3):E2.

    Article  CAS  PubMed  Google Scholar 

  6. Milhorat TH, Kotzen RM, Anzil AP. Stenosis of central canal of spinal cord in man: incidence and pathological findings in 232 autopsy cases. J Neurosurg. 1994;80(4):716–22.

    Article  CAS  PubMed  Google Scholar 

  7. Milhorat TH, Capocelli AL Jr, Anzil AP, Kotzen RM, Milhorat RH. Pathological basis of spinal cord cavitation in syringomyelia: analysis of 105 autopsy cases. J Neurosurg. 1995;82(5):802–12.

    Article  CAS  PubMed  Google Scholar 

  8. Rusbridge C, Greitz D, Iskandar BJ. Syringomyelia: current concepts in pathogenesis, diagnosis, and treatment. J Vet Intern Med. 2006;20(3):469–79.

    Article  PubMed  Google Scholar 

  9. Rusbridge C, Knowler SP, Pieterse L, McFadyen AK. Chiari-like malformation in the griffon Bruxellois. J Small Anim Pract. 2009;50(8):386–93.

    Article  CAS  PubMed  Google Scholar 

  10. Chandler K, Volk H, Rusbridge C, Jeffery N. Syringomyelia in cavalier king Charles spaniels. Vet Rec. 2008;162(10):324.

    Article  PubMed  Google Scholar 

  11. Freeman AC, Platt SR, Kent M, Huguet E, Rusbridge C, Holmes S. Chiari-like malformation and syringomyelia in American Brussels griffon dogs. J Vet Intern Med. 2014;28(5):1551–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Parker JE, Knowler SP, Rusbridge C, Noorman E, Jeffery ND. Prevalence of asymptomatic syringomyelia in cavalier king Charles spaniels. Vet Rec. 2011;168(25):667.

    Article  CAS  PubMed  Google Scholar 

  13. Marino DJ, Loughin CA, Dewey CW, Marino LJ, Sackman JJ, Lesser ML, Akerman MB. Morphometric features of the craniocervical junction region in dogs with suspected Chiari-like malformation determined by combined use of magnetic resonance imaging and computed tomography. Am J Vet Res. 2012;73(1):105–11.

    Article  PubMed  Google Scholar 

  14. Knowler SP, Kiviranta AM, McFadyen AK, Jokinen TS, La Ragione RM, Rusbridge C. Craniometric analysis of the hindbrain and Craniocervical junction of Chihuahua, Affenpinscher and cavalier king Charles spaniel dogs with and without Syringomyelia secondary to Chiari-like malformation. PLoS One. 2017;12(1):e0169898.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Knowler SP, Cross C, Griffiths S, McFadyen AK, Jovanovik J, Tauro A, Kibar Z, Driver CJ, La Ragione RM, Rusbridge C. Use of morphometric mapping to characterise symptomatic Chiari-like malformation, secondary Syringomyelia and associated Brachycephaly in the cavalier king Charles spaniel. PLoS One. 2017;12(1):e0170315.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Driver CJ, Rusbridge C, McGonnell IM, Volk HA. Morphometric assessment of cranial volumes in age-matched cavalier king Charles spaniels with and without syringomyelia. Vet Rec. 2010;167(25):978–9.

    Article  CAS  PubMed  Google Scholar 

  17. Rusbridge C, Knowler SP. Inheritance of occipital bone hypoplasia (Chiari type I malformation) in cavalier king Charles spaniels. J Vet Intern Med. 2004;18(5):673–8.

    Article  PubMed  Google Scholar 

  18. Rusbridge C. Veterinary aspects. In: Flint G, Rusbridge C, editors. Syringomyelia: a disorder of CSF circulation. Heidelberg, New York, Dordecht, London: Springer; 2014. p. 209–30.

    Chapter  Google Scholar 

  19. Rusbridge C. Chiari-like malformation with syringomyelia in the cavalier king Charles spaniel: long-term outcome after surgical management. Veterinary surgery : VS. 2007;36(5):396–405.

    Article  PubMed  Google Scholar 

  20. Elliott NSJ, Bertram CD, Martin BA, Brodbelt AR. Syringomyelia: a review of the biomechanics. J Fluid Struct. 2013;40:1–24.

    Article  Google Scholar 

  21. Alperin N. MR-intracranial compliance and pressure: a method for noninvasive measurement of important neurophysiologic parameters. Methods Enzymol. 2004;386:323–49.

    Article  PubMed  Google Scholar 

  22. Alperin NJ, Lee SH, Loth F, Raksin PB, Lichtor T. MR-intracranial pressure (ICP): a method to measure intracranial elastance and pressure noninvasively by means of MR imaging: baboon and human study. Radiology. 2000;217(3):877–85.

    Article  CAS  PubMed  Google Scholar 

  23. Martin BA, Kalata W, Shaffer N, Fischer P, Luciano M, Loth F: Hydrodynamic and longitudinalimpedance analysis of cerebrospinal fluid dynamics at the craniovertebral junction in type I Chiari malformation. PloS one 2013;8(10):e75335.

  24. Williams B. Simultaneous cerebral and spinal fluid pressure recordings. I. Technique, physiology, and normal results. Acta Neurochir. 1981;58(3–4):167–85.

    Article  CAS  PubMed  Google Scholar 

  25. Williams B. Cerebrospinal fluid pressure changes in response to coughing. Brain. 1976;99(2):331–46.

    Article  CAS  PubMed  Google Scholar 

  26. Martin BA, Loth F. The influence of coughing on cerebrospinal fluid pressure in an in vitro syringomyelia model with spinal subarachnoid space stenosis. Cerebrospinal Fluid Res. 2009;6:17.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Martin BA, Labuda R, Royston TJ, Oshinski JN, Iskandar B, Loth F. Spinal subarachnoid space pressure measurements in an in vitro spinal stenosis model: implications on syringomyelia theories. J Biomech Eng. 2010;132(11):111007.

    Article  PubMed  Google Scholar 

  28. Vavasour IM, Meyers SM, MacMillan EL, Madler B, Li DK, Rauscher A, Vertinsky T, Venu V, MacKay AL, Curt A. Increased spinal cord movements in cervical spondylotic myelopathy. Spine J. 2014;14(10):2344–54.

    Article  PubMed  Google Scholar 

  29. Loderstedt S, Benigni L, Chandler K, Cardwell JM, Rusbridge C, Lamb CR, Volk HA. Distribution of syringomyelia along the entire spinal cord in clinically affected cavalier king Charles spaniels. Vet J. 2011;190(3):359–63.

    Article  PubMed  Google Scholar 

  30. Strahle J, Muraszko KM, Garton HJ, Smith BW, Starr J, Kapurch JR 2nd, Maher CO. Syrinx location and size according to etiology: identification of Chiari-associated syrinx. J Neurosurg Pediatr. 2015;16(1):21–9.

    Article  PubMed  Google Scholar 

  31. Martin BA, Reymond P, Novy J, Baledent O, Stergiopulos N. A coupled hydrodynamic model of the cardiovascular and cerebrospinal fluid system. Am J Phys Heart Circ Phys. 2012;302(7):H1492–509.

    CAS  Google Scholar 

  32. Heil M, Bertram CD. A poroelastic fluid-structure interaction model of syringomyelia. J Fluid Mech. 2016;809:360–89.

    Article  Google Scholar 

  33. Bilston LE, Stoodley MA, Fletcher DF. The influence of the relative timing of arterial and subarachnoid space pulse waves on spinal perivascular cerebrospinal fluid flow as a possible factor in syrinx development. J Neurosurg. 2010;112(4):808–13.

    Article  PubMed  Google Scholar 

  34. Bilston LE, Fletcher DF, Brodbelt AR, Stoodley MA. Arterial pulsation-driven cerebrospinal fluid flow in the perivascular space: a computational model. Comput Methods Biomech Biomed Engin. 2003;6(4):235–41.

    Article  PubMed  Google Scholar 

  35. Bertram CD, Heil M, Poroelastic Fluid A. Structure-interaction model of cerebrospinal fluid dynamics in the cord with Syringomyelia and adjacent subarachnoid-space stenosis. J Biomech Eng. 2017;139(1)

  36. Kim M, Cirovic S. A computational model of the cerebrospinal fluid system incorporating lumped-parameter cranial compartment and one-dimensional distributed spinal compartment. J Biorheol. 2011;25(1):78–87.

    Article  Google Scholar 

  37. Hemley SJ, Biotech B, Tu J, Stoodley MA. Role of the blood-spinal cord-barrier in posttraumatic syringomyelia. Laboratory investigation J Neurosurg-Spine. 2009;11(6):696–704.

    Article  Google Scholar 

  38. Bertram CD, Bilston LE, Stoodley MA. Tensile radial stress in the spinal cord related to arachnoiditis or tethering: a numerical model. Med Biol Eng Comput. 2008;46(7):701–7.

    Article  CAS  PubMed  Google Scholar 

  39. Bertram CD. A numerical investigation of waves propagating in the spinal cord and subarachnoid space in the presence of a syrinx. J Fluid Struct. 2009;25:1189–205.

    Article  Google Scholar 

  40. Elliott NSJ, Lucey AD, Lockerby DA, Brodbelt AR. Fluid-structure interactions in a cylindrical layered wave guide with application in the spinal column to syringomyelia. J Fluid Struct. 2017;70:464–99.

    Article  Google Scholar 

  41. Cirovic S. A coaxial tube model of the cerebrospinal fluid pulse propagation in the spinal column. J Biomech Eng. 2009;131(2):021008.

    Article  PubMed  Google Scholar 

  42. Chang HS, Nakagawa H. Hypothesis on the pathophysiology of syringomyelia based on simulation of cerebrospinal fluid dynamics. J Neurol Neurosurg Psychiatry. 2003;74(3):344–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Carpenter PW, Berkouk K, Lucey AD. Pressure wave propagation in fluid-filled co-axial elastic tubes. Part 2: mechanisms for the pathogenesis of syringomyelia. J Biomech Eng. 2003;125(6):857–63.

    Article  CAS  PubMed  Google Scholar 

  44. Martin BA, Kalata W, Loth F, Royston TJ, Oshinski JN: Syringomyelia hydrodynamics: an in vitro study based on in vivo measurements. J Biomech Eng. 2005;127(7):1110–1120.

  45. Bertram CD. Evaluation by fluid/structure-interaction spinal-cord simulation of the effects of subarachnoid-space stenosis on an adjacent syrinx. J Biomech Eng. 2010;132(6):061009.

    Article  CAS  PubMed  Google Scholar 

  46. Berkouk K, Carpenter PW, Lucey AD. Pressure wave propagation in fluid-filled co-axial elastic tubes. Part 1: basic theory. J Biomech Eng. 2003;125(6):852–6.

    Article  CAS  PubMed  Google Scholar 

  47. Elliott NS, Lockerby DA, Brodbelt AR. A lumped-parameter model of the cerebrospinal system for investigating arterial-driven flow in posttraumatic syringomyelia. Med Eng Phys. 2011;33(7):874–82.

    Article  CAS  PubMed  Google Scholar 

  48. Elliott NS. Syrinx fluid transport: modeling pressure-wave-induced flux across the spinal pial membrane. J Biomech Eng. 2012;134(3):031006.

    Article  CAS  PubMed  Google Scholar 

  49. Cirovic S, Kim M. A one-dimensional model of the spinal cerebrospinal-fluid compartment. J Biomech Eng. 2012;134(2):021005.

    Article  PubMed  Google Scholar 

  50. Clarke EC, Fletcher DF, Stoodley MA, Bilston LE. Computational fluid dynamics modelling of cerebrospinal fluid pressure in Chiari malformation and syringomyelia. J Biomech. 2013;46(11):1801–9.

    Article  PubMed  Google Scholar 

  51. Bertram CD, Brodbelt AR, Stoodley MA. The origins of syringomyelia: numerical models of fluid/structure interactions in the spinal cord. J Biomech Eng-T Asme. 2005;127(7):1099–109.

    Article  CAS  Google Scholar 

  52. Stoverud KH, Alnaes M, Langtangen HP, Haughton V, Mardal KA. Poro-elastic modeling of Syringomyelia - a systematic study of the effects of pia mater, central canal, median fissure, white and gray matter on pressure wave propagation and fluid movement within the cervical spinal cord. Comput Methods Biomech Biomed Engin. 2016;19(6):686–98.

    Article  PubMed  Google Scholar 

  53. Khani M, Xing T, Gibbs C, Oshinski JN, Stewart GR, Zeller JR, Martin BA. Nonuniform moving boundary method for computational fluid dynamics simulation of Intrathecal cerebrospinal flow distribution in a Cynomolgus monkey. J Biomech Eng. 2017;139(8)

  54. Kalata W, Martin BA, Oshinski JN, Jerosch-Herold M, Royston TJ, Loth F. MR measurement of cerebrospinal fluid velocity wave speed in the spinal canal. IEEE Trans Biomed Eng. 2009;56(6):1765–8.

    Article  PubMed  Google Scholar 

  55. Fischbein NJ, Dillon WP, Cobbs C, Weinstein PR. The "presyrinx" state: is there a reversible myelopathic condition that may precede syringomyelia? Neurosurg Focus. 2000;8(3):E4.

    Article  CAS  PubMed  Google Scholar 

  56. Loth F, Yardimci MA, Alperin N. Hydrodynamic modeling of cerebrospinal fluid motion within the spinal cavity. J Biomech Eng. 2001;123(1):71–9.

    CAS  PubMed  Google Scholar 

  57. Heidari Pahlavian S, Yiallourou T, Tubbs RS, Bunck AC, Loth F, Goodin M, Raisee M, Martin BA. The impact of spinal cord nerve roots and denticulate ligaments on cerebrospinal fluid dynamics in the cervical spine. PLoS One. 2014;9(4):e91888.

    Article  PubMed  PubMed Central  Google Scholar 

  58. Schneider CA, Rasband WS, Eliceiri KW. NIH image to ImageJ: 25 years of image analysis. Nat Methods. 2012;9(7):671–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Hezzell MJ, Dennis SG, Humm K, Agee L, Boswood A. Relationships between heart rate and age, bodyweight and breed in 10,849 dogs. J Small Anim Pract. 2013;54(6):318–24.

    Article  CAS  PubMed  Google Scholar 

  60. Cerda-Gonzalez S, Olby NJ, Broadstone R, McCullough S, Osborne JA. Characteristics of cerebrospinal fluid flow in cavalier king Charles spaniels analyzed using phase velocity cine magnetic resonance imaging. Vet Radiol Ultrasound. 2009;50(5):467–76.

    Article  PubMed  Google Scholar 

  61. Yildiz S, Thyagaraj S, Jin N, Zhong X, Heidari Pahlavian S, Martin BA, Loth F, Oshinski J, Sabra KG. Quantifying the influence of respiration and cardiac pulsations on cerebrospinal fluid dynamics using real-time phase-contrast MRI. J Magn Reson Imaging. 2017;46(2):431–9.

    Article  PubMed  Google Scholar 

  62. Martin BA, Sass L, Khani M, Gibbs C, Freeman T, Pluid J, Elliot A, Oshinski JN, Zeller J, Stewart GR et al: Are Monkeys Like Humans? Comparison of Intrathecal CSF Dynamics across Mammalian Species. 4th International Cerebrospinal Fluid Dynamics Symposium 2017 https://csfdynamics.org/atlanta2017/scientific-program/ Accessed 24 July 2017.

  63. Fric R, Lindstrom EK, Ringstad GA, Mardal KA, Eide PK. The association between the pulse pressure gradient at the cranio-cervical junction derived from phase-contrast magnetic resonance imaging and invasively measured pulsatile intracranial pressure in symptomatic patients with Chiari malformation type 1. Acta Neurochir. 2016;158(12):2295–304.

    Article  PubMed  Google Scholar 

  64. Perdiki M, Farooque M, Holtz A, Li GL, Olsson Y. Expression of endothelial barrier antigen immunoreactivity in blood vessels following compression trauma to rat spinal cord. Temporal evolution and relation to the degree of the impact Acta Neuropathol. 1998;96(1):8–12.

    CAS  PubMed  Google Scholar 

  65. Whetstone WD, Hsu JY, Eisenberg M, Werb Z, Noble-Haeusslein LJ. Blood-spinal cord barrier after spinal cord injury: relation to revascularization and wound healing. J Neurosci Res. 2003;74(2):227–39.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We are grateful to the staff at Fitzpatrick Referrals Orthopaedic and Neurology Service for facilitating the MRI imaging and ensuring good patient care and also to the owners and dogs that contributed to our understanding of this condition.

Funding

RL was supported by EPSRC Vacation Bursary in the initial stages of this project (June–August 2014). The bursary was a part of EP/L505092 DTG Training Grant and it covered RL’s salary and consumables.

Fitzpatrick Referrals Ltd., provided support in the form of salaries and materials for authors JJ, and CR. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Availability of data and materials

All data generated or analysed during this study are included in this published article [and its additional files].

Author information

Authors and Affiliations

Authors

Contributions

SC and CR were responsible for the study concept, supervision and funding acquisition. SC and CR and RL were responsible for the methodology. The investigation was performed by RL and SC. Formal analysis was the responsibility of RL and SC. Resources were provided by JJ, CR and HV. SC prepared the original draft which was reviewed and edited by CR, JJ, HV, and RL. SC prepared the figures. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Srdjan Cirovic.

Ethics declarations

Ethics approval

The MRI of the spinal cord of a CKCS was obtained as part of the necessary diagnostic work for an animal undergoing veterinary treatment. Consent was obtained from the owner and the actual identity of the dog remained anonymous.

Consent for publication

Not applicable

Competing interests

RL was as undergraduate student at the University of Surrey where CR and SC are employed. JJ, was and CR is employed by Fitzpatrick Referrals Ltd., Surrey, GU7 Q22. HJV is employed by the Royal Veterinary College, University of London.

The Universities of Surrey and London and Fitzpatrick Referrals did not play a role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript and only provided financial support in the form of authors’ salaries and/or research materials. None of the authors have personal or financial relationships with other people or organizations that might inappropriately influence or bias the content of the paper. There are no patents, products in development, or marketed products to declare.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional file

Additional file 1:

Geometric parameters of the model.ᅟ(XLSX 18 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Cirovic, S., Lloyd, R., Jovanovik, J. et al. Computer simulation of syringomyelia in dogs. BMC Vet Res 14, 82 (2018). https://doi.org/10.1186/s12917-018-1410-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12917-018-1410-7

Keywords