Scaffold Geometry-Imposed Anisotropic Mechanical Loading Guides the Evolution of the Mechanical State of Engineered Cardiovascular Tissues in vitro.

Cardiovascular tissue engineering is a promising approach to develop grafts that, in contrast to current replacement grafts, have the capacity to grow and remodel like native tissues. This approach largely depends on cell-driven tissue growth and remodeling, which are highly complex processes that are difficult to control inside the scaffolds used for tissue engineering. For several tissue engineering approaches, adverse tissue growth and remodeling outcomes were reported, such as aneurysm formation in vascular grafts, and leaflet retraction in heart valve grafts. It is increasingly recognized that the outcome of tissue growth and remodeling, either physiological or pathological, depends at least partly on the establishment of a homeostatic mechanical state, where one or more mechanical quantities in a tissue are maintained in equilibrium. To design long-term functioning tissue engineering strategies, understanding how scaffold parameters such as geometry affect the mechanical state of a construct, and how this state guides tissue growth and remodeling, is therefore crucial. Here, we studied how anisotropic versus isotropic mechanical loading-as imposed by initial scaffold geometry-influences tissue growth, remodeling, and the evolution of the mechanical state and geometry of tissue-engineered cardiovascular constructs in vitro. Using a custom-built bioreactor platform and nondestructive mechanical testing, we monitored the mechanical and geometric changes of elliptical and circular, vascular cell-seeded, polycaprolactone-bisurea scaffolds during 14 days of dynamic loading. The elliptical and circular scaffold geometries were designed using finite element analysis, to induce anisotropic and isotropic dynamic loading, respectively, with similar maximum stretch when cultured in the bioreactor platform. We found that the initial scaffold geometry-induced (an)isotropic loading of the engineered constructs differentially dictated the evolution of their mechanical state and geometry over time, as well as their final structural organization. These findings demonstrate that controlling the initial mechanical state of tissue-engineered constructs via scaffold geometry can be used to influence tissue growth and remodeling and determine tissue outcomes.

Spectroscopic photoacoustic imaging of cartilage damage.

OBJECTIVE: Osteoarthritis (OA) is a chronic joint disease characterized by progressive degradation of cartilage. It affects more than 10% of the people aged over 60 years-old worldwide with a rising prevalence due to the increasingly aging population. OA is a major source of pain, disability, and socioeconomic cost. Currently, the lack of effective diagnosis and affordable imaging options for early detection and monitoring of OA presents the clinic with many challenges. Spectroscopic Photoacoustic (sPA) imaging has the potential to reveal changes in cartilage composition with different degrees of damage, based on optical absorption contrast. DESIGN: In this study, the capabilities of sPA imaging and its potential to characterize cartilage damage were explored. To this end, 15 pieces of cartilage samples from patients undergoing a total joint replacement were collected and were imaged ex vivo with sPA imaging at a wide optical spectral range (between 500 nm and 1,300 nm) to investigate the photoacoustic properties of cartilage tissue. All the PA spectra of the cartilage samples were analyzed and compared to the corresponding histological results. RESULTS: The collagen related PA spectral changes were clearly visible in our imaging data and were related to different degrees of cartilage damage. The results are in good agreement with histology and the current gold standard, i.e., the Mankin score. CONCLUSIONS: This study demonstrates the potential and possible clinical application of sPA imaging in OA.

Unmixing multi-spectral photoacoustic sources in human carotid plaques using non-negative independent component analysis.

Multi-spectral photoacoustic imaging (MSPAI) is promising for morphology assessment of carotid plaques; however, obtaining unique spectral characteristics of chromophores is cumbersome. We used MSPAI and non-negative independent component analysis (ICA) to unmix distinct signal sources in human carotid plaques blindly. The feasibility of the method was demonstrated on a plaque phantom with hemorrhage and cholesterol inclusions, and plaque endarterectomy samples ex vivo. Furthermore, the results were verified with histology using Masson's trichrome staining. Results showed that ICA could separate recent hemorrhages from old hemorrhages. Additionally, the signatures of cholesterol inclusion were also captured for the phantom experiment. Artifacts were successfully removed from signal sources. Histologic examinations showed high resemblance with the unmixed components and confirmed the morphologic distinction between recent and mature hemorrhages. In future pre-clinical studies, unmixing could be used for morphology assessment of intact human plaque samples.

