Development of procedures for the design, optimization and manufacturing of customized orthopaedic and trauma implants: Geometrical/anatomical modelling from 3D medical imaging

Tese de Doutoramento (Programa Doutoral em Engenharia Biomédica)

Assessment of Myocardial Infarction by Cardiac Magnetic Resonance Imaging and Long-Term Mortality

Background: Cardiac magnetic resonance imaging provides detailed anatomical information on infarction. However, few studies have investigated the association of these data with mortality after acute myocardial infarction. Objective: To study the association between data regarding infarct size and anatomy, as obtained from cardiac magnetic resonance imaging after acute myocardial infarction, and long-term mortality. Methods: A total of 1959 reports of “infarct size” were identified in 7119 cardiac magnetic resonance imaging studies, of which 420 had clinical and laboratory confirmation of previous myocardial infarction. The variables studied were the classic risk factors – left ventricular ejection fraction, categorized ventricular function, and location of acute myocardial infarction. Infarct size and acute myocardial infarction extent and transmurality were analyzed alone and together, using the variable named “MET-AMI”. The statistical analysis was carried out using the elastic net regularization, with the Cox model and survival trees. Results: The mean age was 62.3 ± 12 years, and 77.3% were males. During the mean follow-up of 6.4 ± 2.9 years, there were 76 deaths (18.1%). Serum creatinine, diabetes mellitus and previous myocardial infarction were independently associated with mortality. Age was the main explanatory factor. The cardiac magnetic resonance imaging variables independently associated with mortality were transmurality of acute myocardial infarction (p = 0.047), ventricular dysfunction (p = 0.0005) and infarcted size (p = 0.0005); the latter was the main explanatory variable for ischemic heart disease death. The MET-AMI variable was the most strongly associated with risk of ischemic heart disease death (HR: 16.04; 95%CI: 2.64-97.5; p = 0.003). Conclusion: The anatomical data of infarction, obtained from cardiac magnetic resonance imaging after acute myocardial infarction, were independently associated with long-term mortality, especially for ischemic heart disease death.

The Benefits of Prone SPECT Myocardial Perfusion Imaging in Reducing Both Artifact Defects and Patient Radiation Exposure

AbstractBackground:Prone imaging has been demonstrated to minimize diaphragmatic and breast tissue attenuation.Objectives:To determine the role of prone imaging on the reduction of unnecessary rest perfusion studies and coronary angiographies performed, thus decreasing investigation time and radiation exposure.Methods:We examined 139 patients, 120 with an inferior wall and 19 with an anterior wall perfusion defect that might represented attenuation artifact. Post-stress images were acquired in both the supine and prone position. Coronary angiography was used as the “gold standard” for evaluating coronary artery patency. The study was terminated and rest imaging was obviated in the presence of complete improvement of the defect in the prone position. Quantitative interpretation was performed. Results were compared with clinical data and coronary angiographic findings.Results:Prone acquisition correctly revealed defect improvement in 89 patients (89/120) with inferior wall and 12 patients (12/19) with anterior wall attenuation artifact. Quantitative analysis demonstrated statistically significant difference in the mean summed stress scores (SSS) of supine and mean SSS of prone studies in patients with disappearing inferior wall defect in the prone position and patent right coronary artery (true negative results). The mean difference between SSS in supine and in prone position was higher with disappearing than with remaining defects.Conclusion:Technetium-99m (Tc-99m) tetrofosmin myocardial perfusion imaging with the patient in the prone position overcomes soft tissue attenuation; moreover it provides an inexpensive, accurate approach to limit the number of unnecessary rest perfusion studies and coronary angiographies performed.

High-resolution magnetic resonance imaging quantitatively detects individual pancreatic islets.

