Studies on the exaggerated inflammatory response caused by streptococcus suis at systemic and central nervous system levels

Streptococcus suis de type 2 est un microorganisme pathogène d’importance chez le porc. Il est la cause de différentes pathologies ayant comme caractéristique commune la méningite. C’est également un agent émergeant de zoonose : des cas cliniques humains ont récemment été rapportés en Asie. Cependant, la pathogénèse de S. suis n’est pas encore complètement élucidée. Jusqu’à présent, la réponse pro-inflammatoire initiée par S. suis n’a été étudiée qu’in vitro. L’étude du choc septique et de la méningite requiert toujours des modèles expérimentaux appropriés. Au cours de cette étude, nous avons développé un modèle in vivo d’infection chez la souris qui utilise la voie d’inoculation intra-péritonéale. Ce modèle a servi à l’étude de la réponse pro-inflammatoire associée à ce pathogène, tant au niveau systémique qu’au niveau du système nerveux central (SNC). Il nous a également permis de déterminer si la sensibilité aux infections à S. suis pouvait être influencée par des prédispositions génétiques de l’hôte. Le modèle d’infection par S. suis a été mis au point sur des souris de lignée CD1. Les résultats ont démontré une bactériémie élevée pendant les trois jours suivant l’infection. Celle-ci était accompagnée d’une libération rapide et importante de différentes cytokines pro-inflammatoires (TNF-α, IL-6, IL-12p40/p70, IFN-ɣ) et de chémokines (KC, MCP-1 and RANTES), qui ont entraîné un choc septique et la mort de 20 % des animaux. Ensuite, pour confirmer le rôle de l’inflammation sur la mortalité et pour déterminer si les caractéristiques génétiques de l’hôte pouvaient influencer la réponse inflammatoire et l’issue de la maladie, le modèle d’infection a été étendu à deux lignées murines consanguines différentes considérées comme résistante : la lignée C57BL/6 (B6), et sensible : la lignée A/J. Les résultats ont démontré une importante différence de sensibilité entre les souris A/J et les souris B6, avec un taux de mortalité atteignant 100 % à 20 h post-infection (p.i.) pour la première lignée et de seulement 16 % à 36 h p.i. pour la seconde. La quantité de bactéries dans le sang et dans les organes internes était similaire pour les deux lignées. Donc, tout comme dans la lignée CD1, la bactériémie ne semblait pas être liée à la mort des souris. La différence entre les taux de mortalité a été attribuée à un choc septique non contrôlé chez les souris A/J infectées par S. suis. Les souris A/J présentaient des taux exceptionnellement élevés de TNF-α, IL-12p40/p70, IL-1β and IFN- γ, significativement supérieurs à ceux retrouvés dans la lignée B6. Par contre, les niveaux de chémokines étaient similaires entre les lignées, ce qui suggère que leur influence est limitée dans le développement du choc septique dû à S. suis. Les souris B6 avaient une production plus élevée d’IL-10, une cytokine anti-inflammatoire, ce qui suppose que la cascade cytokinaire pro-inflammatoire était mieux contrôlée, entraînant un meilleur taux de survie. Le rôle bénéfique potentiel de l’IL-10 chez les souris infectées par S. suis a été confirmé par deux approches : d’une part en bloquant chez les souris B6 le récepteur cellulaire à l’IL-10 (IL-10R) par un anticorps monoclonal anti-IL-10R de souris et d’autre part en complémentant les souris A/J avec de l’IL-10 de souris recombinante. Les souris B6 ayant reçu le anticorps monoclonal anti-IL-10R avant d’être infectées par S. suis ont développé des signes cliniques aigus similaires à ceux observés chez les souris A/J, avec une mortalité rapide et élevée et des taux de TNF-α plus élevés que les souris infectées non traitées. Chez les souris A/J infectées par S. suis, le traitement avec l’IL-10 de souris recombinante a significativement retardé l’apparition du choc septique. Ces résultats montrent que la survie au choc septique dû à S. suis implique un contrôle très précis des mécanismes pro- et anti-inflammatoires et que la réponse anti-inflammatoire doit être activée simultanément ou très rapidement après le début de la réponse pro-inflammatoire. Grâce à ces expériences, nous avons donc fait un premier pas dans l’identification de gènes associés à la résistance envers S. suis chez l’hôte. Une des réussites les plus importantes du modèle d’infection de la souris décrit dans ce projet est le fait que les souris CD1 ayant survécu à la septicémie présentaient dès 4 jours p.i. des signes cliniques neurologiques clairs et un syndrome vestibulaire relativement similaires à ceux observés lors de méningite à S. suis chez le porc et chez l’homme. L’analyse par hybridation in situ combinée à de l’immunohistochimie des cerveaux des souris CD1 infectées a montré que la réponse inflammatoire du SNC débutait avec une augmentation significative de la transcription du Toll-like receptor (TLR)2 et du CD14 dans les microvaisseaux cérébraux et dans les plexus choroïdes, ce qui suggère que S. suis pourrait se servir de ces structures comme portes d’entrée vers le cerveau. Aussi, le NF-κB (suivi par le système rapporteur de l’activation transcriptionnelle de IκBα), le TNF-α, l’IL-1β et le MCP-1 ont été activés, principalement dans des cellules identifiées comme de la microglie et dans une moindre mesure comme des astrocytes. Cette activation a également été observée dans différentes structures du cerveau, principalement le cortex cérébral, le corps calleux, l’hippocampe, les plexus choroïdes, le thalamus, l’hypothalamus et les méninges. Partout, cette réaction pro-inflammatoire était accompagnée de zones extensives d’inflammation et de nécrose, de démyélinisation sévère et de la présence d’antigènes de S. suis dans la microglie. Nous avons mené ensuite des études in vitro pour mieux comprendre l’interaction entre S. suis et la microglie. Pour cela, nous avons infecté des cellules microgliales de souris avec la souche sauvage virulente (WT) de S. suis, ainsi qu’avec deux mutants isogéniques, un pour la capsule (CPS) et un autre pour la production d’hémolysine (suilysine). Nos résultats ont montré que la capsule était un important mécanisme de résistance à la phagocytose pour S. suis et qu’elle modulait la réponse inflammatoire, en dissimulant les composants pro-inflammatoires de la paroi bactérienne. Par contre, l’absence d’hémolysine, qui est un facteur cytotoxique potentiel, n’a pas eu d’impact majeur sur l’interaction de S. suis avec la microglie. Ces études sur les cellules microgliales ont permis de confirmer les résultats obtenus précédemment in vivo. La souche WT a induit une régulation à la hausse du TLR2 ainsi que la production de plusieurs médiateurs pro-inflammatoires, dont le TNF-α et le MCP-1. S. suis a induit la translocation du NF-kB. Cet effet était plus rapide dans les cellules stimulées par le mutant déficient en CPS, ce qui suggère que les composants de la paroi cellulaire représentent de puissants inducteurs du NF-kB. De plus, la souche S. suis WT a stimulé l’expression de la phosphotyrosine, de la PKC et de différentes cascades liées à l’enzyme mitogen-activated protein kinase (MAPK). Cependant, les cellules microgliales infectées par le mutant déficient en CPS ont montré des profils de phosphorylation plus forts et plus soutenus que celles infectées par le WT. Finalement, la capsule a aussi modulé l’expression de l’oxyde nitrique synthétase inductible (iNOS) induite par S. suis et par la production subséquente d’oxyde nitrique par la microglie. Ceci pourrait être lié in vivo à la neurotoxicité et à la vasodilatation. Nous pensons que ces résultats contribueront à une meilleure compréhension des mécanismes sous-tendant l’induction de l’inflammation par S. suis, ce qui devrait permettre, d’établir éventuellement des stratégies plus efficaces de lutte contre la septicémie et la méningite. Enfin, nous pensons que ce modèle expérimental d’infection chez la souris pourra être utilisé dans l’étude de la pathogénèse d’autres bactéries ayant le SNC pour cible.

