Les mouvements de la prostate relativement aux ganglions pelviens durant un traitement de radiothérapie peut affecter la qualité du traitement et limiter la protection des organes à risque. Ce projet explore trois traitements adaptatifs qui combinent à la fois l’information acquise avec un recalage osseux et un recalage fait sur les marqueurs radio-opaques implantés dans la prostate. L’efficacité de ces traitements selon les mouvements interfractionnaires prostatiques a été évaluée. Cinq patients constituent la base de données de cette étude rétrospective. Chacun possède une tomodensitométrie initiale (CT) en plus de tomodensitométries à faisceau conique (CBCT) acquises lors des jours de traitement. Pour permettre le calcul de dose, un recalage non-rigide du CT a été réalisé sur chacun des CBCT. Les résultats montrent qu’une méthode adaptative peut être choisie pour augmenter le ratio thérapeutique. Ces méthodes semblent être de bons choix pour une stratégie de traitement adaptatif simple et efficace.

OBJECTIF: Déterminer si un régime de Bicalutamide et de Dutasteride est d’efficacité similaire à un régime standard pour la cytoréduction prostatique avant une curiethérapie prostatique par implants permanents et comparer la toxicité de ces deux régimes. MÉTHODE Les patients avec un adénocarcinome de la prostate traités par curiethérapie prostatique ayant un volume prostatique initial de plus de 50 cc ont été randomisés soit à un régime de Bicalutamide 50 mg die, de Dutasteride 0.5 mg die et de tamoxifen 10 mg die pour prévenir la gynécomastie (Groupe expérimental = GE) versus un agoniste de la LHRH (dose de 3 mois) avec du Bicalutamide 50 mg die pour 1 mois (Groupe contrôle = GC). Le volume prostatique total (VPT) a été mesuré à la randomisation (t1) et environ 2.5 mois après le début du traitement de cytoréduction (t2). La curiethérapie par implant permanent d’I-125 à été réalisé environ 3 mois après la randomisation. Le changement de VPT relatif en % à été comparé entre les groupes. Des questionnaires EPIC (Expanded Prostate Index Composite) ont été complétés au t1, t2, ainsi qu’à 6 semaines (t3), 3 mois (t4) et 6 mois (t5) post-curiethérapie. Le profil de changement des scores sexuels et hormonaux du questionnaire EPIC a été évalué de façon qualitative. RÉSULTATS Les premiers 11 patients randomisés dans l’étude ont été analysés (GE n=5, GC ; n=6). Les caractéristiques de base entre les deux groups étaient similaires. Tous les patients présentaient des caractéristiques de cancer de prostate de bas risque (APS< 10, Stade clinique ≤T2a et un score de Gleason ≤6). Le VPT médian au t1 était de 51,00 cc dans le GE et de 55,56 cc dans le GC. Les diminutions médianes entre le volume t1 et le volume t2 en absolu étaient de 18,46 cc dans le GE et de 24,61 cc dans le GC (valeur-p au test exact de Wilcoxon 0,03). Pour ce qui est de la médiane de changement relatif de volume entre t1 et le volume au t2, celles-ci étaient de 36,44 % dans le GE et de 40,83 % dans le GC (valeur-p 0,25). Tous les 11 patients ont pu subir la curiethérapie sans interférence de l’arche pubienne et avec des paramètres dosimétriques acceptables. Les médianes de scores de fonction sexuelle et de sommaire hormonal semblent qualitativement démontrer un meilleur profil au sein du GE. DISCUSSION Cette analyse exploratoire n’a pas la puissance suffisante pour permettre de rejeter ou d’affirmer une hypothèse de non-infériorité du régime expérimental du point de vue de la cytoréduction prostatique. Toutefois, ce régime a permis à tous les patients de subir la curiethérapie sans interférence de l’arche pubienne et avec des paramètres dosimétriques acceptables. Le Bicalutamide et le Dutasteride ont aussi démontré un profil de toxicité favorable. Les résultats de cette analyse précoce justifient donc de continuer à randomiser des patients au sein de cette étude mais ne permettent pas de tirer des conclusions définitives en raison de la petite taille d’échantillon. / OBJECTIVE: To determine if a regimen of Bicalutamide and Dutasteride is of similar efficacy compared to a standard regimen for prostate cytoreduction prior to permanent implant prostate brachytherapy and to compare toxicity for both of these regimens MATERIALS AND METHODS Patients with an adenocarcinoma of the prostate treated with permanent implant brachytherapy having an initial total prostate volume of 50 cc or more were randomized to either Bicalutamide 50 mg daily, Dutasteride 0.5 mg daily and Tamoxifen 10 mg daily (Experimental group = EG) or an LHRH agonist (3 month dose) and Bicalutamide 50 mg daily for one month (Control Group = CG). Total prostate volume (TPV) was evaluated at time of randomization (t1) and approximately 2.5 months after the start of cytoreductive therapy (t2). Permanent implant I-125 brachytherapy was completed 3 months after randomization. Relative TPV change in % and absolute TPV change in cc were compared between groups. Expanded Prostate Index Composite questionnaires were completed by patients at t1, t2, as well as at 6 weeks (t3), 3 months (t4), and 6 months (t5). Profile changes in EPIC sexual and hormonal domains were evaluated qualitatively. RESULTS The first 11 patients registered for the study were analyzed (EG n=5; CG n=6) Baseline parameters were similar between groups. All patients had low risk characteristics (PSA < 10, Clinical stage ≤T2a and Gleason score ≤6). Median baseline TPV was 51.00 cc in EG and 55.56 in CG. Median reduction in absolute TPV was 18.46 cc in EG and 24.61 in CG (Wilcoxon exact test p=0.03). As for median relative change in TPV, this was 36.44 % in EG and 24.61 cc in CG (p=0.25). All 11 patients were able to undergo brachytherapy without pubic arch interference and with acceptable dosimetric parameters. Median EPIC of both sexual and hormonal domains of the EPIC questionnaire seemed to favor the EG over time. DISCUSSION This early exploratory analysis does not have sufficient power to reject or accept the hypothesis of non-inferiority the experimental regimen for prostate volume reduction. Nonetheless, this regimen allowed brachytherapy without pubic arch interference and with acceptable dosimetric parameters. The experimental regimen of Bicalutamide and Dutasteride also seemed to show a favorable toxicity profile. The results of this early analysis thus justify continued randomization for the study but do not allow definitive conclusions due to small sample size.

