A public health risk assessment for yellow fever vaccination: a model exemplified by an outbreak in the state of São Paulo, Brazil

We propose a method to analyse the 2009 outbreak in the region of Botucatu in the state of São Paulo (SP), Brazil, when 28 yellow fever (YF) cases were confirmed, including 11 deaths. At the time of the outbreak, the Secretary of Health of the State of São Paulo vaccinated one million people, causing the death of five individuals, an unprecedented number of YF vaccine-induced fatalities. We apply a mathematical model described previously to optimise the proportion of people who should be vaccinated to minimise the total number of deaths. The model was used to calculate the optimum proportion that should be vaccinated in the remaining, vaccine-free regions of SP, considering the risk of vaccine-induced fatalities and the risk of YF outbreaks in these regions.

Assessment of statistical iterative reconstructions to optimize CT examination protocols for children and young adults

La tomodensitométrie (CT) est une technique d'imagerie dont l'intérêt n'a cessé de croître depuis son apparition dans le début des années 70. Dans le domaine médical, son utilisation est incontournable à tel point que ce système d'imagerie pourrait être amené à devenir victime de son succès si son impact au niveau de l'exposition de la population ne fait pas l'objet d'une attention particulière. Bien évidemment, l'augmentation du nombre d'examens CT a permis d'améliorer la prise en charge des patients ou a rendu certaines procédures moins invasives. Toutefois, pour assurer que le compromis risque - bénéfice soit toujours en faveur du patient, il est nécessaire d'éviter de délivrer des doses non utiles au diagnostic.¦Si cette action est importante chez l'adulte elle doit être une priorité lorsque les examens se font chez l'enfant, en particulier lorsque l'on suit des pathologies qui nécessitent plusieurs examens CT au cours de la vie du patient. En effet, les enfants et jeunes adultes sont plus radiosensibles. De plus, leur espérance de vie étant supérieure à celle de l'adulte, ils présentent un risque accru de développer un cancer radio-induit dont la phase de latence peut être supérieure à vingt ans. Partant du principe que chaque examen radiologique est justifié, il devient dès lors nécessaire d'optimiser les protocoles d'acquisitions pour s'assurer que le patient ne soit pas irradié inutilement. L'avancée technologique au niveau du CT est très rapide et depuis 2009, de nouvelles techniques de reconstructions d'images, dites itératives, ont été introduites afin de réduire la dose et améliorer la qualité d'image.¦Le présent travail a pour objectif de déterminer le potentiel des reconstructions itératives statistiques pour réduire au minimum les doses délivrées lors d'examens CT chez l'enfant et le jeune adulte tout en conservant une qualité d'image permettant le diagnostic, ceci afin de proposer des protocoles optimisés.¦L'optimisation d'un protocole d'examen CT nécessite de pouvoir évaluer la dose délivrée et la qualité d'image utile au diagnostic. Alors que la dose est estimée au moyen d'indices CT (CTDIV0| et DLP), ce travail a la particularité d'utiliser deux approches radicalement différentes pour évaluer la qualité d'image. La première approche dite « physique », se base sur le calcul de métriques physiques (SD, MTF, NPS, etc.) mesurées dans des conditions bien définies, le plus souvent sur fantômes. Bien que cette démarche soit limitée car elle n'intègre pas la perception des radiologues, elle permet de caractériser de manière rapide et simple certaines propriétés d'une image. La seconde approche, dite « clinique », est basée sur l'évaluation de structures anatomiques (critères diagnostiques) présentes sur les images de patients. Des radiologues, impliqués dans l'étape d'évaluation, doivent qualifier la qualité des structures d'un point de vue diagnostique en utilisant une échelle de notation simple. Cette approche, lourde à mettre en place, a l'avantage d'être proche du travail du radiologue et peut être considérée comme méthode de référence.¦Parmi les principaux résultats de ce travail, il a été montré que les algorithmes itératifs statistiques étudiés en clinique (ASIR?, VEO?) ont un important potentiel pour réduire la dose au CT (jusqu'à-90%). Cependant, par leur fonctionnement, ils modifient l'apparence de l'image en entraînant un changement de texture qui pourrait affecter la qualité du diagnostic. En comparant les résultats fournis par les approches « clinique » et « physique », il a été montré que ce changement de texture se traduit par une modification du spectre fréquentiel du bruit dont l'analyse permet d'anticiper ou d'éviter une perte diagnostique. Ce travail montre également que l'intégration de ces nouvelles techniques de reconstruction en clinique ne peut se faire de manière simple sur la base de protocoles utilisant des reconstructions classiques. Les conclusions de ce travail ainsi que les outils développés pourront également guider de futures études dans le domaine de la qualité d'image, comme par exemple, l'analyse de textures ou la modélisation d'observateurs pour le CT.¦-¦Computed tomography (CT) is an imaging technique in which interest has been growing since it first began to be used in the early 1970s. In the clinical environment, this imaging system has emerged as the gold standard modality because of its high sensitivity in producing accurate diagnostic images. However, even if a direct benefit to patient healthcare is attributed to CT, the dramatic increase of the number of CT examinations performed has raised concerns about the potential negative effects of ionizing radiation on the population. To insure a benefit - risk that works in favor of a patient, it is important to balance image quality and dose in order to avoid unnecessary patient exposure.¦If this balance is important for adults, it should be an absolute priority for children undergoing CT examinations, especially for patients suffering from diseases requiring several follow-up examinations over the patient's lifetime. Indeed, children and young adults are more sensitive to ionizing radiation and have an extended life span in comparison to adults. For this population, the risk of developing cancer, whose latency period exceeds 20 years, is significantly higher than for adults. Assuming that each patient examination is justified, it then becomes a priority to optimize CT acquisition protocols in order to minimize the delivered dose to the patient. Over the past few years, CT advances have been developing at a rapid pace. Since 2009, new iterative image reconstruction techniques, called statistical iterative reconstructions, have been introduced in order to decrease patient exposure and improve image quality.¦The goal of the present work was to determine the potential of statistical iterative reconstructions to reduce dose as much as possible without compromising image quality and maintain diagnosis of children and young adult examinations.¦The optimization step requires the evaluation of the delivered dose and image quality useful to perform diagnosis. While the dose is estimated using CT indices (CTDIV0| and DLP), the particularity of this research was to use two radically different approaches to evaluate image quality. The first approach, called the "physical approach", computed physical metrics (SD, MTF, NPS, etc.) measured on phantoms in well-known conditions. Although this technique has some limitations because it does not take radiologist perspective into account, it enables the physical characterization of image properties in a simple and timely way. The second approach, called the "clinical approach", was based on the evaluation of anatomical structures (diagnostic criteria) present on patient images. Radiologists, involved in the assessment step, were asked to score image quality of structures for diagnostic purposes using a simple rating scale. This approach is relatively complicated to implement and also time-consuming. Nevertheless, it has the advantage of being very close to the practice of radiologists and is considered as a reference method.¦Primarily, this work revealed that the statistical iterative reconstructions studied in clinic (ASIR? and VECO have a strong potential to reduce CT dose (up to -90%). However, by their mechanisms, they lead to a modification of the image appearance with a change in image texture which may then effect the quality of the diagnosis. By comparing the results of the "clinical" and "physical" approach, it was showed that a change in texture is related to a modification of the noise spectrum bandwidth. The NPS analysis makes possible to anticipate or avoid a decrease in image quality. This project demonstrated that integrating these new statistical iterative reconstruction techniques can be complex and cannot be made on the basis of protocols using conventional reconstructions. The conclusions of this work and the image quality tools developed will be able to guide future studies in the field of image quality as texture analysis or model observers dedicated to CT.

Aspects financiers liés aux analyses ADN

La surveillance et la maîtrise des coûts liés aux analyses ADN sont très importants, surtout dans un contexte où les ressources financières ne sont pas illimitées. Des décisions quant à la priorisation des prélèvements à analyser s'imposent. La connaissance de l'aspect budgétaire est un des éléments qui peut aider à la prise de décision. Après un bref rappel des bases légales en matière d'analyses ADN en Suisse et de l'évolution législative de ces dernières années, l'article présente différents aspects financiers liés aux analyses ADN. De l'importance des sommes investies dans ces analyses et souvent méconnues dans leur globalité, même au sein des acteurs de la chaîne pénale, aux efforts faits pour optimiser les dépenses liées aux analyses ADN. Des actions sur différents facteurs qui portent sur la communication entre le laboratoire et le mandant des analyses, les changements dans le mode et les stratégies de prélèvements, ou l'utilisation de tests indicatifs moins chers qu'une analyse, nous ont permis d'optimiser nos dépenses pour obtenir plus de résultats positifs.

S-system parameter estimation for noisy metabolic profiles using newton-flow analysis.

Biochemical systems are commonly modelled by systems of ordinary differential equations (ODEs). A particular class of such models called S-systems have recently gained popularity in biochemical system modelling. The parameters of an S-system are usually estimated from time-course profiles. However, finding these estimates is a difficult computational problem. Moreover, although several methods have been recently proposed to solve this problem for ideal profiles, relatively little progress has been reported for noisy profiles. We describe a special feature of a Newton-flow optimisation problem associated with S-system parameter estimation. This enables us to significantly reduce the search space, and also lends itself to parameter estimation for noisy data. We illustrate the applicability of our method by applying it to noisy time-course data synthetically produced from previously published 4- and 30-dimensional S-systems. In addition, we propose an extension of our method that allows the detection of network topologies for small S-systems. We introduce a new method for estimating S-system parameters from time-course profiles. We show that the performance of this method compares favorably with competing methods for ideal profiles, and that it also allows the determination of parameters for noisy profiles.

Study of common database feeding with results coming from different analytical methods in the framework of illicit drugs chemical profiling

Analytical results harmonisation is investigated in this study to provide an alternative to the restrictive approach of analytical methods harmonisation which is recommended nowadays for making possible the exchange of information and then for supporting the fight against illicit drugs trafficking. Indeed, the main goal of this study is to demonstrate that a common database can be fed by a range of different analytical methods, whatever the differences in levels of analytical parameters between these latter ones. For this purpose, a methodology making possible the estimation and even the optimisation of results similarity coming from different analytical methods was then developed. In particular, the possibility to introduce chemical profiles obtained with Fast GC-FID in a GC-MS database is studied in this paper. By the use of the methodology, the similarity of results coming from different analytical methods can be objectively assessed and the utility in practice of database sharing by these methods can be evaluated, depending on profiling purposes (evidential vs. operational perspective tool). This methodology can be regarded as a relevant approach for database feeding by different analytical methods and puts in doubt the necessity to analyse all illicit drugs seizures in one single laboratory or to implement analytical methods harmonisation in each participating laboratory.

Pharmacokinetic and pharmacogenetic factors influencing plasma and cellular antiretroviral drug disposition

Abstract: The improvement in antiretroviral drug therapy has transformed HIV infection into a chronic disease. However, treatment failure and drug toxicity are frequent. Inadequate response to treatment is clearly multifactorial and, therefore, dosage individualisation based on demographic factors, genetic markers and measurement of cellular and plasma drug level may enhance both drug efficacy and tolerability. At present, antiretroviral drugs levels are monitored in plasma, whereas only drugs penetrating into cells are able to exert an antiviral activity, suggesting that cellular drug determination may more confidently reflect drug exposure at the site of pharmacological action. The overall objective of this thesis is to provide a better understanding of the Pharmacokinetic and pharmacogenetic factors influencing the plasma and cellular disposition of antiretroviral drugs. To that endeavour, analytical methods for the measurements of plasma and cellular drug levels have been developed and validated using liquid chromatography methods coupled with ultraviolet and tandem mass spectrometry detection, respectively. Correlations between plasma and cellular exposures were assessed during observational and experimental studies. Cytochrome (CYP) 2B6, efflux transporters (ABCB1, ABCC1, ABCC2 and ABCG2) and orosomucoid (ORM) polymorphisms were determined and were related to plasma and cellular exposures, as well as toxicity of antiretroviral drugs. A Pharmacokinetic population model was developed to characterise inter- and intra-patient variability of atazanavir pharmacokinetics, and to identify covariates influencing drug disposition. In that context, a Pharmacokinetic interaction study between atazanavir and lopinavir, both boosted with ritonavir, has beén conducted to assess the safety and pharmacokinetics of this boosted double-protease inhibitors regimen. Well to moderately-correlated cellular and plasma drug levels are .observed or protease inhibitors, whereas for efavirenz and nevirapine these correlations are weak. Cellular exposure, and CYP2B6 genotype (516G>T) are predictors of efavirenz neuropsychological toxicity. Nevirapine plasma exposure is also influenced by CYPZB6 polymorphism. Nelfinavir cellular exposure appears to be significantly associated only with ABCB1 genotype (3435C>T and intron 26 + 80T>C). Indinavir and lopinavir clearance and lopinavir cellular/plasma exposure ratio are influenced by the concentration of the variant S of ORM, suggesting-a specific binding of these drugs to this variant. Nelfinavir and efavirenz are not influenced by ORM concentration and phenotype. The Pharmacokinetic parameters of atazanavir are adequately described by our population model. The atazanavir-lopinavir interaction study indicates no influence on plasma and cellular atazanavir pharmacokinetics, while limited decrease in lopinavir concentrations was observed after atazanavir addition. The residual variability unexplained by the considered variables suggests that other covariates either uncontrolled at present or remaining to be identified, such as genetic and environmental factors influence antiretroviral drug pharmacokinetics, with substantial impact on treatment efficacy and tolerability. In that context, a comprehensive approach taking into account drug pharmacokinetics and patient genetic background is expected to contribute to increase treatment success, and to reduce the occurrence of adverse drug reactions by stratifying patients in an individualised antiretroviral therapy approach. Résumé Facteurs pharmacocinétiques et pharmacogénétiques influençant l'exposition plasmatique et cellulaire des antirétroviraux Les progrès de la thérapie antirétrovirale ont transformé l'infection par le VIH d'une affection mortelle à une maladie chronique. En dépit de ce succès, l'échec thérapeutique et la toxicité médicamenteuse restent fréquents. Une réponse inadéquate au traitement est clairement multifactorielle et une individualisation de la posologie des médicaments qui se baserait sur les facteurs démographiques et génétiques des patients et sur les taux sanguins des médicaments pourrait améliorer à la fois l'efficacité et la tolérance de la thérapie. Par ailleurs, seules les concentrations plasmatiques sont actuellement considérées pour le suivi thérapeutique des médicaments, alors que les taux cellulaires pourraient mieux refléter l'activité de ses médicaments qui agissent au niveau intracellulaire. L'objectif global de cette thèse était de mieux comprendre les facteurs pharmacocinétiques et pharmacocénétiques influençant l'exposition plasmatique et cellulaire des médicaments antirétroviraux. A cet effet, des méthodes pour quantifier les concentrations plasmatiques et cellulaires des antirétroviraux ont été développées et validées en utilisant la chromatographie liquide couplée à la détection ultraviolette et la spectrométrie de masse en tandem, respectivement. La corrélation entre l'exposition cellulaire et plasmatique de ces médicaments a été étudiée lors d'études observationnelles et expérimentales. Les polymorphismes du cytochrome (CYP) 2B6, ainsi que des transporteurs d'efflux (ABCB1, ABCC1, ABCC2 et ABCG2) et de l'orosomucoïde (ORM) ont été déterminés et corrélés avec l'exposition plasmatique et cellulaire des antirétroviraux, ainsi qu'à leur toxicité. Un modèle de pharmacocinétique de population a été établi afin de caractériser la variabilité inter- et intra-individuelle de l'atazanavir, et d'identifier les covariables pouvant influencer le devenir de ce médicament. Dans ce contexte, une étude d'interaction entre l'atazanavir et le lopinavir a été effectuée afin de déterminer la sécurité et le profil pharmacocinétique de ce régime thérapeutique. Des corrélations modérées à bonnes ont été observées entre les taux cellulaires et plasmatiques des inhibiteurs de protéase, alors que pour l'efavirenz et la névirapine ces corrélations sont faibles. L'exposition cellulaire, ainsi que le génotype du CYP2B6 (516G>T) sont des indices de la toxicité neuropsychologique de l'efavirenz. L'exposition plasmatique de la névirapine est également influencée par le polymorphisme du CYPZB6. L'exposition cellulaire du nelfinavir est significativement associée au génotype du ABCB1 (3435C>T et intron 26 + 80T>C). La clairance de l'indinavir et du lopinavir, ainsi que le rapport entre exposition cellulaire et plasmatique du lopinavir sont influencés par la concentration du variant S de l'ORM, suggérant une liaison spécifique de ces médicaments à ce variant. La clairance du nelfinavir et de l'efavirenz n'est pas influencée ni par la concentration ni par le phénotype de l'ORM. Les paramètres pharmacocinétiques de l'atazanavir ont été décrits de façon adéquate par le modèle de population proposé. De plus, le lopinavir n'influence pas les concentrations plasmatiques et cellulaires de l'atazanavir; alors que celui-ci conduit à une baisse limitée des taux de lopinavir. L'importante variabilité pharmacocinétique des antirétroviraux suggère que d'autres facteurs génétiques et environnementaux -qui restent encore à découvrir- influencent également leur disponibilité. Dans un proche futur, une prise en charge qui tienne. compte de la pharmacocinétique des médicaments et des caractéristiques génétiques du patient devrait permettre d'individualiser le traitement, contribuant certainement à une amélioration de la réponse thérapeutique et à une diminution de la toxicité. Résumé grand public Facteurs pharmacocinétiques et pharmacogénétiques influençant l'exposition plasmatique et cellulaire des antirétroviraux Les progrès effectués dans le traitement de l'infection par le virus de l'immunodéficience humaine acquise (VIH), ont permis de transformer une maladie avec un pronostic sombre, en une maladie chronique traitable avec des médicaments de plus en plus efficaces. Malgré ce succès, de nombreux patients ne répondent pas de façon optimale à leur traitement et/ou souffrent d'effets indésirables médicamenteux entraînant fréquemment une modification de leur thérapie. Actuellement, le suivi de la réponse au traitement s'effectue par la mesure chez les patients de la quantité de virus et du nombre des cellules immunitaires dans le sang, ainsi que par la concentration sanguine des médicaments administrés. Cependant, comme le virus se réplique à l'intérieur de la cellule, la mesure des concentrations médicamenteuses au niveau intracellulaire pourrait mieux refléter l'activité pharmacologique au site d'action. De plus, il a été possible de mettre en évidence la grande variabilité des concentrations plasmatiques de médicaments chez des patients prenant pourtant la même dose de médicament. Comme cette variabilité est notamment due à des facteurs génétiques qui sont susceptibles d'influencer la réponse au traitement antirétroviral, des analyses génétiques ont été également effectuées chez ces patients. Cette thèse a eu pour objectif de mieux comprendre les facteurs pharmacologiques et génétiques influençant l'activité et la toxicité des médicaments antirétroviraux afin de réduire la variabilité de la réponse thérapeutique. A cet effet, une méthode de dosage permettant la quantification des médicaments anti-HIV au niveau intracellulaire a été développée. Par ailleurs, nos études ont également porté .sur les variations génétiques influençant la quantité et l'activité des protéines impliquées dans le métabolisme et dans le transport des médicaments antirétroviraux. Enfin, les conséquences de ces variations sur la réponse clinique et la toxicité du traitement ont été évaluées. Nos études ont mis en évidence des associations significatives entre les variations génétiques considérées et la concentration sanguine, cellulaire et la toxicité de quelques médicaments antirétroviraux. La complémentarité des connaissances pharmacologiques, génétiques et virales pourrait aboutir à une stratégie globale permettant d'individualiser le traitement et la dose administrée, en fonction des caractéristiques propres de chaque patient. Cette approche pourrait contribuer à une optimisation du traitement antirétroviral dans la perspective d'une meilleure- efficacité thérapeutique à long terme et d'une diminution des effets indésirables rencontrés.

The achievable region approach to the optimal control of stochastic systems

The achievable region approach seeks solutions to stochastic optimisation problems by: (i) characterising the space of all possible performances(the achievable region) of the system of interest, and (ii) optimisingthe overall system-wide performance objective over this space. This isradically different from conventional formulations based on dynamicprogramming. The approach is explained with reference to a simpletwo-class queueing system. Powerful new methodologies due to the authorsand co-workers are deployed to analyse a general multiclass queueingsystem with parallel servers and then to develop an approach to optimalload distribution across a network of interconnected stations. Finally,the approach is used for the first time to analyse a class of intensitycontrol problems.

Swiss consensus statement: Recommendations for optimising re-treatment with MabThera (rituximab) in rheumatoid arthritis.

Rituximab is an effective treatment of rheumatoid arthritis (RA), which has been approved for the treatment of moderate to severe disease in patients with an inadequate response to anti-TNF therapies. Rituximab differs from other available biological agents for RA by way of its unique mode of action and unrivalled long dosing interval. The efficacy of rituximab subsides progressively over time and re-therapy is generally required to maintain long term disease control. The timing of re-treatment is currently not well established and varies widely in clinical practice. The present document is a concise recommendation regarding re-treatment with rituximab, based on validated outcomes such as the DAS28 and the EULAR response criteria. The recommendation was established through consensus between practitioners familiar with rituximab therapy in RA. Optimisation of the rituximab re-treatment schedule may improve patient outcomes and balance risks and benefits for the individual patient.

Female labour supply in Spain: The importance of behavioural assumptions and unobserved heterogeneity specification

We estimate four models of female labour supply using a Spanish sampleof married women from 1994, taking into account the complete form of theindividual s budget set. The models differ in the hypotheses relating tothe presence of optimisation errors and/or the way non-workers contributeto the likelihood function. According to the results, the effects of wagesand non-labour income on the labour supply of Spanish married women dependon the specification used. The model which has both preference andoptimisation errors and allows for both voluntarily and involuntarilyunemployed females desiring to participate seems to better fit the evidencefor Spanish married women.

Signalling, mating success and paternal investment in sticklebacks (Gasterosteus aculeatus): A theoretical model

I present an optimisation model that links paternal investment, male display and female choice. Although deviced for sticklebacks, it readily applies to other fish with male guarding behaviour. It relies on a few basic assumptions on the ways hatching success depends on paternal investment and clutch size, and male survival on paternal investment and signaling. Paternal investment is here a state-dependent decision, and signal a condition-dependent handicap by which males inform females of how much they are willing to invest. Series of predictions are derived on female and male breeding strategies, including optimal levels of signaling and paternal investment as functions of clutch size, own condition, and residual reproductive value, as well as alternative strategies such as egg kleptoparasitism. Some predictions already have empirical support, for which the present model provides new interpretations. Other might readily be tested, e.g. by simple clutch-size manipulations.

Extremity exposure in nuclear medicine: preliminary results of a European study.

The Work Package 4 of the ORAMED project, a collaborative project (2008-11) supported by the European Commission within its seventh Framework Programme, is concerned with the optimisation of the extremity dosimetry of medical staff in nuclear medicine. To evaluate the extremity doses and dose distributions across the hands of medical staff working in nuclear medicine departments, an extensive measurement programme has been started in 32 nuclear medicine departments in Europe. This was done using a standard protocol recording all relevant information for radiation exposure, i.e. radiation protection devices and tools. This study shows the preliminary results obtained for this measurement campaign. For diagnostic purposes, the two most-used radionuclides were considered: (99m)Tc and (18)F. For therapeutic treatments, Zevalin(®) and DOTATOC (both labelled with (90)Y) were chosen. Large variations of doses were observed across the hands depending on different parameters. Furthermore, this study highlights the importance of the positioning of the extremity dosemeter for a correct estimate of the maximum skin doses.

Uncertainty Quantification with Support Vector Regression Prediction Models

Uncertainty quantification of petroleum reservoir models is one of the present challenges, which is usually approached with a wide range of geostatistical tools linked with statistical optimisation or/and inference algorithms. The paper considers a data driven approach in modelling uncertainty in spatial predictions. Proposed semi-supervised Support Vector Regression (SVR) model has demonstrated its capability to represent realistic features and describe stochastic variability and non-uniqueness of spatial properties. It is able to capture and preserve key spatial dependencies such as connectivity, which is often difficult to achieve with two-point geostatistical models. Semi-supervised SVR is designed to integrate various kinds of conditioning data and learn dependences from them. A stochastic semi-supervised SVR model is integrated into a Bayesian framework to quantify uncertainty with multiple models fitted to dynamic observations. The developed approach is illustrated with a reservoir case study. The resulting probabilistic production forecasts are described by uncertainty envelopes.

Functional diversity of CD8+ T-cell subsets in chronic infection

Certains pathogènes tels que le virus de l'hépatite C ainsi que le virus de l'immunodéficience humaine sont capables de détourner les mécanismes de défense du système immunitaire adaptatif afin d'établir des infections chroniques chez l'homme. La souche clone-13 du virus de la chorioméningite lymphocytaire est utilisée comme modèle d'études d'infection virale chronique chez la souris. Les raisons qui expliqueraient la persistance de certains virus ne sont pas encore bien définies. Toutefois, il a été montré qu'une exposition prolongée à un environnement inflammatoire ainsi que la présence prolongée de l'antigène sont des facteurs qui vont déclencher un procédé de différentiation particulier des lymphocytes T CD8+. Ces cellules sur-expriment alors des récepteurs inhibiteurs tels que PD-1 tandis que leur capacité à produire des cytokines diminue. Ce procédé de différentiation est appelé « exhaustion » et serait à l'origine de la génération de réponses lymphocytaires non ou peu fonctionnelles entraînant de ce fait la persistance de l'infection. Néanmoins, il a également été démontré que ces lymphocytes maintiennent des fonctions effectrices et qu'ils permettent de limiter la réplication du virus. Afin d'étudier la fonction effectrice résiduelle de ces lymphocytes T, nous avons transféré des cellules provenant de souris infectées chroniquement dans des souris receveuses naïves qui ont à leur tour été infectées avec le virus. Grâce à ces expériences, nous avons démontré que les cellules transférées contiennent des cellules qui sont capables de i) re-proliférer, ii) de protéger les souris contre une infection virale, et de iii) survivre en l'absence d'antigène. Nous avons remarqué que les cellules stimulées de façon chronique maintiennent le même phénotype lorsqu'elles sont transférées dans des souris naïves soumises à une infection virale aiguë. Nous avons de ce fait conclu que les cellules stimulées chroniquement contiennent une sous-population de cellules qui comporte des attributs de cellules T mémoire. D'autre part, nous avons pu identifier le facteur de transcription Tcf-1 comme l'élément essential pour la génération des cellules T ressemblant à des cellules mémoires. D'autre part, nous avons étudié l'impact du niveau de stimulation via le récepteur des cellules T (TCR) sur le phénotype adopté par les lymphocytes T au cours d'une infection chronique. Dans ce but, nous avons généré des souches de virus recombinants qui expriment un épitope modifié de manière à réduire le niveau de stimulation via le TCR. D'autre part, nous avons utilisé un mélange de deux souches virales de manière à moduler spécifiquement la quantité d'un épitope tout en conservant la même charge virale. Nous avons montré que la quantité d'antigène avait plus d'influence sur le phénotype des lymphocytes T que la force d'interaction entre le complexe peptide-CMH et le TCR. De plus, l'apparition de ce phénotype ne semble pas avoir d'impact sur la prolifération des cellules en réponse à une infection primaire ou secondaire. Ainsi, nous proposons un modèle par lequel le procédé d'« exhaustion » des cellules T correspond à une différentiation cellulaire particulière qui est indépendante de la capacité de prolifération des cellules. De manière générale, ces découvertes apportent de nouvelles connaissances sur les sous-catégories de lymphocytes T CD8+ qui sont générés pendant une infection virale chronique. Nous pensons que la réponse effectrice du système immunitaire est maintenue pendant de longues périodes grâce à la présence de cellules par partagent certaines caractéristiques avec des cellules mémoires. L'étude approfondie de ces cellules peut avoir des implications importantes sur l'optimisation des stratégies utilisant l'immunothérapie pour combattre les infections chroniques et cancers.

