Abdominal aortic aneurysms (AAA) are focal enlargements of the aorta that occur preferentially near the aortic bifurcation1. The exact cause of AAA formation is unknown, although many theories suggest that the pathogenesis is multifactorial, with genetic, behavioral, hemodynamic, and environmental factors contributing2,3. While the diagnosis of an abdominal aortic aneurysm can be obtained using non-invasive imaging techniques, the prediction of patient-specific rupture risk is not as precise4,5,6. Surgical repair can reduce the risk of aortic rupture, but operative repair of the aorta carries a high rate of associated morbidity and mortality7. Current surgical practices use the "maximum size criteria", or maximum absolute diameter of the aneurysms, to predict a patient's risk of rupture. Unfortunately, it has been well established that an aneurysm still ruptures below sizes clinically acceptable for surgical repair, meaning that patients with any sized aneurysm carry some risk of rupture8,9,10,11,12,13. Additionally, it is known that historical reports of rupture risk are likely over-estimations of the true rupture risk, meaning many patients are exposed to surgical risk without benefit13. A more accurate assessment of the patient-specific rupture risk is needed to help stratify a patient's risk-benefit ratio of undergoing surgical aneurysm repair.
It has been shown that the spatial stress distribution within an AAA is of critical importance in determining rupture potential and may be a better indicator than maximum diameter14,15,16,17,18. Most of the recent studies that investigate the mechanics of AAA rupture use segmented geometries from X-ray computed tomography (CT) images, and population averaged mechanical properties of aortic tissue measured ex vivo. Finite element (FE) models are then used to predict the vessel wall stresses14,15,16,17,18. However, because the mechanical properties are determined following the tissue excision, it is unclear whether the resulting models accurately depict the resulting in vivo patient-specific stresses. These studies typically assume homogeneous vessel wall material properties and don't account for the highly heterogeneous structure of the aortic wall and thrombus19,20,21,22,23,24,25.
Ultrasound-based elasticity imaging is used clinically to diagnose and monitor a variety of disease pathologies26. This technology provides a non-invasive means to interrogate the physical interactions of soft tissues. Vascular US elasticity imaging has been used as an adjunct imaging modality to clinical US evaluation in the screening and monitoring of AAAs. The combination of these techniques provides both geometric information, such as diameter and length, as well as mechanical data, such as relative stiffness and stiffness variation. While many elasticity imaging techniques require an external load to induce a measurable tissue deformation, the tissue motion to be measured here is induced by changes in the aortic pressure caused by the beating heart. Numerous methods have been published to spatially resolve strain fields in deforming vessels, however, validation studies of these methods have been limited to human patients, animal models, or ex vivo tissue samples27,28,29,30,31,32. To date, few methods allow for creations of custom geometries with spatially varied material properties27,29.
Here we present a method of manufacturing US compatible, tissue-mimicking phantoms that can be tailored to a variety of relevant aortic geometries and material properties for validation of US elastography techniques. Although previous groups have been able to design complex geometry phantoms to mimic AAA geometries using 3D printing technology33,34, printable rubbers are known to have a high attenuation to US and do not have a means to later their material properties. Phantoms are made from polyvinyl alcohol cryogel (PVA-c), which has been previously shown to be ideal for mimicking vascular tissue properties35. These phantoms can be used in US, magnetic resonance, and elastographic imaging36,37,38. The aortic aneurysm geometry was designed similarly to that of the simulation model created by Vorp et al.14. The vessel has a nominal diameter of 22.5 mm and has an aneurysmal bulge that is 64 mm bulge long, 47 mm in diameter and eccentric (β = 0.6)14 to the anterior side of the phantom. The last section mimics the iliac bifurcation with a distal diameter of 15 mm. The phantom was chosen to have a constant thickness of approximately 5 mm. Raghavan et al. reported in a small study that the vessel thickness of AAA ranges from 0.23-4.26 mm, with a median value of 1.48 mm39. A nominal vessel thickness on the larger end of that spectrum was chosen here for manufacturing concerns with the expectation that improved 3D printing techniques will improve the minimum phantom thickness that is able to be molded. Phantom molds were designed in CAD and are 3D printed using commercially available printers and filament.
The molds are injection filled with the PVA-c solution and subjected to a series of freeze/thaw cycles (-20 °C and +20°C) to cross-link the PVA-c polymer and polymerize the gel. The elastic modulus of the PVA-c is controlled by altering the concentration of the PVA-c gel or the number of freeze-thaw cycles. The aneurysmal section of the phantom required loss mold to remove from the inner lumen of the vessel. This was accomplished by the use of a polyvinyl alcohol, 3D printer filament (PVA). Although chemically similar to the PVA-c powder, the PVA filament does not polymerize when frozen and, as such, can be dissolved in water after the PVA-c has been set. Additional sample molds are printed to create tensile testing specimens, in a "dog bone" configuration, with the same PVA-c concentration. These molds undergo the same freeze/thaw cycles and are used for tensile testing to independently measure the elastic modulus of the phantom sections. A background material was manufactured with softer PVA-c, made to simulate tissues of the retroperitoneum40,41. This background phantom was manufactured as a homogeneous axisymmetric cylindrical tube with a 4 cm inner diameter, a 16.5 cm outer diameter, and a length of 16.5 cm. It was made from a 5% PVA solution and subjected to a total of two freeze-thaw cycles.
The final AAA phantoms were placed in the background phantom and connected, via tube fittings and clamps, to a hemodynamic water pump designed to deform the phantoms with physiologic cyclic flows and pressures. The pump speed was set to deliver approximately a 6-7 kPa pressure pulse at a rate of approximately 1 Hz. Ultra sound image sequences of the deforming phantoms were collected, and the pressure normalized strain was calculated to identify differences in the spatially varied mechanical properties. Representative results of the pressure normalized strain images within the vessel region are presented. The increasing regional differences in the normalized strain of the stiffer heterogenous phantoms, relative to the homogeneous phantom, demonstrate the differences in the vessel stiffness and our ability to measure it.