Diastolic right ventricular filling vortex in normal and volume overload states

Functional imaging computational fluid dynamics simulations of right ventricular (RV) inflow fields were obtained by comprehensive software using individual animal-specific dynamic imaging data input from three-dimensional (3-D) real-time echocardiography (RT3D) on a CRAY T-90 supercomputer. Chronically instrumented, lightly sedated awake dogs (n = 7) with normal wall motion (NWM) at control and normal or diastolic paradoxical septal motion (PSM) during RV volume overload were investigated. Up to the E-wave peak, instantaneous inflow streamlines extended from the tricuspid orifice to the RV endocardial surface in an expanding fanlike pattern. During the descending limb of the E-wave, large-scale (macroscopic or global) vortical motions ensued within the filling RV chamber. Both at control and during RV volume overload (with or without PSM), blood streams rolled up from regions near the walls toward the base. The extent and strength of the ring vortex surrounding the main stream were reduced with chamber dilatation. A hypothesis is proposed for a facilitatory role of the diastolic vortex for ventricular filling. The filling vortex supports filling by shunting inflow kinetic energy, which would otherwise contribute to an inflow-impeding convective pressure rise between inflow orifice and the large endocardial surface of the expanding chamber, into the rotational kinetic energy of the vortical motion that is destined to be dissipated as heat. The basic information presented should improve application and interpretation of noninvasive (Doppler color flow mapping, velocity-encoded cine magnetic resonance imaging, etc.) diastolic diagnostic studies and lead to improved understanding and recognition of subtle, flow-associated abnormalities in ventricular dilatation and remodeling.     

RV instantaneous intraventricular diastolic pressure and velocity distributions in normal and volume overload awake dog disease models

Intraventricular diastolic right ventricular (RV) flow field dynamics were studied by functional imaging using three-dimensional (3D) real-time echocardiography with sonomicrometry and computational fluid dynamics in seven awake dogs at control with normal wall motion (NWM) and RV volume overload with diastolic paradoxical septal motion. Burgeoning flow cross section between inflow anulus and chamber walls induces a convective pressure rise, which represents a "convective deceleration load" (CDL). High spatiotemporal resolution dynamic pressure and velocity distributions of the intraventricular RV flow field revealed time-dependent, subtle interactions between intraventricular local acceleration and convective pressure gradients. During the E-wave upstroke, the total pressure gradient along intraventricular flow is the algebraic sum of a pressure decrease contributed by local acceleration and a pressure rise contributed by a convective deceleration that partially counterbalances the local acceleration gradient. This underlies the smallness of early diastolic intraventricular gradients. At peak volumetric inflow, local acceleration vanishes and the total adverse intraventricular gradient is convective. During the E-wave downstroke, the strongly adverse gradient embodies the streamwise pressure augmentations from both local and convective decelerations. It induces flow separation and large-scale vortical motions, stronger in NWM. Their dynamic corollaries on intraventricular pressure and velocity distributions were ascertained. In the NWM pattern, the strong ring-like vortex surrounding the central core encroaches on the area available for flow toward the apex. This results in higher linear velocities later in the downstroke of the E wave than at peak inflow rate. The augmentation of CDL by ventriculoannular disproportion may contribute to E wave and E-to-A ratio depression with chamber dilatation.     

Fluid mechanics of the zebrafish embryonic heart trabeculation

Embryonic heart development is a mechanosensitive process, where specific fluid forces are needed for the correct development, and abnormal mechanical stimuli can lead to malformations. It is thus important to understand the nature of embryonic heart fluid forces. However, the fluid dynamical behaviour close to the embryonic endocardial surface is very sensitive to the geometry and motion dynamics of fine-scale cardiac trabecular surface structures. Here, we conducted image-based computational fluid dynamics (CFD) simulations to quantify the fluid mechanics associated with the zebrafish embryonic heart trabeculae. To capture trabecular geometric and motion details, we used a fish line that expresses fluorescence at the endocardial cell membrane, and high resolution 3D confocal microscopy. Our endocardial wall shear stress (WSS) results were found to exceed those reported in existing literature, which were estimated using myocardial rather than endocardial boundaries. By conducting simulations of single intra-trabecular spaces under varied scenarios, where the translational or deformational motions (caused by contraction) were removed, we found that a squeeze flow effect was responsible for most of the WSS magnitude in the intra-trabecular spaces, rather than the shear interaction with the flow in the main ventricular chamber. We found that trabecular structures were responsible for the high spatial variability of the magnitude and oscillatory nature of WSS, and for reducing the endocardial deformational burden. We further found cells attached to the endocardium within the intra-trabecular spaces, which were likely embryonic hemogenic cells, whose presence increased endocardial WSS. Overall, our results suggested that a complex multi-component consideration of both anatomic features and motion dynamics were needed to quantify the trabeculated embryonic heart fluid mechanics.     