Factor models, VARMA processes and parameter instability with applications in macroeconomics

Avec les avancements de la technologie de l'information, les données temporelles économiques et financières sont de plus en plus disponibles. Par contre, si les techniques standard de l'analyse des séries temporelles sont utilisées, une grande quantité d'information est accompagnée du problème de dimensionnalité. Puisque la majorité des séries d'intérêt sont hautement corrélées, leur dimension peut être réduite en utilisant l'analyse factorielle. Cette technique est de plus en plus populaire en sciences économiques depuis les années 90. Étant donnée la disponibilité des données et des avancements computationnels, plusieurs nouvelles questions se posent. Quels sont les effets et la transmission des chocs structurels dans un environnement riche en données? Est-ce que l'information contenue dans un grand ensemble d'indicateurs économiques peut aider à mieux identifier les chocs de politique monétaire, à l'égard des problèmes rencontrés dans les applications utilisant des modèles standards? Peut-on identifier les chocs financiers et mesurer leurs effets sur l'économie réelle? Peut-on améliorer la méthode factorielle existante et y incorporer une autre technique de réduction de dimension comme l'analyse VARMA? Est-ce que cela produit de meilleures prévisions des grands agrégats macroéconomiques et aide au niveau de l'analyse par fonctions de réponse impulsionnelles? Finalement, est-ce qu'on peut appliquer l'analyse factorielle au niveau des paramètres aléatoires? Par exemple, est-ce qu'il existe seulement un petit nombre de sources de l'instabilité temporelle des coefficients dans les modèles macroéconomiques empiriques? Ma thèse, en utilisant l'analyse factorielle structurelle et la modélisation VARMA, répond à ces questions à travers cinq articles. Les deux premiers chapitres étudient les effets des chocs monétaire et financier dans un environnement riche en données. Le troisième article propose une nouvelle méthode en combinant les modèles à facteurs et VARMA. Cette approche est appliquée dans le quatrième article pour mesurer les effets des chocs de crédit au Canada. La contribution du dernier chapitre est d'imposer la structure à facteurs sur les paramètres variant dans le temps et de montrer qu'il existe un petit nombre de sources de cette instabilité. Le premier article analyse la transmission de la politique monétaire au Canada en utilisant le modèle vectoriel autorégressif augmenté par facteurs (FAVAR). Les études antérieures basées sur les modèles VAR ont trouvé plusieurs anomalies empiriques suite à un choc de la politique monétaire. Nous estimons le modèle FAVAR en utilisant un grand nombre de séries macroéconomiques mensuelles et trimestrielles. Nous trouvons que l'information contenue dans les facteurs est importante pour bien identifier la transmission de la politique monétaire et elle aide à corriger les anomalies empiriques standards. Finalement, le cadre d'analyse FAVAR permet d'obtenir les fonctions de réponse impulsionnelles pour tous les indicateurs dans l'ensemble de données, produisant ainsi l'analyse la plus complète à ce jour des effets de la politique monétaire au Canada. Motivée par la dernière crise économique, la recherche sur le rôle du secteur financier a repris de l'importance. Dans le deuxième article nous examinons les effets et la propagation des chocs de crédit sur l'économie réelle en utilisant un grand ensemble d'indicateurs économiques et financiers dans le cadre d'un modèle à facteurs structurel. Nous trouvons qu'un choc de crédit augmente immédiatement les diffusions de crédit (credit spreads), diminue la valeur des bons de Trésor et cause une récession. Ces chocs ont un effet important sur des mesures d'activité réelle, indices de prix, indicateurs avancés et financiers. Contrairement aux autres études, notre procédure d'identification du choc structurel ne requiert pas de restrictions temporelles entre facteurs financiers et macroéconomiques. De plus, elle donne une interprétation des facteurs sans restreindre l'estimation de ceux-ci. Dans le troisième article nous étudions la relation entre les représentations VARMA et factorielle des processus vectoriels stochastiques, et proposons une nouvelle classe de modèles VARMA augmentés par facteurs (FAVARMA). Notre point de départ est de constater qu'en général les séries multivariées et facteurs associés ne peuvent simultanément suivre un processus VAR d'ordre fini. Nous montrons que le processus dynamique des facteurs, extraits comme combinaison linéaire des variables observées, est en général un VARMA et non pas un VAR comme c'est supposé ailleurs dans la littérature. Deuxièmement, nous montrons que même si les facteurs suivent un VAR d'ordre fini, cela implique une représentation VARMA pour les séries observées. Alors, nous proposons le cadre d'analyse FAVARMA combinant ces deux méthodes de réduction du nombre de paramètres. Le modèle est appliqué dans deux exercices de prévision en utilisant des données américaines et canadiennes de Boivin, Giannoni et Stevanovic (2010, 2009) respectivement. Les résultats montrent que la partie VARMA aide à mieux prévoir les importants agrégats macroéconomiques relativement aux modèles standards. Finalement, nous estimons les effets de choc monétaire en utilisant les données et le schéma d'identification de Bernanke, Boivin et Eliasz (2005). Notre modèle FAVARMA(2,1) avec six facteurs donne les résultats cohérents et précis des effets et de la transmission monétaire aux États-Unis. Contrairement au modèle FAVAR employé dans l'étude ultérieure où 510 coefficients VAR devaient être estimés, nous produisons les résultats semblables avec seulement 84 paramètres du processus dynamique des facteurs. L'objectif du quatrième article est d'identifier et mesurer les effets des chocs de crédit au Canada dans un environnement riche en données et en utilisant le modèle FAVARMA structurel. Dans le cadre théorique de l'accélérateur financier développé par Bernanke, Gertler et Gilchrist (1999), nous approximons la prime de financement extérieur par les credit spreads. D'un côté, nous trouvons qu'une augmentation non-anticipée de la prime de financement extérieur aux États-Unis génère une récession significative et persistante au Canada, accompagnée d'une hausse immédiate des credit spreads et taux d'intérêt canadiens. La composante commune semble capturer les dimensions importantes des fluctuations cycliques de l'économie canadienne. L'analyse par décomposition de la variance révèle que ce choc de crédit a un effet important sur différents secteurs d'activité réelle, indices de prix, indicateurs avancés et credit spreads. De l'autre côté, une hausse inattendue de la prime canadienne de financement extérieur ne cause pas d'effet significatif au Canada. Nous montrons que les effets des chocs de crédit au Canada sont essentiellement causés par les conditions globales, approximées ici par le marché américain. Finalement, étant donnée la procédure d'identification des chocs structurels, nous trouvons des facteurs interprétables économiquement. Le comportement des agents et de l'environnement économiques peut varier à travers le temps (ex. changements de stratégies de la politique monétaire, volatilité de chocs) induisant de l'instabilité des paramètres dans les modèles en forme réduite. Les modèles à paramètres variant dans le temps (TVP) standards supposent traditionnellement les processus stochastiques indépendants pour tous les TVPs. Dans cet article nous montrons que le nombre de sources de variabilité temporelle des coefficients est probablement très petit, et nous produisons la première évidence empirique connue dans les modèles macroéconomiques empiriques. L'approche Factor-TVP, proposée dans Stevanovic (2010), est appliquée dans le cadre d'un modèle VAR standard avec coefficients aléatoires (TVP-VAR). Nous trouvons qu'un seul facteur explique la majorité de la variabilité des coefficients VAR, tandis que les paramètres de la volatilité des chocs varient d'une façon indépendante. Le facteur commun est positivement corrélé avec le taux de chômage. La même analyse est faite avec les données incluant la récente crise financière. La procédure suggère maintenant deux facteurs et le comportement des coefficients présente un changement important depuis 2007. Finalement, la méthode est appliquée à un modèle TVP-FAVAR. Nous trouvons que seulement 5 facteurs dynamiques gouvernent l'instabilité temporelle dans presque 700 coefficients.