Cette étude porte sur le développement d’un outil de contrôle qualité basé sur l’expérience, dérivé du concept de frontière stochastique en économie et s’appuyant sur des connaissances géométriques spécifiques au patient pour améliorer la qualité des traitements de curiethérapie à haut débit de dose pour le cancer de la prostate. Cent plans cliniques de curiethérapie à haut débit de dose de la prostate ont été utilisés dans cette étude, dans laquelle l’échographie transrectale était la seule modalité d’imagerie. Une fraction unique de 15 Gy a était prescrite à tous ces patients. Un algorithme de recuit simulé de planification inverse a été appliqué pour réaliser tous les plans et Oncentra Prostate a été employé comme système d’imagerie et de planification du traitement en temps réel. Les recommandations relatives aux paramètres de dose de la société américaine de curiethérapie pour la cible et les organes à risque ont été suivies. Les relations entre les paramètres géométriques et les paramètres dosimétriques d’intérêt sont examinées. Les paramètres géométriques sont liés aux dimensions anatomiques des patients et ceux associés aux cathéters. Pour déterminer les paramètres géométriques dominants dans un modèle de frontière stochastique donné, les relations monotones entre les paramètres géométriqueset les paramètres dosimétriques d’intérêt sont mesurées avec une approche non paramétrique, à savoir le coefficient de corrélation de Spearman. Ensuite, une recherche de force brute est effectuée pour un modèle donné dans lequel différents modèles, incluant toutes les combinaisons possibles des paramètres géométriques dominantes, sont optimisés. L’optimisation est accomplie en utilisant une méthode de vraisemblance maximale implémentée dans le progiciel de calcul statistique R, avec son algorithme de recuit simulée généralisée. Le test du rapport de vraisemblance et sa valeur-p correspondante sont utilisés pour comparer la signification statistique de l’ajout de nouveaux paramètres géométriques aux modèles. Un modèle de production pour la cible et un modèle de coût pour chacun des organes à risque sont développés pour le traitement par curiethérapie à haut débit de dose guidé par l’échographie transrectale. De plus, pour valider si chacun des modèles développés est universel, nous l’appliquons à une autre catégorie de traitement de la curiethérapie à haut débit de dose, dans laquelle la tomodensitométrie était utilisée comme modalité d’imagerie plutôt que de l’échographie transrectale. Ainsi, une nouvelle cohorte de cent plans cliniques curiethérapie à haut débit de dose guidés par la tomodensitométrie est prise en compte. Un modèle de frontière stochastique de production pour la cible et trois modèles de coût pour les organes à risque basés sur la tomodensitométrie sont développés. Enfin, les modèles intégrés de la tomodensitométrie et de l’échographie transrectale sont comparés. / This thesis focuses on developing an experience-based quality control (QC) tool, derived from the concept of stochastic frontier (SF) analysis in economics and based on patient-specific geometric knowledge to improve the quality of the high-dose-rate brachytherapy (HDR-BT) treatment for prostate cancer. One hundred clinical HDR prostate BT plans, using the transrectal ultrasound (TRUS) asthe only imaging modality, all treated with a single fraction of 15 Gy, and made using Inverse PlanningSimulated Annealing (IPSA) algorithm, are studied. Also, Oncentra Prostate system is employed as the real-time 3D prostate imaging and treatment planning system. American Brachytherapy Society dose parameter recommendations for target and organs at risk (OARs) were followed. Relationships between all the different geometric parameters (GPs) and the four dosimetric parameters (DPs) V100 of the prostate, V75 of the bladder and rectum, and D10 of the urethra were examined. Geometricinformation of the patients and catheters are considered as different GPs. To find the dominant GPs in a given SF model, monotonic relationships between the GPs and DPs of interest are measured using a nonparametric approach: the Spearman correlation coefficient. Then, to determine the optimal SF model for each of the target production SF, and the OARs cost SF models, brute-force searches are performed. Different SF models including all the possible combinations of the dominant GPs in the SF model under study are optimized. Optimization is done using a maximum likelihood method implemented in the statistical computing package R, along with its Generalized Simulated Annealing algorithm. The likelihood ratio test and its corresponding p-value are used to compare the statistical significance of adding new GPs to SF models. A production SF (PSF) model for the target, and a costSF (CSF) model for each of the bladder, rectum, and urethra are developed for TRUS-guided HDR-BTtreatment. The difference between the dose value of a plan obtained by IPSA and the one predicted by an SF model is explored. Additionally, to verify if each of the models developed for the TRUS-guided category of the HDR-BT treatment for prostate is universal, we apply it on another category of HDR-BT treatment, in which computed tomography (CT) was used as the imaging modality. So, a different cohort of one hundred clinical CT-guided HDR-BT plans is taken into consideration. A target production SF and three OARs cost SF models are developed for the CT-based plans. Subsequently, the built-in SF models for the TRUS-based and CT-based plans are compared.