OBJECTIVE-We studied whether manganese-enhanced high-field magnetic resonance (MR) imaging (MEHFMRI) could quantitatively detect individual islets in situ and in vivo and evaluate changes in a model of experimental diabetes.RESEARCH DESIGN AND METHODS-Whole pancreata from untreated (n = 3), MnCl(2) and glucose-injected mice (n = 6), and mice injected with either streptozotocin (STZ; n = 4) or citrate buffer (n = 4) were imaged ex vivo for unambiguous evaluation of islets. Exteriorized pancreata of MnCl(2) and glucose-injected mice (n = 6) were imaged in vivo to directly visualize the gland and minimize movements. In all cases, MR images were acquired in a 14.1 Testa scanner and correlated with the corresponding (immuno)histological sections.RESULTS-In ex vivo experiments, MEHFMRI distinguished different pancreatic tissues and evaluated the relative abundance of islets in the pancreata of normoglycemic mice. MEHFMRI also detected a significant decrease in the numerical and volume density of islets in STZ-injected mice. However, in the latter measurements the loss of beta-cells was undervalued under the conditions tested. The experiments on the externalized pancreata confirmed that MEHFMRI could visualize native individual islets in living, anesthetized mice.CONCLUSIONS-Data show that MEHFMRI quantitatively visualizes individual islets in the intact mouse pancreas, both ex vivo and in vivo. Diabetes 60:2853-2860, 2011

Improved temporal resolution for functional studies with reduced number of segments with three-dimensional echo planar imaging.

PURPOSE: To introduce a new k-space traversal strategy for segmented three-dimensional echo planar imaging (3D EPI) that encodes two partitions per radiofrequency excitation, effectively reducing the number excitations used to acquire a 3D EPI dataset by half. METHODS: The strategy was evaluated in the context of functional MRI applications for: image quality compared with segmented 3D EPI, temporal signal-to-noise ratio (tSNR) (the ability to detect resting state networks compared with multislice two-dimensional (2D) EPI and segmented 3D EPI, and temporal resolution (the ability to separate cardiac- and respiration-related fluctuations from the desired blood oxygen level-dependent signal of interest). RESULTS: Whole brain images with a nominal voxel size of 2 mm isotropic could be acquired with a temporal resolution under half a second using traditional parallel imaging acceleration up to 4× in the partition-encode direction and using novel data acquisition speed-up of 2× with a 32-channel coil. With 8× data acquisition speed-up in the partition-encode direction, 3D reduced excitations (RE)-EPI produced acceptable image quality without introduction of noticeable additional artifacts. Due to increased tSNR and better characterization of physiological fluctuations, the new strategy allowed detection of more resting state networks compared with multislice 2D-EPI and segmented 3D EPI. CONCLUSION: 3D RE-EPI resulted in significant increases in temporal resolution for whole brain acquisitions and in improved physiological noise characterization compared with 2D-EPI and segmented 3D EPI. Magn Reson Med 72:786-792, 2014. © 2013 Wiley Periodicals, Inc.

New Horizons in Forensic Imaging: Post Mortem CT Angiography - Benefits and Practical Application.

Puropse/Aim: To learn about the developement of post mortem CT angiography, its indications, benefits, pitfalls and practical application. Content Organization: A. Developement of post mortem CT angiography B. Technical prerequisites C. Practical application of post mortem CT angiography (preparation of the body, injection of contrast agent, examination protocol) D. Indications and benefits (including a comparison with conventional autopsy) E. Interpretation of imaging data (with case demonstrations) F. Artifacts, pitfalls and limitations G. Current and potential future use. Summary: This exhibit demonstrates the developement, application and interpretation of post mortem CT angiography. Teaching points: 1. post mortem CT angiography is feasible and useful for identification of the cause of death 2. depending on the indication it can be superior to autopsy 3. limitations and artifacts need to be known for interpreta

Imaging of experience-dependent structural plasticity in the mouse neocortex in vivo

The functionality of adult neocortical circuits can be altered by novel experiences or learning. This functional plasticity appears to rely on changes in the strength of neuronal connections that were established during development. Here we will describe some of our studies in which we have addressed whether structural changes, including the remodeling of axons and dendrites with synapse formation and elimination, could underlie experience-dependent plasticity in the adult neocortex. Using 2-photon laser-scanning microscopes and transgenic mice expressing GFP in a subset of pyramidal cells, we have observed that a small subset of dendritic spines continuously appear and disappear on a daily basis, whereas the majority of spines persists for months. Axonal boutons from different neuronal classes displayed similar behavior, although the extent of remodeling varied. Under baseline conditions, new spines in the barrel cortex were mostly transient and rarely survived for more than a week. However, when every other whisker was trimmed, the generation and loss of persistent spines was enhanced. Ultrastructural reconstruction of previously imaged spines and boutons showed that new spines slowly form synapses. New spines persisting for a few days always had synapses, whereas very young spines often lacked synapses. New synapses were predominantly found on large, multi-synapse boutons, suggesting that spine growth is followed by synapse formation, preferentially on existing boutons. Altogether our data indicate that novel sensory experience drives the stabilization of new spines on subclasses of cortical neurons and promotes the formation of new synapses. These synaptic changes likely underlie experience-dependent functional remodeling of specific neocortical circuits.

Cost evaluation of cardiac magnetic resonance imaging versus coronary angiography for the work-up of coronary artery disease

Background: CMR has recently emerged as a robust and reliable technique to assess coronary artery disease (CAD). A negative perfusion CMR test predicts low event rates of 0.3-0.5%/year. Invasive coronary angiography (CA) remains the "gold standard" for the evaluation of CAD in many countries.Objective: Assessing the costs of the two strategies in the European CMR registry for the work-up of known or suspected CAD from a health care payer perspective. Strategy 1) a CA to all patients or 2) a CA only to patients who are diagnosed positive for ischemia in a prior CMR.Method and results: Using data of the European CMR registry (20 hospitals, 11'040 consecutive patients) we calculated the proportion of patients who were diagnosed positive (20.6%), uncertain (6.5%), and negative (72.9%) after the CMR test in patients with known or suspected CAD (n=2'717). No other medical test was performed to patients who were negative for ischemia. Positive diagnosed patients had a coronary angiography. Those with uncertain diagnosis had additional tests (84.7%: stress echocardiography, 13.1%: CCT, 2.3% SPECT), these costs were added to the CMR strategy costs. Information from costs for tests in Germany and Switzerland were used. A sensibility analysis was performed for inpatient CA. For costs see figure. Results - costs.Discussion: The CMR strategy costs less than the CA strategy for the health insurance systems both, in Germany and Switzerland. While lower in costs, the CMR strategy is a non-invasive one, does not expose to radiation, and yields additional information on cardiac function, viability, valves, and great vessels. Developing the use of CMR instead of CA might imply some reduction in costs together with superior patient safety and comfort, and a better utilization of resources at the hospital level. Document introduit le : 01.12.2011

