Skip to main content

Ex-vivo quantification of ovine pia arachnoid complex biomechanical properties under uniaxial tension

Abstract

Background

The pia arachnoid complex (PAC) is a cerebrospinal fluid-filled tissue conglomerate that surrounds the brain and spinal cord. Pia mater adheres directly to the surface of the brain while the arachnoid mater adheres to the deep surface of the dura mater. Collagen fibers, known as subarachnoid trabeculae (SAT) fibers, and microvascular structure lie intermediately to the pia and arachnoid meninges. Due to its structural role, alterations to the biomechanical properties of the PAC may change surface stress loading in traumatic brain injury (TBI) caused by sub-concussive hits. The aim of this study was to quantify the mechanical and morphological properties of ovine PAC.

Methods

Ovine brain samples (n = 10) were removed from the skull and tissue was harvested within 30 min post-mortem. To access the PAC, ovine skulls were split medially from the occipital region down the nasal bone on the superior and inferior aspects of the skull. A template was used to remove arachnoid samples from the left and right sides of the frontal and occipital regions of the brain. 10 ex-vivo samples were tested with uniaxial tension at 2 mm s−1, average strain rate of 0.59 s−1, until failure at < 5 h post extraction. The force and displacement data were acquired at 100 Hz. PAC tissue collagen fiber microstructure was characterized using second-harmonic generation (SHG) imaging on a subset of n = 4 stained tissue samples. To differentiate transverse blood vessels from SAT by visualization of cell nuclei and endothelial cells, samples were stained with DAPI and anti-von Willebrand Factor, respectively. The Mooney-Rivlin model for average stress–strain curve fit was used to model PAC material properties.

Results

The elastic modulus, ultimate stress, and ultimate strain were found to be 7.7 ± 3.0, 2.7 ± 0.76 MPa, and 0.60 ± 0.13, respectively. No statistical significance was found across brain dissection locations in terms of biomechanical properties. SHG images were post-processed to obtain average SAT fiber intersection density, concentration, porosity, tortuosity, segment length, orientation, radial counts, and diameter as 0.23, 26.14, 73.86%, 1.07 ± 0.28, 17.33 ± 15.25 µm, 84.66 ± 49.18°, 8.15%, 3.46 ± 1.62 µm, respectively.

Conclusion

For the sizes, strain, and strain rates tested, our results suggest that ovine PAC mechanical behavior is isotropic, and that the Mooney-Rivlin model is an appropriate curve-fitting constitutive equation for obtaining material parameters of PAC tissues.

Introduction

The pia-arachnoid complex (PAC), commonly known as the leptomeninges of the brain, provides mechanical stability for the brain and spinal cord in all vertebrate species [1, 2]. The PAC lies deep to the dura mater and superficial to the brain, and is comprised of the pia mater, arachnoid mater and the subarachnoid space (SAS). The SAS is filled with cerebrospinal fluid (CSF) and contains subarachnoid trabeculae (SAT) which span the intermediate area of the meninges (Fig. 1). SAT are collagen fibers that develop as columns or sheets spanning the SAS, with lateral and transverse orientations relative to the pia and arachnoid mater [1, 3]. The unique PAC morphology suggests that one of its potential functions is to protect the brain from traumatic brain injury (TBI) through energy distribution and absorption [4].

Fig. 1
figure1

Anatomical drawing of pia arachnoid complex showing the a dura mater, b arachnoid mater, c subarachnoid space containing the arachnoid trabeculae fibers and blood vessels, d pia mater, and e brain tissue

Existing studies have identified the significance of the physiological and biomechanical properties of PAC in central nervous system (CNS) pathologies, such as TBI, spinal canal stenosis, and Space Flight Associated Neuro-Ocular syndrome (also called SANS), as well as in applications such as CSF-based drug delivery [4,5,6,7]. Neurosurgical dissection of spinal arachnoiditis, or inflammation of the PAC surrounding the spinal cord, is a procedure believed to reduce CSF flow obstruction by removing SAT fibers [7,8,9,10,11,12]. Gottschalk et al. found in-vivo flow alterations due to spinal arachnoiditis [13]. The relevance of SAT fibers has also been investigated in the context of their potential role in CSF drug delivery and solute transport [14,15,16]. Additionally, Killer et al. showed a potential for increased intraocular pressure caused by SAT fibers and by extension this could potentially affect astronauts with space flight associated neuro-ocular syndrome [17]. Taken together, available data implicates the PAC and SAT in CNS pathologies. Further investigations of the PAC are needed to improve both treatment options and preventative care.

Increased knowledge of PAC biomechanical properties may improve understanding of CNS pathology following trauma and enhance potential post-trauma treatments. Currently, a limited understanding of the mechanical behavior of the PAC and SAT hinders the development of accurate pathophysiological models of the CNS. A common method of modeling mechanical responses in the cranial cavity during TBI is finite element analysis. In order to increase model accuracy, biomechanical properties and microanatomy of the PAC need to be quantified. Previous studies have excluded or used assumed PAC biomechanical properties in computational fluid dynamic-models to investigate CSF dynamics [14, 15, 18]. SAT fibers have been found to potentially increase CSF mixing and pressure gradients throughout the CNS [14, 16, 19], and prior studies have quantified bovine PAC properties under uniaxial tension, traction, and shear [20,21,22]. However, more data on PAC across species are needed to show inter-species differentiation and elucidate how PAC biomechanical properties may affect transient stress loading on the brain during TBI.

