Nonrigid Motion Compensation of Free Breathing Acquired Myocardial Perfusion Data

In this work, we present a novel method to compensate the movement in images acquired during free breathing using first-pass gadolinium enhanced, myocardial perfusion magnetic resonance imaging (MRI). First, we use independent component analysis (ICA) to identify the optimal number of independent components (ICs) that separate the breathing motion from the intensity change induced by the contrast agent. Then, synthetic images are created by recombining the ICs, but other then in previously published work (Milles et al. 2008), we omit the component related to motion, and therefore, the resulting reference image series is free of motion. Motion compensation is then achieved by using a multi-pass non-rigid image registration scheme. We tested our method on 15 distinct image series (5 patients) consisting of 58 images each and we validated our method by comparing manually tracked intensity profiles of the myocardial sections to automatically generated ones before and after registration. The average correlation to the manually obtained curves before registration 0:89 0:11 was increased to 0:98 0:02

Automatic motion compensation of free breathing acquired myocardial perfusion data by using independent component analysis

Images acquired during free breathing using first-pass gadolinium-enhanced myocardial perfusion magnetic resonance imaging (MRI) exhibit a quasiperiodic motion pattern that needs to be compensated for if a further automatic analysis of the perfusion is to be executed. In this work, we present a method to compensate this movement by combining independent component analysis (ICA) and image registration: First, we use ICA and a time?frequency analysis to identify the motion and separate it from the intensity change induced by the contrast agent. Then, synthetic reference images are created by recombining all the independent components but the one related to the motion. Therefore, the resulting image series does not exhibit motion and its images have intensities similar to those of their original counterparts. Motion compensation is then achieved by using a multi-pass image registration procedure. We tested our method on 39 image series acquired from 13 patients, covering the basal, mid and apical areas of the left heart ventricle and consisting of 58 perfusion images each. We validated our method by comparing manually tracked intensity profiles of the myocardial sections to automatically generated ones before and after registration of 13 patient data sets (39 distinct slices). We compared linear, non-linear, and combined ICA based registration approaches and previously published motion compensation schemes. Considering run-time and accuracy, a two-step ICA based motion compensation scheme that first optimizes a translation and then for non-linear transformation performed best and achieves registration of the whole series in 32 ± 12 s on a recent workstation. The proposed scheme improves the Pearsons correlation coefficient between manually and automatically obtained time?intensity curves from .84 ± .19 before registration to .96 ± .06 after registration

Numerical study of ultrasound induced non-linear shape and size bubble oscillations in viscoelastic media

The theoretical study of forced bubble oscillations is motivated by the importance of cavitation bubbles and oscillating encapsulated microbubbles (i.e. contrast agents) in medical sciences. In more details,theoretical studies on bubble dynamics addressing the sound-bubble interaction phenomenon provide the basis for understanding the dynamics of contrast agent microbubbles used in medical diagnosis and of non-linearly oscillating cavitation bubbles in the case of high-intensity ultrasound therapy. Moreover, the inclusion of viscoelasticity is of vital importance for an accurate theoretical analysis since most biological tissues and fluids exhibit non-Newtonian behavior.

Nanopartículas superparamagnéticas encapsuladas com polímeros para tratamento de câncer por hipertermia

O câncer é uma das maiores causas de mortalidade no Brasil e no mundo, com potencial de crescimento nas próximas décadas. Um tipo de tratamento promissor é a hipertermia magnética, procedimento no qual as células tumorais morrem pelo efeito do calor gerado por partículas magnéticas após a aplicação de campo magnético alternado em frequências adequadas. Tais partículas também são capazes de atuar como agentes de contraste para imageamento por ressonância magnética, um poderoso método de diagnóstico para identificação de células neoplásicas, formando a combinação conhecida como theranostics (terapia e diagnóstico). Neste trabalho foram sintetizadas nanopartículas de óxido de ferro por método de coprecipitação com posterior encapsulação por técnica de nano spray drying, visando sua aplicação no tratamento de câncer por hipertermia e como agente de contraste para imageamento por ressonância magnética. Para a encapsulação foram utilizadas matrizes poliméricas de Maltodextrina com Polissorbato 80, Pluronic F68, Eudragit® S100 e PCL com Pluronic F68, escolhidos com o intuito de formar partículas que dispersem bem em meio aquoso e que consigam atingir alvo tumoral após administração no corpo do paciente. Parâmetros de secagem pelo equipamento Nano Spray Dryer, como temperatura, solvente e concentração de reagentes, foram avaliados. As partículas formadas foram caracterizadas por Microscopia Eletrônica de Varredura, Difração de Raios-X, Análise Termogravimétrica, Espalhamento de Luz Dinâmico, Espectroscopia de Infravermelho, magnetismo quanto a magnetização de saturação e temperatura, citotoxicidade e potencial de aquecimento. Tais procedimentos indicaram que o método de coprecipitação produziu nanopartículas de magnetita de tamanho em torno 20 nm, superparamagnéticas a temperatura ambiente, sem potencial citotóxico. A técnica de nano spray drying foi eficiente para a formação de partículas com tamanho em torno de 1 μm, também superparamagnéticas, biocompatíveis e com propriedades magnéticas adequadas e para aplicações pretendidas. Destaca-se a amostra com Pluronic, OF-10/15-1P, que apresentou magnetização de saturação de 68,7 emu/g e interação específica com células tumorais.

Synthesis of gadolinium-labeled shell-crosslinked nanoparticles for magnetic resonance imaging applications