The borehole-fluid effect in electrical resistivity imaging

Fluid that fills boreholes in crosswell electrical resistivity investigations provides the necessary electrical contact between the electrodes and the rock formation but it is also the source of image artifacts in standard inversions that do not account for the effects of the boreholes. The image distortions can be severe for large resistivity contrasts between the rock formation and borehole fluid and for large borehole diameters. We have carried out 3D finite-element modeling using an unstructured-grid approach to quantify the magnitude of borehole effects for different resistivity contrasts, borehole diameters, and electrode configurations. Relatively common resistivity contrasts of 100:1 and borehole diameters of 10 and 20 cm yielded, for a bipole length of 5 m, apparent resistivity underestimates of approximately 12% and 32% when using AB-MN configurations and apparent resistivity overestimates of approximately 24% and 95% when using AM-BN configurations. Effects are generally more severe at shorter bipole spacings. We report the results obtained by either including or ignoring the boreholes in inversions of 3D field data from a test site in Switzerland, where approximately 10,000 crosswell resistivity-tomography measurements were made across six acquisition planes among four boreholes. Inversions of raw data that ignored the boreholes filled with low-resistivity fluid paradoxically produced high-resistivity artifacts around the boreholes. Including correction factors based on the modeling results fora ID model with and without the boreholes did not markedly improve the images. The only satisfactory approach was to use a 3D inversion code that explicitly incorporated the boreholes in the actual inversion. This new approach yielded an electrical resistivity image that was devoid of artifacts around the boreholes and that correlated well with coincident crosswell radar images.