The objective of the present study was to quantify nonlinear elastic properties under uniaxial tension and derive the appropriate material parameters of fresh ovine PAC, using the Mooney-Rivlin curve-fit model. Additionally, we aimed to characterize SAT structure and orientation through second-harmonic generation (SHG) imaging, with an automated collagen fiber structural analysis. Results of this study improve the understanding of the structure and biomechanical function of the PAC across species, which will advance finite element models and potentially lead to improved understanding and treatment of PAC-related pathologies.

Methods

PAC tissue collection

Ten ovine brains from approximately 1-year-old animals were harvested from USDA inspected sheep at the University of Idaho Vandal Meats facility. The animals were euthanized using a captive bolt and immediately exsanguinated and decapitated. The captive bolt was applied on the forehead, slightly above a line drawn between the eyes. All samples were visually inspected for brain tissue damage and prepared for brain removal at < 10 min (min) post-mortem. Brain samples were extracted within 30 min post-mortem to minimize potential PAC breakdown. To access the PAC, the lower mandible was removed by making clean proximal cuts through the mastoid muscle, connective tissues, and temporomandibular joints on either side of the mouth (Fig. 2a). A strip of skin ~ 6 cm wide was removed from the medial portion of the skull to expose a clear strip of bone to split. The skull was then flipped upside down to expose the inferior portion of the skull. A steel wedge and hammer were used to split through the hard and soft palate of the mouth. The head was flipped again and split from the nasal bone along the exposed strip of bone to the superior occipital region. The skull was flipped upside down again to split the inferior portion of the occipital region, and all connective tissue was cut. During this procedure, the skull was examined to ensure that the PAC tissue remained intact. The skull was then separated carefully on the medial cut and the brain was removed in its entirety.

Fig. 2
figure2

Procedure for dissection of the ovine PAC tissue. a Initial midline incision is made on superior and inferior aspect of the skull with lower mandible removes. Cuts are made sequentially (1–4). b Location of dissection of occipital and frontal PAC tissue sample removal and visual image of dissected ovine brain

The brain was kept moist post-harvest with artificial CSF (EcoCyte Bioscience, Austin, Texas USA) while four samples were removed from the left and right occipital and frontal lobes (Fig. 2b). Each sample was measured with a clear rectangular template measuring 25 mm long, 10 mm wide, and 3 mm thick. A sample thickness of approximately 3 mm was chosen to include the PAC layer and underlying brain tissue. Extraction locations varied slightly so that each specimen would contain at least one sulcus or gyrus, and to avoid any obvious physical damage from the captive bolt traumatic injury site. Samples were visually inspected for tissue damage and either submerged in artificial CSF for biomechanical testing or fixed in 4% paraformaldehyde overnight for imaging. Fixed samples were washed 3 times for 5 min each in phosphate buffered saline (PBS) and stored in PBS at 4 °C and imaged at a later time. Biomechanical tests were performed within 5 h post-mortem, using samples submerged in artificial CSF immediately following dissection.

PAC biomechanical testing preparation

Immediately prior to biomechanical testing, individual brain samples were placed in a 49-mm diameter petri dish with the arachnoid mater facing down and continually moistened with artificial CSF. Room temperature (22 °C) was maintained throughout uniaxial tension tests. Each sample cut to ~ 15 mm long × 5 mm wide and placed on a c-shaped template with a 10-mm long window (Fig. 3a). Each end was secured in this c-clamp using cyanoacrylate glue. The soft brain matter was gently teased away using a scalpel and flat-head forceps until the translucent PAC was isolated. The dissected tissue was visually confirmed to be the PAC using a dissecting microscope, and tissue integrity was assessed before the sample was considered testable.

Fig. 3
figure3

a PAC with brain matter removed and placed in C-shaped template used for fixture into the load cell. b Bioreactor apparatus with mounted sample within grips and orientation of load cell unixial tension

PAC thickness measurement

To avoid potentially damaging testable PAC samples, a separate set of brain samples (n = 10) were used for an average PAC thickness measurement. A caliper (Model 293-340-30 Mitutoyo, Japan), equipped with a ratchet thimble, was used to measure the thickness of each PAC sample three times to obtain an average sample thickness. The sample was gently placed flat on the lower caliper clamp. The caliper was slowly tightened until the rachet thimble turned over and a visual examination showed the caliper clamps were in full contact of both side of the PAC. The average PAC thickness was calculated and used for subsequent calculations related to biomechanical tests.

Uniaxial tensile testing

A custom small-scale uniaxial tensile load frame with a 150 g maximum capacity load cell (Model 31, Honeywell, Columbus, OH) was used to evaluate biomechanical properties of each PAC sample [23]. PAC samples were mounted and secured into custom soft tissue grips, and the paper c-clamp frame was cut to ensure that only the PAC tissue was loaded. Once mounted, width and length of the sample were measured using ImageJ (NIH, Bethesda, MD). PAC samples were pulled-to-failure in tension at a constant rate of 2 mm s−1 (Fig. 3b) while a LabVIEW program (National Instruments, Austin TX) recorded force and displacement data at 100 Hz. A total of 22 ovine PAC samples were tested under uniaxial tension from the frontal and occipital lobes. Samples were discarded if they failed at the grip (n = 6), a large blood vessel was present within the sample (n = 2), the cross-section was non-uniform (n = 1), or the increasing stress–strain curve was non-monotonic (n = 3). All 10 samples included in the analysis failed at the mid-substance and had a uniform cross-section.

Mooney-Rivlin curve fitting and material parameter estimation

The force–displacement data collected from the experiments were analyzed using a custom MATLAB script (MathWorks, Natick, MA). A Mooney-Rivlin curve fit model (Eq. 1) for uniaxial tension was used to fit the obtained average stress–strain curve (Fig. 4), with stretch ratio \(\lambda =\varepsilon +1\). Mooney-Rivlin constants,\({C}_{10},{C}_{01 },{C}_{20}\), were estimated by a nonlinear least-squares fit using the MATLAB curve fit tool (Vers. R2018B).