Shell-crosslinked knedel-like nanoparticles (SCKs; knedel is a Polish term for dumplings) were derivatized with gadolinium Shell chelates and studied as robust magnetic-resonance-imaging-active structures with hydrodynamic diameters of 40 +/- 3 nm. SCKs possessing an amphiphilic core-shell morphology were produced from the aqueous assembly of diblock copolymers of poly(acrylic acid) (PAA) and poly(methyl acrylate) (PMA), PAA(52)-b-PMA(128), and subsequent covalent crosslinking by amidation upon reaction with 2,2'-(ethylenedioxy)bis(ethylamine) throughout the shell layer. The properties of these materials, including non-toxicity towards mammalian cells, non-immunogenicity within mice, and capability for polyvalent targeting, make them ideal candidates for utilization within biological systems. The synthesis of SCKs derivatized with Gd-III and designed for potential use as a unique nanometer-scale contrast agent for MRI applications is described herein. Utilization of an amino-functionalized diethylenetriaminepentaacetic acid-Gd analogue allowed for direct covalent conjugation throughout the hydrophilic shell layer of the SCKs and served to increase the rotational correlation lifetime of the Gd. In addition, the highly hydrated nature of the shell layer in which the Gd was located allowed for rapid water exchange; thus, the resulting material demonstrated large ionic relaxivities (39 s(-1) mM(-1)) in an applied magnetic field of 0.47 T at 40 degrees C and, as a result of the large loading capacity of the material, also demonstrated high molecular relaxivities (20 000 s(-1) mM(-1)).

MRI pressure and stress measurement in novel homogeneous soft solids

This paper presents MRI measurements of a novel semi solid MR contrast agent to pressure. The agent is comprised of potassium chloride cross linked carageenan gum at a concentration of 2% w/v, with micron size lipid coated bubbles of air at a concentration of 3% v/v. The choice for an optimum suspending medium, the methods of production and the preliminary MRI results are presented herein. The carageenan gum is shown to be ideally elastic for compressions relating to volume changes less than 15%, in contrast to the inelastic gellan gum also tested. Although slightly lower than that of gellan gum, carageenan has a water diffusion coefficient of 1.72×10-9 m2.s-1 indicating its suitability to this purpose. RARE imaging is performed whilst simultaneously compressing test and control samples and a maximum sensitivity of 1.6% MR signal change per % volume change is found which is shown to be independent of proton density variations due to the presence of microbubbles and compression. This contrast agent could prove useful for numerous applications, and particularly in chemical engineering. More generally the method allows the user to non-invasively image with MRI any process that causes, within the solid, local changes either in bubble size or bubble shape. © 2008 American Institute of Physics.

Radiotherapy Treatment Assessment using DCE-MRI

Abstract

The goal of modern radiotherapy is to precisely deliver a prescribed radiation dose to delineated target volumes that contain a significant amount of tumor cells while sparing the surrounding healthy tissues/organs. Precise delineation of treatment and avoidance volumes is the key for the precision radiation therapy. In recent years, considerable clinical and research efforts have been devoted to integrate MRI into radiotherapy workflow motivated by the superior soft tissue contrast and functional imaging possibility. Dynamic contrast-enhanced MRI (DCE-MRI) is a noninvasive technique that measures properties of tissue microvasculature. Its sensitivity to radiation-induced vascular pharmacokinetic (PK) changes has been preliminary demonstrated. In spite of its great potential, two major challenges have limited DCE-MRI’s clinical application in radiotherapy assessment: the technical limitations of accurate DCE-MRI imaging implementation and the need of novel DCE-MRI data analysis methods for richer functional heterogeneity information.

This study aims at improving current DCE-MRI techniques and developing new DCE-MRI analysis methods for particular radiotherapy assessment. Thus, the study is naturally divided into two parts. The first part focuses on DCE-MRI temporal resolution as one of the key DCE-MRI technical factors, and some improvements regarding DCE-MRI temporal resolution are proposed; the second part explores the potential value of image heterogeneity analysis and multiple PK model combination for therapeutic response assessment, and several novel DCE-MRI data analysis methods are developed.

I. Improvement of DCE-MRI temporal resolution. First, the feasibility of improving DCE-MRI temporal resolution via image undersampling was studied. Specifically, a novel MR image iterative reconstruction algorithm was studied for DCE-MRI reconstruction. This algorithm was built on the recently developed compress sensing (CS) theory. By utilizing a limited k-space acquisition with shorter imaging time, images can be reconstructed in an iterative fashion under the regularization of a newly proposed total generalized variation (TGV) penalty term. In the retrospective study of brain radiosurgery patient DCE-MRI scans under IRB-approval, the clinically obtained image data was selected as reference data, and the simulated accelerated k-space acquisition was generated via undersampling the reference image full k-space with designed sampling grids. Two undersampling strategies were proposed: 1) a radial multi-ray grid with a special angular distribution was adopted to sample each slice of the full k-space; 2) a Cartesian random sampling grid series with spatiotemporal constraints from adjacent frames was adopted to sample the dynamic k-space series at a slice location. Two sets of PK parameters’ maps were generated from the undersampled data and from the fully-sampled data, respectively. Multiple quantitative measurements and statistical studies were performed to evaluate the accuracy of PK maps generated from the undersampled data in reference to the PK maps generated from the fully-sampled data. Results showed that at a simulated acceleration factor of four, PK maps could be faithfully calculated from the DCE images that were reconstructed using undersampled data, and no statistically significant differences were found between the regional PK mean values from undersampled and fully-sampled data sets. DCE-MRI acceleration using the investigated image reconstruction method has been suggested as feasible and promising.