A cell therapy approach for erythropoietin delivery using encapsulated human primary fibroblasts

Summary Cell therapy has emerged as a strategy for the treatment of various human diseases. Cells can be transplanted considering their morphological and functional properties to restore a tissue damage, as represented by blood transfusion, bone marrow or pancreatic islet cells transplantation. With the advent of the gene therapy, cells also were used as biological supports for the production of therapeutic molecules that can act either locally or at distance. This strategy represents the basis of ex vivo gene therapy characterized by the removal of cells from an organism, their genetic modification and their implantation into the same or another individual in a physiologically suitable location. The tissue or biological function damage dictates the type of cells chosen for implantation and the required function of the implanted cells. The general aim of this work was to develop an ex vivo gene therapy approach for the secretion of erythropoietin (Epo) in patients suffering from Epo-responsive anemia, thus extending to humans, studies previously performed with mouse cells transplanted in mice and rats. Considering the potential clinical application, allogeneic primary human cells were chosen for practical and safety reasons. In contrast to autologous cells, the use of allogeneic cells allows to characterize a cell lineage that can be further transplanted in many individuals. Furthermore allogeneic cells avoid the potential risk of zoonosis encountered with xenogeneic cells. Accordingly, the immune reaction against this allogeneic source was prevented by cell macro- encapsulation that prevents cell-to-cell contact with the host immune system and allows to easy retrieve the implanted device. The first step consisted in testing the survival of various human primary cells that were encapsulated and implanted for one month in the subcutaneous tissue of immunocompetent and naturally or therapeutically immunodepressed mice, assuming that xenogeneic applications constitute a stringent and representative screening before human transplantation. A fibroblast lineage from the foreskin of a young donor, DARC 3.1 cells, showed the highest mean survival score. We have then performed studies to optimize the manufacturing procedures of the encapsulation device for successful engraftment. The development of calcifications on the polyvinyl alcohol (PVA) matrix serving as a scaffold for enclosed cells into the hollow fiber devices was reported after one month in vivo. Various parameters, including matrix rinsing solutions, batches of PVA and cell lineages were assessed for their respective role in the development of the phenomenon. We observed that the calcifications could be totally prevented by using ultra-pure sterile water instead of phosphate buffer saline solution in the rinsing procedure of the PVA matrix. Moreover, a higher lactate dehydrogenase activity of the cells was found to decrease calcium depositions due to more acidic microenvironment, inhibiting the calcium precipitation. After the selection of the appropriate cell lineage and the optimization of encapsulation conditions, a retroviral-based approach was applied to DARC 3.1 fibroblasts for the transduction of the human Epo cDNA. Various modifications of the retroviral vector and the infection conditions were performed to obtain clinically relevant levels of human Epo. The insertion of a post-transcriptional regulatory element from the woodchuck hepatitis virus as well as of a Kozak consensus sequence led to a 7.5-fold increase in transgene expression. Human Epo production was further optimized by increasing the multiplicity of infection and by selecting high producer cells allowing to reach 200 IU hEpo/10E6 cells /day. These modified cells were encapsulated and implanted in vivo in the same conditions as previously described. All the mouse strains showed a sustained increase in their hematocrit and a high proportion of viable cells were observed after retrieval of the capsules. Finally, in the perspective of human application, a syngeneic model using encapsulated murine myoblasts transplanted in mice was realized to investigate the roles of both the host immune response and the cells metabolic requirements. Various loading densities and anti-inflammatory as well as immunosuppressive drugs were studied. The results showed that an immune process is responsible of cell death in capsules loaded at high cell density. A supporting matrix of PVA was shown to limit the cell density and to avoid early metabolic cell death, preventing therefore the immune reaction. This study has led to the development of encapsulated cells of human origin producing clinically relevant amounts of human EPO. This work resulted also to the optimization of cell encapsulation technical parameters allowing to begin a clinical application in end-stage renal failure patients. Résumé La thérapie cellulaire s'est imposée comme une stratégie de traitement potentiel pour diverses maladies. Si l'on considère leur morphologie et leur fonction, les cellules peuvent être transplantées dans le but de remplacer une perte tissulaire comme c'est le cas pour les transfusions sanguines ou les greffes de moelle osseuse ou de cellules pancréatiques. Avec le développement de la thérapie génique, les cellules sont également devenues des supports biologiques pour la production de molécules thérapeutiques. Cette stratégie représente le fondement de la thérapie génique ex vivo, caractérisée par le prélèvement de cellules d'un organisme, leur modification génétique et leur implantation dans le même individu ou dans un autre organisme. Le choix du type de cellule et la fonction qu'elle doit remplir pour un traitement spécifique dépend du tissu ou de la fonction biologique atteintes. Le but général de ce travail est de développer .une approche par thérapie génique ex vivo de sécrétion d'érythropoïétine (Epo) chez des patients souffrant d'anémie, prolongeant ainsi des travaux réalisés avec des cellules murines implantées chez des souris et des rats. Dans cette perpective, notre choix s'est porté sur des cellules humaines primaires allogéniques. En effet, contrairement aux cellules autologues, une caractérisation unique de cellules allogéniques peut déboucher sur de nombreuses applications. Par ailleurs, l'emploi de cellules allogéniques permet d'éviter les riques de zoonose que l'on peut rencontrer avec des cellules xénogéniques. Afin de protéger les cellules allogéniques soumises à une réaction immunitaire, leur confinement dans des macro-capsules cylindriques avant leur implantation permet d'éviter leur contact avec les cellules immunitaires de l'hôte, et de les retrouver sans difficulté en cas d'intolérance ou d'effet secondaire. Dans un premier temps, nous avons évalué la survie de différentes lignées cellulaires humaines primaires, une fois encapsulées et implantées dans le tissu sous-cutané de souris, soit immunocompétentes, soit immunodéprimées naturellement ou par l'intermédiaire d'un immunosuppresseur. Ce modèle in vivo correspond à des conditions xénogéniques et représente par conséquent un environnement de loin plus hostile pour les cellules qu'une transplantation allogénique. Une lignée fibroblastique issue du prépuce d'un jeune enfant, nommée DARC 3 .1, a montré une remarquable résistance avec un score de survie moyen le plus élevé parmi les lignées testées. Par la suite, nous nous sommes intéressés aux paramètres intervenant dans la réalisation du système d'implantation afin d'optimaliser les conditions pour une meilleure adaptation des cellules à ce nouvel environnement. En effet, en raison de l'apparition, après un mois in vivo, de calcifications au niveau de la matrice de polyvinyl alcohol (PVA) servant de support aux cellules encapsulées, différents paramètres ont été étudiés, tels que les procédures de fabrication, les lots de PVA ou encore les lignées cellulaires encapsulées, afin de mettre en évidence leur rôle respectif dans la survenue de ce processus. Nous avons montré que l'apparition des calcifications peut être totalement prévenue par l'utilisation d'eau pure au lieu de tampon phosphaté lors du rinçage des matrices de PVA. De plus, nous avons observe qu'un taux de lactate déshydrogénase cellulaire élevé était corrélé avec une diminution des dépôts de calcium au sein de la matrice en raison d'un micro-environnement plus acide inhibant la précipitation du calcium. Après sélection de la lignée cellulaire appropriée et de l'optimisation des conditions d'encapsulation, une modification génétique des fibroblastes DARC 3.1 a été réalisée par une approche rétrovirale, permettant l'insertion de l'ADN du gène de l'Epo dans le génome cellulaire. Diverses modifications, tant au niveau génétique qu'au niveau des conditions d'infection, ont été entreprises afin d'obtenir des taux de sécrétion d'Epo cliniquement appropriés. L'insertion dans la séquence d'ADN d'un élément de régulation post¬transcriptionnelle dérivé du virus de l'hépatite du rongeur (« woodchuck ») ainsi que d'une séquence consensus appelée « Kozak » ont abouti à une augmentation de sécrétion d'Epo 7.5 fois plus importante. De même, l'optimisation de la multiplicité d'infection et la sélection plus drastique des cellules hautement productrices ont permis finalement d'obtenir une sécrétion correspondant à 200 IU d'Epo/10E6 cells/jour. Ces cellules génétiquement modifiées ont été encapsulées et implantées in vivo dans les mêmes conditions que celles décrites plus haut. Toutes les souris transplantées ont montré une augmentation significative de leur hématocrite et une proportion importante de cellules présentait une survie conservée au moment de l'explantation des capsules. Finalement, dans la perspective d'une application humaine, un modèle syngénique a été proposé, basé sur l'implantation de myoblastes murins encapsulés dans des souris, afin d'investiguer les rôles respectifs de la réponse immunitaire du receveur et des besoins métaboliques cellulaires sur leur survie à long terme. Les cellules ont été encapsulées à différentes densités et les animaux transplantés se sont vus administrer des injections de molécules anti-inflammatoires ou immunosuppressives. Les résultats ont démontré qu'une réaction immunologique péri-capsulaire était à la base du rejet cellulaire dans le cas de capsules à haute densité cellulaire. Une matrice de PVA peut limiter cette densité et éviter une mort cellulaire précoce due à une insuffisance métabolique et par conséquent prévenir la réaction immunitaire. Ce travail a permis le développement de cellules encapsulées d'origine humaine sécrétant des taux d'Epo humaine adaptés à des traitements cliniques. De pair avec l'optimalisation des paramètres d'encapsulation, ces résultats ont abouti à l'initiation d'une application clinique destinée à des patients en insuffisance rénale terminale.