Le niveau réduit d’AMPc dans la surexpression des protéines G(alpha)i et la prolifération accrue des cellules du muscle lisse vasculaire des rats spontanément hypertendus

Nous avons précédemment montré que les cellules musculaires lisses vasculaires(CMLV) des rats spontanément hypertendus (SHR) présentent une expression augmentée des protéines G inhibitrices (Gi) et une prolifération cellulaire accrue par rapport aux CMLV des rats Wystar-Kyoto (WKY). Le niveau d'AMPc s’est également avéré plus faible dans les CMLV de SHR. La présente étude a donc été entreprise afin d'examiner la contribution de la diminution du niveau intracellulaire d'AMPc à l’augmentation de l'expression des protéines Gi et à la prolifération accrue des CMLV de SHR et de continuer à explorer les mécanismes moléculaires sous-jacents responsables de cette réponse. Les CMLV de SHR ont montré par rapport aux CMLV des WKY une expression accrue de Giα-2 et Giα-3 qui a été diminué d'une manière dépendante de concentration par le dbcAMP, un analogue d'AMPc perméable à la membrane cellulaire. En outre, les fonctions augmentées des protéines Gi comme démontrées par l'amplification de l’inhibition de l'adénylate cyclase par les hormones inhibitrices et l'activité forskoline (FSK)-stimulée de l’adénylate cyclase par une faible concentration de GTPγS dans les CMLV de SHR ont également été restaurées aux niveaux de WKY par le dbcAMP. La prolifération accrue des CMLV de SHR a également été atténuée par le dbcAMP et la forskoline, un activateur de l'adénylate cyclase. De plus, dbcAMP a restauré la production augmentée d'anion superoxyde (O2-), l'activité de la NAD(P)H oxydase et l’expression accrue des protéines Nox 4 et p47phox observée dans les CMLV de SHR jusqu’au niveau contrôle. Par ailleurs, la phosphorylation accrue des PDGF-R, EGF-R, c-Src et ERK1/2 énoncée par les CMLV de SHR a également été diminuée par le dbcAMP d'une manière dépendante de concentration. Ces résultats suggèrent que le niveau réduit d'AMPc intracellulaire montré par les CMLV de SHR contribue à l'expression accrue des protéines Gi et à l’hyperprolifération cellulaire à travers l’augmentation du stress oxydatif, la transactivation des EGF-R, PDGF-R et la voie de signalisation des MAP kinases.

The role of Microsomal prostaglandin synthase-1 (mPGES-1) and Ephrin B2 in Scleroderma