L’objectif général de cette thèse est d’utiliser les connaissances en physique de la radiation, en programmation informatique et en équipement informatique à la haute pointe de la technologie pour améliorer les traitements du cancer. En particulier, l’élaboration d’un plan de traitement en radiothérapie peut être complexe et dépendant de l’utilisateur. Cette thèse a pour objectif de simplifier la planification de traitement actuelle en curiethérapie de la prostate à haut débit de dose (HDR). Ce projet a débuté à partir d’un algorithme de planification inverse largement utilisé, la planification de traitement inverse par recuit simulé (IPSA). Pour aboutir à un algorithme de planification inverse ultra-rapide et automatisé, trois algorithmes d’optimisation multicritères (MCO) ont été mis en oeuvre. Suite à la génération d’une banque de plans de traitement ayant divers compromis avec les algorithmes MCO, un plan de qualité a été automatiquement sélectionné. Dans la première étude, un algorithme MCO a été introduit pour explorer les frontières de Pareto en curiethérapie HDR. L’algorithme s’inspire de la fonctionnalité MCO intégrée au système Raystation (RaySearch Laboratories, Stockholm, Suède). Pour chaque cas, 300 plans de traitement ont été générés en série pour obtenir une approximation uniforme de la frontière de Pareto. Chaque plan optimal de Pareto a été calculé avec IPSA et chaque nouveau plan a été ajouté à la portion de la frontière de Pareto où la distance entre sa limite supérieure et sa limite inférieure était la plus grande. Dans une étude complémentaire, ou dans la seconde étude, un algorithme MCO basé sur la connaissance (kMCO) a été mis en oeuvre pour réduire le temps de calcul de l’algorithme MCO. Pour ce faire, deux stratégies ont été mises en oeuvre : une prédiction de l’espace des solutions cliniquement acceptables à partir de modèles de régression et d’un calcul parallèle des plans de traitement avec deux processeurs à six coeurs. En conséquence, une banque de plans de traitement de petite taille (14) a été générée et un plan a été sélectionné en tant que plan kMCO. L’efficacité de la planification et de la performance dosimétrique ont été comparées entre les plans approuvés par le médecin et les plans kMCO pour 236 cas. La troisième et dernière étude de cette thèse a été réalisée en coopération avec Cédric Bélanger. Un algorithme MCO (gMCO) basé sur l’utilisation d’un environnement de développement compatible avec les cartes graphiques a été mis en oeuvre pour accélérer davantage le calcul. De plus, un algorithme d’optimisation quasi-Newton a été implémenté pour remplacer le recuit simulé dans la première et la deuxième étude. De cette manière, un millier de plans de traitement avec divers compromis et équivalents à ceux générés par IPSA ont été calculés en parallèle. Parmi la banque de plans de traitement généré par l’agorithme gMCO, un plan a été sélectionné (plan gMCO). Le temps de planification et les résultats dosimétriques ont été comparés entre les plans approuvés par le médecin et les plans gMCO pour 457 cas. Une comparaison à grande échelle avec les plans approuvés par les radio-oncologues montre que notre dernier algorithme MCO (gMCO) peut améliorer l’efficacité de la planification du traitement (de quelques minutes à 9:4 s) ainsi que la qualité dosimétrique des plans de traitements (des plans passant de 92:6% à 99:8% selon les critères dosimétriques du groupe de traitement oncologique par radiation (RTOG)). Avec trois algorithmes MCO mis en oeuvre, cette thèse représente un effort soutenu pour développer un algorithme de planification inverse ultra-rapide, automatique et robuste en curiethérapie HDR. / The overall purpose of this thesis is to use the knowledge of radiation physics, computer programming and computing hardware to improve cancer treatments. In particular, designing a treatment plan in radiation therapy can be complex and user-dependent, and this thesis aims to simplify current treatment planning in high dose rate (HDR) prostate brachytherapy. This project was started from a widely used inverse planning algorithm, Inverse Planning Simulated Annealing (IPSA). In order to eventually lead to an ultra-fast and automatic inverse planning algorithm, three multi-criteria optimization (MCO) algorithms were implemented. With MCO algorithms, a desirable plan was selected after computing a set of treatment plans with various trade-offs. In the first study, an MCO algorithm was introduced to explore the Pareto surfaces in HDR brachytherapy. The algorithm was inspired by the MCO feature integrated in the Raystation system (RaySearch Laboratories, Stockholm, Sweden). For each case, 300 treatment plans were serially generated to obtain a uniform approximation of the Pareto surface. Each Pareto optimal plan was computed with IPSA, and each new plan was added to the Pareto surface portion where the distance between its upper boundary and its lower boundary was the largest. In a companion study, or the second study, a knowledge-based MCO (kMCO) algorithm was implemented to shorten the computation time of the MCO algorithm. To achieve this, two strategies were implemented: a prediction of clinical relevant solution space with previous knowledge, and a parallel computation of treatment plans with two six-core CPUs. As a result, a small size (14) plan dataset was created, and one plan was selected as the kMCO plan. The planning efficiency and the dosimetric performance were compared between the physician-approved plans and the kMCO plans for 236 cases. The third and final study of this thesis was conducted in cooperation with Cédric Bélanger. A graphics processing units (GPU) based MCO (gMCO) algorithm was implemented to further speed up the computation. Furthermore, a quasi-Newton optimization engine was implemented to replace simulated annealing in the first and the second study. In this way, one thousand IPSA equivalent treatment plans with various trade-offs were computed in parallel. One plan was selected as the gMCO plan from the calculated plan dataset. The planning time and the dosimetric results were compared between the physician-approved plans and the gMCO plans for 457 cases. A large-scale comparison against the physician-approved plans shows that our latest MCO algorithm (gMCO) can result in an improved treatment planning efficiency (from minutes to 9:4 s) as well as an improved treatment plan dosimetric quality (Radiation Therapy Oncology Group (RTOG) acceptance rate from 92.6% to 99.8%). With three implemented MCO algorithms, this thesis represents a sustained effort to develop an ultra-fast, automatic and robust inverse planning algorithm in HDR brachytherapy.

