Moving Domain Computational Fluid Dynamics to Interface with an Embryonic Model of Cardiac Morphogenesis

Peristaltic contraction of the embryonic heart tube produces time- and spatial-varying wall shear stress (WSS) and pressure gradients (∇P) across the atrioventricular (AV) canal. Zebrafish (Danio rerio) are a genetically tractable system to investigate cardiac morphogenesis. The use of Tg(fli1a:EGFP)y1 transgenic embryos allowed for delineation and two-dimensional reconstruction of the endocardium. This time-varying wall motion was then prescribed in a two-dimensional moving domain computational fluid dynamics (CFD) model, providing new insights into spatial and temporal variations in WSS and ∇P during cardiac development. The CFD simulations were validated with particle image velocimetry (PIV) across the atrioventricular (AV) canal, revealing an increase in both velocities and heart rates, but a decrease in the duration of atrial systole from early to later stages. At 20-30 hours post fertilization (hpf), simulation results revealed bidirectional WSS across the AV canal in the heart tube in response to peristaltic motion of the wall. At 40-50 hpf, the tube structure undergoes cardiac looping, accompanied by a nearly 3-fold increase in WSS magnitude. At 110-120 hpf, distinct AV valve, atrium, ventricle, and bulbus arteriosus form, accompanied by incremental increases in both WSS magnitude and ∇P, but a decrease in bi-directional flow. Laminar flow develops across the AV canal at 20-30 hpf, and persists at 110-120 hpf. Reynolds numbers at the AV canal increase from 0.07±0.03 at 20-30 hpf to 0.23±0.07 at 110-120 hpf (p< 0.05, n=6), whereas Womersley numbers remain relatively unchanged from 0.11 to 0.13. Our moving domain simulations highlights hemodynamic changes in relation to cardiac morphogenesis; thereby, providing a 2-D quantitative approach to complement imaging analysis. © 2013 Lee et al

Imaging of the Lamina Cribrosa using Swept-Source Optical Coherence Tomography.

The lamina cribrosa (LC) is the presumed site of axonal injury in glaucoma. Its deformation has been suggested to contribute to optic neuropathy by impeding axoplasmic flow within the optic nerve fibers, leading to apoptosis of retinal ganglion cells. To visualize the LC in vivo, optical coherence tomography (OCT) has been applied. Spectral domain (SD)-OCT, used in conjunction with recently introduced enhanced depth imaging (EDI)-OCT, has improved visualization of deeper ocular layers, but in many individuals it is still limited by inadequate resolution, poor image contrast and insufficient depth penetrance. The posterior laminar surface especially is not viewed clearly using these methods. New generation high-penetration (HP)-OCTs, also known as swept-source (SS)-OCT, are capable to evaluate the choroid in vivo to a remarkable level of detail. SS-OCTs use a longer wavelength (1,050 nm instead of 840 nm) compared to the conventional techniques. We review current knowledge of the LC, findings from trials that use SD-OCT and EDI-OCT, and our experience with a prototype SS-OCT to visualize the LC in its entirety. Key Points What is known? •     The LC is the presumed site of axonal injury in glaucoma •     Compared to spectral domain-OCT, enhanced depth imaging-OCT improves imaging of the LC •     Even so, currently used SD-OCT techniques are restricted by poor wavelength penetrance of the deeper ocular layers What our findings add? •    SS-OCT may be a superior imaging modality for deep ocular structures •    Prior studies used SS-OCT to evaluate choroidal thickness in both healthy and 'normal tension glaucoma' eyes •    SS-OCT enables good evaluation of three-dimension (3D) lamina cribrosa morphology. How to cite this article: Nuyen B, Mansouri K, Weinreb RN. Imaging of the Lamina Cribrosa using Swept-Source Optical Coherence Tomography. J Current Glau Prac 2012;6(3): 113-119

Modeling the variability of shapes of a human placenta

While it is well-understood what a normal human placenta should look like, a deviation from the norm can take many possible shapes. In this paper we propose a mechanism for this variability based on the change in the structure of the vascular tree

Three-dimensional instabilities and transition of steady and pulsatile axisymmetric stenotic flows

Published versio

Inverse modelling of image-based patient-specific blood vessels : zero-pressure geometry and in vivo stress incorporation

In vivo visualization of cardiovascular structures is possible using medical images. However, one has to realize that the resulting 3D geometries correspond to in vivo conditions. This entails an internal stress state to be present in the in vivo measured geometry of e.g. a blood vessel due to the presence of the blood pressure. In order to correct for this in vivo stress, this paper presents an inverse method to restore the original zero-pressure geometry of a structure, and to recover the in vivo stress field of the final, loaded structure. The proposed backward displacement method is able to solve the inverse problem iteratively using fixed point iterations, but can be significantly accelerated by a quasi-Newton technique in which a least-squares model is used to approximate the inverse of the Jacobian. The here proposed backward displacement method allows for a straightforward implementation of the algorithm in combination with existing structural solvers, even if the structural solver is a black box, as only an update of the coordinates of the mesh needs to be performed