Fig. 4
figure4

Stress–strain curves of the ovine PAC specimens when subjected to uniaxial tension. Average curve shown in dark blue

$${\sigma }_{Uniaxial}^{Eng}=2{C}_{10}\left(\lambda -\frac{1}{{\lambda }^{2}}\right)+2{C}_{01}\left(\lambda -\frac{1}{{\lambda }^{3}}\right)+4{C}_{20}\left(\lambda -\frac{1}{{\lambda }^{2}}\right)\left({\lambda }^{2}+\frac{2}{\lambda }-3\right)$$
(1)

Young’s modulus (E) was calculated by tensile stress (\(\sigma\)) over engineering extensional strain (\(\varepsilon\)) over the linear region (Eq. 2). Tensile stress (\(\sigma\)) is the force exerted on the sample F divided by the cross-sectional area of the sample A. Engineering strain (ε) is the change in sample length ΔL divided by the initial length L0 of the sample (Eq. 3). The strain rate was calculated as the strain over the time the force was applied. The stretch ratio used was \(\lambda\) =\(\varepsilon +1\). MATLAB was used to calculate Young’s modulus and ultimate stress and strain were ascertained as the highest values before PAC rupture.

$$E=\frac{\sigma }{\varepsilon }=\frac{F{L}_{0}}{A\Delta L}$$
(2)
$$\varepsilon =\frac{\Delta L}{{L}_{0}}$$
(3)

Second Harmonic Generation Imaging (SHG) and immunofluorescence

SHG imaging was used to visualize the collagen fiber microstructure in PAC samples. Images were taken at 20 × magnification and 860 nm excitation wavelength using a confocal/multiphoton microscope (Fluoview 1000, Olympus, Valley Center, PA) (Fig. 5a). To visualize microvascular structure and cell nuclei, the samples were fixed in 4% paraformaldehyde for 24 h and washed 3 × with PBS + 0.1% Triton X-100 (PBST). The samples were then blocked with 10% bovine serum albumin in PBST for 1-h. To stain blood vessels, samples were incubated with anti-von Willebrand Factor (ab6994, Abcam, Cambridge MA) primary antibody at 1:1000 dilution ratio in 5% bovine serum albumin in PBST on a rocker overnight at 4 °C (Fig. 5b). Samples were then washed 3 × for 5 min in PBST and incubated with anti-goat Alexa Fluoro 488 secondary antibody (1:2000 dilution ratio) for 1 h at room temperature. Cell nuclei were co-stained with 4,6-diamidino-2-phenylindole (DAPI) (Invitrogen). Samples were washed 3 times for 5 min each in PBST and imaged with the confocal microscope (Fig. 5c).

Fig. 5
figure5

a Second Harmonic Generation (SHG) imaging of the PAC at 860 nm. b Anti-Von Willebrand Factor applied at a concentration by volume of 1:1000 to view the membrane lining of endothelial cells in blood vessels. c DAPI applied at concentration by volume of 1:1000 to view cell nuclei

An automated image-based collagen fiber detection method for use with multi-photon microscopy [24] was used to derive fiber orientation and structure in multiple planes (Fig. 6). Slice thickness for z-direction stacks was chosen to be 2.32 µm. This method quantified fiber intersection density, concentration, porosity, tortuosity, segment length, orientation, radial counts, and diameter. These parameters resulted in imaged samples (n = 4) having approximately 20 to 50 slices.

Fig. 6
figure6

Single image of the SHG image stack from fiber analysis code developed by Koch et al. Blue lines depict fiber analysis result (skeleton) used to compute fiber-related parameters (see methods) [24]

Statistical analysis

Descriptive statistics were obtained for stress–strain parameters resulting from different PAC locations (occipital and frontal, and left and right hemispheres). To determine if regional differences existed between tissue sample locations, a two-sample independent t-test with unequal variance was applied using MATLAB. The parameters tested were Young’s Modulus, ultimate stress, and ultimate strain. Differences were considered significant at p-value < 0.05.

Results

PAC mechanical properties

Testing was performed less than five hours post-mortem with an average strain rate of 0.59 ± 0.12 s−1. A two-sample independent t-test applied to the occipital and frontal lobes, and left and right hemispheres showed no statistically significant differences between regions for Young’s modulus, ultimate stress, or ultimate strain (p > 0.05) (Table 1). Thus, all samples were pooled to find average length, width, thickness, Young’s modulus, ultimate stress, and ultimate strain of 4.06 ± 0.61, 3.57 ± 0.90 mm, 70. ± 30. µm, 7.7 ± 3.0, 2.7 ± 0.13 MPa, and 0.60 ± 0.14, respectively. The isotropic material parameters in the Mooney-Rivlin curve fit, C10, C01, and C20 (Eq. 1) were calculated to be 1, −1.004 and 0.629 MPa, respectively (SSE = 0.6552, R2 = 0.985), with 95% confidence bounds.

Table 1 Review of previously published literature documenting PAC mechanical properties and comparison of results to the present study

SHG imaging

SHG imaging allowed for visualization of collagen microstructure in terms of sheets, pillars, and blood vessel walls in a subset of PAC samples (n = 4, Fig. 5a). Anti-von Willebrand Factor and DAPI staining confirmed the presence of endothelial cells and cell nuclei throughout the PAC samples (Fig. 5b, c). This process helped differentiate SAT from transverse blood vessels. SAT fibers showed straight and crimp-like morphology (Fig. 5a and Fig. 6). Using fiber tracing and analysis software (Fig. 6), the average SAT fiber intersection density, concentration, porosity, tortuosity, segment length, orientation, radial counts, and diameter were found to be 0.23, 26.14, 73.86%, 1.07 ± 0.28, 17.33 ± 15.25 µm, 84.66 ± 49.18°, 8.15%, and 3.46 ± 1.62 µm, respectively.

