Automated 3D mapping & shape analysis of the lateral ventricles via fluid registration of multiple surface-based atlases

We developed and validated a new method to create automated 3D parametric surface models of the lateral ventricles, designed for monitoring degenerative disease effects in clinical neuroscience studies and drug trials. First we used a set of parameterized surfaces to represent the ventricles in a manually labeled set of 9 subjects' MRIs (atlases). We fluidly registered each of these atlases and mesh models to a set of MRIs from 12 Alzheimer's disease (AD) patients and 14 matched healthy elderly subjects, and we averaged the resulting meshes for each of these images. Validation experiments on expert segmentations showed that (1) the Hausdorff labeling error rapidly decreased, and (2) the power to detect disease-related alterations monotonically improved as the number of atlases, N, was increased from 1 to 9. We then combined the segmentations with a radial mapping approach to localize ventricular shape differences in patients. In surface-based statistical maps, we detected more widespread and intense anatomical deficits as we increased the number of atlases, and we formulated a statistical stopping criterion to determine the optimal value of N. Anterior horn anomalies in Alzheimer's patients were only detected with the multi-atlas segmentation, which clearly outperformed the standard single-atlas approach.

Ultrastructural alterations of choroid plexuses of lateral ventricles of rats (rattus norvegicus) submitted to experimental chronic alcoholism

Adult male rats (Wistar lineage) were alcoholized with sugar cane liquor diluted at 30° GL during 300 days and sacrificed every 60 days in 5 stages. Samples of choroid plexuses of lateral ventricles were collected and examined at transmission electronic microscope to detect possible ultrastructural alterations and to raise possible pathological correlations. Gradual changes were observed in these animals during all the experiment: dilatation and enlargement of cisternae of Golgi complex, dilatation of RER, presence of digestive vacuoles and a large amount of pinocytic vesicles as well as vesicles with electronlucent content throughout cytoplasm, as well as an enlargement of intercellular space between basolateral interdigitation of the cells and of the connective tissue. The changes observed in the epithelium and connective tissue of choroid plexuses specially in 240 and 300 days of treatment are presumably due to a disturbance in hydroelectrolitic homeostasis, contributing to several morpho-functional disturbs of central nervous system. No changes were observed in the control group animals.

Computational fluid dynamics of a centrifugal heart pump as biventricular assist device(BVAD)

The results of a recent study have shown that there is a severe shortage of donor hearts to meet the demand of patients suffering from acute heart failures, and patients who received a left ventricular assist device (LVAD) have extended lives. However, some of them develop right heart failure syndrome, and these patients required a right ventricular assist device (RVAD). Hence, current research focus is in the development of a bi-ventricular assist device (Bi-VAD). Computational Fluid Dynamics (CFD) is useful for estimating blood damage for design of a Bi-VAD centrifugal heart pump to meet the demand of the left and right ventricles of a normal hearts with a flow rate of 5 lit/min and the supply pressure of 100 mmHg for the left ventricle and 20 mmHg for the right ventricle. Numerical studies have been conducted to predict pressure, flow rate, the velocity profiles, and streamlines in a continuous flow Bi-VAD using. Based on the predictions of numerical simulations, only few flow regions in the Bi-VAD exhibited signs of velocity profiles and stagnation points, thereby signifying potentially low levels of thrombogenesis.

Passive control of a bi-ventricular assist device : an experimental and numerical investigation

