La leptospirose féline : sondage sérologique et de PCR urinaire chez des chats sains et des chats atteints de maladie rénale

La leptospirose est une zoonose à distribution mondiale dont la prévalence chez le chat varie géographiquement de 4.8% à 35%. Bien que l’exposition féline à Leptospira spp. soit rapportée dans des études sérologiques, les conséquences cliniques de cette maladie chez le chat sont peu connues. Le but principal de cette étude était de comparer le statut sérologique et de porteur (PCR urinaire) de Leptospira spp. entre des chats sains et des chats atteints de maladie rénale (insulte rénale aigue et maladie rénale chronique de stades IIb, III et IV). Une étude préliminaire pour valider la sensibilité et la spécificité analytiques de la PCR de Leptospira spp. réalisée par le Laboratoire de Diagnostic Moléculaire de la FMV sur l’urine de chat a été effectuée. La validation in vitro a démontré que la technique de PCR est efficace pour déterminer la présence de leptospires pathogènes dans l’urine du chat. Dans le cadre de l’étude principale, 251 chats ont été recrutés entre janvier 2010 et mars 2012,. De ceux-ci, 240 ont été inclus et divisés en 2 groupes (chats sains (C=125) et chats atteints de maladie rénale (MR=115) en se basant sur un examen physique ainsi que sur des résultats d’hématologie, de biochimie et d’analyse d’urine. Tous les chats recrutés ont également été examinés sérologiquement par test de micro-agglutination pour la présence d’anticorps contre Leptospira spp. (résultat considéré positif si ≥1 :100) et par PCR pour la présence de Leptospira spp. dans l’urine. Le pourcentage prédit de séropositivité pour Leptospira spp. était significativement plus élevé chez les chats atteints de maladie rénale (13,7%) que chez les chats sains (5%) (p=0,02). Les sérovars impliqués étaient Pomona (n=16), Bratislava (n=8) et Grippotyphosa (n=1). De plus, les chats séropositifs pour Pomona présentaient des titres significativement plus élevés que pour les autres sérovars (p=0,04). L’excrétion de Leptospira spp. a été confirmée par PCR dans l’urine de huit chats. Des 26 chats séropositifs, quatre (C=2, MR=2) se sont également révélés PCR positifs. La prévalence a été plus élevée chez les chats du groupe MR (5.3%; 6/113) lorsque comparée à celle du groupe C (1.6%; 2/125), mais cette différence ne s’est pas révélée statistiquement significative (C=0,9% , MR= 5,5% ; p = 0,09). L’âge, le sexe et le milieu de vie (urbain versus rural) n’ont pas influencé le statut sérologique ou d’excrétion pour Leptospira spp. Le pourcentage prédit de séropositivité était significativement plus élevée chez les chasseurs (p<0.01) et pendant les mois de juin à août (p=0.02). La présence d’un autre chat à la maison a également significativement augmenté ce pourcentage (p<0.01), mais la présence d’un chien ne l’a pas influencé. Lors de l’évaluation du PCR par le modèle GGE, seules les variables « contact avec raton laveur » et « contact avec mouffettes » sont ressorties statistiquement significatives (p≤0.03). Le rôle que joue Leptospira spp. comme agent étiologique de maladie rénale chez le chat demeure incertain. Toutefois, la différence significative de statut sérologique entre les chats sains et les chats atteints de maladie rénale suggère que la leptospirose pourrait être une cause sous-diagnostiquée de maladie rénale chez cette espèce. Dans cette étude, plusieurs porteurs asymptomatiques ont été identifiés, ce qui suggère que l’espèce féline puisse être un acteur sous-estimé dans la transmission de la bactérie aux humains.

La régulation du gène CYP19A1 dans les cellules de granulosa bovine in vitro

L’oestradiol joue un rôle important dans la reproduction en général, particulièrement dans la croissance folliculaire chez la vache. La production de l’œstradiol nécessite l’expression du gène CYP19A1 suite à la stimulation des cellules de granulosa par l’hormone folliculostimulante (FSH) ou le facteur de croissance insulinique de type 1 (IGF-1). Chez la vache, il existe six promoteurs (1.1 ; 1.2 ; 1.3 ; 1.4 ; 1.5 et 2) qui dirigent la transcription du gène CYP19A1 dans les cellules de la granulosa. Le principal promoteur qui dirige la transcription au niveau de l’ovaire (cellules de granulosa) est le promoteur 2 (P2). Cependant, l’effet de la FSH et de l’IGF-1 sur l’activation de ces promoteurs d’aromatase demeure mal connu. De plus, la demi-vie du transcrit CYP19A1 est très courte avec une région 3’UTR relativement longue. L’analyse de la séquence 3’UTR montre la présence des motifs ARE (séquence riche en AU), des études antérieur montrent que ces séquences impliquent dans la régulation de la stabilité ou la dégradation de l’ARNm, ce qui est fort probable que la courte demi-vie de l’ARNm CYP19A1 est sous le contrôle post-transcriptionel. L’objectif de la thèse visait à étudier la régulation de l’expression du gène CYP19A1 chez la vache. Il y a deux thèmes soit étude de la régulation transcriptionnelle ciblant le promoteur et soit étude de la régulation post-transcriptionnelle impliquant la région 3’non traduite (3’UTR). Le premier objectif vise à étudier la régulation transcriptionnelle du gène CYP19A1. Nous avons étudié l'activité du promoteur ovarien bovin dans deux modèles de cellules de la granulosa, les cellules lutéinisées et nonlutéinisées in vitro, suite à une stimulation des cellules par la FSH ou IGF-1. Nous avons également évalué la voie de signalisation impliquée dans la régulation des différents promoteurs en utilisant un RT-PCR et un gène rapporteur (les différents promoteurs d’aromatase ont été insérés dans le vecteur pGL3promoter en amont du gène exprimant la luciférase). Les résultats de RT-PCR démontrent que la FSH et l’IGF-1 augmentent les concentrations d’ARNm provenant des deux promoteurs 2 et 1.1 dans les cellules de la granulosa non lutéinisées. Des expériences subséquentes ont montré que la FSH stimule le promoteur 2 via la voie PKA tandis que l'IGF-1 stimule le promoteur 2 via la voie PKC. La FSH et l’IGF-1 stimulent l’expression du promoteur 1.1 via la voie PI3K. L’analyse de l’activité luciférase démontre que dans les cellules de granulosa lutéinisées, la FSH stimule le promoteur 1.1 de façon dose dépendante et ne semble y avoir aucun effet significatif sur le promoteur 2. Nous avons donc comparé l’activité du promoteur PII/P2 humain, du rat, de la chèvre et de la vache dans les cellules de granulosa bovine lutéinisées. Le résultat le plus significatif est que le promoteur 2 bovine (et caprine) dépend de plusieurs facteurs de transcription (NR5A2, FOXL2) comparé au promoteur PII humain et celui du promoteur proximal du rat qui dépendent principalement de l'AMPc. En effet, nos résultats ont démontré une expression raisonnablement robuste du P2 bovine lorsque les cellules sont traitées à la forskoline, NR5A2 et FOXL2. Le facteur FOXL2 semble déterminer l'activité du promoteur 2 chez le ruminant. Le deuxième objectif vise à étudier la régulation post-transcriptionnelle du gène CYP19A1. Pour ce faire, nous avons déterminé la séquence minimale de l'ARNm CYP19A1 requise pour la régulation de sa demi-vie. Différents séquences de la région 3’UTR ont été insérés dans le vecteur pGL3promoter en aval du gène exprimant la luciférase ou soit dans le vecteur pGEMTeasy. Le vecteur pGL3promoter a été transfecté dans les cellules de granulosa lutéinisées pour évaluer l'impact de la séquence 3'UTR sur l'expression du gène rapporteur de la luciférase, alors que le vecteur pGEMTeasy a été utilisé pour la transcription in vitro afin de générer de l’ARNm. Ce dernier sera utilisé en réaction croisée au UV avec des extraits protéiques pour démontrer l’association du complexe ARNm/protéine. L’analyse de l’activité luciférase a permis d’identifier une séquence de 200 pb située entre 926 et 1134 pb de la région 3'UTR de l’ARNm CYP19A1 qui a réduit significativement l’activité de la luciférase. Selon les analyses de la réaction croisée au UV, une ou plusieurs protéines de 66 et 80 kDA se lient spécifiquement à la séquence de 200 pb qui réduit l’activité de luciférase. Cette protéine s'exprime dans les cellules de granulosa, mais n’a pas été détectée dans d'autres tissus comme le foie et le cœur. Par ailleurs, l’utilisation du gène rapporteur sensible à la FSH a suscité l’intérêt d'une compagnie pharmaceutique qui vend de l’equine chorionic gonadotropin (eCG) pour lui permettre de distinguer facilement l’eCG ayant une forte activité FSH et donc, avoir un produit commercial plus efficace et de meilleure qualité. Dans cette étude, nous avons développé un système de bioessai à la FSH basé sur la transfection des cellules avec un récepteur à la FSH et un gène rapporteur colorimètrique qui permet d’estimer l’activité de la FSH dans le sérum de la jument et qui pourrait être applicable au niveau de la ferme/industrie.

Étude du goulot d’étranglement dans la transmission mère-enfant du virus de l’hépatite C.

La transmission mère-enfant (TME) du virus de l’hépatite C (VHC) est la première cause d’acquisition de l’infection chez les enfants des pays développés. Celle-ci prend place dans <10% des cas. Toutefois, dans le cas d’une coinfection maternelle avec le virus de l’immunodéficience de type 1 (VIH-1), ce taux est accru alors qu’il n’existe aucune intervention préventive de la TME du VHC. Le VHC arbore une diversité importante qui est le résultat d’une réplication exempte de mécanisme de correction. Il est donc retrouvé chez son hôte sous la forme d’un spectre de virions génétiquement apparentés mais différents qu’on appelle quasiespèce. Lorsque le VHC est transmis entre adultes, seulement un nombre limité de variantes sont responsables de l’infection, c’est ce qu’on appelle un goulot d’étranglement génétique. L’existence d’un tel profil de transmission lors de la TME du VHC restait, jusqu’à maintenant, à confirmer. En se basant sur la détection par RT-PCR de la virémie à la naissance, la TME du VHC est réputée prendre place in utero et peripartum, une dynamique de transmission qui reste à démontrer. Ici, nous rapportons une analyse longitudinale de la TME du VHC par séquençage de nouvelle génération chez 5 paires mère-enfant dont 3 mères sont également coinfectées avec le VIH-1. L’analyse de l’identité des variantes virales basée sur la séquence nucléotidique des régions hypervariables 1-2 de la glycoprotéine E2 (positions 1491-1787 de l’isolat H77) révèle qu’un nombre limité de variantes virales sont transmises de la mère à l’enfant lorsque la mère est seulement infectée par le VHC (n = 1-4 variantes transmises). Dans le cas de la coinfection maternelle avec le VIH-1, ce nombre est toutefois drastiquement plus important (n = 111-118). La détection de variantes retrouvées chez la mère au deuxième trimestre et l’enfant mais non détectées subséquemment chez la mère témoigne que la TME du VHC peut prendre place aussi tôt que lors du deuxième trimestre de grossesse. Finalement, nous montrons que la dynamique d’infection chez l’enfant implique une augmentation transitoire de la virémie concomitante avec une perte de diversité de la quasiespèce. Dans l’ensemble ces résultats sont les premiers à démontrer directement l’existence d’un goulot d’étranglement lors de la TME du VHC. Celui-ci serait moins restringent dans le cas de la coinfection maternelle avec le VIH-1. Cette transmission peut prendre place aussi tôt que lors du deuxième trimestre de grossesse et il semblerait qu’un spectre limité de variantes soit responsable pour l’établissement de l’essentiel de la production virale chez le jeune enfant.

Emergence of a new type of porcine circovirus (PCV) in swine: a type 1 and type 2 PCV recombinant.

In late September 2008, tissue samples from piglets experiencing an acute outbreak of porcine reproductive and respiratory syndrome (PRRS) were submitted to the Veterinary diagnostic service of the University of Montreal. Several diagnostic assays were performed including a multiplex real-time quantitative PCR assay (mrtqPCR) for the detection and differentiation of porcine circovirus (PCV) type 2a and 2b genotypes in the lung and lymph nodes. The pig samples were found to be positive for PCV2a using the mrtqPCR but odd results were obtained. The Ct values obtained with mrtqPCR probes targeting the ORF1 and ORF2 of PCV2 were not as expected which suggested the presence of genomic variations in the PCV2 viral genome. Ultimately, a total of three diagnostic cases with mrtqPCR unusual results were investigated. After virus isolation and sequence analyses, a new type of PCV was identified in those three cases. Based on sequence analyses, this new PCV genome contains the ORF1 of PCV1 and the ORF2 of PCV2a and its entire viral genome nucleotide identity compared to PCV1, PCV2a and 2b are 86.4%, 88.7% and 86.5%, respectively. It is proposed to name this new PCV by taking into account the nomenclature of Segales et al. (2008) and by indicating the origin of the ORF1 at first and the origin of the ORF2 in second. Consequently, the name proposed for this new PCV is PCV1/2a. The prevalence of PCV1/2a seems to be very low in Quebec, Canada (2.5% of PCV positive cases), and its origin is now in debate.

La génomique évolutive mitochondriale révèle des échanges génétiques et la ségrégation chez les Gloméromycètes

Les champignons mycorhiziens à arbuscules (CMA) sont des organismes microscopiques du sol qui jouent un rôle crucial dans les écosystèmes naturels et que l’on retrouve dans tous les habitats de la planète. Ils vivent en relation symbiotique avec la vaste majorité des plantes terrestres. Ils sont des biotrophes obligatoires, c'est-à-dire qu'ils ne peuvent croître qu'en présence d'une plante hôte. Cette symbiose permet entre autres à la plante d'acquérir des nutriments supplémentaires, en particulier du phosphore et du nitrate. Malgré le fait que cette symbiose apporte des services importants aux écosystèmes, la richesse des espèces, la structure des communautés, ainsi que la diversité fonctionnelle des CMA sont mal connues et l'approfondissement des connaissances dans ces domaines dépend d’outils de diagnostic moléculaire. Cependant, la présence de polymorphisme nucléaire intra-isolat combiné à un manque de données génomiques dans différents groupes phylogénétique de ces champignons complique le développement de marqueurs moléculaires et la détermination de l'affiliation évolutive à hauts niveaux de résolution (c.a.d. entre espèces génétiquement similaires et/ou isolats de la même espèce). . Pour ces raisons, il semble une bonne alternative d’utiliser un système génétique différent en ciblant le génome mitochondrial, qui a été démontré homogène au sein d'un même isolat de CMA. Cependant, étant donné le mode de vie particulier de ces organismes, une meilleure compréhension des processus évolutifs mitochondriaux est nécessaire afin de valoriser l'utilisation de tels marqueurs dans des études de diversité et en génétique des populations. En ce sens, mon projet de doctorat consistait à investiguerétudier: i) les vecteurs de divergences inter-isolats et -espèces génétiquement rapprochéesphylogénétiquement apparentées, ii) la plasticité des génomes mitochondriaux, iii) l'héritabilité mitochondriale et les mécanismes potentiels de ségrégation, ainsi que iv) la diversité mitochondriale intra-isolat in situ. À l'aide de la génomique mitochondriale comparative, en utilisant le séquençage nouvelle génération, on a démontré la présence de variation génétique substantielle inter-isolats et -espèces, engendrées par l'invasion d'éléments mobiles dans les génomes mitochondriaux des CMA, donnant lieu à une évolution moléculaire rapide des régions intergéniques. Cette variation permettait de développer des marqueurs spécifiques à des isolats de la même espèce. Ensuite, à l'aide d'une approche analytique par réseaux de gènes sur des éléments mobiles, on a été en mesure de démontrer des évènements de recombinaisons homologues entre des haplotypes mitochondriaux distincts, menant à des réarrangements génomiques. Cela a permis d'ouvrir les perspectives sur la dynamique mitochondriale et l'hétéroplasmie dans un même isolatsuggère une coexistence de différents haplotypes mitochondriaux dans les populations naturelles et que les cultures monosporales pourraient induirent une sous-estimation de la diversité allélique mitochondriale. Cette apparente contradiction avec l'homogénéité mitochondriale intra-isolat généralement observée, a amené à investiguer étudier les échanges génétiques à l'aide de croisements d'isolats génétiquement distincts. Malgré l'observation de quelques spores filles hétéroplasmiques, l'homoplasmie était le statut par défaut dans toutes les cultures monosporales, avec un biais en faveur de l'un des haplotypes parentaux. Ces résultats suggèrent que la ségrégation opère durant la formation de la spore et/ou le développement de la coloniedu mycélium. De plus, ils supportent la présence d'une machinerie protéique de ségrégation mitochondriale chez les CMAAMF, où l'ensemble des gènes impliqués dans ce mécanisme ont été retrouvé et sont orthologues aux autres champignons. Finalement, on est revenue aux sources avecon a étudié le polymorphisme mitochondrial intra-isolat à l'aide d'une approche conventionnelle de PCR en utilisant une Taq polymérase de haute fidélité, suivie de clonage et de séquençage Sanger, sur deux isolats de R. irregularis. Cela a permis l'observation d'hétéroplasmie in situ, ainsi que la co-expression de variantes de variantes de protéines'ARNm dans une souche in vitro. Les résultats suggèrent que d'autres études basées sur le séquençage nouvelle génération aurait potentiellement ignorée cette variation, offrant ainsi plusieurs nouveaux arguments permettant de considérer les CMA comme des organismes possédant une population de génomes mitochondriaux et nucléaires distincts.