La sclérodermie (sclérose systémique, ScS) est une maladie auto-immune du tissu conjonctif caractérisée  par  l’épaississement  de  la  peau,  l’apparition  spontanée de lésions cicatricielles, des maladies   des   vaisseaux   sanguins,   divers   degrés   d’inflammation,   en   association   avec   un   système   immunitaire hyperactif. La pathogénèse exacte de cette maladie est inconnue et aucun traitement approprié   n’est   disponible.   La fibrose est un élément distinctif de la maladie de ScS et est considérée   résulter   d’une   incapacité   à   mettre   fin   de   façon   appropriée   à   la   réponse   normale   de   réparation   des   plaies.   L’analyse   histologique   du   stade   initial   de   la   ScS   révèle   une   infiltration   périvasculaire de cellules mononucléaires dans le derme, associée à une synthèse accrue de collagène dans les fibroblastes environnants. Ainsi, la compréhension des moyens de contrôler le stade inflammatoire de la ScS pourrait être bénéfique pour contrôler la progression de la maladie peu après son apparition. La mPGES-1 est une enzyme inductible qui agit en aval de la cyclo- oxygénase (COX) pour catalyser spécifiquement la conversion de la prostaglandine (PG) H2 en PGE2. La mPGES-1  joue  un  rôle  clé  dans  l’inflammation,  la  douleur  et  l’arthrite;;  toutefois,  le   rôle de la mPGES-1 dans les mécanismes de fibrose, spécifiquement en rapport avec la ScS humaine, est inconnu. Mon laboratoire a précédemment montré que les souris à mPGES-1 nulle sont résistantes à la fibrose   cutanée   induite   par   la   bléomycine,   à   l’inflammation,   à   l’épaississement  cutané,  à  la  production  de  collagène  et  à  la  formation  de  myofibroblastes.  Sur  la   base  de  ces  résultats,  j’ai  formulé  l’hypothèse  que  l’inhibition pharmacologique de la mPGES-1 régulera à la baisse la production de médiateurs pro-inflammatoires et pro-fibreux au cours de la maladie   de   ScS.   Afin   d’explorer   le   rôle   de   la   mPGES-1   dans   l’inflammation   et   la   fibrose   associées  à  la  maladie  de  ScS,  j’ai  d’abord  examiné  l’expression  de  la  mPGES-1 dans la peau normale comparativement à des biopsies de peau extraites de patients atteints de ScS. Mes résultats ont montré que la mPGES-1 est nettement élevée dans la peau de patients atteints de ScS en comparaison avec la peau humaine normale. De plus, les niveaux de PGE2 dérivés de la mPGES-1 étaient également significativement plus élevés dans les fibroblastes cutanés isolés de patients  atteints  de  ScS  comparativement  aux  fibroblastes  isolés  de  témoins  sains.  J’ai  également   étudié  l’effet  de  l’inhibition pharmacologique de la mPGES-1  sur  l’expression  de  marqueurs  pro- fibreux.   Mes   études   ont   montré   que   l’expression   de   médiateurs   pro-fibreux clés (α-SMA, endothéline-1, collagène de type 1 et facteur de croissance du tissu conjonctif (FCTC)) est élevée dans les fibroblastes cutanés ScS en comparaison avec les fibroblastes cutanés normaux. Un traitement avec un inhibiteur de la mPGES-1 a eu pour effet de réduire significativement l’expression  de  l’α-SMA,  de  l’endothéline-1, du collagène de type 1 mais pas du FCTC dans les fibroblastes  ScS,  sans  effet  significatif  sur  les  fibroblastes  normaux.  J’ai  en  outre  examiné  l’effet   de  l’inhibition  de  la  mPGES-1 sur des cytokines pro-inflammatoires clés impliquées dans la pathologie de la ScS, incluant IL-6, IL-8 et MCP-1.  L’inhibition  pharmacologique  de  la  mPGES- 1 a eu pour effet de réduire significativement les niveaux de production de cytokines pro- inflammatoires IL6, IL8 et MCP-1 dans les fibroblastes avec lésion ScS comparativement à des fibroblastes non traités. De plus, les patients atteints de ScS ont présenté des niveaux plus élevés de p-AKT, de p-FAK et de p-SMAD3 en comparaison avec les fibroblastes cutanés normaux. L’inhibiteur  de  la  mPGES-1 a pu réguler à la baisse cette expression accrue de p-AKT et de p- FAK, mais pas de p-SMAD3,  dans  les  fibroblastes  ScS.  Ces  résultats  ont  suggéré  que  l’inhibition   de la mPGES-1 pourrait être une méthode viable pour réduire le développement de sclérose cutanée et constituent une cible thérapeutique potentielle pour contrôler les mécanismes fibreux et inflammatoires associés à la pathophysiologie de la maladie de ScS. L’un   des   autres   processus   critiques   reliés   à   l’évolution de la réponse fibreuse associée à la maladie de ScS est la différenciation des fibroblastes en des cellules activées spécialisées iii iv appelées myofibroblastes, responsables de déclencher une signalisation adhésive excessive et le dépôt excessif de matrice extracellulaire,   conduisant   à   la   destruction   de   l’architecture   de   l’organe.   Ainsi,   l’identification   des   facteurs   endogènes   qui   initient/   favorisent   la   différenciation   fibroblaste-myofibroblaste peut mener à des stratégies thérapeutiques prometteuses pour contrôler  l’excès  de  signalisation  adhésive  et  de  fibrose  associé  à  la  maladie  de  ScS.  Des  études   antérieures  dans  le  domaine  de  la  biologie  du  cancer  ont  suggéré  que  l’éphrine  B2,  une  protéine   transmembranaire appartenant à la famille des éphrines, est impliquée dans la signalisation adhésive   et   le   remodelage   extracellulaire.   Cependant,   son   rôle   dans   la   fibrose   n’a   jamais   été   exploré.   Dans   la   deuxième   partie   de   mon   étude,   j’ai   donc   étudié   le   rôle   de   l’éphrine   B2   dans   la   fibrose.   Mes   études   montrent   que   l’expression   de   l’éphrine   B2   est   significativement   augmentée   dans la peau humaine ScS comparativement à la peau normale. Plus important encore, le traitement in vitro de   fibroblastes   de   la   peau   humaine   normale   avec   de   l’éphrine   B2   recombinante est capable de transformer des fibroblastes en cellules myofibroblastiques manifestant toutes les caractéristiques myofibroblastiques typiques, incluant la formation accrue de  fibres  de  tension,  des  adhérences  focales,  l’activation  accrue  de  la  FAK,  un  accroissement  de   l’expression  et  de  la  migration  de  fibroblastes  et  de  leur  adhérence  à  la  fibronectine  à  la  fois  chez   les   fibroblastes   cutanés   normaux   et   ScS.   En   outre,   j’ai   traité   des   souris   avec   de   l’éphrine   B2   recombinante et montré que ces souris ont développé une fibrose cutanée significative associée à une épaisseur dermique et à une synthèse de collagène augmentées, une teneur en hydroxyproline (teneur en collagène) accrue et un nombre accru de myofibroblastes exprimant de   l’α-SMA, une activation augmentée de la FAK et de marqueurs pro-fibreux incluant le collagène de type 1 et le FCTC. Dans  l’ensemble,  mes  études  ont  identifié  deux  médiateurs  endogènes  cruciaux  impliqués  dans  la   propagation  de  l’inflammation  et  de  la  fibrose  associées  à  la  maladie  de  ScS.  L’inhibition  de  la   mPGES-1   pourrait   représenter   une   bonne   stratégie   alternative   pour   contrer   l’inflammation   et   la   fibrose au moins durant les stades précoces de la maladie de ScS. De plus, une signalisation excessive   de   l’éphrine B2 favorise la signalisation adhésive et fibreuse en déclenchant la différenciation   de   fibroblastes   en   myofibroblastes   par   l’activation   de   la   voie   de   signalisation   de   la  FAK.  Ainsi,  l’inhibition  d’éphrine  B2  bloquera  la  formation  de  fibroblastes-myofibroblastes et régulera à la baisse la fibrose associée à la maladie de ScS. En somme, la mPGES-1  et  l’éphrine   B2 semblent toutes deux des cibles attrayantes pour le traitement de la ScS et des troubles fibreux qui y sont reliés.