Second, for high temporal resolution DCE-MRI, a new PK model fitting method was developed to solve PK parameters for better calculation accuracy and efficiency. This method is based on a derivative-based deformation of the commonly used Tofts PK model, which is presented as an integrative expression. This method also includes an advanced Kolmogorov-Zurbenko (KZ) filter to remove the potential noise effect in data and solve the PK parameter as a linear problem in matrix format. In the computer simulation study, PK parameters representing typical intracranial values were selected as references to simulated DCE-MRI data for different temporal resolution and different data noise level. Results showed that at both high temporal resolutions (<1s) and clinically feasible temporal resolution (~5s), this new method was able to calculate PK parameters more accurate than the current calculation methods at clinically relevant noise levels; at high temporal resolutions, the calculation efficiency of this new method was superior to current methods in an order of 102. In a retrospective of clinical brain DCE-MRI scans, the PK maps derived from the proposed method were comparable with the results from current methods. Based on these results, it can be concluded that this new method can be used for accurate and efficient PK model fitting for high temporal resolution DCE-MRI.

II. Development of DCE-MRI analysis methods for therapeutic response assessment. This part aims at methodology developments in two approaches. The first one is to develop model-free analysis method for DCE-MRI functional heterogeneity evaluation. This approach is inspired by the rationale that radiotherapy-induced functional change could be heterogeneous across the treatment area. The first effort was spent on a translational investigation of classic fractal dimension theory for DCE-MRI therapeutic response assessment. In a small-animal anti-angiogenesis drug therapy experiment, the randomly assigned treatment/control groups received multiple fraction treatments with one pre-treatment and multiple post-treatment high spatiotemporal DCE-MRI scans. In the post-treatment scan two weeks after the start, the investigated Rényi dimensions of the classic PK rate constant map demonstrated significant differences between the treatment and the control groups; when Rényi dimensions were adopted for treatment/control group classification, the achieved accuracy was higher than the accuracy from using conventional PK parameter statistics. Following this pilot work, two novel texture analysis methods were proposed. First, a new technique called Gray Level Local Power Matrix (GLLPM) was developed. It intends to solve the lack of temporal information and poor calculation efficiency of the commonly used Gray Level Co-Occurrence Matrix (GLCOM) techniques. In the same small animal experiment, the dynamic curves of Haralick texture features derived from the GLLPM had an overall better performance than the corresponding curves derived from current GLCOM techniques in treatment/control separation and classification. The second developed method is dynamic Fractal Signature Dissimilarity (FSD) analysis. Inspired by the classic fractal dimension theory, this method measures the dynamics of tumor heterogeneity during the contrast agent uptake in a quantitative fashion on DCE images. In the small animal experiment mentioned before, the selected parameters from dynamic FSD analysis showed significant differences between treatment/control groups as early as after 1 treatment fraction; in contrast, metrics from conventional PK analysis showed significant differences only after 3 treatment fractions. When using dynamic FSD parameters, the treatment/control group classification after 1st treatment fraction was improved than using conventional PK statistics. These results suggest the promising application of this novel method for capturing early therapeutic response.

The second approach of developing novel DCE-MRI methods is to combine PK information from multiple PK models. Currently, the classic Tofts model or its alternative version has been widely adopted for DCE-MRI analysis as a gold-standard approach for therapeutic response assessment. Previously, a shutter-speed (SS) model was proposed to incorporate transcytolemmal water exchange effect into contrast agent concentration quantification. In spite of richer biological assumption, its application in therapeutic response assessment is limited. It might be intriguing to combine the information from the SS model and from the classic Tofts model to explore potential new biological information for treatment assessment. The feasibility of this idea was investigated in the same small animal experiment. The SS model was compared against the Tofts model for therapeutic response assessment using PK parameter regional mean value comparison. Based on the modeled transcytolemmal water exchange rate, a biological subvolume was proposed and was automatically identified using histogram analysis. Within the biological subvolume, the PK rate constant derived from the SS model were proved to be superior to the one from Tofts model in treatment/control separation and classification. Furthermore, novel biomarkers were designed to integrate PK rate constants from these two models. When being evaluated in the biological subvolume, this biomarker was able to reflect significant treatment/control difference in both post-treatment evaluation. These results confirm the potential value of SS model as well as its combination with Tofts model for therapeutic response assessment.

In summary, this study addressed two problems of DCE-MRI application in radiotherapy assessment. In the first part, a method of accelerating DCE-MRI acquisition for better temporal resolution was investigated, and a novel PK model fitting algorithm was proposed for high temporal resolution DCE-MRI. In the second part, two model-free texture analysis methods and a multiple-model analysis method were developed for DCE-MRI therapeutic response assessment. The presented works could benefit the future DCE-MRI routine clinical application in radiotherapy assessment.

Radiotherapy Assessment Using Diffusion Weighted MRI

Purpose: There are two goals of this study. The first goal of this study is to investigate the feasibility of using classic textural feature extraction in radiotherapy response assessment among a unique cohort of early stage breast cancer patients who received the single-dose preoperative radiotherapy. The second goal of this study is to investigate the clinical feasibility of using classic texture features as potential biomarkers which are supplementary to regional apparent diffusion coefficient in gynecological cancer radiotherapy response assessment.

Methods and Materials: For the breast cancer study, 15 patients with early stage breast cancer were enrolled in this retrospective study. Each patient received a single-fraction radiation treatment, and DWI and DCE-MRI scans were conducted before and after the radiotherapy. DWI scans were acquired using a spin-echo EPI sequence with diffusion weighting factors of b = 0 and b = 500 mm2/s, and the apparent diffusion coefficient (ADC) maps were calculated. DCE-MRI scans were acquired using a T1-weighted 3D SPGR sequence with a temporal resolution of about 1 minute. The contrast agent (CA) was intravenously injected with a 0.1 mmol/kg bodyweight dose at 2 ml/s. Two parameters, volume transfer constant (Ktrans) and kep were analyzed using the two-compartment Tofts pharmacokinetic model. For pharmacokinetic parametric maps and ADC maps, 33 textural features were generated from the clinical target volume (CTV) in a 3D fashion using the classic gray level co-occurrence matrix (GLCOM) and gray level run length matrix (GLRLM). Wilcoxon signed-rank test was used to determine the significance of each texture feature’s change after the radiotherapy. The significance was set to 0.05 with Bonferroni correction.