Jusqu'à ce jour, très peu de travaux empiriques ont évalué les aspects psychologiques associés au cancer de la prostate. Cette étude transversale vise donc à comparer la prévalence des difficultés psychologiques (i.e., anxiété et dépression), psychophysiologiques (i.e., insomnie, fatigue) et des difficultés sexuelles en fonction de trois modalités de traitement pour le cancer de la prostate (radiothérapie, curiethérapie et prostatectomie radicale) et du temps écoulé depuis le diagnostic. L'échantillon est composé de 861 hommes ayant complété une batterie de questionnaires. Les résultats démontrent que les difficultés sexuelles constituent le problème le plus fréquent suivies de l'insomnie et de l'anxiété. De plus, les symptômes psychologiques et psychophysiologiques restent relativement stables dans le temps alors que les symptômes physiques tendent à se résorber. Enfin, les patients ayant reçu de la radiothérapie ont un risque plus élevé de souffrir d'humeur dépressive, de détresse psychologique globale et de fatigue comparativement aux patients n'ayant pas reçu ce traitement.

Les calculs de dose en curiethérapie à bas débit reposent depuis 1995 sur un formalisme qui considère le milieu irradié comme étant homogène et constitué d’eau afin de mieux accommoder les temps de calcul clinique. Ce travail présente une amélioration de ce formalisme dosimétrique issu du protocole TG-43 de l’AAPM (American Association of Physicists in Medicine), utilisé actuellement dans les systèmes de planification de traitement (SPT) clinique. Avec l’ajout d’un tracé radiologique pour tenir compte des hétérogénéités et en séparant la dose primaire de la dose diffusée, il est possible de raffiner les calculs dosimétriques. Or, cette modification au calcul le rendant plus complexe était coûteuse en temps d’exécution jusqu’à très récemment. Elle se voit matérialisée dans ce travail avec l’arrivée récente de la technologie GPGPU, les calculs scientifiques sur périphériques graphiques. En exploitant le parallélisme des calculs de dose en curiethérapie à bas débit pour les implants permanents, l’algorithme baptisé TG-43-RT pour TG-43 avec tracé de rayons (Ray-Tracing en anglais) permet d’obtenir des facteurs d’accélération de l’ordre de 103 par rapport au CPU. Seulement 0.5 s par source est requis dans une géométrie de 1003 voxels, et les dosimétries d’anatomies hétérogènes sont sensiblement améliorées. Le TG-43-RT corrige les dépôts de dose en aval des hétérogénéités et réduit les effets indésirables de l’atténuation inter-sources. Des écarts de dose de plus de 80% au-delà d’une calcification avaient de quoi influencer les HDV, or, le TG-43-RT ramène ces déviations en-deçà de l’incertitude de 8.7% (2σ) concédée par les SPT, tel qu’établi dans le TG-138. Les simulations Monte Carlo ont été employées pour servir de référence absolue à la méthode développée et quantifier les améliorations dosimétriques relatives à l’actuel calcul de dose. L’algorithme a été testé d’un point de vue clinique dans un fantôme de prostate avec des sources d’125I et dans un fantôme de sein avec des sources de 103Pd pour unifier les corrections apportées. / Brachytherapy dose calculations have been relying since 1995 on a formalism that considers the whole geometry as a homogeneous water tank. This gives the opportunity to compute dose distributions within a reasonable clinical timeframe, but with considerable approximations that can influence or even change the treatment. This work presents an upgraded version of this formalism derived from the TG-43 protocol (AAPM), currently used in actual treatment planning systems. The new algorithm includes a primary and scatter dose separation using ray-tracing operation to account for heterogeneities through the medium. These modifications to the dose calculation are very power consuming and too long for the clinical needs when executed on modern CPUs. With the GPGPU technology, a GPU-driven algorithm allows a complex handling of the anatomic heterogeneities in the dose calculation and keeps execution times below 0.5 s/seed. The algorithm was named TG-43-RT for TG-43 with Ray-Tracing and accelerations factors of three order of magnitude were obtained over a regular CPU implementation. The TG-43-RT algorithm adequatly corrects dose deposition along heterogeneities and it reduces interseed attenuation effects. Dose is scored in the medium instead of water, which also applies a severe correction for high heterogeneous medium like the breast. Deviations of more than 80% in dose deposition were obtained over calcification inside the prostate when compared to Monte Carlo simulations. This kind of deviation influences DVH shape. However, the TG-43-RT was keeping those deviations within the 8.7% uncertainty range (2σ) associated with a regular TG-43 evaluation as stated in the TG-138. Monte Carlo simulations were used as a gold standard to get an absolute dose calculation reference. The algorithm was also tested in a prostate phantom with 125I seeds and in a breast phantom with 103Pd seeds to mimic real anatomic geometries and unify the applied corrections.