Initial scaffold thickness affects the emergence of a geometrical and mechanical equilibrium in engineered cardiovascular tissues.

In situ cardiovascular tissue-engineering can potentially address the shortcomings of the current replacement therapies, in particular, their inability to grow and remodel. In native tissues, it is widely accepted that physiological growth and remodelling occur to maintain a homeostatic mechanical state to conserve its function, regardless of changes in the mechanical environment. A similar homeostatic state should be reached for tissue-engineered (TE) prostheses to ensure proper functioning. For in situ tissue-engineering approaches obtaining such a state greatly relies on the initial scaffold design parameters. In this study, it is investigated if the simple scaffold design parameter initial thickness, influences the emergence of a mechanical and geometrical equilibrium state in in vitro TE constructs, which resemble thin cardiovascular tissues such as heart valves and arteries. Towards this end, two sample groups with different initial thicknesses of myofibroblast-seeded polycaprolactone-bisurea constructs were cultured for three weeks under dynamic loading conditions, while tracking geometrical and mechanical changes temporally using non-destructive ultrasound imaging. A mechanical equilibrium was reached in both groups, although at different magnitudes of the investigated mechanical quantities. Interestingly, a geometrically stable state was only established in the thicker constructs, while the thinner constructs' length continuously increased. This demonstrates that reaching geometrical and mechanical stability in TE constructs is highly dependent on functional scaffold design.

Perfusion dynamics assessment with Power Doppler ultrasound in skeletal muscle during maximal and submaximal cycling exercise.

PURPOSE: Assessment of limitations in the perfusion dynamics of skeletal muscle may provide insight in the pathophysiology of exercise intolerance in, e.g., heart failure patients. Power doppler ultrasound (PDUS) has been recognized as a sensitive tool for the detection of muscle blood flow. In this volunteer study (N = 30), a method is demonstrated for perfusion measurements in the vastus lateralis muscle, with PDUS, during standardized cycling exercise protocols, and the test-retest reliability has been investigated. METHODS: Fixation of the ultrasound probe on the upper leg allowed for continuous PDUS measurements. Cycling exercise protocols included a submaximal and an incremental exercise to maximal power. The relative perfused area (RPA) was determined as a measure of perfusion. Absolute and relative reliability of RPA amplitude and kinetic parameters during exercise (onset, slope, maximum value) and recovery (overshoot, decay time constants) were investigated. RESULTS: A RPA increase during exercise followed by a signal recovery was measured in all volunteers. Amplitudes and kinetic parameters during exercise and recovery showed poor to good relative reliability (ICC ranging from 0.2-0.8), and poor to moderate absolute reliability (coefficient of variation (CV) range 18-60%). CONCLUSIONS: A method has been demonstrated which allows for continuous (Power Doppler) ultrasonography and assessment of perfusion dynamics in skeletal muscle during exercise. The reliability of the RPA amplitudes and kinetics ranges from poor to good, while the reliability of the RPA increase in submaximal cycling (ICC = 0.8, CV = 18%) is promising for non-invasive clinical assessment of the muscle perfusion response to daily exercise.

Three-dimensional ultrasound strain imaging of skeletal muscles.