For the gynecological cancer study, 12 female patients with gynecologic cancer treated with fractionated external beam radiotherapy (EBRT) combined with high dose rate (HDR) intracavitary brachytherapy were studied. Each patient first received EBRT treatment followed by five fractions of HDR treatment. Before EBRT and before each fraction of brachytherapy, Diffusion Weighted MRI (DWI-MRI) and CT scans were acquired. DWI scans were acquired in sagittal plane utilizing a spin-echo echo-planar imaging sequence with weighting factors of b = 500 s/mm2 and b = 1000 s/mm2, one set of images of b = 0 s/mm2 were also acquired. ADC maps were calculated using linear least-square fitting method. Distributed diffusion coefficient (DDC) maps and stretching parameter α were calculated. For ADC and DDC maps, 33 classic texture features were generated utilizing the classic gray level run length matrix (GLRLM) and gray level co-occurrence matrix (GLCOM) from high-risk clinical target volume (HR-CTV). Wilcoxon signed-rank statistics test was applied to determine the significance of each feature’s numerical value change after radiotherapy. Significance level was set to 0.05 with multi-comparison correction if applicable.

Results: For the breast cancer study, regarding ADC maps calculated from DWI-MRI, 24 out of 33 CTV features changed significantly after the radiotherapy. For DCE-MRI pharmacokinetic parameters, all 33 CTV features of Ktrans and 33 features of kep changed significantly.

For the gynecological cancer study, regarding ADC maps, 28 out of 33 HR-CTV texture features showed significant changes after the EBRT treatment. 28 out of 33 HR-CTV texture features indicated significant changes after HDR treatments. The texture features that indicated significant changes after HDR treatments are the same as those after EBRT treatment. 28 out of 33 HR-CTV texture features showed significant changes after whole radiotherapy treatment process. The texture features that indicated significant changes for the whole treatment process are the same as those after HDR treatments.

Conclusion: Initial results indicate that certain classic texture features are sensitive to radiation-induced changes. Classic texture features with significant numerical changes can be used in monitoring radiotherapy effect. This might suggest that certain texture features might be used as biomarkers which are supplementary to ADC and DDC for assessment of radiotherapy response in breast cancer and gynecological cancer.

Theranostic Fibers for Simultaneous Imaging and Drug Delivery

New methods for creating theranostic systems with simultaneous encapsulation of therapeutic, diagnostic, and targeting agents are much sought after. This work reports for the first time the use of coaxial electrospinning to prepare such systems in the form of core–shell fibers. Eudragit S100 was used to form the shell of the fibers, while the core comprised poly(ethylene oxide) loaded with the magnetic resonance contrast agent Gd(DTPA) (Gd(III) diethylenetriaminepentaacetate hydrate) and indomethacin as a model therapeutic agent. The fibers had linear cylindrical morphologies with clear core–shell structures, as demonstrated by electron microscopy. X-ray diffraction and differential scanning calorimetry proved that both indomethacin and Gd(DTPA) were present in the fibers in the amorphous physical form. This is thought to be a result of intermolecular interactions between the different components, the presence of which was suggested by infrared spectroscopy. In vitro dissolution tests indicated that the fibers could provide targeted release of the active ingredients through a combined mechanism of erosion and diffusion. The proton relaxivities for Gd(DTPA) released from the fibers into tris buffer increased (r1 = 4.79–9.75 s–1 mM–1; r2 = 7.98–14.22 s–1 mM–1) compared with fresh Gd(DTPA) (r1 = 4.13 s–1 mM–1 and r2 = 4.40 s–1 mM–1), which proved that electrospinning has not diminished the contrast properties of the complex. The new systems reported herein thus offer a new platform for delivering therapeutic and imaging agents simultaneously to the colon.

Suivi par élastographie ultrasonore après réparation endovasculaire d’anévrisme aorto-iliaque : étude de faisabilité in vivo