Effects of wall motion and compliance on flow patterns in the ascending aorta

Helical flows have been observed in the ascending aorta in vivo, and geometric curvature has been suggested to be a major contributing factor. We employed magnetic resonance imaging (MRI) and velocity mapping to develop a computational model to examine the effects of curvature and also wall compliance and movement upon flow patterns. In the computational model, MRI-derived geometry and velocities were imposed as boundary conditions, which included both radial expansion-contraction and translational motion of the wall. The computed results were in agreement with the MR data only when full wall motion was included in the model, suggesting that the flow patterns observed in the ascending aorta arise not only from geometric curvature of the arch but also from the motion of the aorta resulting from its attachment to the beating heart.     

Right ventricular function and three-dimensional modeling using computer-aided design

Right ventricular (RV) volumetric and morphological analysis is complicated by the trabeculations and geometric configuration of the RV chamber. To improve RV analysis, custom computer-aided design programs were employed to obtain RV volumes and three-dimensional models from biplane ventriculograms. Biplane RV ventriculograms were analyzed from 14 anesthetized dogs and 22 RV casts. Computed volumes were highly correlated with reference RV volumes (r = 0.98, n = 36, P less than 0.01) with a range of 5-73 ml. Three-dimensional wire-frame and solid models constructed from the ventriculographic images provided excellent detail and a new perspective in chamber shape. This modeling technique was then used to examine RV volumes, geometric conformation, and regional shortening in 10 pigs during inotropic stimulation and preload reduction. Changes in RV volumes, ejection fraction, and regional motion were detected as well as alterations in chamber conformation. In summary 1) computer-aided design offers an accurate and simplified means to compute RV volumes using basic microcomputer equipment, and 2) three-dimensional reconstruction provided a unique view of RV geometry and a means to examine regional RV function.     

Accuracy and reproducibility of CFD predicted wall shear stress using 3D ultrasound images

Computational fluid dynamics (CFD) flow simulation techniques have the potential to enhance our understanding of how haemodynamic factors are involved in atherosclerosis. Recently, 3D ultrasound has emerged as an alternative to other 3D imaging techniques, such as magnetic resonance angiography (MRA). The method can be used to generate realistic vascular geometry suitable for CFD simulations. In order to assess accuracy and reproducibility of the procedure from image acquisition to reconstruction to CFD simulation, a human carotid artery bifurcation phantom was scanned three times using 3D ultrasound. The geometry was reconstructed and flow simulations were carried out on the three sets as well as on a model generated using computer aided design (CAD) from the geometric information given by the manufacturer. It was found that the three reconstructed sets showed good reproducibility as well as satisfactory quantitative agreement with the CAD model. Analyzing two selected locations probably representing the 'worst cases,' accuracy comparing ultrasound and CAD reconstructed models was estimated to be between 7.2% and 7.7% of the maximum instantaneous WSS and reproducibility comparing the three scans to be between 8.2% and 10.7% of their average maximum.     

Quantitative detection of cartilage surfaces and ligament geometry of the wrist using an imaging cryomicrotome system

Biomechanical models may aid in improving diagnosis and treatment of wrist joint disorders. As input, geometrical information is required for model development. Previous studies acquired some elements of the average wrist joint geometry. However, there is a close geometric functional match between articulating surfaces and ligament geometry. Therefore, biomechanical models need to be fed with the geometric data of individual joints. This study is aimed at acquiring geometric data of cartilage surfaces and ligaments from individual wrist joints by using a cryomicrotome imaging system and the evaluation of inter- and intra-observer variability of the data.The 3D geometry of 30 cartilage surfaces and 15 ligaments in three cadaver wrists was manually detected and quantitatively reconstructed. The inter- and intra-observer variability of the cartilage surface detection was 0.14 and 0.19 mm, respectively. For the position of the radius attachment of the dorsal radiocarpal ligament (DRC), the observer variations were 0.12 and 0.65 mm, for intra-/inter-observer, respectively. For the DRC attachment on the triquetrum, the observer variations were 0.22 and 1.19 mm.Anatomic reconstruction from 3D cryomicrotome images offer a method to obtain unique geometry data of the entire wrist joint for modeling purposes.     