Discussion

This study quantified the mechanical properties of the PAC, with the goal of improving simulations of TBI and other CNS pathologies. The protocol for material testing and the use of the Mooney-Rivlin model was similar to previously published studies. However, to our knowledge, this study represents the first biomechanical and morphological quantification of ovine PAC tissue under uniaxial tension. Methods for brain extraction/dissection and uniaxial tension experiments were closely followed according to previous studies. Fresh ovine PAC samples (n = 10) were kept on ice and submerged in artificial CSF between dissection and testing locations. Samples were removed from ice, taken out of their respective containers and kept moist at room temperature, and tested within 5 h post-extraction. SHG imaging and fiber analysis of the collagen microstructure within the PAC indicated an isotropic collagen structure (for example shown by a wide range of variation in fiber orientation: 84.66 ± 49.18°). Taken together, these findings provide normative values for PAC material properties and structure that can be used for modeling and understanding disease states.

Comparison of biomechanical results to previous studies

PAC biomechanical properties under uniaxial tension have been found to vary from ~ 7 to 65 MPa across a range of species and testing protocols (Table 1 and Fig. 7) [22, 25, 26]. In our study, Young’s modulus was found to be at the lower range of those values at 7.7 ± 3.0 MPa. Prior work using a testing methodology most comparable to our study identified the Young’s modulus of bovine PAC to be 7.52 MPa, at a strain rate of 0.5 s−1 [22]. Another study in bovine spinal cord estimated the spinal PAC to be 17 MPa at a strain rate of 0.05 s−1 [25]. Here, biomechanical differences across the four tested regions of the brain were found to be similar, with no statistically significant differences detected in Young’s modulus, ultimate stress and strain, or thickness. These findings are consistent with prior results in bovine models [21, 25].

Fig. 7
figure7

Relation of Young’s Modulus and strain rate computed for previous studies published in the literature based on isolated PAC in the spine and brain across bovine, ovine, porcine, and rabbit

To our knowledge, human PAC biomechanical properties have not been documented, and therefore a direct comparison with our results is not possible. While studies using human tissues are limited, the volume fraction in post-mortem human SAT was found to be 5 to 10% greater in frontal regions of the brain compared to other regions [27]. In principle, these differences in volume fraction could lead to relatively small regional differences in mechanical properties, but these relatively small regional differences were not possible to confirm in the present study. Future work may confirm or deny the relationship between region-dependent SAT volume fraction and stress-loading during TBI.

Ours is also one of the first studies to measure brain PAC directly in “fresh” unfrozen samples. The measured PAC thickness averaged 70 ± 30 µm and was measured using a micrometer with a thimble compression spring. A prior study used a cryostat sliced sample embedded in cow cheek muscles and an optical microscope to measure average bovine brain PAC thickness as 23.6 ± 5.8 µm [22]. While existing data in brain PAC are limited, prior studies have found spinal PAC thickness to vary from roughly 10–300 µm, depending on the measurement method [21, 25, 26, 28]. Using a camera image obtained from the side of the sample, female ovine spinal cord PAC was found to have a thickness of 200 ± 40 µm [28]. Finally, the thickness of rabbit spinal cord PAC has been estimated at 12 ± 3 µm [26]. Collectively, variations in thickness across studies likely result from the different methodologies, species, and regions analyzed. These variations underscore the importance of generating unique PAC material parameters for simulations based on the location of the simulated trauma.

PAC material model

Similar to previous studies investigating CNS tissue biomechanics [29,30,31] a Mooney-Rivlin model was applied to estimate material constants that were based on experimentally generated stress–strain curves. The present study found the Ovine PAC Mooney-Rivlin material parameters C10, C01, and C20 (Eq. 1) to be 1, −1.004 and 0.629 MPa, respectively. The parameters C10, C01, and C20 are material constants related to the distortional response of the tissue, which are determined according to curve-fitting accuracy requirements in this study. However, the constitutive relationship described by Eq. 1 is phenomenological. The use of a phenomenological constitutive relationship has the characteristic that its material constants do not bear a clear physiological meaning. Jin et al. also used the Mooney-Rivlin model based on Bovine PAC under uniaxial tension and found constants to be C31 = 1.3 (kPa), C32 = 56 (kPa), C41 = 27.5 (kPa), C42 = 7.9 (kPa), C51 = 16.55 (MPa), C52 = 8.45 (MPa) [32]. Other studies have used similar methods regarding the mathematical model for biological tissues [30, 33,34,35]. The Mooney-Rivlin model was also chosen since previous research [36] has found the relation between shear stress and shear strain to be linear for brain tissue. However, this assumption may not be correct for PAC tissue that is separate from the brain parenchyma. A constituent-based stress–strain relationship could help account for the main wall constituents and their structural properties [37], as described by the image-based fiber detection method in our study, giving thus a clear physiological meaning to the constants defining the elastic properties of the tissue.

Structural morphology

This study utilized a novel method for SAT morphology quantification. SHG imaging with an 860 nm wavelength was used to visualize and analyze the type I collagen SAT structure. As the majority of cranial SAT are type I collagen [38], SHG is an effective tool for producing high fidelity images of the SAT structure. Additionally, using these images, SAT fiber intersection density, concentration, porosity, tortuosity, segment length, orientation, radial counts, and diameter were characterized.