For the last two decades heart disease has been the highest single cause of death for the human population. With an alarming number of patients requiring heart transplant, and donations not able to satisfy the demand, treatment looks to mechanical alternatives. Rotary Ventricular Assist Devices, VADs, are miniature pumps which can be implanted alongside the heart to assist its pumping function. These constant flow devices are smaller, more efficient and promise a longer operational life than more traditional pulsatile VADs. The development of rotary VADs has focused on single pumps assisting the left ventricle only to supply blood for the body. In many patients however, failure of both ventricles demands that an additional pulsatile device be used to support the failing right ventricle. This condition renders them hospital bound while they wait for an unlikely heart donation. Reported attempts to use two rotary pumps to support both ventricles concurrently have warned of inherent haemodynamic instability. Poor balancing of the pumps’ flow rates quickly leads to vascular congestion increasing the risk of oedema and ventricular ‘suckdown’ occluding the inlet to the pump. This thesis introduces a novel Bi-Ventricular Assist Device (BiVAD) configuration where the pump outputs are passively balanced by vascular pressure. The BiVAD consists of two rotary pumps straddling the mechanical passive controller. Fluctuations in vascular pressure induce small deflections within both pumps adjusting their outputs allowing them to maintain arterial pressure. To optimise the passive controller’s interaction with the circulation, the controller’s dynamic response is optimised with a spring, mass, damper arrangement. This two part study presents a comprehensive assessment of the prototype’s ‘viability’ as a support device. Its ‘viability’ was considered based on its sensitivity to pathogenic haemodynamics and the ability of the passive response to maintain healthy circulation. The first part of the study is an experimental investigation where a prototype device was designed and built, and then tested in a pulsatile mock circulation loop. The BiVAD was subjected to a range of haemodynamic imbalances as well as a dynamic analysis to assess the functionality of the mechanical damper. The second part introduces the development of a numerical program to simulate human circulation supported by the passively controlled BiVAD. Both investigations showed that the prototype was able to mimic the native baroreceptor response. Simulating hypertension, poor flow balancing and subsequent ventricular failure during BiVAD support allowed the passive controller’s response to be assessed. Triggered by the resulting pressure imbalance, the controller responded by passively adjusting the VAD outputs in order to maintain healthy arterial pressures. This baroreceptor-like response demonstrated the inherent stability of the auto regulating BiVAD prototype. Simulating pulmonary hypertension in the more observable numerical model, however, revealed a serious issue with the passive response. The subsequent decrease in venous return into the left heart went unnoticed by the passive controller. Meanwhile the coupled nature of the passive response not only decreased RVAD output to reduce pulmonary arterial pressure, but it also increased LVAD output. Consequently, the LVAD increased fluid evacuation from the left ventricle, LV, and so actually accelerated the onset of LV collapse. It was concluded that despite the inherently stable baroreceptor-like response of the passive controller, its lack of sensitivity to venous return made it unviable in its present configuration. The study revealed a number of other important findings. Perhaps the most significant was that the reduced pulse experienced during constant flow support unbalanced the ratio of effective resistances of both vascular circuits. Even during steady rotary support therefore, the resulting ventricle volume imbalance increased the likelihood of suckdown. Additionally, mechanical damping of the passive controller’s response successfully filtered out pressure fluctuations from residual ventricular function. Finally, the importance of recognising inertial contributions to blood flow in the atria and ventricles in a numerical simulation were highlighted. This thesis documents the first attempt to create a fully auto regulated rotary cardiac assist device. Initial results encourage development of an inlet configuration sensitive to low flow such as collapsible inlet cannulae. Combining this with the existing baroreceptor-like response of the passive controller will render a highly stable passively controlled BiVAD configuration. The prototype controller’s passive interaction with the vasculature is a significant step towards a highly stable new generation of artificial heart.

Simulation and development of a mock circulation loop with variable compliance