Développement des marqueurs moléculaires spécifiques pour l'identification et la quantification de deux espèces de champignons mycorhiziens arbusculaires en utilisant la PCR en temps réel

Les champignons mycorhizien à arbuscules (CMA) sont des organismes pouvant établir des symbioses avec 80% des plantes terrestres. Les avantages d'une telle symbiose sont de plus en plus caractérisés et exploités en agriculture. Par contre, jusqu'à maintenant, il n'existe aucun outil permettant à la fois l'identification et la quantification de ces champignons dans le sol de façon fiable et rapide. Un tel outil permettrait, entre autres, de mieux comprendre les dynamiques des populations des endomycorhizes dans le sol. Pour les producteurs d'inoculum mycorhiziens, cela permettrait également d'établir un suivi de leurs produits en champs et d'avoir un contrôle de qualité de plus sur leurs inoculants. C'est ce que nous avons tenté de développer au sein du laboratoire du Dr. Hijri. Depuis environ une trentaine d'années, des outils d'identification et/ou de quantification ont été développés en utilisant les profiles d'acides gras, les isozymes, les anticorps et finalement l'ADN nucléaire. À ce jour, ces méthodes d’identification et de quantification sont soit coûteuses, soit imprécises. Qui plus est, aucune méthode ne permet à la fois la quantification et l’identification de souches particulières de CMA. L’ADN mitochondrial ne présente pas le même polymorphisme de séquence que celui qui rend l’ADN nucléaire impropre à la quantification. C'est pourquoi nous avons analysé les séquences d’ADN mitochondrial et sélectionné les régions caractéristiques de deux espèces de champignons mycorhiziens arbusculaires (CMA). C’est à partir de ces régions que nous avons développé des marqueurs moléculaires sous forme de sondes et d’amorces TaqMan permettant de quantifier le nombre de mitochondries de chacune de ces espèces dans un échantillon d’ADN. Nous avons ensuite tenté de déterminer une unité de quantification des CMA, soit un nombre de mitochondries par spore. C’est alors que nous avons réalisé que la méthode de préparation des échantillons de spores ainsi que la méthode d’extraction d’ADN avaient des effets significatifs sur l’unité de quantification de base. Nous avons donc optimisé ces protocoles, avant d’en e tester l’application sur des échantillons de sol et de racines ayant été inoculés avec chacune des deux espèces cibles. À ce stade, cet outil est toujours semi-quantificatif, mais il permet 9 l’identification précise de deux espèces de CMA compétentes dans des milieux saturés en phosphore inorganique. Ces résultats , en plus d’être prometteurs, ont permis d’augmenter les connaissances méthodologiques reliées à la quantification des CMA dans le sol, et suggèrent qu’à cause de leurs morphologies différentes, l’élaboration d’un protocole de quantification standardisé pour toutes les espèces de CMA demeure un objectif complexe, qui demande de nouvelles études in vivo.