A three-dimensional multiscale model for the prediction of thrombus growth under flow with single-platelet resolution

Modeling thrombus growth in pathological flows allows evaluation of risk under patient-specific pharmacological, hematological, and hemodynamical conditions. We have developed a 3D multiscale framework for the prediction of thrombus growth under flow on a spatially resolved surface presenting collagen and tissue factor (TF). The multiscale framework is composed of four coupled modules: a Neural Network (NN) that accounts for platelet signaling, a Lattice Kinetic Monte Carlo (LKMC) simulation for tracking platelet positions, a Finite Volume Method (FVM) simulator for solving convection-diffusion-reaction equations describing agonist release and transport, and a Lattice Boltzmann (LB) flow solver for computing the blood flow field over the growing thrombus. A reduced model of the coagulation cascade was embedded into the framework to account for TF-driven thrombin production. The 3D model was first tested against in vitro microfluidics experiments of whole blood perfusion with various antiplatelet agents targeting COX-1, P₂Y₁, or the IP receptor. The model was able to accurately capture the evolution and morphology of the growing thrombus. Certain problems of 2D models for thrombus growth (artifactual dendritic growth) were naturally avoided with realistic trajectories of platelets in 3D flow. The generalizability of the 3D multiscale solver enabled simulations of important clinical situations, such as cylindrical blood vessels and acute flow narrowing (stenosis). Enhanced platelet-platelet bonding at pathologically high shear rates (e.g., von Willebrand factor unfolding) was required for accurately describing thrombus growth in stenotic flows. Overall, the approach allows consideration of patient-specific platelet signaling and vascular geometry for the prediction of thrombotic episodes.     

A Novel Application of Musculoskeletal Ultrasound Imaging

Ultrasound is an attractive modality for imaging muscle and tendon motion during dynamic tasks and can provide a complementary methodological approach for biomechanical studies in a clinical or laboratory setting. Towards this goal, methods for quantification of muscle kinematics from ultrasound imagery are being developed based on image processing. The temporal resolution of these methods is typically not sufficient for highly dynamic tasks, such as drop-landing. We propose a new approach that utilizes a Doppler method for quantifying muscle kinematics. We have developed a novel vector tissue Doppler imaging (vTDI) technique that can be used to measure musculoskeletal contraction velocity, strain and strain rate with sub-millisecond temporal resolution during dynamic activities using ultrasound. The goal of this preliminary study was to investigate the repeatability and potential applicability of the vTDI technique in measuring musculoskeletal velocities during a drop-landing task, in healthy subjects. The vTDI measurements can be performed concurrently with other biomechanical techniques, such as 3D motion capture for joint kinematics and kinetics, electromyography for timing of muscle activation and force plates for ground reaction force. Integration of these complementary techniques could lead to a better understanding of dynamic muscle function and dysfunction underlying the pathogenesis and pathophysiology of musculoskeletal disorders.     

Numerical simulation of blood flow in the left ventricle and aortic sinus using magnetic resonance imaging and computational fluid dynamics

Understanding cardiac blood flow patterns has many applications in analysing haemodynamics and for the clinical assessment of heart function. In this study, numerical simulations of blood flow in a patient-specific anatomical model of the left ventricle (LV) and the aortic sinus are presented. The realistic 3D geometry of both LV and aortic sinus is extracted from the processing of magnetic resonance imaging (MRI). Furthermore, motion of inner walls of LV and aortic sinus is obtained from cine-MR image analysis and is used as a constraint to a numerical computational fluid dynamics (CFD) model based on the moving boundary approach. Arbitrary Lagrangian–Eulerian finite element method formulation is used for the numerical solution of the transient dynamic equations of the fluid domain. Simulation results include detailed flow characteristics such as velocity, pressure and wall shear stress for the whole domain. The aortic outflow is compared with data obtained by phase-contrast MRI. Good agreement was found between simulation results and these measurements.     

Noninvasive reconstruction of the three-dimensional ventricular activation sequence during pacing and ventricular tachycardia in the canine heart