L’immunoprotéasome : producteur de peptides-CMH I et régulateur de l’expression génique

Le système ubiquitine-protéasome est le principal mécanisme par lequel les protéines intracellulaires sont dégradées. Le protéasome dit constitutif (PC) est donc essentiel à l’homéostasie mais aussi à la régulation de la majorité des processus cellulaires importants. La découverte d’un deuxième type de protéasome, appelé immunoprotéasome (IP), soulève toutefois de nouvelles questions. Pourquoi existe-t-il plus d’un type de protéasome ? L’IP a-t-il des rôles redondants ou complémentaires avec le PC ? L’IP étant présent principalement dans les cellules immunitaires ou stimulées par des cytokines, plusieurs groupes ont tenté de définir son rôle dans la réponse immunitaire. Or, l’implication de son homologue constitutif dans un éventail de processus non spécifiquement immunitaires nous laisse croire que l’IP pourrait lui aussi avoir un impact beaucoup plus large. L’objectif de cette thèse était donc de caractériser certains rôles cellulaires de l’IP dans les cellules dendritiques. Nous avons d’abord étudié l’impact global de l’IP sur la présentation antigénique de classe I. Ce faisant, nous avons pu déterminer ses deux contributions principales, soit l’augmentation drastique du nombre et de la diversité des peptides présentés sur les complexes majeurs d’histocompatibilité de classe I. Les différences de clivage entre le PC et l’IP pourraient expliquer en partie cette diversité du répertoire peptidique, notamment par l’affinité apparente de l’IP pour les régions protéiques non structurées. Dans un deuxième temps, nous avons dévoilé un nouveau rôle de l’IP sur un processus dépassant le cadre immunitaire : la transcription. Nous avons découvert que l’IP modifie l’abondance des ARNm en agissant principalement au niveau de leur synthèse. L’impact de l’IP sur le transcriptome est majeur et serait dû en partie à une dégradation différente de facteurs de transcription des familles IRF, STAT et NF-kB. Les cellules dendritiques IP-déficientes activent moins efficacement les lymphocytes T CD8+ et nous croyons que cette défaillance est causée (du moins en partie) par la perturbation transcriptomique provoquée par l’absence d’IP. Il importe donc de comprendre les différents rôles moléculaires de l’IP afin de mieux définir sa contribution globale au fonctionnement de la cellule et comprendre l’avantage évolutif, au niveau de l’organisme, procuré par une telle plasticité du système ubiquitine-protéasome.

The specific in vivo role of PPARgamma and its downstream signaling pathway in the pathophysiology of Osteoarthritis