Les maladies cardiovasculaires sont la première cause de mortalité dans le monde et les anévrismes de l’aorte abdominale (AAAs) font partie de ce lot déplorable. Un anévrisme est la dilatation d’une artère pouvant conduire à la mort. Une rupture d’AAA s’avère fatale près de 80% du temps. Un moyen de traiter les AAAs est l’insertion d’une endoprothèse (SG) dans l’aorte, communément appelée la réparation endovasculaire (EVAR), afin de réduire la pression exercée par le flux sanguin sur la paroi. L’efficacité de ce traitement est compromise par la survenue d’endofuites (flux sanguins entre la prothèse et le sac anévrismal) pouvant conduire à la rupture de l’anévrisme. Ces flux sanguins peuvent survenir à n’importe quel moment après le traitement EVAR. Une surveillance par tomodensitométrie (CT-scan) annuelle est donc requise, augmentant ainsi le coût du suivi post-EVAR et exposant le patient à la radiation ionisante et aux complications des contrastes iodés. L’endotension est le concept de dilatation de l’anévrisme sans la présence d’une endofuite apparente au CT-scan. Après le traitement EVAR, le sang dans le sac anévrismal coagule pour former un thrombus frais, qui deviendra progressivement un thrombus plus fibreux et plus organisé, donnant lieu à un rétrécissement de l’anévrisme. Il y a très peu de données dans la littérature pour étudier ce processus temporel et la relation entre le thrombus frais et l’endotension. L’étalon d’or du suivi post-EVAR, le CT-scan, ne peut pas détecter la présence de thrombus frais. Il y a donc un besoin d’investir dans une technique sécuritaire et moins coûteuse pour le suivi d’AAAs après EVAR. Une méthode récente, l’élastographie dynamique, mesure l’élasticité des tissus en temps réel. Le principe de cette technique repose sur la génération d’ondes de cisaillement et l’étude de leur propagation afin de remonter aux propriétés mécaniques du milieu étudié. Cette thèse vise l’application de l’élastographie dynamique pour la détection des endofuites ainsi que de la caractérisation mécanique des tissus du sac anévrismal après le traitement EVAR. Ce projet dévoile le potentiel de l’élastographie afin de réduire les dangers de la radiation, de l’utilisation d’agent de contraste ainsi que des coûts du post-EVAR des AAAs. L’élastographie dynamique utilisant le « Shear Wave Imaging » (SWI) est prometteuse. Cette modalité pourrait complémenter l’échographie-Doppler (DUS) déjà utilisée pour le suivi d’examen post-EVAR. Le SWI a le potentiel de fournir des informations sur l’organisation fibreuse du thrombus ainsi que sur la détection d’endofuites. Tout d’abord, le premier objectif de cette thèse consistait à tester le SWI sur des AAAs dans des modèles canins pour la détection d’endofuites et la caractérisation du thrombus. Des SGs furent implantées dans un groupe de 18 chiens avec un anévrisme créé au moyen de la veine jugulaire. 4 anévrismes avaient une endofuite de type I, 13 avaient une endofuite de type II et un anévrisme n’avait pas d’endofuite. Des examens échographiques, DUS et SWI ont été réalisés à l’implantation, puis 1 semaine, 1 mois, 3 mois et 6 mois après le traitement EVAR. Une angiographie, un CT-scan et des coupes macroscopiques ont été produits au sacrifice. Les régions d’endofuites, de thrombus frais et de thrombus organisé furent identifiées et segmentées. Les valeurs de rigidité données par le SWI des différentes régions furent comparées. Celles-ci furent différentes de façon significative (P < 0.001). Également, le SWI a pu détecter la présence d’endofuites où le CT-scan (1) et le DUS (3) ont échoué. Dans la continuité de ces travaux, le deuxième objectif de ce projet fut de caractériser l’évolution du thrombus dans le temps, de même que l’évolution des endofuites après embolisation dans des modèles canins. Dix-huit anévrismes furent créés dans les artères iliaques de neuf modèles canins, suivis d’une endofuite de type I après EVAR. Deux gels embolisants (Chitosan (Chi) ou Chitosan-Sodium-Tetradecyl-Sulfate (Chi-STS)) furent injectés dans le sac anévrismal pour promouvoir la guérison. Des examens échographiques, DUS et SWI ont été effectués à l’implantation et après 1 semaine, 1 mois, 3 mois et 6 mois. Une angiographie, un CT-scan et un examen histologique ont été réalisés au sacrifice afin d’évaluer la présence, le type et la grosseur de l’endofuite. Les valeurs du module d’élasticité des régions d’intérêts ont été identifiées et segmentées sur les données pathologiques. Les régions d’endofuites et de thrombus frais furent différentes de façon significative comparativement aux autres régions (P < 0.001). Les valeurs d’élasticité du thrombus frais à 1 semaine et à 3 mois indiquent que le SWI peut évaluer la maturation du thrombus, de même que caractériser l’évolution et la dégradation des gels embolisants dans le temps. Le SWI a pu détecter des endofuites où le DUS a échoué (2) et, contrairement au CT-scan, détecter la présence de thrombus frais. Finalement, la dernière étape du projet doctoral consistait à appliquer le SWI dans une phase clinique, avec des patients humains ayant déjà un AAA, pour la détection d’endofuite et la caractérisation de l’élasticité des tissus. 25 patients furent sélectionnés pour participer à l’étude. Une comparaison d’imagerie a été produite entre le SWI, le CT-scan et le DUS. Les valeurs de rigidité données par le SWI des différentes régions (endofuite, thrombus) furent identifiées et segmentées. Celles-ci étaient distinctes de façon significative (P < 0.001). Le SWI a détecté 5 endofuites sur 6 (sensibilité de 83.3%) et a eu 6 faux positifs (spécificité de 76%). Le SWI a pu détecter la présence d’endofuites où le CT-scan (2) ainsi que le DUS (2) ont échoué. Il n’y avait pas de différence statistique notable entre la rigidité du thrombus pour un AAA avec endofuite et un AAA sans endofuite. Aucune corrélation n’a pu être établie de façon significative entre les diamètres des AAAs ainsi que leurs variations et l’élasticité du thrombus. Le SWI a le potentiel de détecter les endofuites et caractériser le thrombus selon leurs propriétés mécaniques. Cette technique pourrait être combinée au suivi des AAAs post-EVAR, complémentant ainsi l’imagerie DUS et réduisant le coût et l’exposition à la radiation ionisante et aux agents de contrastes néphrotoxiques.

Applications of functional nucleic acids in imaging, sensing and drug delivery and investigations of DNAzyme-metal interactions