Single-beat imaging of myocardial activation promises to aid in both cardiovascular research and clinical medicine. In the present study we validate a three-dimensional (3D) cardiac electrical imaging (3DCEI) technique with the aid of simultaneous 3D intracardiac mapping to assess its capability to localize endocardial and epicardial initiation sites and image global activation sequences during pacing and ventricular tachycardia (VT) in the canine heart. Body surface potentials were measured simultaneously with bipolar electrical recordings in a closed-chest condition in healthy canines. Computed tomography images were obtained after the mapping study to construct realistic geometry models. Data analysis was performed on paced rhythms and VTs induced by norepinephrine (NE). The noninvasively reconstructed activation sequence was in good agreement with the simultaneous measurements from 3D cardiac mapping with a correlation coefficient of 0.74 ± 0.06, a relative error of 0.29 ± 0.05, and a root mean square error of 9 ± 3 ms averaged over 460 paced beats and 96 ectopic beats including premature ventricular complexes, couplets, and nonsustained monomorphic VTs and polymorphic VTs. Endocardial and epicardial origins of paced beats were successfully predicted in 72% and 86% of cases, respectively, during left ventricular pacing. The NE-induced ectopic beats initiated in the subendocardium by a focal mechanism. Sites of initial activation were estimated to be ～7 mm from the measured initiation sites for both the paced beats and ectopic beats. For the polymorphic VTs, beat-to-beat dynamic shifts of initiation site and activation pattern were characterized by the reconstruction. The present results suggest that 3DCEI can noninvasively image the 3D activation sequence and localize the origin of activation of paced beats and NE-induced VTs in the canine heart with good accuracy. This 3DCEI technique offers the potential to aid interventional therapeutic procedures for treating ventricular arrhythmias arising from epicardial or endocardial sites and to noninvasively assess the mechanisms of these arrhythmias.     

Optimisation of a Stirred Bioreactor through the Use of a Novel Holographic Correlation Velocimetry Flow Measurement Technique

We describe a method for measuring three dimensional (3D) velocity fields of a fluid at high speed, by combining a correlation-based approach with in-line holography. While this method utilizes tracer particles contained within the flow, our method does not require the holographic reconstruction of 3D images. The direct flow reconstruction approach developed here allows for measurements at seeding densities in excess of the allowable levels for techniques based on image or particle reconstruction, thus making it suited for biological flow measurement, such as the flow in bioreactor. We outline the theory behind our method, which we term Holographic Correlation Velocimetry (HCV), and subsequently apply it to both synthetic and laboratory data. Moreover, because the system is based on in-line holography, it is very efficient with regard to the use of light, as it does not rely on side scattering. This efficiency could be utilized to create a very high quality system at a modest cost. Alternatively, this efficiency makes the system appropriate for high-speed flows and low exposure times, which is essential for imaging dynamic systems.     

Altered Right Ventricular Kinetic Energy Work Density and Viscous Energy Dissipation in Patients with Pulmonary Arterial Hypertension: A Pilot Study Using 4D Flow MRI

Introduction

Right ventricular (RV) function has increasingly being recognized as an important predictor for morbidity and mortality in patients with pulmonary arterial hypertension (PAH). The increased RV after-load increase RV work in PAH. We used time-resolved 3D phase contrast MRI (4D flow MRI) to derive RV kinetic energy (KE) work density and energy loss in the pulmonary artery (PA) to better characterize RV work in PAH patients.

Methods

4D flow and standard cardiac cine images were obtained in ten functional class I/II patients with PAH and nine healthy subjects. For each individual, we calculated the RV KE work density and the amount of viscous dissipation in the PA.

Results

PAH patients had alterations in flow patterns in both the RV and the PA compared to healthy subjects. PAH subjects had significantly higher RV KE work density than healthy subjects (94.7±33.7 mJ/mL vs. 61.7±14.8 mJ/mL, p = 0.007) as well as a much greater percent PA energy loss (21.1±6.4% vs. 2.2±1.3%, p = 0.0001) throughout the cardiac cycle. RV KE work density and percent PA energy loss had mild and moderate correlations with RV ejection fraction.

Conclusion

This study has quantified two kinetic energy metrics to assess RV function using 4D flow. RV KE work density and PA viscous energy loss not only distinguished healthy subjects from patients, but also provided distinction amongst PAH patients. These metrics hold promise as imaging markers for RV function.     

Quantification of left and right ventricular kinetic energy using four-dimensional intracardiac magnetic resonance imaging flow measurements