Foraminiferal and Pteropod Records in the Late Quaternary Sequence off Northern Kerala and Their Oceanographic and Climatic Significance

In the present investigation, an attempt is made to study late Quaternary foraminiferal and pteropod records of the shelf of northern Kerala and to evaluate their potentiality in paleocenographic and paleoclimatic reconstruction. The study gives details of sediment cores, general characteristics of foraminifera and pteropod species recorded from the examined samples and their systematic classification, spatial distribution of Recent foraminifera and pteropods and their response to varying bathymetry, nature of substrate, organic matter content in sediment and hydrography across the shelf. An attempt is also made to establish an integrated chronostratigraphy for the examined core sections. An effort is also made to identify microfaunal criteria useful in biostratigraphic division in shallow marine core sections. An attempt is made to infer various factors responsible for the change in microfaunal assemblage. Reconstruction of sea level changes during the last 36,000 years was attempted based on the pteropod record. The study reveals a bathymetric control on benthic/planktic (BF/PF) foraminiferal and pteropods/planktic foraminiferal (Pt/PF) abundance ratio. Bathymetric distribution pattern of BF/PF ratio is opposite to the (Pt/PF) ratio with decreasing trend of former from the shore across the shelf. Quantitative benthic foraminiferal record in the surficial sediments reveals a positive correlation between the diversity and bathymetry. R-mode cluster analysis performed on 30n significant Recent benthic foraminiferal, determines three major assemblage.