Imaging molasse and quaternary sediments in Lake Geneva, Switzerland, with 3-D high-resolution seismic reflection methods: A case study

A high-resolution three-dimensional (3-D) seismic reflection survey was conducted in Lake Geneva, near the city of Lausanne, Switzerland, as part of a project for developing such seismic techniques. Using a single 48-channel streamer, the 3-D site with an area of 1200 m x 600 m was surveyed in 10 days. A variety of complex geologic structures (e.g. thrusts, folds, channel-fill) up to similar to150 m below the water bottom were obtained with a 15 in.(3) water gun. The 3-D data allowed the construction of an accurate velocity model and the distinction of five major seismic facies within the Lower Freshwater Molasse (Aquitanian) and the Quaternary sedimentary units. Additionally, the Plateau Molasse (PM) and Subalpine Molasse (SM) erosional surface, "La Paudeze" thrust fault (PM-SM boundary) and the thickness of Quaternary sediments were accurately delineated in 3-D.

Imaging EEG synchronization in schizophrenia patients with S-estimator

Recently a new measure of the cooperative behavior of simultaneous time series was introduced (Carmeli et al. NeuroImage 2005). This measure called S-estimator is defined from the embedding dimension in a state space. S-estimator quantifies the amount of synchronization within a data set by comparing the actual dimensionality of the set with the expected full dimensionality of the asynchronous set. It has the advantage of being a multivariate measure over traditionally used in systems neuroscience bivariate measures of synchronization. Multivariate measures of synchronization are of particular interest for applications in the field of modern multichannel EEG research, since they easily allow mapping of local and/or regional synchronization and are compatible with other imaging techniques. We applied Sestimator to the analysis of EEG synchronization in schizophrenia patients vs. matched controls. The whole-head mapping with S-estimator revealed a specific pattern of local synchronization in schizophrenia patients. The differences in the landscape of synchronization included decreased local synchronization in the territories over occipital and midline areas and increased synchronization over temporal areas. In frontal areas, the S-estimator revealed a tendency for an asymmetry: decreased S-values over the left hemisphere were adjacent to increased values over the right hemisphere. Separate calculations showed reproducibility of this pattern across the main EEG frequency bands. The maintenance of the same synchronization landscape across EEG frequencies probably implies the structural changes in the cortical circuitry of schizophrenia patients. These changes are regionally specific and suggest that schizophrenia is a misconnectivity rather than hypo- or hyper-connectivity disorder.

Orthogonal organic and aqueous-based washes of tissue sections to enhance protein sensitivity by MALDI imaging mass spectrometry.

Imaging mass spectrometry (IMS) is an emergent and innovative approach for measuring the composition, abundance and regioselectivity of molecules within an investigated area of fixed dimension. Although providing unprecedented molecular information compared with conventional MS techniques, enhancement of protein signature by IMS is still necessary and challenging. This paper demonstrates the combination of conventional organic washes with an optimized aqueous-based buffer for tissue section preparation before matrix-assisted laser desorption/ionization (MALDI) IMS of proteins. Based on a 500 mM ammonium formate in water-acetonitrile (9:1; v/v, 0.1% trifluororacetic acid, 0.1% Triton) solution, this buffer wash has shown to significantly enhance protein signature by profiling and IMS (~fourfold) when used after organic washes (70% EtOH followed by 90% EtOH), improving the quality and number of ion images obtained from mouse kidney and a 14-day mouse fetus whole-body tissue sections, while maintaining a similar reproducibility with conventional tissue rinsing. Even if some protein losses were observed, the data mining has demonstrated that it was primarily low abundant signals and that the number of new peaks found is greater with the described procedure. The proposed buffer has thus demonstrated to be of high efficiency for tissue section preparation providing novel and complementary information for direct on-tissue MALDI analysis compared with solely conventional organic rinsing.