SAT fiber morphology has traditionally been studied with scanning electron microscopy, transmission electron microscopy, and brightfield microscopy. Scanning electron microscopy and transmission electron microscopy showed that the structural morphology of human (post-mortem) SAT consisted of pillars, columns, sheets, branched fibrils, and other complex structures with fibril diameters that varied from 0.5–3 µm [3]. Optical coherence tomography was utilized to show that SAT fiber diameter varies from 19.2–45.5 µm [27]. Fiber diameter in the human (post-mortem) bulbar subarachnoid space was also found to be between 0.2–1 µm [39]. The present study used SHG imaging in concert with a program developed by Koch et al. [24] and found fiber diameter of ovine brain PAC to be 3.5 ± 1.6 µm. Fibril diameter variation is likely due to different species or varied anatomical sampling locations. Few studies have explored the linkage between morphology and biomechanical properties. Fabris et al., performed atomic force microscopy indentation and immunofluorescent staining on ex-vivo rat PAC and found significantly different stiffness with respect to reflecting trabecular density [40]. They also correlated stiffness with increased vascularization and vimentin density. The linkage between the image-based morphological characteristics and their specific impact on biomechanical properties should be further explored.

Biological tissues typically have a non-linear stress–strain behavior, which is characterized by a toe region at low strains. Collagen tissues have a wavy structure (crimp) that is thought to produce this behavior due to its unfolding or rotation of fibers parallel to the stretch direction.[41,42,43]. Based on the SHG imaging of the collagen structure in the PAC, which was performed at the zero-stretch state of the tissue, it is possible that the non-linear stress–strain behavior we observed is a result of a re-orientation of the underlying collagen fiber network along the stretch direction. Additional studies are needed to determine how stretch affects the collagen fiber orientation and structural parameters in the subarachnoid space. To further understand possible location variances in this matrix, histological analysis of the PAC and the collagen crimp structure is needed.

Relevance of results to TBI modeling

In TBI, stress loads on the brain may be dependent on PAC and SAT morphology, similar to the significant biomechanical changes observed in spinal cord tissue properties when the PAC is removed [28]. These morphologies may change after multiple sub-concussive or concussive impacts, altering the impact stress loads on the brain from future head injuries and ultimately making the subject potentially prone to TBI. The removal or absence of pia has been shown, through simulations, to have an important effect on pressure distribution in the spine [44, 45]. The correlation between how TBI affects morphological changes and how those changes impact biomechanical properties is still largely unknown, but merits further investigation. Accurate representation of this tissue in simulations of TBI is critical for protective devices to be developed in high-risk situations such as impact sports. The provided biomechanical and morphological parameter results give information about a new animal model that can be used for finite element analysis experiments and interspecies comparisons.

Limitations

Variability across studies may be dependent on extraction and isolation methods [46] since PAC is sensitive to temperature and moisture content, and known to degrade and change rapidly post-mortem [28]. The present study employed methods analogous to prior work [22]. To avoid potential structural alteration with the testable tissue samples, thickness measurements were averaged from a separate set of samples (n = 10).

The resting state of the PAC is pre-stressed in situ on the brain. When removed from the surface of the brain, the tissue shrinks in size to a zero-stress condition. Similar to previous studies, a zero-stress condition was assumed. In the present study, the load frame coordinates were lowered by 2 mm after the PAC sample was secured in the soft grips to eliminate pre-stress conditions. An average strain rate of 0.59 s−1 was applied because this was the maximum strain rate that could be delivered by our apparatus. Specimen width was found to change by as much as 63% between its pre-stressed condition situated on the brain, to its unstressed condition in the C-clamp pre-test. To remain consistent with previous studies, room temperature (22 °C) was maintained throughout uniaxial tension tests [21, 22, 47]. Dehydration was shown to increase stiffness and strength in tissues made of collagen [48]. Other biological tissues show biomechanical alterations due to dehydration including the cellular level [49,50,51,52]. Anti-drying methods include saline solution spray [53, 54], silicone adhesive coating [55], silicon oil [56], or submersion in a saline bath [57]. The present study kept the samples hydrated with artificial CSF. Hartmann et al. developed the only known technique for in-vivo human SAT visualization and reported the depths of different neural structures [58]. This technique could be modified to obtain SAT morphological parameters.

In addition to the potential of tissue degradation, this study had a relatively low sample size. This was due to the number of brains and PAC samples that experienced visible tissue damage. Future studies may develop a gentler dissection method to protect the tissues. Additionally, average sample thickness was quantified using samples from multiple areas of the brain. While tests showed no significant difference in the mechanical properties of the four sample locations, future studies could control thickness based on the area (occipital, frontal and etcetera). Finally, this study did not control for age or sex. With higher sample sizes it would be possible to determine variances that occur with age or sex in addition to location.

Importantly, this and many previous studies were carried out using animal tissues. As our work potentially indicates that mechanical properties of the PAC vary based on species, quantifying the difference in mechanical properties relative to human tissue is vital for generating appropriate mechanical parameters for simulations of human TBI. Additional work could be carried out in concussed and non-concussed human samples of individual SAT fibers. These tests could be done along multiple axis with multiple strain rates. Future studies could be also be conducted with an alternative apparatus to test for higher strain rates that would be more applicable in TBI. An increased sample size could improve statistical significance of stress strain relationship between brain location and species.

Conclusion

This study represents the first biomechanical characterization of ovine PAC tissue. Methods were successfully developed to rapidly harvest and test samples across four extraction locations. Results are consistent and reproducible with previous studies of PAC biomechanical properties. This can be used to narrow hypotheses made about the stress–strain curve of PAC across different species and, in the future, humans. Further development of the understanding of PAC properties can be used to establish boundary conditions and nonhomogeneous material properties for numerical modeling of TBI under variable stress loads.