An Alternative Approach To Computer System Security Monitoring And Enhancement Through System Call Sequence Analysis

Modern computer systems are plagued with stability and security problems: applications lose data, web servers are hacked, and systems crash under heavy load. Many of these problems or anomalies arise from rare program behavior caused by attacks or errors. A substantial percentage of the web-based attacks are due to buffer overflows. Many methods have been devised to detect and prevent anomalous situations that arise from buffer overflows. The current state-of-art of anomaly detection systems is relatively primitive and mainly depend on static code checking to take care of buffer overflow attacks. For protection, Stack Guards and I-leap Guards are also used in wide varieties.This dissertation proposes an anomaly detection system, based on frequencies of system calls in the system call trace. System call traces represented as frequency sequences are profiled using sequence sets. A sequence set is identified by the starting sequence and frequencies of specific system calls. The deviations of the current input sequence from the corresponding normal profile in the frequency pattern of system calls is computed and expressed as an anomaly score. A simple Bayesian model is used for an accurate detection.Experimental results are reported which show that frequency of system calls represented using sequence sets, captures the normal behavior of programs under normal conditions of usage. This captured behavior allows the system to detect anomalies with a low rate of false positives. Data are presented which show that Bayesian Network on frequency variations responds effectively to induced buffer overflows. It can also help administrators to detect deviations in program flow introduced due to errors.

Mean Squared Residue Based Biclustering Algorithms for the Analysis of Gene Expression Data