Functional nucleic acids (FNA), including nucleic acids catalysts (ribozymes and DNAzymes) and ligands (aptamers), have been discovered in nature or isolated in a laboratory through a process called in vitro selection. They are nucleic acids with functions similar to protein enzymes or antibodies. They have been developed into sensors with high sensitivity and selectivity; it is realized by converting the reaction catalyzed by a DNAzyme/ribozyme or the binding event of an aptamer to a fluorescent, colorimetric or electrochemical signal. While a number of studies have been reported for in vitro sensing using DNAzymes or aptamers, there are few reports on in vivo sensing or imaging. MRI is a non-invasive imaging technique; smart MRI contrast agents were synthesized for molecular imaging purposes. However, their rational design remains a challenge due to the difficulty to predict molecular interactions. Chapter 2 focuses on rational design of smart T1-weighted MRI contrast agents with high specificity based on DNAzymes and aptamers. It was realized by changing the molecular weight of the gadolinium conjugated DNA strand with the analytes, which lead to analyte-specific water proton relaxation responses and contrast changes on an MRI image. The designs are general; the high selectivity of FNA was retained. Most FNA-based fluorescent sensors require covalent labeling of fluorophore/quencher to FNAs, which incurrs extra expenses and could interfere the function of FNAs. Chapter 3 describes a new sensor design avoiding the covalent labeling of fluorophore and quencher. The fluorescence of malachite green (MG) was regulated by the presence of adenosine. Conjugate of aptamers of MG and adenosine and a bridge strand were annealed in a solution containing MG. The MG aptamer did not bind MG because of its hybridization to the bridge strand, resulting in low fluorescence signal of MG. The hybridization was weakened in the presence of adenosine, leading to the binding of MG to its aptamer and a fluorescence increase. The sensor has comparable detection limit (20 micromolar) and specificity to its labeled derivatives. Enzymatic activity of most DNAzymes requires metal cations. The research on the metal-DNAzyme interaction is of interest and challenge to scientists because of the lack of structural information. Chapters 4 presents the research on the characterization of the interaction between a Cu2+-dependent DNAzyme and Cu2+. Electron paramagnetic resonance (EPR) and UV-Vis spectroscopy were used to probe the binding of Cu2+ to the DNAzyme; circular dichroism was used to probe the conformational change of the DNAzyme induced by Cu2+. It was proposed that the conformational change by the Cu2+ binding is important for the activity of the DNAzyme. Chapter 5 reports the dependence of the activity of 8-17 DNAzyme on the presence of both Pb2+ and other metal cations including Zn2+, Cd2+ and Mg2+. It was discovered that presence of those metal cations can be cooperative or inhibitive to 8-17 activity. It is hypothesized that the 8-17 DNAzyme had multiple binding sites for metal cations based on the results. Cisplatin is effective killing tumor cells, but with significant side effects, which can be minimized by its targeted delivery. Chapter 6 focuses on the effort to functionalize liposomes encapsulating cisplatin by an aptamer that selectively bind nucleolin, an overexpressed protein by breast cancer cells. The study proved the selective cytotoxicity to breast cancer cells of the aptamer-functionalized liposome.

Oxidative stress and its haemodynamic consequences in chronic kidney disease

Chronic kidney disease (CKD) is associated with increased cardiovascular risk in comparison with the general population. This can be observed even in the early stages of CKD, and rises in proportion to the degree of renal impairment. Not only is cardiovascular disease (CVD) more prevalent in CKD, but its nature differs too, with an excess of morbidity and mortality associated with congestive cardiac failure, arrhythmia and sudden death, as well as the accelerated atherosclerosis which is also observed. Conventional cardiovascular risk factors such as hypertension, dyslipidaemia, obesity, glycaemia and smoking, are highly prevalent amongst patients with CKD, although in many of these examples the interaction between risk factor and disease differs from that which exists in normal renal function. Nevertheless, the extent of CVD cannot be fully explained by these conventional risk factors, and non-conventional factors specific to CKD are now recognised to contribute to the burden of CVD. Oxidative stress is a state characterised by excessive production of reactive oxygen species (ROS) and other radical species, a reduction in the capacity of antioxidant systems, and disturbance in normal redox homeostasis with depletion of protective vascular signalling molecules such as nitric oxide (NO). This results in oxidative damage to macromolecules such as lipids, proteins and DNA which can alter their functionality. Moreover, many enzymes are sensitive to redox regulation such that oxidative modification to cysteine thiol groups results in activation of signalling cascades which result in adverse cardiovascular effects such as vascular and endothelial dysfunction. Endothelial dysfunction and oxidative stress are present in association with many conventional cardiovascular risk factors, and can be observed even prior to the development of overt, clinical, vascular pathology, suggesting that these phenomena represent the earliest stages of CVD. In the presence of CKD, there is increased ROS production due to upregulated NADPH oxidase (NOX), increase in a circulating asymmetric dimethylarginine (ADMA), uncoupling of endothelial nitric oxide synthase (eNOS) as well as other mechanisms. There is also depletion in exogenous antioxidants such as ascorbic acid and tocopherol, and a reduction in activity of endogenous antioxidant systems regulated by the master gene regulator Nrf-2. In previous studies, circulating markers of oxidative stress have been shown to be increased in CKD, together with a reduction in endothelial function in a stepwise fashion relating to the severity of renal impairment. Not only is CVD linked to oxidative stress, but the progression of CKD itself is also in part dependent on redox sensitive mechanisms. For example, administration of the ROS scavenger tempol attenuates renal injury and reduces renal fibrosis seen on biopsy in a mouse model of CKD, whilst conversely, supplementation with the NOS inhibitor L-NAME causes proteinuria and renal impairment. Previous human studies examining the effect of antioxidant administration on vascular and renal function have been conflicting however. The work contained in this thesis therefore examines the effect of antioxidant administration on vascular and endothelial function in CKD. Firstly, 30 patients with CKD stages 3 – 5, and 20 matched hypertensive controls were recruited. Participants with CKD had lower ascorbic acid, higher TAP and ADMA, together with higher augmentation index and pulse wave velocity. There was no difference in baseline flow mediated dilatation (FMD) between groups. Intravenous ascorbic acid increased TAP and O2-, and reduced central BP and augmentation index in both groups, and lowered ADMA in the CKD group only. No effect on FMD was observed. The effects of ascorbic acid on kidney function was then investigated, however this was hindered by the inherent drawbacks of existing methods of non-invasively measuring kidney function. Arterial spin labelling MRI is an emerging imaging technique which allows measurement of renal perfusion without administration of an exogenous contrast agent. The technique relies upon application of an inversion pulse to blood within the vasculature proximal to the kidneys, which magnetically labels protons allowing measurement upon transit to the kidney. At the outset of this project local experience using ASL MRI was limited and there ensued a prolonged pre-clinical phase of testing with the aim of optimising imaging strategy. A study was then designed to investigate the repeatability of ASL MRI in a group of 12 healthy volunteers with normal renal function. The measured T1 longitudinal relaxation times and ASL MRI perfusion values were in keeping with those found in the literature; T1 time was 1376 ms in the cortex and 1491 ms in the whole kidney ROI, whilst perfusion was 321 mL/min/100g in the cortex, and 228 mL/min/100g in the whole kidney ROI. There was good reproducibility demonstrated on Bland Altman analysis, with a CVws was 9.2% for cortical perfusion and 7.1% for whole kidney perfusion. Subsequently, in a study of 17 patients with CKD and 24 healthy volunteers, the effects of ascorbic acid on renal perfusion was investigated. Although no change in renal perfusion was found following ascorbic acid, it was found that ASL MRI demonstrated significant differences between those with normal renal function and participants with CKD stages 3 – 5, with increased cortical and whole kidney T1, and reduced cortical and whole kidney perfusion. Interestingly, absolute perfusion showed a weak but significant correlation with progression of kidney disease over the preceding year. Ascorbic acid was therefore shown to have a significant effect on vascular biology both in CKD and in those with normal renal function, and to reduce ADMA only in patients with CKD. ASL MRI has shown promise as a non-invasive investigation of renal function and as a biomarker to identify individuals at high risk of progressive renal impairment.