We aimed to quantify kinetic energy (KE) during the entire cardiac cycle of the left ventricle (LV) and right ventricle (RV) using four-dimensional phase-contrast magnetic resonance imaging (MRI). KE was quantified in healthy volunteers (n = 9) using an in-house developed software. Mean KE through the cardiac cycle of the LV and the RV were highly correlated (r(2) = 0.96). Mean KE was related to end-diastolic volume (r(2) = 0.66 for LV and r(2) = 0.74 for RV), end-systolic volume (r(2) = 0.59 and 0.68), and stroke volume (r(2) = 0.55 and 0.60), but not to ejection fraction (r(2) < 0.01, P = not significant for both). Three KE peaks were found in both ventricles, in systole, early diastole, and late diastole. In systole, peak KE in the LV was lower (4.9 ± 0.4 mJ, P = 0.004) compared with the RV (7.5 ± 0.8 mJ). In contrast, KE during early diastole was higher in the LV (6.0 ± 0.6 mJ, P = 0.004) compared with the RV (3.6 ± 0.4 mJ). The late diastolic peaks were smaller than the systolic and early diastolic peaks (1.3 ± 0.2 and 1.2 ± 0.2 mJ). Modeling estimated the proportion of KE to total external work, which comprised ～0.3% of LV external work and 3% of RV energy at rest and 3 vs. 24% during peak exercise. The higher early diastolic KE in the LV indicates that LV filling is more dependent on ventricular suction compared with the RV. RV early diastolic filling, on the other hand, may be caused to a higher degree of the return of the atrioventricular plane toward the base of the heart. The difference in ventricular geometry with a longer outflow tract in the RV compared with the LV explains the higher systolic KE in the RV.     

3D Fluid-Structure Interaction Canine Heart Model with Patch to Quantify Mechanical Conditions for Optimal Myocardium Stem Cell Growth and Tissue Regeneration

Right ventricular (RV) dysfunction is a common cause of heart failure in patients with congenital heart defects and often leads to impaired functional capacity and premature death. Myocardial tissue regeneration techniques are being developed for the potential that viable myocardium may be regenerated to replace scar tissues in the heart or used as patch material in heart surgery. 3D computational RV/LV/Patch models with fluid-structure interactions (FSI) were constructed based on data from a healthy dog heart to obtain local fluid dynamics and structural stress/strain information and identify optimal conditions under which tissue regeneration techniques could achieve best outcome. RV/LV/Patch geometry and blood pressure data were obtained from a dog following established procedures. Four FSI models were used to quantify the influence of different patch materials (Dacron scaffold, treated pericardium) on local environment around the patch area, especially focusing on the thickness and stiffness of the patch. Our results indicated that changes in patch stiffness had little impact on the ejection fraction of the right ventricle because the total patch area was small. However, patch stiffness had huge impact on local RV maximum principal stress (Stress-P1) and strain (Strain-P1) around the patch area. Compared to the no-patch model, patch models had increased Stress-P1 and decreased Strain-P1 values in the patch area. Softer patches were associated with greater stress/strain variations. Thinner patch led to complex local flow environment which may have impact on myocytes seeding and RV remodeling. Our multi-physics RV/LV/Patch FSI model can serve as a useful tool to investigate cellular biology and tissue regeneration under localized flow and structural stress environment.     

Measurement of hemodynamics in human carotid artery using ultrasound and computational fluid dynamics.

The objective of the study was to investigate the feasibility of using computational fluid dynamic modeling (CFD) with noninvasive ultrasound measurements to determine time-variant three-dimensional (3D) carotid arterial hemodynamics in humans in vivo. The effects of hyperoxia and hypoxic hypercapnia on carotid artery local hemodynamics were examined by use of this approach. Six normotensive volunteers followed a double-blind randomized crossover design. Blood pressure, heart rate, and carotid blood flow were measured while subjects breathed normal air, a mixture of 5% CO(2) and 15% O(2) (hypoxic hypercapnia), and 100% O(2) (hyperoxia). Carotid artery geometry was reconstructed on the basis of B-mode ultrasound images by using purpose-built image processing software. Time-variant 3D carotid hemodynamics were estimated by using finite volume-based CFD. Systemic blood pressure was not significantly affected by hyperoxia or hypoxic hypercapnia, but heart rate decreased significantly with hyperoxia. There was an increase in diastolic flow velocity in the external carotid artery after hypoxic hypercapnia, but otherwise carotid blood flow velocities did not change significantly. Compared with normal air, hyperoxic conditions were associated with a decrease in the width of the region of flow separation in the external carotid artery. During hyperoxia, there was also an increase in the minimum and a decrease in maximum shear stress in the bifurcation and hence a reduction in cyclic variation in shear stress. Hypoxic hypercapnia was associated with a reduced duration of flow separation in the external carotid artery and an increase in the minimum shear stress without affecting the cyclic variation in shear stress. This study demonstrates the feasibility of using noninvasive ultrasound techniques in conjunction with CFD to describe time-variant 3D hemodynamics in the human carotid arterial bifurcation in vivo.     