L'arthrose est une maladie articulaire dégénérative, avec une pathogenèse inconnue. Des études récentes suggèrent que l'activation du facteur de transcription du récepteur activateur de la prolifération des peroxysomes (PPAR) gamma est une cible thérapeutique pour ce maladie. Les agonistes du PPARγ inhibent l'inflammation et réduisent la synthèse des produits de dégradation du cartilage in vitro et in vivo. Cependant, des études utilisant des agonistes du PPARγ n’élucident pas les effets exacts médiés par ce gène complexe. En effet, certains de ces agonistes ont la capacité de régulariser d'autres voies de signalisation indépendantes de PPARγ, ainsi entraînant des effets secondaires graves. Afin d'obtenir une efficacité thérapeutique avec potentiellement moins de problèmes de sécurité, il est donc essentiel d'élucider, in vivo, le rôle exact de PPARγ dans la physiopathologie OA. Mon projet de thèse permettra de déterminer, pour la première fois, le rôle spécifique de PPARγ in vivo dans la physiopathologie OA. Les souris utilisées pour l’étude avaient une délétion conditionnelle du gène PPARγ dans le cartilage. Ces dernières ont été générées en employant le système LoxP/Cre. Pour tester cette hypothèse, j'ai généré deux types de souris avec une délétion au PPARγ, (a) une suppression du gène PPARγ spécifiquement dans le cartilage germinale pour l'étude de l'arthrose liée au développement et à l'âge et (b) la suppression inductible du gène PPARγ spécifiquement dans le cartilage chez la souris adulte pour les études OA. L’étude précédente dans notre laboratoire, utilisant ces souris ayant une délétion au gène PPARγ germinales, montre que ces souris présentent des anomalies du développement du cartilage. J'ai également exploré si ces souris qui présentent des défauts précoces du développement ont toutes les modifications phénotypiques dans le cartilage au cours du vieillissement. Mes résultats ont montré que les souris adultes, ayant une délétion au gène PPARγ, ont présenter un phénotype de l'arthrose spontanée associée à une dégradation du cartilage, l’hypocellularité, la fibrose synoviale. Cette étude a montré que PPARγ est un régulateur essentiel pour le cartilage, et c’est le manque (l’absence) de ce dernier qui conduit à un phénotype de l'arthrose spontanée accélérée (American Journal of Pathologie). A partir de ce but de l'étude, on n’a pas pu vérifier si ces souris présentaient l’OA spontanée en raison des défauts de développement ou à la suite de la délétion du gène PPARγ. Pour contourner les défauts de développement, j'ai généré des souris ayant une délétion du gène PPARγ spécifiquement dans le cartilage inductible avec le système Col2rTACre. Ces souris ont été soumises à modèle de la chirurgie OA (DMM: déstabilisation du ménisque médial) et les résultats révèlent que les souris PPARγ KO ont une dégradation accélérée du cartilage, une hypocellularité, une fibrose synoviale et une augmentation de l'expression des marqueurs cataboliques et des marqueurs inflammatoire. La perte de PPAR dans le cartilage articulaire est un évènement critique qui initie la dégradation de cartilage dans OA. Les études récentes suggèrent que le procès d’autophagie, une forme de survie cellulaire programmée, est altéré pendant l’OA et peut contribuer vers une protection diminuée des cellules, résultant la dégradation du cartilage. J’ai donc exploré le rôle de PPARγ dans la protection des cellules en déterminant l’effet de manque de PPARγ dans le cartilage par l’expression de mTOR (régulateur négatif principal d’autophagie) et les gènes d’autophagie durant OA. Mes résultats ont montré que les souris KO PPARγ présentent également une augmentation sur l'expression de mTOR et une diminution sur l’expression des marqueurs autophagiques en comparaison avec les chondrocytes articulaires isolés des souris contrôles OA. J'ai suggéré l'hypothèse que PPARγ contrôle la régulation de la signalisation de mTOR/autophagie, et finalement la mort des chondrocytes et l’expression des facteurs cataboliques et les facteurs inflammatoire. Pour tester cette hypothèse, j’ai fait la transfection des chondrocytes arthrosiques PPARγ-KO avec le vecteur d’expression de PPARγ pour déterminer si la restauration de l'expression de PPARγ peut sauver le phénotype des cellules PPARγ-KO OA. J'ai observé que la restauration de l'expression de PPARγ dans les cellules PPARγ-KO en présence du vecteur d'expression PPARγ, a pu considérablement régulariser négativement l'expression de mTOR et mettre en règle positivement l'expression des gènes autophagiques ainsi que le sauvetage significative de l'expression du collagène de type II et l’aggrecan et de baisser de manière significative l'expression de marqueurs cataboliques critiques et des marqueurs inflammatoires. Pour prouver que l’augmentation de la signalisation de mTOR et la diminution de l'autophagie est responsable du phénotype OA accélérée observée dans les souris PPARγ KO in vivo, j'ai généré les souris doubles KO PPARγ- mTOR inductible spécifique du cartilage en utilisant le système Col2 - rtTA -Cre et soumis ces souris à DMM modèle de l'arthrose. Mes résultants démontrent que les souris avec PPARγ- mTOR doubles KO ont été significativement protégés contre les OA DMM induites associées à une protection significative contre la destruction du cartilage, la perte de protéoglycanes et la perte de chondro-cellularité par rapport aux souris témoins. Considérant que mTOR est un répresseur majeur de l'autophagie, j'ai trouvé que l'expression de deux marqueurs de l'autophagie critiques (ULK1 et LC3B) a été significativement plus élevée dans les chondrocytes extraits les souris doubles KO PPARγ-mTOR par rapport aux souris témoins. En plus, les études de sauvetage in vitro en utilisant le vecteur d'expression PPAR et les études in vivo utilisant les souris doubles KO PPARγ- mTOR montrent que PPARγ est impliqué dans la régulation de la protéine signalant de mTOR/autophagie dans le cartilage articulaire. Ces résultats contournent PPARγ et sa signalisation en aval de mTOR/autophagie en tant que cibles thérapeutiques potentielles pour le traitement de l'arthrose.

L’impact de l’inhibition de la iNOS et de la kininase I dans les effets délétères du récepteur B1 des kinines dans le diabète de type 2

Les kinines sont des peptides vasoactifs et des neuromédiateurs centraux impliqués dans un bon nombre de processus biologiques et inflammatoires. Elles agissent sur deux types de récepteurs (R) couplés aux protéines G, le RB2 constitutif et le RB1 qui est induit par le stress oxydatif et les cytokines pro-inflammatoires via le facteur de transcription nucléaire, le NF-kB. Le RB1 est un puissant activateur de la iNOS et il augmente son expression chez le rat insulino-résistant. Dans ce modèle de rats soumis à une diète riche en D-glucose, un traitement d’une semaine avec un antagoniste non peptidique du RB1, le SSR240612, renverse la plupart des complications diabétiques. Ces travaux nous mènent à émettre l’hypothèse que la iNOS contribue aux effets délétères du RB1 chez le rat insulino-résistant. Nous avons donc évalué les effets d’un traitement prolongé d’une semaine soit avec le 1400W (1 mg/kg x 2 fois/jour), un inhibiteur sélectif de la iNOS, ou avec le Mergetpa (1mg/kg x 2 fois/jour), un inhibiteur non sélectif de la carboxypeptidase M (CPM) qui supprime la formation d’agonistes du RB1. Ces deux traitements devraient reproduire les effets bénéfiques de l’antagoniste du RB1 (SSR240612). En effet, le 1400W et le Mergetpa corrigent l’hyperglycémie, la résistance à l’insuline (l’indice HOMA), l’allodynie au froid et l’expression de plusieurs marqueurs de l’inflammation (Cox-2, iNOS, RB1, IL-1, anion superoxyde). Ces résultats confirment la contribution de la iNOS dans les effets délétères du RB1 chez le rat insulino-résistant. Dans un autre volet, ce mémoire vise à mieux comprendre l’impact de l’inhibition du RB1 par le SSR240612 (10 μg/g/jour) combiné ou pas avec le Pioglitazone (1.6 mg/g/jour) (un anti-diabétique de la famille des thiazolidinediones, le TZD) dans un modèle de diabète de type 2 associé à l’obésité chez la souris C57BL/6J soumise à une diète riche vi en gras pendant vingt semaines. Un traitement pendant deux semaines avec le TZD corrige l’intolérance au glucose et réduit les taux plasmatiques d’insuline alors que le SSR n’a pas d’effet. Les traitements combinés du TZD avec le SSR corrigent davantage la perte de la sensibilité à l’insuline et réduisent les taux plasmatiques de leptine. Les résultats obtenus suggèrent que le SSR n’apporte pas l’effet bénéfique souhaité, dans ce modèle avancé de diabète de type 2, contrairement au modèle des rats insulino-résistants (pré-diabétiques).

The role of transcription factor Nrf2 in osteoarthritis