A gamma camera count rate saturation correction method for whole-body planar imaging.

Whole-body (WB) planar imaging has long been one of the staple methods of dosimetry, and its quantification has been formalized by the MIRD Committee in pamphlet no 16. One of the issues not specifically addressed in the formalism occurs when the count rates reaching the detector are sufficiently high to result in camera count saturation. Camera dead-time effects have been extensively studied, but all of the developed correction methods assume static acquisitions. However, during WB planar (sweep) imaging, a variable amount of imaged activity exists in the detector's field of view as a function of time and therefore the camera saturation is time dependent. A new time-dependent algorithm was developed to correct for dead-time effects during WB planar acquisitions that accounts for relative motion between detector heads and imaged object. Static camera dead-time parameters were acquired by imaging decaying activity in a phantom and obtaining a saturation curve. Using these parameters, an iterative algorithm akin to Newton's method was developed, which takes into account the variable count rate seen by the detector as a function of time. The algorithm was tested on simulated data as well as on a whole-body scan of high activity Samarium-153 in an ellipsoid phantom. A complete set of parameters from unsaturated phantom data necessary for count rate to activity conversion was also obtained, including build-up and attenuation coefficients, in order to convert corrected count rate values to activity. The algorithm proved successful in accounting for motion- and time-dependent saturation effects in both the simulated and measured data and converged to any desired degree of precision. The clearance half-life calculated from the ellipsoid phantom data was calculated to be 45.1 h after dead-time correction and 51.4 h with no correction; the physical decay half-life of Samarium-153 is 46.3 h. Accurate WB planar dosimetry of high activities relies on successfully compensating for camera saturation which takes into account the variable activity in the field of view, i.e. time-dependent dead-time effects. The algorithm presented here accomplishes this task.

Confirmation of functional zones within the human subthalamic nucleus: patterns of connectivity and sub-parcellation using diffusion weighted imaging.

The subthalamic nucleus (STN) is a small, glutamatergic nucleus situated in the diencephalon. A critical component of normal motor function, it has become a key target for deep brain stimulation in the treatment of Parkinson's disease. Animal studies have demonstrated the existence of three functional sub-zones but these have never been shown conclusively in humans. In this work, a data driven method with diffusion weighted imaging demonstrated that three distinct clusters exist within the human STN based on brain connectivity profiles. The STN was successfully sub-parcellated into these regions, demonstrating good correspondence with that described in the animal literature. The local connectivity of each sub-region supported the hypothesis of bilateral limbic, associative and motor regions occupying the anterior, mid and posterior portions of the nucleus respectively. This study is the first to achieve in-vivo, non-invasive anatomical parcellation of the human STN into three anatomical zones within normal diagnostic scan times, which has important future implications for deep brain stimulation surgery.

Research development of terahertz imaging systems

The electromagnetic radiation at a terahertz frequencies (from 0.1 THz to 10 THz) is situated in the frequency band comprised between the optical band and the radio band. The interest of the scientific community in this frequency band has grown up due to its large capabilities to develop innovative imaging systems. The terahertz waves are able to generate extremely short pulses that achieve good spatial resolution, good penetration capabilities and allow to identify microscopic structures using spectral analysis. The work carried out during the period of the grant has been based on the developement of system working at the aforementioned frequency band. The main system is based on a total power radiometer working at 0.1 THz to perform security imaging. Moreover, the development of this system has been useful to gain knowledge in the behavior of the component systems at this frequency band. Moreover, a vectorial network analyzer has been used to characterize materials and perform active raster imaging. A materials measurement system has been designed and used to measure material properties as permittivity, losses and water concentration. Finally, the design of a terahertz time-domain spectrometer (THz-TDS) system has been started. This system will allow to perform tomographic measurement with very high penetration resolutions while allowing the spectral characterization of the sample material. The application range of this kind of system is very wide: from the identification of cancerous tissues of a skin to the characterization of the thickness of a painted surface of a car.