Imaging of thoracic aortic injury

International audience

Suivi par élastographie ultrasonore après réparation endovasculaire d’anévrisme aorto-iliaque : étude de faisabilité in vivo

Les maladies cardiovasculaires sont la première cause de mortalité dans le monde et les anévrismes de l’aorte abdominale (AAAs) font partie de ce lot déplorable. Un anévrisme est la dilatation d’une artère pouvant conduire à la mort. Une rupture d’AAA s’avère fatale près de 80% du temps. Un moyen de traiter les AAAs est l’insertion d’une endoprothèse (SG) dans l’aorte, communément appelée la réparation endovasculaire (EVAR), afin de réduire la pression exercée par le flux sanguin sur la paroi. L’efficacité de ce traitement est compromise par la survenue d’endofuites (flux sanguins entre la prothèse et le sac anévrismal) pouvant conduire à la rupture de l’anévrisme. Ces flux sanguins peuvent survenir à n’importe quel moment après le traitement EVAR. Une surveillance par tomodensitométrie (CT-scan) annuelle est donc requise, augmentant ainsi le coût du suivi post-EVAR et exposant le patient à la radiation ionisante et aux complications des contrastes iodés. L’endotension est le concept de dilatation de l’anévrisme sans la présence d’une endofuite apparente au CT-scan. Après le traitement EVAR, le sang dans le sac anévrismal coagule pour former un thrombus frais, qui deviendra progressivement un thrombus plus fibreux et plus organisé, donnant lieu à un rétrécissement de l’anévrisme. Il y a très peu de données dans la littérature pour étudier ce processus temporel et la relation entre le thrombus frais et l’endotension. L’étalon d’or du suivi post-EVAR, le CT-scan, ne peut pas détecter la présence de thrombus frais. Il y a donc un besoin d’investir dans une technique sécuritaire et moins coûteuse pour le suivi d’AAAs après EVAR. Une méthode récente, l’élastographie dynamique, mesure l’élasticité des tissus en temps réel. Le principe de cette technique repose sur la génération d’ondes de cisaillement et l’étude de leur propagation afin de remonter aux propriétés mécaniques du milieu étudié. Cette thèse vise l’application de l’élastographie dynamique pour la détection des endofuites ainsi que de la caractérisation mécanique des tissus du sac anévrismal après le traitement EVAR. Ce projet dévoile le potentiel de l’élastographie afin de réduire les dangers de la radiation, de l’utilisation d’agent de contraste ainsi que des coûts du post-EVAR des AAAs. L’élastographie dynamique utilisant le « Shear Wave Imaging » (SWI) est prometteuse. Cette modalité pourrait complémenter l’échographie-Doppler (DUS) déjà utilisée pour le suivi d’examen post-EVAR. Le SWI a le potentiel de fournir des informations sur l’organisation fibreuse du thrombus ainsi que sur la détection d’endofuites. Tout d’abord, le premier objectif de cette thèse consistait à tester le SWI sur des AAAs dans des modèles canins pour la détection d’endofuites et la caractérisation du thrombus. Des SGs furent implantées dans un groupe de 18 chiens avec un anévrisme créé au moyen de la veine jugulaire. 4 anévrismes avaient une endofuite de type I, 13 avaient une endofuite de type II et un anévrisme n’avait pas d’endofuite. Des examens échographiques, DUS et SWI ont été réalisés à l’implantation, puis 1 semaine, 1 mois, 3 mois et 6 mois après le traitement EVAR. Une angiographie, un CT-scan et des coupes macroscopiques ont été produits au sacrifice. Les régions d’endofuites, de thrombus frais et de thrombus organisé furent identifiées et segmentées. Les valeurs de rigidité données par le SWI des différentes régions furent comparées. Celles-ci furent différentes de façon significative (P < 0.001). Également, le SWI a pu détecter la présence d’endofuites où le CT-scan (1) et le DUS (3) ont échoué. Dans la continuité de ces travaux, le deuxième objectif de ce projet fut de caractériser l’évolution du thrombus dans le temps, de même que l’évolution des endofuites après embolisation dans des modèles canins. Dix-huit anévrismes furent créés dans les artères iliaques de neuf modèles canins, suivis d’une endofuite de type I après EVAR. Deux gels embolisants (Chitosan (Chi) ou Chitosan-Sodium-Tetradecyl-Sulfate (Chi-STS)) furent injectés dans le sac anévrismal pour promouvoir la guérison. Des examens échographiques, DUS et SWI ont été effectués à l’implantation et après 1 semaine, 1 mois, 3 mois et 6 mois. Une angiographie, un CT-scan et un examen histologique ont été réalisés au sacrifice afin d’évaluer la présence, le type et la grosseur de l’endofuite. Les valeurs du module d’élasticité des régions d’intérêts ont été identifiées et segmentées sur les données pathologiques. Les régions d’endofuites et de thrombus frais furent différentes de façon significative comparativement aux autres régions (P < 0.001). Les valeurs d’élasticité du thrombus frais à 1 semaine et à 3 mois indiquent que le SWI peut évaluer la maturation du thrombus, de même que caractériser l’évolution et la dégradation des gels embolisants dans le temps. Le SWI a pu détecter des endofuites où le DUS a échoué (2) et, contrairement au CT-scan, détecter la présence de thrombus frais. Finalement, la dernière étape du projet doctoral consistait à appliquer le SWI dans une phase clinique, avec des patients humains ayant déjà un AAA, pour la détection d’endofuite et la caractérisation de l’élasticité des tissus. 25 patients furent sélectionnés pour participer à l’étude. Une comparaison d’imagerie a été produite entre le SWI, le CT-scan et le DUS. Les valeurs de rigidité données par le SWI des différentes régions (endofuite, thrombus) furent identifiées et segmentées. Celles-ci étaient distinctes de façon significative (P < 0.001). Le SWI a détecté 5 endofuites sur 6 (sensibilité de 83.3%) et a eu 6 faux positifs (spécificité de 76%). Le SWI a pu détecter la présence d’endofuites où le CT-scan (2) ainsi que le DUS (2) ont échoué. Il n’y avait pas de différence statistique notable entre la rigidité du thrombus pour un AAA avec endofuite et un AAA sans endofuite. Aucune corrélation n’a pu être établie de façon significative entre les diamètres des AAAs ainsi que leurs variations et l’élasticité du thrombus. Le SWI a le potentiel de détecter les endofuites et caractériser le thrombus selon leurs propriétés mécaniques. Cette technique pourrait être combinée au suivi des AAAs post-EVAR, complémentant ainsi l’imagerie DUS et réduisant le coût et l’exposition à la radiation ionisante et aux agents de contrastes néphrotoxiques.