Combined pulmonary stenosis and insufficiency preserves myocardial contractility in the developing heart of growing swine at midterm follow-up.

This study was conducted to determine the effects of chronic combined pulmonary stenosis and pulmonary insufficiency (PSPI) on right (RV) and left ventricular (LV) function in young, growing swine. Six pigs with combined PSPI were studied, and data were compared with previously published data of animals with isolated pulmonary insufficiency and controls. Indexes of systolic function (stroke volume, ejection fraction, and cardiac functional reserve), myocardial contractility (slope of the end-systolic pressure-volume and change in pressure over time-end-diastolic volume relationship), and diastolic compliance were assessed within 2 days of intervention and 3 mo later. Magnetic resonance imaging was used to quantify pulmonary insufficiency and ventricular volumes. The conductance catheter was used to obtain indexes of the cardiac functional reserve, diastolic compliance, and myocardial contractility from pressure-volume relations acquired at rest and under dobutamine infusion. In the PSPI group, the pulmonary regurgitant fraction was 34.3 +/- 5.8%, the pressure gradient across the site of pulmonary stenosis was 20.9 +/- 20 mmHg, and the average RV peak systolic pressure was 70% systemic at 12 wk follow-up. Biventricular resting cardiac outputs and cardiac functional reserves were significantly limited (P < 0.05), LV diastolic compliance significantly decreased (P < 0.05), but RV myocardial contractility significantly enhanced (P < 0.05) compared with control animals at 3-mo follow-up. In the young, developing heart, chronic combined PSPI impairs biventricular systolic pump function and diastolic compliance but preserves RV myocardial contractility.     

Optimization of shunt placement for the Norwood surgery using multi-domain modeling

An idealized systemic-to-pulmonary shunt anatomy is parameterized and coupled to a closed loop, lumped parameter network (LPN) in a multidomain model of the Norwood surgical anatomy. The LPN approach is essential for obtaining information on global changes in cardiac output and oxygen delivery resulting from changes in local geometry and physiology. The LPN is fully coupled to a custom 3D finite element solver using a semi-implicit approach to model the heart and downstream circulation. This closed loop multidomain model is then integrated with a fully automated derivative-free optimization algorithm to obtain optimal shunt geometries with variable parameters of shunt diameter, anastomosis location, and angles. Three objective functions: (1) systemic; (2) coronary; and (3) combined systemic and coronary oxygen deliveries are maximized. Results show that a smaller shunt diameter with a distal shunt-brachiocephalic anastomosis is optimal for systemic oxygen delivery, whereas a more proximal anastomosis is optimal for coronary oxygen delivery and a shunt between these two anatomies is optimal for both systemic and coronary oxygen deliveries. Results are used to quantify the origin of blood flow going through the shunt and its relationship with shunt geometry. Results show that coronary artery flow is directly related to shunt position.     

From Dynamic Live Cell Imaging to 3D Ultrastructure: Novel Integrated Methods for High Pressure Freezing and Correlative Light-Electron Microscopy

Background

In cell biology, the study of proteins and organelles requires the combination of different imaging approaches, from live recordings with light microscopy (LM) to electron microscopy (EM).

Methodology

To correlate dynamic events in adherent cells with both ultrastructural and 3D information, we developed a method for cultured cells that combines confocal time-lapse images of GFP-tagged proteins with electron microscopy. With laser micro-patterned culture substrate, we created coordinates that were conserved at every step of the sample preparation and visualization processes. Specifically designed for cryo-fixation, this method allowed a fast freezing of dynamic events within seconds and their ultrastructural characterization. We provide examples of the dynamic oligomerization of GFP-tagged myotubularin (MTM1) phosphoinositides phosphatase induced by osmotic stress, and of the ultrastructure of membrane tubules dependent on amphiphysin 2 (BIN1) expression.

Conclusion

Accessible and versatile, we show that this approach is efficient to routinely correlate functional and dynamic LM with high resolution morphology by EM, with immuno-EM labeling, with 3D reconstruction using serial immuno-EM or tomography, and with scanning-EM.