Availability of data and material

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Abbreviations

CSF:

Cerebrospinal Fluid

CNS:

Central Nervous System

DAPI:

4,6-Diamidino-2-phenylindole

PAC:

Pia-Arachnoid Complex

PBS:

Phosphate Buffered Saline

PBST:

Phosphate Buffered Saline with Triton

SAT:

Subarachnoid Trabeculae

SHG:

Second Harmonic Generation

TBI:

Traumatic Brain Injury

References

  1. 1.

    Mortazavi MM, Quadri SA, Khan MA, Gustin A, Suriya SS, Hassanzadeh T, et al. Subarachnoid trabeculae: a comprehensive review of their embryology, histology, morphology, and surgical significance. World Neurosurg. 2018;111:279–90.

    PubMed  Article  PubMed Central  Google Scholar 

  2. 2.

    Moore KL, Dalley AF. Clinically oriented anatomy. 4th ed. Philadelphia: Lippincott Williams & Wilkins; 1999.

    Google Scholar 

  3. 3.

    Saboori P, Sadegh A. Histology and morphology of the brain subarachnoid trabeculae. Anat Res Int. 2015;2015:279814.

    PubMed  PubMed Central  Google Scholar 

  4. 4.

    Scott GG, Coats B. Microstructural characterization of the pia-arachnoid complex using optical coherence tomography. IEEE Trans Med Imaging. 2015;34:1452–9.

    PubMed  Article  PubMed Central  Google Scholar 

  5. 5.

    Margulies S, Prange M, Myers BS, Maltese MR, Ji S, Ning X, et al. Shaken baby syndrome: a flawed biomechanical analysis. Forensic Sci Int. 2006;164:278–9.

    PubMed  Article  PubMed Central  Google Scholar 

  6. 6.

    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.

    PubMed  PubMed Central  Article  Google Scholar 

  7. 7.

    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:701–7.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  8. 8.

    Davidoff CL, Liu S, Wong JHY, Koustais S, Rogers JM, Stoodley MA. Treatment of syringomyelia in patients with arachnoiditis at the craniocervical junction. World Neurosurg. 2017;107:565–73.

    PubMed  Article  PubMed Central  Google Scholar 

  9. 9.

    Koyanagi I, Iwasaki Y, Hida K, Houkin K. Clinical features and pathomechanisms of syringomyelia associated with spinal arachnoiditis. Surg Neurol. 2005;63:350–5.

    PubMed  Article  PubMed Central  Google Scholar 

  10. 10.

    Cheng S, Stoodley MA, Wong J, Hemley S, Fletcher DF, Bilston LE. The presence of arachnoiditis affects the characteristics of CSF flow in the spinal subarachnoid space: a modelling study. J Biomech. 2012;45:1186–91.

    PubMed  Article  PubMed Central  Google Scholar 

  11. 11.

    Klekamp J, Iaconetta G, Batzdorf U, Samii M. Syringomyelia associated with foramen magnum arachnoiditis. J Neurosurg. 2002;97:317–22.

    PubMed  PubMed Central  Google Scholar 

  12. 12.

    Bilston LE, Fletcher DF, Stoodley MA. Focal spinal arachnoiditis increases subarachnoid space pressure: a computational study. Clin Biomech. 2006;21:579–84.

    CAS  Article  Google Scholar 

  13. 13.

    Gottschalk A, Schmitz B, Mauer UM, Bornstedt A, Steinhoff S, Danz B, et al. Dynamic visualization of arachnoid adhesions in a patient with idiopathic syringomyelia using high-resolution cine magnetic resonance imaging at 3T. J Magn Reson Imaging. 2010;32:218–22.

    PubMed  Article  PubMed Central  Google Scholar 

  14. 14.

    Tangen KM, Hsu Y, Zhu DC, Linninger AA. CNS wide simulation of flow resistance and drug transport due to spinal microanatomy. J Biomech. 2015;48:2144–54.

    PubMed  Article  PubMed Central  Google Scholar 

  15. 15.

    Stockman HW. Effect of anatomical fine structure on the flow of cerebrospinal fluid in the spinal subarachnoid space. J Biomech Eng. 2006;128:106–14.

    PubMed  Article  PubMed Central  Google Scholar 

  16. 16.

    Stockman HW. Effect of anatomical fine structure on the dispersion of solutes in the spinal subarachnoid space. J Biomech Eng. 2007;129:666–75.

    PubMed  Article  PubMed Central  Google Scholar 

  17. 17.

    Marshall-Goebel K, Terlevic R, Gerlach DA, Kuehn S, Mulder E, Rittweger J. Lower body negative pressure reduces optic nerve sheath diameter during head-down tilt. J Appl Physiol. 2017;1985(123):1139–44.

    Article  CAS  Google Scholar 

  18. 18.

    Linge SO, Haughton V, Lovgren AE, Mardal KA, Langtangen HP. CSF flow dynamics at the craniovertebral junction studied with an idealized model of the subarachnoid space and computational flow analysis. AJNR Am J Neuroradiol. 2010;31:185–92.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  19. 19.

    Gupta S, Soellinger M, Grzybowski DM, Boesiger P, Biddiscombe J, Poulikakos D, et al. Cerebrospinal fluid dynamics in the human cranial subarachnoid space: an overlooked mediator of cerebral disease. I. Computational model. J R Soc Interface. 2010;7:1195–204.

    PubMed  PubMed Central  Article  Google Scholar 

  20. 20.

    Jin X, Yang KH, King AI. Mechanical properties of bovine pia-arachnoid complex in shear. J Biomech. 2011;44:467–74.

    PubMed  Article  PubMed Central  Google Scholar 

  21. 21.

    Jin X, Ma C, Zhang L, Yang KH, King AI, Dong G, et al. Biomechanical response of the bovine pia-arachnoid complex to normal traction loading at varying strain rates. Stapp Car Crash J. 2007;51:115–26.

    PubMed  PubMed Central  Google Scholar 

  22. 22.

    Jin X, Lee JB, Leung LY, Zhang L, Yang KH, King AI. Biomechanical response of the bovine pia-arachnoid complex to tensile loading at varying strain-rates. Stapp Car Crash J. 2006;50:637–49.

    PubMed  PubMed Central  Google Scholar 

  23. 23.

    Raveling AR, Theodossiou SK, Schiele NR. A 3D printed mechanical bioreactor for investigating mechanobiology and soft tissue mechanics. MethodsX. 2018;5:924–32.

    PubMed  PubMed Central  Article  Google Scholar 

  24. 24.

    Koch RG, Tsamis A, D’Amore A, Wagner WR, Watkins SC, Gleason TG, et al. A custom image-based analysis tool for quantifying elastin and collagen micro-architecture in the wall of the human aorta from multi-photon microscopy. J Biomech. 2014;47:935–43.

    PubMed  PubMed Central  Article  Google Scholar 

  25. 25.

    Kimpara H, Nakahira Y, Iwamoto M, Miki K, Ichihara K, Kawano S, et al. Investigation of anteroposterior head-neck responses during severe frontal impacts using a brain-spinal cord complex FE model. Stapp Car Crash J. 2006;50:509–44.

    PubMed  PubMed Central  Google Scholar 

  26. 26.

    Ozawa H, Matsumoto T, Ohashi T, Sato M, Kokubun S. Mechanical properties and function of the spinal pia mater. J Neurosurg Spine. 2004;1:122–7.

    PubMed  Article  PubMed Central  Google Scholar 

  27. 27.

    Benko N, Luke E, Alsanea Y, Coats B. Spatial distribution of human arachnoid trabeculae. J Anat. 2020;237:275.

    PubMed  Article  PubMed Central  Google Scholar 

  28. 28.

    Ramo NL, Troyer KL, Puttlitz CM. Viscoelasticity of spinal cord and meningeal tissues. Acta Biomater. 2018;75:253–62.

    PubMed  Article  PubMed Central  Google Scholar 

  29. 29.

    Jin X, Mao H, Yang KH, King AI. Constitutive modeling of pia-arachnoid complex. Ann Biomed Eng. 2014;42:812–21.

    PubMed  Article  PubMed Central  Google Scholar 

  30. 30.

    Rashid B, Destrade M, Gilchrist MD. Mechanical characterization of brain tissue in simple shear at dynamic strain rates. J Mech Behav Biomed Mater. 2013;28:71–85.

    PubMed  Article  PubMed Central  Google Scholar 

  31. 31.

    Hosseini Farid M, Eslaminejad A, Ziejewski M, Karami G. A study on the effects of strain rates on characteristics of brain tissue. In ASME 2017 International Mechanical Engineering Congress and Exposition, 2017.

  32. 32.

    Jin X. Biomechanical response and constitutive modeling of bovine pia-arachnoid complex. 3380487 Ph.D., Ann Arbor: Wayne State University, 2009.

  33. 33.

    Leclerc GE, Debernard L, Foucart F, Robert L, Pelletier KM, Charleux F, et al. Characterization of a hyper-viscoelastic phantom mimicking biological soft tissue using an abdominal pneumatic driver with magnetic resonance elastography (MRE). J Biomech. 2012;45:952–7.

    PubMed  PubMed Central  Article  Google Scholar 

  34. 34.

    Mendis KK, Stalnaker RL, Advani SH. A constitutive relationship for large deformation finite element modeling of brain tissue. J Biomech Eng. 1995;117:279–85.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  35. 35.

    Hu T, Desai JP. Characterization of soft-tissue material properties: large deformation analysis. Berlin, Heidelberg: Springer; 2004. p. 28–37.

    Google Scholar 

  36. 36.

    Saccomandi G, Vergori L. Generalised Mooney-Rivlin models for brain tissue: a theoretical perspective. Int J Non Linear Mech. 2019;109:9–14.

    Article  Google Scholar 

  37. 37.

    Tsamis A, Stergiopulos N. Arterial remodeling in response to hypertension using a constituent-based model. Am J Physiol Heart Circ Physiol. 2007;293:H3130–9.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  38. 38.

    Haines DE, Mihailoff GA. Fundamental neuroscience for basic and clinical applications (Fifth edition). 2018. https://SK8ES4MC2L.search.serialssolutions.com/?sid=sersol&SS_jc=TC0001882702&title=Fundamental%20neuroscience%20for%20basic%20and%20clinical%20applications

  39. 39.

    Killer HE, Laeng HR, Flammer J, Groscurth P. Architecture of arachnoid trabeculae, pillars, and septa in the subarachnoid space of the human optic nerve: anatomy and clinical considerations. Br J Ophthalmol. 2003;87:777–81.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  40. 40.

    Fabris G, Zeynep MS, Kurt M. Micromechanical heterogeneity of the rat pia-arachnoid complex. Acta Biomater. 2019;100:29–37.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  41. 41.

    Guerin HA, Elliott DM. The role of fiber-matrix interactions in a nonlinear fiber-reinforced strain energy model of tendon. J Biomech Eng. 2005;127:345–50.

    PubMed  Article  PubMed Central  Google Scholar 

  42. 42.

    Viidik A, Danielson CC, Oxlund H. On fundamental and phenomenological models, structure and mechanical properties of collagen, elastin and glycosaminoglycan complexes. Biorheology. 1982;19:437–51.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  43. 43.

    Fratzl P, Misof K, Zizak I, Rapp G, Amenitsch H, Bernstorff S. Fibrillar structure and mechanical properties of collagen. J Struct Biol. 1998;122:119–22.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  44. 44.

    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:686–98.

    PubMed  Article  PubMed Central  Google Scholar 

  45. 45.

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

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  46. 46.

    Garo A, Hrapko M, van Dommelen JA, Peters GW. Towards a reliable characterisation of the mechanical behaviour of brain tissue: the effects of post-mortem time and sample preparation. Biorheology. 2007;44:51–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  47. 47.

    Aimedieu P, Grebe R. Tensile strength of cranial pia mater: preliminary results. J Neurosurg. 2004;100:111–4.

    PubMed  Article  PubMed Central  Google Scholar 

  48. 48.

    Betsch DF, Baer E. Structure and mechanical properties of rat tail tendon. Biorheology. 1980;17:83–94.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  49. 49.

    Nyman JS, Roy A, Shen X, Acuna RL, Tyler JH, Wang X. The influence of water removal on the strength and toughness of cortical bone. J Biomech. 2006;39:931–8.

    PubMed  PubMed Central  Article  Google Scholar 

  50. 50.

    Thornton GM, Shrive NG, Frank CB. Altering ligament water content affects ligament pre-stress and creep behaviour. J Orthop Res. 2001;19:845–51.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  51. 51.

    Jayasuriya AC, Scheinbeim JI, Lubkin V, Bennett G, Kramer P. Piezoelectric and mechanical properties in bovine cornea. J Biomed Mater Res A. 2003;66:260–5.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  52. 52.

    Mendelson NH, Thwaites JJ. Cell wall mechanical properties as measured with bacterial thread made from Bacillus subtilis. J Bacteriol. 1989;171:1055–62.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  53. 53.

    Darvish KK, Crandall JR. Nonlinear viscoelastic effects in oscillatory shear deformation of brain tissue. Med Eng Phys. 2001;23:633–45.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  54. 54.

    Takhounts EG, Crandall JR, Darvish K. On the importance of nonlinearity of brain tissue under large deformations. Stapp Car Crash J. 2003;47:79–92.

    PubMed  PubMed Central  Google Scholar 

  55. 55.

    Fallenstein GT, Hulce VD, Melvin JW. Dynamic mechanical properties of human brain tissue. J Biomech. 1969;2:217–26.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  56. 56.

    Shen F, Tay TE, Li JZ, Nigen S, Lee PV, Chan HK. Modified Bilston nonlinear viscoelastic model for finite element head injury studies. J Biomech Eng. 2006;128:797–801.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  57. 57.

    Tamura A, Omori K, Miki K, Lee JB, Yang KH, King AI. Mechanical characterization of porcine abdominal organs. Stapp Car Crash J. 2002;46:55–69.

    PubMed  PubMed Central  Google Scholar 

  58. 58.

    Hartmann K, Stein KP, Neyazi B, Sandalcioglu IE. First in vivo visualization of the human subarachnoid space and brain cortex via optical coherence tomography. Ther Adv Neurol Disord. 2019;12:1756286419843040.

    PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Goutham Burla1, Email: gburla@uidaho.edu. Thanks to Claire L. Majors1 for artwork in Figs. 1 and 2, Email: clairemajors11@gmail.com. James Nasados, Email: jnasados@uidaho.edu