Heart disease is attributed as the highest cause of death in the world. Although this could be alleviated by heart transplantation, there is a chronic shortage of donor hearts and so mechanical solutions are being considered. Currently, many Ventricular Assist Devices (VADs) are being developed worldwide in an effort to increase life expectancy and quality of life for end stage heart failure patients. Current pre-clinical testing methods for VADs involve laboratory testing using Mock Circulation Loops (MCLs), and in vivo testing in animal models. The research and development of highly accurate MCLs is vital to the continuous improvement of VAD performance. The first objective of this study was to develop and validate a mathematical model of a MCL. This model could then be used in the design and construction of a variable compliance chamber to improve the performance of an existing MCL as well as form the basis for a new miniaturised MCL. An extensive review of literature was carried out on MCLs and mathematical modelling of their function. A mathematical model of a MCL was then created in the MATLAB/SIMULINK environment. This model included variable features such as resistance, fluid inertia and volumes (resulting from the pipe lengths and diameters); compliance of Windkessel chambers, atria and ventricles; density of both fluid and compressed air applied to the system; gravitational effects on vertical columns of fluid; and accurately modelled actuators controlling the ventricle contraction. This model was then validated using the physical properties and pressure and flow traces produced from a previously developed MCL. A variable compliance chamber was designed to reproduce parameters determined by the mathematical model. The function of the variability was achieved by controlling the transmural pressure across a diaphragm to alter the compliance of the system. An initial prototype was tested in a previously developed MCL, and a variable level of arterial compliance was successfully produced; however, the complete range of compliance values required for accurate physiological representation was not able to be produced with this initial design. The mathematical model was then used to design a smaller physical mock circulation loop, with the tubing sizes adjusted to produce accurate pressure and flow traces whilst having an appropriate frequency response characteristic. The development of the mathematical model greatly assisted the general design of an in vitro cardiovascular device test rig, while the variable compliance chamber allowed simple and real-time manipulation of MCL compliance to allow accurate transition between a variety of physiological conditions. The newly developed MCL produced an accurate design of a mechanical representation of the human circulatory system for in vitro cardiovascular device testing and education purposes. The continued improvement of VAD test rigs is essential if VAD design is to improve, and hence improve quality of life and life expectancy for heart failure patients.

P2 purinergic receptor mRNA in rat and human sinoatrial node and other heart regions.

It is known that adenosine 5'-triphosphate (ATP) is a cotransmitter in the heart. Additionally, ATP is released from ischemic and hypoxic myocytes. Therefore, cardiac-derived sources of ATP have the potential to modify cardiac function. ATP activates P2X(1-7) and P2Y(1-14) receptors; however, the presence of P2X and P2Y receptor subtypes in strategic cardiac locations such as the sinoatrial node has not been determined. An understanding of P2X and P2Y receptor localization would facilitate investigation of purine receptor function in the heart. Therefore, we used quantitative PCR and in situ hybridization to measure the expression of mRNA of all known purine receptors in rat left ventricle, right atrium and sinoatrial node (SAN), and human right atrium and SAN. Expression of mRNA for all the cloned P2 receptors was observed in the ventricles, atria, and SAN of the rat. However, their abundance varied in different regions of the heart. P2X(5) was the most abundant of the P2X receptors in all three regions of the rat heart. In rat left ventricle, P2Y(1), P2Y(2), and P2Y(14) mRNA levels were highest for P2Y receptors, while in right atrium and SAN, P2Y(2) and P2Y(14) levels were highest, respectively. We extended these studies to investigate P2X(4) receptor mRNA in heart from rats with coronary artery ligation-induced heart failure. P2X(4) receptor mRNA was upregulated by 93% in SAN (P < 0.05), while a trend towards an increase was also observed in the right atrium and left ventricle (not significant). Thus, P2X(4)-mediated effects might be modulated in heart failure. mRNA for P2X(4-7) and P2Y(1,2,4,6,12-14), but not P2X(2,3) and P2Y(11), was detected in human right atrium and SAN. In addition, mRNA for P2X(1) was detected in human SAN but not human right atrium. In human right atrium and SAN, P2X(4) and P2X(7) mRNA was the highest for P2X receptors. P2Y(1) and P2Y(2) mRNA were the most abundant for P2Y receptors in the right atrium, while P2Y(1), P2Y(2), and P2Y(14) were the most abundant P2Y receptor subtypes in human SAN. This study shows a widespread distribution of P2 receptor mRNA in rat heart tissues but a more restricted presence and distribution of P2 receptor mRNA in human atrium and SAN. This study provides further direction for the elucidation of P2 receptor modulation of heart rate and contractility.