Computational Biology is the research are that contributes to the analysis of biological data through the development of algorithms which will address significant research problems.The data from molecular biology includes DNA,RNA ,Protein and Gene expression data.Gene Expression Data provides the expression level of genes under different conditions.Gene expression is the process of transcribing the DNA sequence of a gene into mRNA sequences which in turn are later translated into proteins.The number of copies of mRNA produced is called the expression level of a gene.Gene expression data is organized in the form of a matrix. Rows in the matrix represent genes and columns in the matrix represent experimental conditions.Experimental conditions can be different tissue types or time points.Entries in the gene expression matrix are real values.Through the analysis of gene expression data it is possible to determine the behavioral patterns of genes such as similarity of their behavior,nature of their interaction,their respective contribution to the same pathways and so on. Similar expression patterns are exhibited by the genes participating in the same biological process.These patterns have immense relevance and application in bioinformatics and clinical research.Theses patterns are used in the medical domain for aid in more accurate diagnosis,prognosis,treatment planning.drug discovery and protein network analysis.To identify various patterns from gene expression data,data mining techniques are essential.Clustering is an important data mining technique for the analysis of gene expression data.To overcome the problems associated with clustering,biclustering is introduced.Biclustering refers to simultaneous clustering of both rows and columns of a data matrix. Clustering is a global whereas biclustering is a local model.Discovering local expression patterns is essential for identfying many genetic pathways that are not apparent otherwise.It is therefore necessary to move beyond the clustering paradigm towards developing approaches which are capable of discovering local patterns in gene expression data.A biclusters is a submatrix of the gene expression data matrix.The rows and columns in the submatrix need not be contiguous as in the gene expression data matrix.Biclusters are not disjoint.Computation of biclusters is costly because one will have to consider all the combinations of columans and rows in order to find out all the biclusters.The search space for the biclustering problem is 2 m+n where m and n are the number of genes and conditions respectively.Usually m+n is more than 3000.The biclustering problem is NP-hard.Biclustering is a powerful analytical tool for the biologist.The research reported in this thesis addresses the problem of biclustering.Ten algorithms are developed for the identification of coherent biclusters from gene expression data.All these algorithms are making use of a measure called mean squared residue to search for biclusters.The objective here is to identify the biclusters of maximum size with the mean squared residue lower than a given threshold. All these algorithms begin the search from tightly coregulated submatrices called the seeds.These seeds are generated by K-Means clustering algorithm.The algorithms developed can be classified as constraint based,greedy and metaheuristic.Constarint based algorithms uses one or more of the various constaints namely the MSR threshold and the MSR difference threshold.The greedy approach makes a locally optimal choice at each stage with the objective of finding the global optimum.In metaheuristic approaches particle Swarm Optimization(PSO) and variants of Greedy Randomized Adaptive Search Procedure(GRASP) are used for the identification of biclusters.These algorithms are implemented on the Yeast and Lymphoma datasets.Biologically relevant and statistically significant biclusters are identified by all these algorithms which are validated by Gene Ontology database.All these algorithms are compared with some other biclustering algorithms.Algorithms developed in this work overcome some of the problems associated with the already existing algorithms.With the help of some of the algorithms which are developed in this work biclusters with very high row variance,which is higher than the row variance of any other algorithm using mean squared residue, are identified from both Yeast and Lymphoma data sets.Such biclusters which make significant change in the expression level are highly relevant biologically.

β-glucosidases from Aspergillus unguis NII 08123: Molecular characterization, properties and applications