Introduction: L'arthrose est caractérisée par une destruction progressive du cartilage, une inflammation synoviale, et un remodelage de l’os sous-chondral avec une production excessive des médiateurs inflammatoires et cataboliques. Nous avons démontré que le niveau du 4-hydroxynonénal (4-HNE), un produit de la peroxydation lipidique, est augmenté dans le cartilage humain arthrosique sans qu’on sache le mécanisme exacte impliqué dans l’augmentation de cette molécule. Des données de la littérature indiquent que l’accumulation du HNE est contrôlée par l’action de la glutathione S-transférase A4-4 (GSTA4-4), une enzyme impliquée dans la détoxification du HNE. Au niveau transcriptionel, l’expression de cette enzyme est régulée par la transactivation du facteur de transcription Nrf2. Objectif: L’objectif de cette étude vise à démontrer que l’augmentation du HNE dans le cartilage arthrosique est attribuée, en partie, à l’altération de l’expression de la GSTA4-4 et de Nrf2. Méthode: Le niveau d’expression de la GSTA4-4 et de Nrf2 a été mesurée par Western blot et par PCR en temps réel dans le cartilage humain arthrosique et dans le cartilage provenant des souris atteintes d’arthrose. Pour démontrer le rôle du Nrf2 dans l’arthrose, les chondrocytes humains arthrosiques ont été traités par l’interleukine 1beta (IL-1β) ou par le H2O2 en présence ou en absence des activateurs du Nrf2 tels que le Protandim®, AI, et du 6-Gingérol. Par ailleurs, les chondrocytes ont été transfectés par un vecteur d’expression de Nrf2 puis traités par l’IL-β. En utilisant le modèle d’arthrose chez la souris, les animaux ont été traités par voie orale de 10 mg/kg/jour de Protandim® pendant 8 semaines. Résultats: Nous avons observé une diminution significative de l’expression de la GSTA4-4 et de Nrf2 dans le cartilage humain et murin arthrosique. L'activation de Nrf2 bloque la stimulation de la métalloprotéinase-13 (MMP-13), la prostaglandine E2 (PGE2) et de l'oxyde nitrique (NO) par l’IL-1β. En outre, nous avons montré que l'activation Nrf2 protège les cellules contre la mort cellulaire induite par H2O2. Fait intéressant, l'administration orale de Protandim® réduit la production du HNE par l'intermédiaire de l’activation de la GSTA4. Nous avons démontré que le niveau d’expression de la GSTA4-4 et de Nrf2 diminue dans le cartilage provenant des patients et des souris atteints d’arthrose. De plus, la surexpression de ce facteur nucléaire Nrf2 empêche la production du HNE et la MMP-13 et l’inactivation de la GSTA4-4. Dans notre modèle expérimental d’arthrose induite par déstabilisation du ménisque médial chez la souris, nous avons trouvé que l'administration orale de Protandim® à 10 mg / kg / jour réduit les lésions du cartilage. Conclusion: Cette étude est de la première pour démontrer le rôle physiopathologique du Nrf2 in vitro et in vivo. Nos résultats démontrent que l’activation du Nrf2 est essentielle afin de maintenir l’expression de la GSTA4-4 et de réduire le niveau du HNE. Le fait que les activateurs du Nrf2 abolissent la production de la HNE et aussi un certain nombre de facteurs connus pour être impliqués dans la pathogenèse de l’arthrose les rend des agents cliniquement utiles pour la prévention de la maladie.

Dose/frequency: A critical factor in the administration of glucan as immunostimulant to Indian white shrimp Fenneropenaeus indicus

The immunostimulatory effect of an alkali insoluble glucan extracted from marine yeast isolate Candida sake S165 was tested in Fenneropenaeus indicus. Post larvae (PL) of F. indicus, fed glucan incorporated diet at varying concentrations (0.05, 0.1, 0.2, 0.3, 0.4 g glucan/100 g feed) for 21 days were challenged orally with white spot syndrome virus (WSSV). Maximum survival was observed in PL fed the 0.2% glucan incorporated diet. Subsequently the feed incorporated with 0.2% glucan was fed to F. indicus post larvae at different feeding intervals, i.e. daily, once every two days, once every five days, once every seven days and once every ten days. After 40 days, the prawns were challenged orally with WSSV and post challenge survival was recorded. Shrimp feed containing 0.2% glucan when administered once every seven days gave maximum survival. This was supported by haematological data obtained from adult F. indicus, i.e. total haemocyte count, phenoloxidase activity and nitroblue tetrazolium reduction (NBT). The present observation confirms the importance of dose and frequency of administration of immunostimulants in shrimp health management

Anti-lipopolysaccharide factor and crustin-III, the anti-white spot virus peptides in Penaeus monodon: Control of viral infection by up-regulation

White Spot Syndrome Virus (WSSV) is the most devastating disease affecting shrimp culture around the world. Though, considerable progress has been made in the detection and molecular characterization of WSSV in recent years, information pertaining to immune gene expression in shrimps with respect to WSSV infection remains limited. In this context, the present study was undertaken to understand the differential expression of antimicrobial peptide (AMP) genes in the haemocytes of Penaeus monodon in response to WSSV infection on a time-course basis employing semi-quantitative RT-PCR. The present work analyzes the expression profile of six AMP genes (ALF, crustin-1, crustin-2, crustin-3, penaeidin-3 and penaeidin-5), eight WSSV genes (DNA polymerase, endonuclease, immediate early gene, latency related gene, protein kinase, ribonucleotide reductase, thymidine kinase and VP28) and three control genes (18S rRNA, β-actin and ELF) in P. monodon in response to WSSV challenge. Penaeidins were found to be up-regulated during early hours of infection and crustin-3 during late period of infection. However, ALF was found to be up-regulated early to late period of WSSV infection. The present study suggests that AMPs viz. ALF and crustin-3 play an important role in antiviral defense in shrimps. WSSV gene transcripts were detected post-challenge day 1 itself and increased considerably day 5 onwards. Evaluation of the control genes confirmed ELF as the most reliable control gene followed by 18S rRNA and β-actin for gene expression studies in shrimps. This study indicated the role of AMPs in the protection of shrimps against viral infection and their possible control through the up-regulation of AMPs

A new isoform of anti-lipopolysaccharide factor identified from the blue swimmer crab, Portunus pelagicus: Molecular characteristics and phylogeny