Replication of the Frank-Starling response in a Mock Circulation loop

Mock circulation loops (MCLs) are used to evaluate cardiovascular devices prior to in-vivo trials; however they lack the vital autoregulatory responses that occur in humans. This study aimed to develop and implement a left and right ventricular Frank-Starling response in a MCL. A proportional controller based on ventricular end diastolic volume was used to control the driving pressure of the MCL’s pneumatically operated ventricles. Ventricular pressure-volume loops and end systolic pressure-volume relationships were produced for a variety of healthy and pathological conditions and compared with human data to validate the simulated Frank-Starling response. The non-linear Frank-Starling response produced in this study successfully altered left and right ventricular contractility with changing preload and was validated with previously reported data. This improvement to an already detailed MCL has resulted in a test rig capable of further refining cardiovascular devices and reducing the number of in-vivo trials.

Phenomenological modeling of cell-to-cell and beat-to-beat variability in isolated Guinea Pig ventricular myocytes

Experimental action potential (AP) recordings in isolated ventricular myoctes display significant temporal beat-to-beat variability in morphology and duration. Furthermore, significant cell-to-cell differences in AP also exist even for isolated cells originating from the same region of the same heart. However, current mathematical models of ventricular AP fail to replicate the temporal and cell-to-cell variability in AP observed experimentally. In this study, we propose a novel mathematical framework for the development of phenomenological AP models capable of capturing cell-to-cell and temporal variabilty in cardiac APs. A novel stochastic phenomenological model of the AP is developed, based on the deterministic Bueno-Orovio/Fentonmodel. Experimental recordings of AP are fit to the model to produce AP models of individual cells from the apex and the base of the guinea-pig ventricles. Our results show that the phenomenological model is able to capture the considerable differences in AP recorded from isolated cells originating from the location. We demonstrate the closeness of fit to the available experimental data which may be achieved using a phenomenological model, and also demonstrate the ability of the stochastic form of the model to capture the observed beat-to-beat variablity in action potential duration.

Volume estimation from multiplanar 2D ultrasound images using a remote electromagnetic position and orientation sensor

A system is described for calculating volume from a sequence of multiplanar 2D ultrasound images. Ultrasound images are captured using a video digitising card (Hauppauge Win/TV card) installed in a personal computer, and regions of interest transformed into 3D space using position and orientation data obtained from an electromagnetic device (Polbemus, Fastrak). The accuracy of the system was assessed by scanning 10 water filled balloons (13-141 ml), 10 kidneys (147  200 ml) and 16 fetal livers (8  37 ml) in water using an Acuson 128XP/10 (5 MHz curvilinear probe). Volume was calculated using the ellipsoid, planimetry, tetrahedral and ray tracing methods and compared with the actual volume measured by weighing (balloons) and water displacement (kidneys and livers). The mean percentage error for the ray tracing method was 0.9 ± 2.4%, 2.7 ± 2.3%, 6.6 ± 5.4% for balloons, kidneys and livers, respectively. So far the system has been used clinically to scan fetal livers and lungs, neonate brain ventricles and adult prostate glands.

Inflow cannula design for biventricular assist devices