In this study, a multi-dimensional strain estimation method is presented to assess local relative deformation in three orthogonal directions in 3D space of skeletal muscles during voluntary contractions. A rigid translation and compressive deformation of a block phantom, that mimics muscle contraction, is used as experimental validation of the 3D technique and to compare its performance with respect to a 2D based technique. Axial, lateral and (in case of 3D) elevational displacements are estimated using a cross-correlation based displacement estimation algorithm. After transformation of the displacements to a Cartesian coordinate system, strain is derived using a least-squares strain estimator. The performance of both methods is compared by calculating the root-mean-squared error of the estimated displacements with the calculated theoretical displacements of the phantom experiments. We observe that the 3D technique delivers more accurate displacement estimations compared to the 2D technique, especially in the translation experiment where out-of-plane motion hampers the 2D technique. In vivo application of the 3D technique in the musculus vastus intermedius shows good resemblance between measured strain and the force pattern. Similarity of the strain curves of repetitive measurements indicates the reproducibility of voluntary contractions. These results indicate that 3D ultrasound is a valuable imaging tool to quantify complex tissue motion, especially when there is motion in three directions, which results in out-of-plane errors for 2D techniques.

Patient Specific Wall Stress Analysis and Mechanical Characterization of Abdominal Aortic Aneurysms Using 4D Ultrasound.

OBJECTIVES: The aim of this study was to perform wall stress analysis (WSA) using 4D ultrasound (US) in 40 patients with an abdominal aortic aneurysm (AAA). The geometries and wall stress results were compared with computed tomography (CT) in seven patients. Additionally, the WSA models were calibrated using 4D motion estimation, resulting in patient specific material parameters that were compared among patients. METHODS: 4D-US images were acquired for 40 patients (AAA diameter 27-52 mm). Patient specific AAA geometries and wall motion were extracted from the 4D-US. WSA was performed and corresponding patient specific material properties were derived. For seven patients, CT data were available and analyzed for geometry and wall stress comparison. RESULTS: The 4D-US based 99th percentile wall stress ranged from 198 to 390 kPa. Regression analysis showed no significant relation between wall stress and diameter of the AAA. The similarity indices between US and CT were very good and ranged between 0.90 and 0.96, and the 25th, 50th, 75th, and 95th percentile wall stresses of the US and CT data were in agreement. The characterized patient specific shear modulus had a median of 1.1 MPa (interquartile range, 0.7-1.4 MPa). Based on the maximum AAA diameter, the AAAs were divided in a small, medium, and large diameter groups. The largest AAAs revealed an increased wall stiffness compared with the smallest AAAs. CONCLUSIONS: 4D ultrasound is applicable for wall stress analysis of AAAs, and offers the opportunity to perform wall stress analysis over time, also for AAAs who do not qualify for a CT or magnetic resonance imaging. Moreover, the patient specific material properties can be determined, which could possibly improve risk assessment.

Full 2D displacement vector and strain tensor estimation for superficial tissue using beam-steered ultrasound imaging.

Ultrasound strain imaging is used to measure local tissue deformations. Usually, only strains along the ultrasound beam are estimated, because those estimates are most precise, due to the availability of phase information. For estimating strain in other directions we propose to steer the ultrasound beam at an angle, which allows estimating different projections of the 2D strain tensor, while phase information remains available. This study investigates beam steering at maximally three different angles to determine the full 2D strain tensor. The method was tested on simulated and experimental data of an inclusion phantom and a vessel phantom. The combination of data from a non-steered acquisition and acquisitions at a large positive and an equally large but negative steering angle enabled the most precise estimation of the strain components. The method outperforms conventional methods that do not use beam steering.

Cardiac biplane strain imaging: initial in vivo experience.