Lignocellulosic biomass is probably the best alternative resource for biofuel production and it is composed mainly of cellulose, hemicelluloses and lignin. Cellulose is the most abundant among the three and conversion of cellulose to glucose is catalyzed by the enzyme cellulase. Cellulases are groups of enzymes act synergistically upon cellulose to produce glucose and comprise of endoglucanase, cellobiohydrolase and β-glucosidase. β -glucosidase assumes great importance due to the fact that it is the rate limiting enzyme. Endoglucanases (EG) produces nicks in the cellulose polymer exposing reducing and non reducing ends, cellobiohydrolases (CBH) acts upon the reducing or non reducing ends to liberate cellobiose units, and β - glucosidases (BGL) cleaves the cellobiose to liberate glucose completing the hydrolysis. . β -glucosidases undergo feedback inhibition by their own product- β glucose, and cellobiose which is their substrate. Few filamentous fungi produce glucose tolerant β - glucosidases which can overcome this inhibition by tolerating the product concentration to a particular threshold. The present study had targeted a filamentous fungus producing glucose tolerant β - glucosidase which was identified by morphological as well as molecular method. The fungus showed 99% similarity to Aspergillus unguis strain which comes under the Aspergillus nidulans group where most of the glucose tolerant β -glucosidase belongs. The culture was designated the strain number NII 08123 and was deposited in the NII culture collection at CSIR-NIIST. β -glucosidase multiplicity is a common occurrence in fungal world and in A.unguis this was demonstrated using zymogram analysis. A total 5 extracellular isoforms were detected in fungus and the expression levels of these five isoforms varied based on the carbon source available in the medium. Three of these 5 isoforms were expressed in higher levels as identified by the increased fluorescence (due to larger amounts of MUG breakdown by enzyme action) and was speculated to contribute significantly to the total _- β glucosidase activity. These isoforms were named as BGL 1, BGL3 and BGL 5. Among the three, BGL5 was demonstrated to be the glucose tolerant β -glucosidase and this was a low molecular weight protein. Major fraction was a high molecular weight protein but with lesser tolerance to glucose. BGL 3 was between the two in both activity and glucose tolerance.121 Glucose tolerant .β -glucosidase was purified and characterized and kinetic analysis showed that the glucose inhibition constant (Ki) of the protein is 800mM and Km and Vmax of the enzyme was found to be 4.854 mM and 2.946 mol min-1mg protein-1respectively. The optimumtemperature was 60°C and pH 6.0. The molecular weight of the purified protein was ~10kDa in both SDS as well as Native PAGE indicating that the glucose tolerant BGL is a monomeric protein.The major β -glucosidase, BGL1 had a pH and temperature optima of 5.0 and 60 °C respectively. The apparent molecular weight of the Native protein is 240kDa. The Vmax and Km was 78.8 mol min-1mg protein-1 and 0.326mM respectively. Degenerate primers were designed for glycosyl hydrolase families 1, 3 and 5 and the BGL genes were amplified from genomic DNA of Aspergillus unguis. The sequence analyses performed on the amplicons results confirmed the presence of all the three genes. Amplicon with a size of ~500bp was sequenced and which matched to a GH1 –BGL from Aspergillus oryzae. GH3 degenerate primers producing amplicons were sequenced and the sequences matched to β - glucosidase of GH3 family from Aspergillus nidulans and Aspergillus acculateus. GH5 degenerate primers also gave amplification and sequencing results indicated the presence of GH5 family BGL gene in the Aspergillus unguis genomic DNA.From the partial gene sequencing results, specific as well as degenerate primers were designed for TAIL PCR. Sequencing results of the 1.0 Kb amplicon matched Aspergillus nidulans β -glucosidase gene which belongs to the GH1 family. The sequence mainly covered the N-Terminal region of the matching peptide. All the three BGL proteins ie. BGL1, BGL3 and BGL5 were purified by chromatography an electro elution from Native PAGE gels and were subjected to MALDI-TOF mass spectrometric analysis. The results showed that BGL1 peptide mass matched to . β -glucosidase-I of Aspergillus flavus which is a 92kDa protein with 69% protein coverage. The glucose tolerant β -glucosidase BGL5 mass matched to the catalytic C-terminal domain of β -glucosidase-F from Emericella nidulans, but the protein coverage was very low compared to the size of the Emericella nidulans protein. While comparing the size of BGL5 from Aspergillus unguis, the protein sequence coverage is more than 80%. BGL F is a glycosyl hydrolase family 3 protein.The properties of BGL5 seem to be very unique, in that it is a GH3 β -glucosidase with a very low molecular weight of ~10kDa and at the same time having catalytic activity and glucose 122 tolerance which is as yet un-described in GH β -glucosidases. The occurrence of a fully functional 10kDA protein with glucose tolerant BGL activity has tremendous implications both from the points of understanding the structure function relationships as well as for applications of BGL enzymes. BGL-3 showed similarity to BGL1 of Aspergillus aculateus which was another GH3 β -glucosidase. It may be noted that though PCR could detect GH1, GH3 and GH5 β-glucosidases in the fungus, the major isoforms BGL1 BGL3 and BGL5 were all GH3 family enzymes. This would imply that β-glucosidases belonging to other families may also co-exist in the fungus and the other minor isoforms detected in zymograms may account for them. In biomass hydrolysis, GT-BGL containing BGL enzyme was supplemented to cellulase and the performances of blends were compared with a cocktail where commercial β- glucosidase was supplemented to the biomass hydrolyzing enzyme preparation. The cocktail supplemented with A unguis BGL preparation yielded 555mg/g sugar in 12h compared to the commercial enzyme preparation which gave only 333mg/g in the same period and the maximum sugar yield of 858 mg/g was attained in 36h by the cocktail containing A. unguis BGL. While the commercial enzyme achieved almost similar sugar yield in 24h, there was rapid drop in sugar concentration after that, indicating probably the conversion of glucose back to di-or oligosaccharides by the transglycosylation activity of the BGl in that preparation. Compared this, the A.unguis enzyme containing preparation supported peak yields for longer duration (upto 48h) which is important for biomass conversion to other products since the hydrolysate has to undergo certain unit operations before it goes into the next stage ie – fermentation in any bioprocesses for production of either fuels or chemicals.. Most importantly the Aspergillus unguis BGL preparation yields approximately 1.6 fold increase in the sugar release compared to the commercial BGL within 12h of time interval and 2.25 fold increase in the sugar release compared to the control ie. Cellulase without BGL supplementation. The current study therefore leads to the identification of a potent new isolate producing glucose tolerant β - glucosidase. The organism identified as Aspergillus unguis comes under the Aspergillus nidulans group where most of the GT-BGL producers belong and the detailed studies showed that the glucose tolerant β -glucosidase was a very low molecular weight protein which probably belongs to the glycosyl hydrolase family 3. Inhibition kinetic studies helped to understand the Ki and it is the second highest among the nidulans group of Aspergilli. This has promoted us for a detailed study regarding the mechanism of glucose tolerance. The proteomic 123 analyses clearly indicate the presence of GH3 catalytic domain in the protein. Since the size of the protein is very low and still its active and showed glucose tolerance it is speculated that this could be an entirely new protein or the modification of the existing β -glucosidase with only the catalytic domain present in it. Hydrolysis experiments also qualify this BGL, a suitable candidate for the enzyme cocktail development for biomass hydrolysis

Characterization and Pathogenicity of Vibrio cholerae and Vibrio vulnificus from Marine environments