Cardiovascular diseases are a leading cause of death throughout the developed world. With the demand for donor hearts far exceeding the supply, a bridge-to-transplant or permanent solution is required. This is currently achieved with ventricular assist devices (VADs), which can be used to assist the left ventricle (LVAD), right ventricle (RVAD), or both ventricles simultaneously (BiVAD). Earlier generation VADs were large, volume-displacement devices designed for temporary support until a donor heart was found. The latest generation of VADs use rotary blood pump technology which improves device lifetime and the quality of life for end stage heart failure patients. VADs are connected to the heart and greater vessels of the patient through specially designed tubes called cannulae. The inflow cannulae, which supply blood to the VAD, are usually attached to the left atrium or ventricle for LVAD support, and the right atrium or ventricle for RVAD support. Few studies have characterized the haemodynamic difference between the two cannulation sites, particularly with respect to rotary RVAD support. Inflow cannulae are usually made of metal or a semi-rigid polymer to prevent collapse with negative pressures. However suction, and subsequent collapse, of the cannulated heart chamber can be a frequent occurrence, particularly with the relatively preload insensitive rotary blood pumps. Suction events may be associated with endocardial damage, pump flow stoppages and ventricular arrhythmias. While several VAD control strategies are under development, these usually rely on potentially inaccurate sensors or somewhat unreliable inferred data to estimate preload. Fixation of the inflow cannula is usually achieved through suturing the cannula, often via a felt sewing ring, to the cannulated chamber. This technique extends the time on cardiopulmonary bypass which is associated with several postoperative complications. The overall objective of this thesis was to improve the placement and design of rotary LVAD and RVAD inflow cannulae to achieve enhanced haemodynamic performance, reduced incidence of suction events, reduced levels of postoperative bleeding and a faster implantation procedure. Specific objectives were: * in-vitro evaluation of LVAD and RVAD inflow cannula placement, * design and in-vitro evaluation of a passive mechanism to reduce the potential for heart chamber suction, * design and in-vitro evaluation of a novel suture-less cannula fixation device. In order to complete in-vitro evaluation of VAD inflow cannulae, a mock circulation loop (MCL) was developed to accurately replicate the haemodynamics in the human systemic and pulmonary circulations. Validation of the MCL’s haemodynamic performance, including the form and magnitude of pressure, flow and volume traces was completed through comparisons of patient data and the literature. The MCL was capable of reproducing almost any healthy or pathological condition, and provided a useful tool to evaluate VAD cannulation and other cardiovascular devices. The MCL was used to evaluate inflow cannula placement for rotary VAD support. Left and right atrial and ventricular cannulation sites were evaluated under conditions of mild and severe heart failure. With a view to long term LVAD support in the severe left heart failure condition, left ventricular inflow cannulation was preferred due to improved LVAD efficiency and reduced potential for thrombus formation. In the mild left heart failure condition, left atrial cannulation was preferred to provide an improved platform for myocardial recovery. Similar trends were observed with RVAD support, however to a lesser degree due to a smaller difference in right atrial and ventricular pressures. A compliant inflow cannula to prevent suction events was then developed and evaluated in the MCL. As rotary LVAD or RVAD preload was reduced, suction events occurred in all instances with a rigid inflow cannula. Addition of the compliant segment eliminated suction events in all instances. This was due to passive restriction of the compliant segment as preload dropped, thus increasing the VAD circuit resistance and decreasing the VAD flow rate. Therefore, the compliant inflow cannula acted as a passive flow control / anti-suction system in LVAD and RVAD support. A novel suture-less inflow cannula fixation device was then developed to reduce implantation time and postoperative bleeding. The fixation device was evaluated for LVAD and RVAD support in cadaveric animal and human hearts attached to a MCL. LVAD inflow cannulation was achieved in under two minutes with the suture-less fixation device. No leakage through the suture-less fixation device – myocardial interface was noted. Continued development and in-vivo evaluation of this device may result in an improved inflow cannulation technique with the potential for off-bypass insertion. Continued development of this research, in particular the compliant inflow cannula and suture-less inflow cannulation device, will result in improved postoperative outcomes, life span and quality of life for end-stage heart failure patients.

Passive control of a biventricular assist device with compliant inflow cannulae