Implantation of 3D-Printed Patient-Specific Aneurysm Models into Cadaveric Specimens: A New Training Paradigm to Allow for Improvements in Cerebrovascular Surgery and Research.

AimTo evaluate the feasibility of implanting 3D-printed brain aneurysm model in human cadavers and to assess their utility in neurosurgical research, complex case management/planning, and operative training.MethodsTwo 3D-printed aneurysm models, basilar apex and middle cerebral artery, were generated and implanted in four cadaveric specimens. The aneurysms were implanted at the same anatomical region as the modeled patient. Pterional and orbitozygomatic approaches were done on each specimen. The aneurysm implant, manipulation capabilities, and surgical clipping were evaluated.ResultsThe 3D aneurysm models were successfully implanted to the cadaveric specimens' arterial circulation in all cases. The features of the neck in terms of flexibility and its relationship with other arterial branches allowed for the practice of surgical maneuvering characteristic to aneurysm clipping. Furthermore, the relationship of the aneurysm dome with the surrounding structures allowed for better understanding of the aneurysmal local mass effect. Noticeably, all of these observations were done in a realistic environment provided by our customized embalming model for neurosurgical simulation.Conclusion3D aneurysms models implanted in cadaveric specimens may represent an untapped training method for replicating clip technique; for practicing certain approaches to aneurysms specific to a particular patient; and for improving neurosurgical research

MODELLING THE INFLUENCE OF NUCLEUS ELASTICITY ON CELL INVASION IN FIBER NETWORKS AND MICROCHANNELS

Cell migration in highly constrained extracellular matrices is exploited in scaffold-based tissue engineering and is fundamental in a wide variety of physiological and pathological phenomena, among others in cancer invasion and development. Research into the critical processes involved in cell migration has mainly focused on cell adhesion and proteolytic degradation of the external environment. However, rising evidence has recently shown that a number of cell-derived biophysical and mechanical parameters, among others nucleus stiffness and cell deformability, plays a major role in cell motility, especially in the ameboid-like migration mode in 3D confined tissue structures. We here present an extended cellular Potts model (CPM) first used to simulate a micro-fabricated migration chip, which tests the active invasive behavior of cancer cells into narrow channels. As distinct features of our approach, cells are modeled as compartmentalized discrete objects, differentiated in the nucleus and in the cytosolic region, while the migration chamber is composed of channels of different widths. We find that cell motile phenotype and velocity in open spaces (i.e., 2D flat surfaces or large channels) are not significantly influenced by cell elastic properties. On the contrary, the migratory behavior of cells within subcellular and subnuclear structures strongly relies on the deformability of the cytosol and of the nuclear cluster, respectively. Further, we characterize two migration dynamics: a stepwise way, characterized by fluctuations in cell length, within channels smaller than nucleus dimensions and a smooth sliding (i.e., maintaining constant cell length) behavior within channels larger than the nuclear cluster. These resulting observations are then extended looking at cell migration in an artificial fiber network, which mimics cell invasion in a 3D extracellular matrix. In particular, in this case, we analyze the effect of variations in elasticity of the nucleus on cell movement. In order to summarize, with our simulated migration assays, we demonstrate that the dimensionality of the environment strongly affects the migration phenotype and we suggest that the cytoskeletal and nuclear elastic characteristics correlate with the tumor cell's invasive potentia

Comprehensive 4D velocity mapping of the heart and great vessels by cardiovascular magnetic resonance

<p>Abstract</p> <p>Background</p> <p>Phase contrast cardiovascular magnetic resonance (CMR) is able to measure all three directional components of the velocities of blood flow relative to the three spatial dimensions and the time course of the heart cycle. In this article, methods used for the acquisition, visualization, and quantification of such datasets are reviewed and illustrated.</p> <p>Methods</p> <p>Currently, the acquisition of 3D cine (4D) phase contrast velocity data, synchronized relative to both cardiac and respiratory movements takes about ten minutes or more, even when using parallel imaging and optimized pulse sequence design. The large resulting datasets need appropriate post processing for the visualization of multidirectional flow, for example as vector fields, pathlines or streamlines, or for retrospective volumetric quantification.</p> <p>Applications</p> <p>Multidirectional velocity acquisitions have provided 3D visualization of large scale flow features of the healthy heart and great vessels, and have shown altered patterns of flow in abnormal chambers and vessels. Clinically relevant examples include retrograde streams in atheromatous descending aortas as potential thrombo-embolic pathways in patients with cryptogenic stroke and marked variations of flow visualized in common aortic pathologies. Compared to standard clinical tools, 4D velocity mapping offers the potential for retrospective quantification of flow and other hemodynamic parameters.</p> <p>Conclusions</p> <p>Multidirectional, 3D cine velocity acquisitions are contributing to the understanding of normal and pathologically altered blood flow features. Although more rapid and user-friendly strategies for acquisition and analysis may be needed before 4D velocity acquisitions come to be adopted in routine clinical CMR, their capacity to measure multidirectional flows throughout a study volume has contributed novel insights into cardiovascular fluid dynamics in health and disease.</p