Modélisation pharmacocinétique en imagerie par résonance magnétique et en tomographie d’émission par positrons appliquée à un modèle de glioblastome chez le rat

Résumé : En imagerie médicale, il est courant d’associer plusieurs modalités afin de tirer profit des renseignements complémentaires qu’elles fournissent. Par exemple, la tomographie d’émission par positrons (TEP) peut être combinée à l’imagerie par résonance magnétique (IRM) pour obtenir à la fois des renseignements sur les processus biologiques et sur l’anatomie du sujet. Le but de ce projet est d’explorer les synergies entre l’IRM et la TEP dans le cadre d’analyses pharmacocinétiques. Plus spécifiquement, d’exploiter la haute résolution spatiale et les renseignements sur la perfusion et la perméabilité vasculaire fournis par l’IRM dynamique avec agent de contraste afin de mieux évaluer ces mêmes paramètres pour un radiotraceur TEP injecté peu de temps après. L’évaluation précise des paramètres de perfusion du radiotraceur devrait permettre de mieux quantifier le métabolisme et de distinguer l’accumulation spécifique et non spécifique. Les travaux ont porté sur deux radiotraceurs de TEP (18F-fluorodésoxyglucose [FDG] et 18F-fluoroéthyle-tyrosine [FET]) ainsi que sur un agent de contraste d’IRM (acide gadopentétique [Gd DTPA]) dans un modèle de glioblastome chez le rat. Les images ont été acquises séquentiellement, en IRM, puis en TEP, et des prélèvements sanguins ont été effectués afin d’obtenir une fonction d’entrée artérielle (AIF) pour chaque molécule. Par la suite, les images obtenues avec chaque modalité ont été recalées et l’analyse pharmacocinétique a été effectuée par régions d’intérêt (ROI) et par voxel. Pour le FDG, un modèle irréversible à 3 compartiments (2 tissus) a été utilisé conformément à la littérature. Pour la FET, il a été déterminé qu’un modèle irréversible à 2 tissus pouvait être appliqué au cerveau et à la tumeur, alors qu’un modèle réversible à 2 tissus convenait aux muscles. La possibilité d’effectuer une conversion d’AIF (sanguine ou dérivée de l’image) entre le Gd DTPA et la FET, ou vice versa, a aussi été étudiée et s’est avérée faisable dans le cas des AIF sanguines obtenues à partir de l’artère caudale, comme c’est le cas pour le FDG. Finalement, l’analyse pharmacocinétique combinée IRM et TEP a relevé un lien entre la perfusion du Gd-DTPA et du FDG, ou de la FET, pour les muscles, mais elle a démontré des disparités importantes dans la tumeur. Ces résultats soulignent la complexité du microenvironnement tumoral (p. ex. coexistence de divers modes de transport pour une même molécule) et les nombreux défis rencontrées lors de sa caractérisation chez le petit animal.