The genus Vibrioof the family Vibrionaceae are Gram negative, oxidasepositive, rod- or curved- rodshaped facultative anaerobes, widespread in marine and estuarine environments. Vibrio species are opportunistic human pathogens responsible for diarrhoeal disease, gastroenteritis, septicaemia and wound infections and are also pathogens of aquatic organisms, causing infections to crustaceans, bivalves and fishes. In the present study, marine environmental samples like seafood and water and sediment samples from aquafarms and mangroves were screened for the presence of Vibrio species. Of the134 isolates obtained from the various samples, 45 were segregated to the genus Vibrio on the basis of phenotypic characterization.like Gram staining, oxidase test, MoF test and salinity tolerance. Partial 16S rDNA sequence analysis was utilized for species level identification of the isolates and the strains were identified as V. cholerae(N=21), V. vulnificus(N=18), V. parahaemolyticus(N=3), V. alginolyticus (N=2) and V. azureus (N=1). The genetic relatedness and variations among the 45 Vibrio isolates were elucidated based on 16S rDNA sequences. Phenotypic characterization of the isolates was based on their response to 12 biochemical tests namely Voges-Proskauers’s (VP test), arginine dihydrolase , tolerance to 3% NaCl test, ONPG test that detects β-galactosidase activity, and tests for utilization of citrate, ornithine, mannitol, arabinose, sucrose, glucose, salicin and cellobiose. The isolates exhibited diverse biochemical patterns, some specific for the species and others indicative of their environmental source.Antibiogram for the isolates was determined subsequent to testing their susceptibility to 12 antibiotics by the disc diffusion method. Varying degrees of resistance to gentamycin (2.22%), ampicillin(62.22%), nalidixic acid (4.44%), vancomycin (86.66), cefixime (17.77%), rifampicin (20%), tetracycline (42.22%) and chloramphenicol (2.22%) was exhibited. All the isolates were susceptible to streptomycin, co-trimoxazole, trimethoprim and azithromycin. Isolates from all the three marine environments exhibited multiple antibiotic resistance, with high MAR index value. The molecular typing methods such as ERIC PCR and BOX PCR revealed intraspecies relatedness and genetic heterogeneity within the environmental isolatesof V. cholerae and V. vulnificus. The 21 strains of V. choleraewere serogroupedas non O1/ non O139 by screening for the presence O1rfb and O139 rfb marker genes by PCR. The virulence/virulence associated genes namely ctxA, ctxB, ace, VPI, hlyA, ompU, rtxA, toxR, zot, nagst, tcpA, nin and nanwere screened in V. cholerae and V. vulnificusstrains.The V. vulnificusstrains were also screened for three species specific genes viz., cps, vvhand viu. In V. cholerae strains, the virulence associated genes like VPI, hlyA, rtxA, ompU and toxR were confirmed by PCR. All the isolates, except for strain BTOS6, harbored at least one or a combination of the tested genes and V. choleraestrain BTPR5 isolated from prawn hosted the highest number of virulence associated genes. Among the V. vulnificusstrains, only 3 virulence genes, VPI, toxR and cps, were confirmed out of the 16 tested and only 7 of the isolates had these genes in one or more combinations. Strain BTPS6 from aquafarm and strain BTVE4 from mangrove samples yielded positive amplification for the three genes. The toxRgene from 9 strains of V. choleraeand 3 strains of V. vulnificus were cloned and sequenced for phylogenetic analysis based on nucleotide and the amino acid sequences. Multiple sequence alignment of the nucleotide sequences and amino acid sequences of the environmental strains of V. choleraerevealed that the toxRgene in the environmental strains are 100% homologous to themselves and to the V. choleraetoxR gene sequence available in the Genbank database. The 3 strains of V. vulnificus displayed high nucleotide and amino acid sequence similarity among themselves and to the sequences of V. cholerae and V. harveyi obtained from the GenBank database, but exhibited only 72% homology to the sequences of its close relative V. vulnificus. Structure prediction of the ToxR protein of Vibrio cholerae strain BTMA5 was by PHYRE2 software. The deduced amino acid sequence showed maximum resemblance with the structure of DNA-binding domain of response regulator2 from Escherichia coli k-12 Template based homology modelling in PHYRE2 successfully modelled the predicted protein and its secondary structure based on protein data bank (PDB) template c3zq7A. The pathogenicity studies were performed using the nematode Caenorhabditiselegansas a model system. The assessment of pathogenicity of environmental strain of V. choleraewas conducted with E. coli strain OP50 as the food source in control plates, environmental V. cholerae strain BTOS6, negative for all tested virulence genes, to check for the suitability of Vibrio sp. as a food source for the nematode;V. cholerae Co 366 ElTor, a clinical pathogenic strain and V. cholerae strain BTPR5 from seafood (Prawn) and positive for the tested virulence genes like VPI, hlyA, ompU,rtxA and toxR. It was found that V. cholerae strain BTOS6 could serve as a food source in place of E. coli strain OP50 but behavioral aberrations like sluggish movement and lawn avoidance and morphological abnormalities like pharyngeal and intestinal distensions and bagging were exhibited by the worms fed on V. cholerae Co 366 ElTor strain and environmental BTPR5 indicating their pathogenicity to the nematode. Assessment of pathogenicity of the environmental strains of V. vulnificus was performed with V. vulnificus strain BTPS6 which tested positive for 3 virulence genes, namely, cps, toxRand VPI, and V. vulnificus strain BTMM7 that did not possess any of the tested virulence genes. A reduction was observed in the life span of worms fed on environmental strain of V. vulnificusBTMM7 rather than on the ordinary laboratory food source, E. coli OP50. Behavioral abnormalities like sluggish movement, lawn avoidance and bagging were also observed in the worms fed with strain BTPS6, but the pharynx and the intestine were intact. The presence of multi drug resistant environmental Vibrio strainsthat constitute a major reservoir of diverse virulence genes are to be dealt with caution as they play a decisive role in pathogenicity and horizontal gene transfer in the marine environments.

Differential induction, isolation, physicochemical and molecular characterization of temperate phages of environmental Vibrio cholerae as evidence of phage mediated Horizontal Gene Transfer