Funding

This work was supported by an Institutional Development Award (IDeA) from the National Institute of General Medical Sciences (NIGMS) of the National Institutes of health (NIH) under Grant #P20GM1033408 (Martin & Schiele) and #4U54GM104944-04TBD (Martin), NIH National Institute of Neurological Disorders and Stroke under Grant #R01NS111283, American Heart Association under Grant [19TPA34860008 to BCWT] (Tanner), and National Science Foundation under Grant [1656450 to BCWT] (Tanner).

Author information

Affiliations

Authors

Contributions

Study conception and design: BAM, NRS, GKM, GP, MM. Acquisition of data: BAM, GCN, GB, GKM, SKT, NRS, CLM. Analysis and interpretation of data: BAM, GCN, GB, NRS, SKT, GKM, AT. Drafting of manuscript: BAM, GCN, GB, SKT, GKM. Critical revision: BAM, MM, BT, NRS, AT, GP. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Bryn A. Martin.

Ethics declarations

Ethics approval and consent to participate

This study was approved as exempt by the local governing IACUC review board.

Consent for publication

N/A

Competing interests

BAM is a full-time employee of Alcyone Therapeutics and has received grant support from Minnetronix Neuro, Biogen, Genentech, Voyager Therapeutics, and Alcyone Therapeutics. BAM is a scientific advisory board member for Alcyone Lifesciences, International Society for Hydrocephalus and CSF Disorders, and the Chiari and Syringomyelia Foundation and has served as a consultant to Flux Neuroscience, Genentech, Roche, Biogen, Praxis Medicines, InviCRO, SwanBio Therapeutics, Cerebral Therapeutics, Minnetronix Neuro, and CereVasc.

Additional information

Publisher's Note

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

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Conley Natividad, G., Theodossiou, S.K., Schiele, N.R. et al. Ex-vivo quantification of ovine pia arachnoid complex biomechanical properties under uniaxial tension. Fluids Barriers CNS 17, 68 (2020). https://doi.org/10.1186/s12987-020-00229-w

Download citation