In this study, first we propose a biplane strain imaging method using a commercial ultrasound system, yielding estimation of the strain in three orthogonal directions. Secondly, an animal model of a child's heart was introduced that is suitable to simulate congenital heart disease and was used to test the method in vivo. The proposed approach can serve as a framework to monitor the development of cardiac hypertrophy and fibrosis. A 2D strain estimation technique using radio frequency (RF) ultrasound data was applied. Biplane image acquisition was performed at a relatively low frame rate (<100 Hz) using a commercial platform with an RF interface. For testing the method in vivo, biplane image sequences of the heart were recorded during the cardiac cycle in four dogs with an aortic stenosis. Initial results reveal the feasibility of measuring large radial, circumferential and longitudinal cumulative strain (up to 70%) at a frame rate of 100 Hz. Mean radial strain curves of a manually segmented region-of-interest in the infero-lateral wall show excellent correlation between the measured strain curves acquired in two perpendicular planes. Furthermore, the results show the feasibility and reproducibility of assessing radial, circumferential and longitudinal strains simultaneously. In this preliminary study, three beagles developed an elevated pressure gradient over the aortic valve (Deltap: 100-200 mmHg) and myocardial hypertrophy. One dog did not develop any sign of hypertrophy (Deltap = 20 mmHg). Initial strain (rate) results showed that the maximum strain (rate) decreased with increasing valvular stenosis (-50%), which is in accordance with previous studies. Histological findings corroborated these results and showed an increase in fibrotic tissue for the hearts with larger pressure gradients (100, 200 mmHg), as well as lower strain and strain rate values.

In vivo validation of cardiac output assessment in non-standard 3D echocardiographic images.

Automatic segmentation of the endocardial surface in three-dimensional (3D) echocardiographic images is an important tool to assess left ventricular (LV) geometry and cardiac output (CO). The presence of speckle noise as well as the nonisotropic characteristics of the myocardium impose strong demands on the segmentation algorithm. In the analysis of normal heart geometries of standardized (apical) views, it is advantageous to incorporate a priori knowledge about the shape and appearance of the heart. In contrast, when analyzing abnormal heart geometries, for example in children with congenital malformations, this a priori knowledge about the shape and anatomy of the LV might induce erroneous segmentation results. This study describes a fully automated segmentation method for the analysis of non-standard echocardiographic images, without making strong assumptions on the shape and appearance of the heart. The method was validated in vivo in a piglet model. Real-time 3D echocardiographic image sequences of five piglets were acquired in radiofrequency (rf) format. These ECG-gated full volume images were acquired intra-operatively in a non-standard view. Cardiac blood flow was measured simultaneously by an ultrasound transit time flow probe positioned around the common pulmonary artery. Three-dimensional adaptive filtering using the characteristics of speckle was performed on the demodulated rf data to reduce the influence of speckle noise and to optimize the distinction between blood and myocardium. A gradient-based 3D deformable simplex mesh was then used to segment the endocardial surface. A gradient and a speed force were included as external forces of the model. To balance data fitting and mesh regularity, one fixed set of weighting parameters of internal, gradient and speed forces was used for all data sets. End-diastolic and end-systolic volumes were computed from the segmented endocardial surface. The cardiac output derived from this automatic segmentation was validated quantitatively by comparing it with the CO values measured from the volume flow in the pulmonary artery. Relative bias varied between 0 and -17%, where the nominal accuracy of the flow meter is in the order of 10%. Assuming the CO measurements from the flow probe as a gold standard, excellent correlation (r = 0.99) was observed with the CO estimates obtained from image segmentation.

Identifiability analysis of the standard pharmacokinetic models in DCE MR imaging of tumours.

The usage of dynamic contrast-enhanced MRI (DCE-MRI) as a clinical tool is still widely assessed. Application of the standard pharmacokinetic models to obtain physiologically relevant parameter values using DCE-MRI in tumours is not trivial, when the temporal resolution is low. Mathematical analysis and analysis by simulation of the identifiability for the generalized and extended Kety models was executed. Parameter estimation was executed using synthetic data sets and maximum likelihood estimation (MLE). The influence of temporal resolution was examined. The generalized and extended Kety model showed a large bias in the parameter estimates (10-120%) for sampling times >4 s, although the estimated variance was relatively low (<1%). This was in accordance with the generated contour plots of the hyperplane of the MLE cost-function. The influence of measurement noise on the input and output turned out to be less significant than the temporal resolution.