Anti-lipopolysaccharide factors are small proteins that bind and neutralize lipopolysaccharide and exhibit potent antimicrobial activities. This study presents the molecular characterization and phylogenetic analysis of the first ALF isoform (Pp-ALF1; JQ745295) identified from the hemocytes of Portunus pelagicus. The full length cDNA of Pp-ALF1 consisted of 880 base pairs encoding 293 amino acids with an ORF of 123 amino acids and contains a putative signal peptide of 24 amino acids. Pp-ALF1 possessed a predicted molecular weight (MW) of 13.86 kDa and theoretical isoelectric point (pI) of 8.49. Two highly conserved cysteine residues and putative LPS binding domain were observed in Pp-ALF1. Peptide model of Pp-ALF1 consisted of two α-helices crowded against a four-strand β-sheet. Comparison of amino acid sequences and neighbor joining tree showed that Pp-ALF1 has a maximum similarity (46%) to ALF present in Portunus trituberculatus followed by 39% similarity to ALF of Eriocheir sinensis and 38% similarity to ALFs of Scylla paramamosain and Scylla serrata. Pp-ALF1 is found to be a new isoform of ALF family and its characteristic similarity with other known ALFs signifies its role in protection against invading pathogens.

Length-weight relationship, relative condition factor (Kn) and morphometry of Arius subrostratus (Valenciennes, 1840) from a coastal wetland in Kerala

The length – weight relationship and relative condition factor of the shovel nose catfish, Arius subrostratus (Valenciennes, 1840) from Champakkara backwater were studied by examination of 392 specimens collected during June to September 2008. These fishes ranged from 6 to 29 cm in total length and 5.6 to 218 g in weight. The relation between the total length and weight of Arius subrostratus is described as Log W = -1.530+2.6224 log L for males, Log W = - 2.131 + 3.0914 log L for females and Log W = - 1.742 + 2.8067 log L for sexes combined. The mean relative condition factor (Kn) values ranged from 0.75 to 1.07 for males, 0.944 to 1.407 for females and 0.96 to 1.196 for combined sexes. The length-weight relationship and relative condition factor showed that the well-being of A. subrostratus is good. The morphometric measurements of various body parts and meristic counts were recorded. The morphometric measurements were found to be non-linear and there is no significant difference observed between the two sexes. From the present investigation, the fin formula can be written as D: I, 7; P: I, 12; A: 17 – 20; C: 26 – 32. There is no change in meristic counts with the increase in body length.

Volume isotope shifts in low lying transitions of Au I

Electronic factors in the volume isotope shift have been calculated in an ab initio way with the relativistic Dirac-Fock method for a number of different optical single/and two-photon transitions in Au I. The agreement with a semi-empirical method is within 10% for the resonance transition. For this one and a few other transitions the effect of core excitation has been analyzed with the Multi-configuration Dirac-Fock method as well, and it was found to reduce the electronic factor in the order of 5 %.

A Factor analysis of hidrochemical composition of Llobregat river basin

Hydrogeological research usually includes some statistical studies devised to elucidate mean background state, characterise relationships among different hydrochemical parameters, and show the influence of human activities. These goals are achieved either by means of a statistical approach or by mixing models between end-members. Compositional data analysis has proved to be effective with the first approach, but there is no commonly accepted solution to the end-member problem in a compositional framework. We present here a possible solution based on factor analysis of compositions illustrated with a case study. We find two factors on the compositional bi-plot fitting two non-centered orthogonal axes to the most representative variables. Each one of these axes defines a subcomposition, grouping those variables that lay nearest to it. With each subcomposition a log-contrast is computed and rewritten as an equilibrium equation. These two factors can be interpreted as the isometric log-ratio coordinates (ilr) of three hidden components, that can be plotted in a ternary diagram. These hidden components might be interpreted as end-members. We have analysed 14 molarities in 31 sampling stations all along the Llobregat River and its tributaries, with a monthly measure during two years. We have obtained a bi-plot with a 57% of explained total variance, from which we have extracted two factors: factor G, reflecting geological background enhanced by potash mining; and factor A, essentially controlled by urban and/or farming wastewater. Graphical representation of these two factors allows us to identify three extreme samples, corresponding to pristine waters, potash mining influence and urban sewage influence. To confirm this, we have available analysis of diffused and widespread point sources identified in the area: springs, potash mining lixiviates, sewage, and fertilisers. Each one of these sources shows a clear link with one of the extreme samples, except fertilisers due to the heterogeneity of their composition. This approach is a useful tool to distinguish end-members, and characterise them, an issue generally difficult to solve. It is worth note that the end-member composition cannot be fully estimated but only characterised through log-ratio relationships among components. Moreover, the influence of each endmember in a given sample must be evaluated in relative terms of the other samples. These limitations are intrinsic to the relative nature of compositional data

Robust Factor Analysis for Compositional Data

Factor analysis as frequent technique for multivariate data inspection is widely used also for compositional data analysis. The usual way is to use a centered logratio (clr) transformation to obtain the random vector y of dimension D. The factor model is then y = Λf + e (1) with the factors f of dimension k < D, the error term e, and the loadings matrix Λ. Using the usual model assumptions (see, e.g., Basilevsky, 1994), the factor analysis model (1) can be written as Cov(y) = ΛΛT + ψ (2) where ψ = Cov(e) has a diagonal form. The diagonal elements of ψ as well as the loadings matrix Λ are estimated from an estimation of Cov(y). Given observed clr transformed data Y as realizations of the random vector y. Outliers or deviations from the idealized model assumptions of factor analysis can severely effect the parameter estimation. As a way out, robust estimation of the covariance matrix of Y will lead to robust estimates of Λ and ψ in (2), see Pison et al. (2003). Well known robust covariance estimators with good statistical properties, like the MCD or the S-estimators (see, e.g. Maronna et al., 2006), rely on a full-rank data matrix Y which is not the case for clr transformed data (see, e.g., Aitchison, 1986). The isometric logratio (ilr) transformation (Egozcue et al., 2003) solves this singularity problem. The data matrix Y is transformed to a matrix Z by using an orthonormal basis of lower dimension. Using the ilr transformed data, a robust covariance matrix C(Z) can be estimated. The result can be back-transformed to the clr space by C(Y ) = V C(Z)V T where the matrix V with orthonormal columns comes from the relation between the clr and the ilr transformation. Now the parameters in the model (2) can be estimated (Basilevsky, 1994) and the results have a direct interpretation since the links to the original variables are still preserved. The above procedure will be applied to data from geochemistry. Our special interest is on comparing the results with those of Reimann et al. (2002) for the Kola project data