The resurgence of the enteric pathogen Vibrio cholerae, the causative organism of epidemic cholera, remains a major health problem in many developing countries like India. The southern Indian state of Kerala is endemic to cholera. The outbreaks of cholera follow a seasonal pattern in regions of endemicity. Marine aquaculture settings and mangrove environments of Kerala serve as reservoirs for V. cholerae. The non-O1/non-O139 environmental isolates of V. cholerae with incomplete ‘virulence casette’ are to be dealt with caution as they constitute a major reservoir of diverse virulence genes in the marine environment and play a crucial role in pathogenicity and horizontal gene transfer. The genes coding cholera toxin are borne on, and can be infectiously transmitted by CTXΦ, a filamentous lysogenic vibriophages. Temperate phages can provide crucial virulence and fitness factors affecting cell metabolism, bacterial adhesion, colonization, immunity, antibiotic resistance and serum resistance. The present study was an attempt to screen the marine environments like aquafarms and mangroves of coastal areas of Alappuzha and Cochin, Kerala for the presence of lysogenic V. cholerae, to study their pathogenicity and also gene transfer potential. Phenotypic and molecular methods were used for identification of isolates as V. cholerae. The thirty one isolates which were Gram negative, oxidase positive, fermentative, with or without gas production on MOF media and which showed yellow coloured colonies on TCBS (Thiosulfate Citrate Bile salt Sucrose) agar were segregated as vibrios. Twenty two environmental V. cholerae strains of both O1 and non- O1/non-O139 serogroups on induction with mitomycin C showed the presence of lysogenic phages. They produced characteristic turbid plaques in double agar overlay assay using the indicator strain V. cholerae El Tor MAK 757. PCR based molecular typing with primers targeting specific conserved sequences in the bacterial genome, demonstrated genetic diversity among these lysogen containing non-O1 V. cholerae . Polymerase chain reaction was also employed as a rapid screening method to verify the presence of 9 virulence genes namely, ctxA, ctxB, ace, hlyA, toxR, zot,tcpA, ninT and nanH, using gene specific primers. The presence of tcpA gene in ALPVC3 was alarming, as it indicates the possibility of an epidemic by accepting the cholera. Differential induction studies used ΦALPVC3, ΦALPVC11, ΦALPVC12 and ΦEKM14, underlining the possibility of prophage induction in natural ecosystems, due to abiotic factors like antibiotics, pollutants, temperature and UV. The efficiency of induction of prophages varied considerably in response to the different induction agents. The growth curve of lysogenic V. cholerae used in the study drastically varied in the presence of strong prophage inducers like antibiotics and UV. Bacterial cell lysis was directly proportional to increase in phage number due to induction. Morphological characterization of vibriophages by Transmission Electron Microscopy revealed hexagonal heads for all the four phages. Vibriophage ΦALPVC3 exhibited isometric and contractile tails characteristic of family Myoviridae, while phages ΦALPVC11 and ΦALPVC12 demonstrated the typical hexagonal head and non-contractile tail of family Siphoviridae. ΦEKM14, the podophage was distinguished by short non-contractile tail and icosahedral head. This work demonstrated that environmental parameters can influence the viability and cell adsorption rates of V. cholerae phages. Adsorption studies showed 100% adsorption of ΦALPVC3 ΦALPVC11, ΦALPVC12 and ΦEKM14 after 25, 30, 40 and 35 minutes respectively. Exposure to high temperatures ranging from 50ºC to 100ºC drastically reduced phage viability. The optimum concentration of NaCl required for survival of vibriophages except ΦEKM14 was 0.5 M and that for ΦEKM14 was 1M NaCl. Survival of phage particles was maximum at pH 7-8. V. cholerae is assumed to have existed long before their human host and so the pathogenic clones may have evolved from aquatic forms which later colonized the human intestine by progressive acquisition of genes. This is supported by the fact that the vast majority of V. cholerae strains are still part of the natural aquatic environment. CTXΦ has played a critical role in the evolution of the pathogenicity of V. cholerae as it can transmit the ctxAB gene. The unusual transformation of V. cholerae strains associated with epidemics and the emergence of V. cholera O139 demonstrates the evolutionary success of the organism in attaining greater fitness. Genetic changes in pathogenic V. cholerae constitute a natural process for developing immunity within an endemically infected population. The alternative hosts and lysogenic environmental V. cholerae strains may potentially act as cofactors in promoting cholera phage ‘‘blooms’’ within aquatic environments, thereby influencing transmission of phage sensitive, pathogenic V. cholerae strains by aquatic vehicles. Differential induction of the phages is a clear indication of the impact of environmental pollution and global changes on phage induction. The development of molecular biology techniques offered an accessible gateway for investigating the molecular events leading to genetic diversity in the marine environment. Using nucleic acids as targets, the methods of fingerprinting like ERIC PCR and BOX PCR, revealed that the marine environment harbours potentially pathogenic group of bacteria with genetic diversity. The distribution of virulence associated genes in the environmental isolates of V. cholerae provides tangible material for further investigation. Nucleotide and protein sequence analysis alongwith protein structure prediction aids in better understanding of the variation inalleles of same gene in different ecological niche and its impact on the protein structure for attaining greater fitness of pathogens. The evidences of the co-evolution of virulence genes in toxigenic V. cholerae O1 from different lineages of environmental non-O1 strains is alarming. Transduction studies would indicate that the phenomenon of acquisition of these virulence genes by lateral gene transfer, although rare, is not quite uncommon amongst non-O1/non-O139 V. cholerae and it has a key role in diversification. All these considerations justify the need for an integrated approach towards the development of an effective surveillance system to monitor evolution of V. cholerae strains with epidemic potential. Results presented in this study, if considered together with the mechanism proposed as above, would strongly suggest that the bacteriophage also intervenes as a variable in shaping the cholera bacterium, which cannot be ignored and hinting at imminent future epidemics.

Optimization-Based Robotic Manipulation for Safe Interaction with Human Operators

This thesis investigates a method for human-robot interaction (HRI) in order to uphold productivity of industrial robots like minimization of the shortest operation time, while ensuring human safety like collision avoidance. For solving such problems an online motion planning approach for robotic manipulators with HRI has been proposed. The approach is based on model predictive control (MPC) with embedded mixed integer programming. The planning strategies of the robotic manipulators mainly considered in the thesis are directly performed in the workspace for easy obstacle representation. The non-convex optimization problem is approximated by a mixed-integer program (MIP). It is further effectively reformulated such that the number of binary variables and the number of feasible integer solutions are drastically decreased. Safety-relevant regions, which are potentially occupied by the human operators, can be generated online by a proposed method based on hidden Markov models. In contrast to previous approaches, which derive predictions based on probability density functions in the form of single points, such as most likely or expected human positions, the proposed method computes safety-relevant subsets of the workspace as a region which is possibly occupied by the human at future instances of time. The method is further enhanced by combining reachability analysis to increase the prediction accuracy. These safety-relevant regions can subsequently serve as safety constraints when the motion is planned by optimization. This way one arrives at motion plans that are safe, i.e. plans that avoid collision with a probability not less than a predefined threshold. The developed methods have been successfully applied to a developed demonstrator, where an industrial robot works in the same space as a human operator. The task of the industrial robot is to drive its end-effector according to a nominal sequence of grippingmotion-releasing operations while no collision with a human arm occurs.

Estimating the motion of an underwater robot from a monocular image sequence

When underwater vehicles perform navigation close to the ocean floor, computer vision techniques can be applied to obtain quite accurate motion estimates. The most crucial step in the vision-based estimation of the vehicle motion consists on detecting matchings between image pairs. Here we propose the extensive use of texture analysis as a tool to ameliorate the correspondence problem in underwater images. Once a robust set of correspondences has been found, the three-dimensional motion of the vehicle can be computed with respect to the bed of the sea. Finally, motion estimates allow the construction of a map that could aid to the navigation of the robot

Detection of matchings in a sequence of underwater images through texture analysis

This paper presents an approach to ameliorate the reliability of the correspondence points relating two consecutive images of a sequence. The images are especially difficult to handle, since they have been acquired by a camera looking at the sea floor while carried by an underwater robot. Underwater images are usually difficult to process due to light absorption, changing image radiance and lack of well-defined features. A new approach based on gray-level region matching and selective texture analysis significantly improves the matching reliability