Rotary ventricular assist device (VAD) support of the cardiovascular system is susceptible to suction events due to the limited preload sensitivity of these devices. This may be of particular concern with rotary biventricular support (BiVAD) where the native, flow-balancing Starling response is diminished in both ventricles. The reliability of sensor and sensor-less based control systems which aim to control VAD flow based on preload have limitations and thus an alternative solution is desired. This study introduces a compliant inflow cannula (CIC) which could improve the preload sensitivity of a rotary VAD by passively altering VAD flow depending on preload. To evaluate the design, both the CIC and a standard rigid inflow cannula were inserted into a mock circulation loop to enable biventricular heart failure support using configurations of atrial and ventricular inflow, and arterial outflow cannulation. A range of left (LVAD) and right VAD (RVAD) rotational speeds were tested as well as step changes in systemic/pulmonary vascular resistance to alter relative preloads, with resulting flow rates recorded. Simulated suction events were observed, particularly at higher VAD speeds, during support with the rigid inflow cannula, while the CIC prevented suction events under all circumstances. The compliant section passively restricted its internal diameter as preload was reduced, which increased the VAD circuit resistance and thus reduced VAD flow. Therefore, a compliant inflow cannula could potentially be used as a passive control system to prevent suction events in rotary left, right and biventricular support.

Anatomic fitting of total artificial hearts for in vivo evaluation

Successful anatomic fitting of a total artificial heart (TAH) is vital to achieve optimal pump hemodynamics after device implantation. Although many anatomic fitting studies have been completed in humans prior to clinical trials, few reports exist that detail the experience in animals for in vivo device evaluation. Optimal hemodynamics are crucial throughout the in vivo phase to direct design iterations and ultimately validate device performance prior to pivotal human trials. In vivo evaluation in a sheep model allows a realistically sized representation of a smaller patient, for which smaller third-generation TAHs have the potential to treat. Our study aimed to assess the anatomic fit of a single device rotary TAH in sheep prior to animal trials and to use the data to develop a threedimensional, computer-aided design (CAD)-operated anatomic fitting tool for future TAH development. Following excision of the native ventricles above the atrio-ventricular groove, a prototype TAH was inserted within the chest cavity of six sheep (28–40 kg).Adjustable rods representing inlet and outlet conduits were oriented toward the center of each atrial chamber and the great vessels, with conduit lengths and angles recorded for future analysis. A threedimensional, CAD-operated anatomic fitting tool was then developed, based on the results of this study, and used to determine the inflow and outflow conduit orientation of the TAH. The mean diameters of the sheep left atrium, right atrium, aorta, and pulmonary artery were 39, 33, 12, and 11 mm, respectively. The center-to-center distance and outer-edge-to-outer-edge distance between the atria, found to be 39 ± 9 mm and 72 ± 17 mm in this study, were identified as the most critical geometries for successful TAH connection. This geometric constraint restricts the maximum separation allowable between left and right inlet ports of a TAH to ensure successful alignment within the available atrial circumference.

Quantitative genetic modeling of lateral ventricular shape and volume using multi-atlas fluid image alignment in twins

Despite substantial progress in measuring the 3D profile of anatomical variations in the human brain, their genetic and environmental causes remain enigmatic. We developed an automated system to identify and map genetic and environmental effects on brain structure in large brain MRI databases . We applied our multi-template segmentation approach ("Multi-Atlas Fluid Image Alignment") to fluidly propagate hand-labeled parameterized surface meshes into 116 scans of twins (60 identical, 56 fraternal), labeling the lateral ventricles. Mesh surfaces were averaged within subjects to minimize segmentation error. We fitted quantitative genetic models at each of 30,000 surface points to measure the proportion of shape variance attributable to (1) genetic differences among subjects, (2) environmental influences unique to each individual, and (3) shared environmental effects. Surface-based statistical maps revealed 3D heritability patterns, and their significance, with and without adjustments for global brain scale. These maps visualized detailed profiles of environmental versus genetic influences on the brain, extending genetic models to spatially detailed, automatically computed, 3D maps.