Time-lapse and probabilistic inversion strategies for plane-wave electromagnetic data.

L'utilisation efficace des systèmes géothermaux, la séquestration du CO2 pour limiter le changement climatique et la prévention de l'intrusion d'eau salée dans les aquifères costaux ne sont que quelques exemples qui démontrent notre besoin en technologies nouvelles pour suivre l'évolution des processus souterrains à partir de la surface. Un défi majeur est d'assurer la caractérisation et l'optimisation des performances de ces technologies à différentes échelles spatiales et temporelles. Les méthodes électromagnétiques (EM) d'ondes planes sont sensibles à la conductivité électrique du sous-sol et, par conséquent, à la conductivité électrique des fluides saturant la roche, à la présence de fractures connectées, à la température et aux matériaux géologiques. Ces méthodes sont régies par des équations valides sur de larges gammes de fréquences, permettant détudier de manières analogues des processus allant de quelques mètres sous la surface jusqu'à plusieurs kilomètres de profondeur. Néanmoins, ces méthodes sont soumises à une perte de résolution avec la profondeur à cause des propriétés diffusives du champ électromagnétique. Pour cette raison, l'estimation des modèles du sous-sol par ces méthodes doit prendre en compte des informations a priori afin de contraindre les modèles autant que possible et de permettre la quantification des incertitudes de ces modèles de façon appropriée. Dans la présente thèse, je développe des approches permettant la caractérisation statique et dynamique du sous-sol à l'aide d'ondes EM planes. Dans une première partie, je présente une approche déterministe permettant de réaliser des inversions répétées dans le temps (time-lapse) de données d'ondes EM planes en deux dimensions. Cette stratégie est basée sur l'incorporation dans l'algorithme d'informations a priori en fonction des changements du modèle de conductivité électrique attendus. Ceci est réalisé en intégrant une régularisation stochastique et des contraintes flexibles par rapport à la gamme des changements attendus en utilisant les multiplicateurs de Lagrange. J'utilise des normes différentes de la norme l2 pour contraindre la structure du modèle et obtenir des transitions abruptes entre les régions du model qui subissent des changements dans le temps et celles qui n'en subissent pas. Aussi, j'incorpore une stratégie afin d'éliminer les erreurs systématiques de données time-lapse. Ce travail a mis en évidence l'amélioration de la caractérisation des changements temporels par rapport aux approches classiques qui réalisent des inversions indépendantes à chaque pas de temps et comparent les modèles. Dans la seconde partie de cette thèse, j'adopte un formalisme bayésien et je teste la possibilité de quantifier les incertitudes sur les paramètres du modèle dans l'inversion d'ondes EM planes. Pour ce faire, je présente une stratégie d'inversion probabiliste basée sur des pixels à deux dimensions pour des inversions de données d'ondes EM planes et de tomographies de résistivité électrique (ERT) séparées et jointes. Je compare les incertitudes des paramètres du modèle en considérant différents types d'information a priori sur la structure du modèle et différentes fonctions de vraisemblance pour décrire les erreurs sur les données. Les résultats indiquent que la régularisation du modèle est nécessaire lorsqu'on a à faire à un large nombre de paramètres car cela permet d'accélérer la convergence des chaînes et d'obtenir des modèles plus réalistes. Cependent, ces contraintes mènent à des incertitudes d'estimations plus faibles, ce qui implique des distributions a posteriori qui ne contiennent pas le vrai modèledans les régions ou` la méthode présente une sensibilité limitée. Cette situation peut être améliorée en combinant des méthodes d'ondes EM planes avec d'autres méthodes complémentaires telles que l'ERT. De plus, je montre que le poids de régularisation des paramètres et l'écart-type des erreurs sur les données peuvent être retrouvés par une inversion probabiliste. Finalement, j'évalue la possibilité de caractériser une distribution tridimensionnelle d'un panache de traceur salin injecté dans le sous-sol en réalisant une inversion probabiliste time-lapse tridimensionnelle d'ondes EM planes. Etant donné que les inversions probabilistes sont très coûteuses en temps de calcul lorsque l'espace des paramètres présente une grande dimension, je propose une stratégie de réduction du modèle ou` les coefficients de décomposition des moments de Legendre du panache de traceur injecté ainsi que sa position sont estimés. Pour ce faire, un modèle de résistivité de base est nécessaire. Il peut être obtenu avant l'expérience time-lapse. Un test synthétique montre que la méthodologie marche bien quand le modèle de résistivité de base est caractérisé correctement. Cette méthodologie est aussi appliquée à un test de trac¸age par injection d'une solution saline et d'acides réalisé dans un système géothermal en Australie, puis comparée à une inversion time-lapse tridimensionnelle réalisée selon une approche déterministe. L'inversion probabiliste permet de mieux contraindre le panache du traceur salin gr^ace à la grande quantité d'informations a priori incluse dans l'algorithme. Néanmoins, les changements de conductivités nécessaires pour expliquer les changements observés dans les données sont plus grands que ce qu'expliquent notre connaissance actuelle des phénomenès physiques. Ce problème peut être lié à la qualité limitée du modèle de résistivité de base utilisé, indiquant ainsi que des efforts plus grands devront être fournis dans le futur pour obtenir des modèles de base de bonne qualité avant de réaliser des expériences dynamiques. Les études décrites dans cette thèse montrent que les méthodes d'ondes EM planes sont très utiles pour caractériser et suivre les variations temporelles du sous-sol sur de larges échelles. Les présentes approches améliorent l'évaluation des modèles obtenus, autant en termes d'incorporation d'informations a priori, qu'en termes de quantification d'incertitudes a posteriori. De plus, les stratégies développées peuvent être appliquées à d'autres méthodes géophysiques, et offrent une grande flexibilité pour l'incorporation d'informations additionnelles lorsqu'elles sont disponibles. -- The efficient use of geothermal systems, the sequestration of CO2 to mitigate climate change, and the prevention of seawater intrusion in coastal aquifers are only some examples that demonstrate the need for novel technologies to monitor subsurface processes from the surface. A main challenge is to assure optimal performance of such technologies at different temporal and spatial scales. Plane-wave electromagnetic (EM) methods are sensitive to subsurface electrical conductivity and consequently to fluid conductivity, fracture connectivity, temperature, and rock mineralogy. These methods have governing equations that are the same over a large range of frequencies, thus allowing to study in an analogous manner processes on scales ranging from few meters close to the surface down to several hundreds of kilometers depth. Unfortunately, they suffer from a significant resolution loss with depth due to the diffusive nature of the electromagnetic fields. Therefore, estimations of subsurface models that use these methods should incorporate a priori information to better constrain the models, and provide appropriate measures of model uncertainty. During my thesis, I have developed approaches to improve the static and dynamic characterization of the subsurface with plane-wave EM methods. In the first part of this thesis, I present a two-dimensional deterministic approach to perform time-lapse inversion of plane-wave EM data. The strategy is based on the incorporation of prior information into the inversion algorithm regarding the expected temporal changes in electrical conductivity. This is done by incorporating a flexible stochastic regularization and constraints regarding the expected ranges of the changes by using Lagrange multipliers. I use non-l2 norms to penalize the model update in order to obtain sharp transitions between regions that experience temporal changes and regions that do not. I also incorporate a time-lapse differencing strategy to remove systematic errors in the time-lapse inversion. This work presents improvements in the characterization of temporal changes with respect to the classical approach of performing separate inversions and computing differences between the models. In the second part of this thesis, I adopt a Bayesian framework and use Markov chain Monte Carlo (MCMC) simulations to quantify model parameter uncertainty in plane-wave EM inversion. For this purpose, I present a two-dimensional pixel-based probabilistic inversion strategy for separate and joint inversions of plane-wave EM and electrical resistivity tomography (ERT) data. I compare the uncertainties of the model parameters when considering different types of prior information on the model structure and different likelihood functions to describe the data errors. The results indicate that model regularization is necessary when dealing with a large number of model parameters because it helps to accelerate the convergence of the chains and leads to more realistic models. These constraints also lead to smaller uncertainty estimates, which imply posterior distributions that do not include the true underlying model in regions where the method has limited sensitivity. This situation can be improved by combining planewave EM methods with complimentary geophysical methods such as ERT. In addition, I show that an appropriate regularization weight and the standard deviation of the data errors can be retrieved by the MCMC inversion. Finally, I evaluate the possibility of characterizing the three-dimensional distribution of an injected water plume by performing three-dimensional time-lapse MCMC inversion of planewave EM data. Since MCMC inversion involves a significant computational burden in high parameter dimensions, I propose a model reduction strategy where the coefficients of a Legendre moment decomposition of the injected water plume and its location are estimated. For this purpose, a base resistivity model is needed which is obtained prior to the time-lapse experiment. A synthetic test shows that the methodology works well when the base resistivity model is correctly characterized. The methodology is also applied to an injection experiment performed in a geothermal system in Australia, and compared to a three-dimensional time-lapse inversion performed within a deterministic framework. The MCMC inversion better constrains the water plumes due to the larger amount of prior information that is included in the algorithm. The conductivity changes needed to explain the time-lapse data are much larger than what is physically possible based on present day understandings. This issue may be related to the base resistivity model used, therefore indicating that more efforts should be given to obtain high-quality base models prior to dynamic experiments. The studies described herein give clear evidence that plane-wave EM methods are useful to characterize and monitor the subsurface at a wide range of scales. The presented approaches contribute to an improved appraisal of the obtained models, both in terms of the incorporation of prior information in the algorithms and the posterior uncertainty quantification. In addition, the developed strategies can be applied to other geophysical methods, and offer great flexibility to incorporate additional information when available.