Effect of human TGF-beta on the gene expression profile of Schistosoma mansoni adult worms

Schistosoma mansoni is responsible for schistosomiasis, a parasitic disease that affects 200 million people worldwide. Molecular mechanisms of host-parasite interaction are complex and involve a crosstalk between host signals and parasite receptors. TGF-beta signaling pathway has been shown to play an important role in S. mansoni development and embryogenesis. In particular human (h) TGF-beta has been shown to bind to a S. mansoni receptor, transduce a signal that regulates the expression of a schistosome target gene. Here we describe 381 parasite genes whose expression levels are affected by in vitro treatment with hTGF-beta. Among these differentially expressed genes we highlight genes related to morphology, development and cell cycle that could be players of cytokine effects on the parasite. We confirm by qPCR the expression changes detected with microarrays for 5 out of 7 selected genes. We also highlight a set of non-coding RNAs transcribed from the same loci of protein-coding genes that are differentially expressed upon hTCF-beta treatment. These datasets offer potential targets to be explored in order to understand the molecular mechanisms behind the possible role of hTGF-beta effects on parasite biology. (C) 2012 Elsevier B.V. All rights reserved.

Haplotypes of susceptibility to chronic periodontitis in the Interleukin 8 gene do not influence protein level in the gingival crevicular fluid

Objective: Previously, we identified that the ATC/TTC haplotype formed by polymorphisms in the Interleukin-(IL)8 gene conferred susceptibility to chronic periodontitis (CP). The aim of the study was to investigate whether the IL8 haplotype ATC/TTC was associated with the volume of gingival crevicular fluid (GCF), the concentration of interleukin IL-8 in the GCF, as well as periodontal conditions in patients with CP in comparison to controls without CP. Methods: Seventy-nine individuals (CP: n = 41, controls: n = 38) were grouped according to the presence (susceptible for CP) or absence (not susceptible for CP) of the IL8 ATC/TTC haplotype. After periodontal clinical evaluation, they were subdivided by the presence or absence of CP. GCF was collected from each patient and the IL-8 levels were determined by ELISA. The GCF volume of each subject was measured by means of a calibrated electronic device. Comparisons of means between carriers and non-carriers of the ATC/TTC haplotype were evaluated using the Mann-Whitney test. Linear regression and stepwise linear regression analysis were used to analyse the association of the GCF volume with potential covariates and their contribution for the phenotype. Results: We did not find significant differences of both periodontal conditions and IL-8 concentration in the GCF of patients with the presence or absence of the IL8 ATC/TTC haplotype. However, the GCF volume was significantly higher amongst the patients affected by CP that are absent for the IL8 ATC/TTC haplotype. In addition, linear regression analysis showed a statistically significant association between GCF volume and CP, IL8 haplotype ATC/TTC and IL-8 concentration. Conclusions: The IL8 haplotype of susceptibility to CP was neither associated with IL-8 cytokine levels nor with clinical periodontal parameters. Also, CP, IL8 haplotype and IL-8 concentration showed a positive association with the GCF volume levels in the studied patients. (c) 2012 Published by Elsevier Ltd.

Impact of complex NOTCH1 mutations on survival in paediatric T-cell leukaemia

Background: Molecular alterations occur frequently in T-ALL and the potential impact of those abnormalities on outcome is still controversial. The current study aimed to test whether NOTCH1 mutations and additional molecular abnormalities would impact T-ALL outcome in a series of 138 T-ALL paediatric cases. Methods: T-ALL subtypes, status of SIL-TAL1 fusion, ectopic expression of TLX3, and mutations in FBXW7, KRAS, PTEN and NOTCH1 were assessed as overall survival (OS) and event-free survival (EFS) prognostic factors. OS and EFS were determined using the Kaplan-Meier method and compared using the log-rank test. Results: The frequencies of mutations were 43.5% for NOTCH1, while FBXW7, KRAS and PTEN exhibited frequencies of 19.1%, 9.5% and 9.4%, respectively. In 78.3% of cases, the coexistence of NOTCH1 mutations and other molecular alterations was observed. In multivariate analysis no statistical association was revealed between NOTCH1 mutations and any other variable analyzed. The mean length of the follow-up was 68.4 months and the OS was 50.7%. SIL-TAL1 was identified as an adverse prognostic factor. NOTCH1 mutation status was not associated with outcome, while the presence of NOTCH1 complex mutations (indels) were associated with a longer overall survival (p = 0.031) than point mutations. Conclusion: NOTCH1 mutations alone or in combination with FBXW7 did not impact T-ALL prognosis. Nevertheless, complex NOTCH1 mutations appear to have a positive impact on OS and the SIL-TAL1 fusion was validated as a negative prognostic marker in our series of T-ALL.

Development of a recombinant fusion protein based on the dynein light chain LC8 for non-viral gene delivery

The low efficiency of gene transfer is a recurrent problem in DNA vaccine development and gene therapy studies using non-viral vectors such as plasmid DNA (pDNA). This is mainly due to the fact that during their traffic to the target cell's nuclei, plasmid vectors must overcome a series of physical, enzymatic and diffusional barriers. The main objective of this work is the development of recombinant proteins specifically designed for pDNA delivery, which take advantage of molecular motors like dynein, for the transport of cargos from the periphery to the centrosome of mammalian cells. A DNA binding sequence was fused to the N-terminus of the recombinant human dynein light chain LC8. Expression studies indicated that the fusion protein was correctly expressed in soluble form using E. coli BL21(DE3) strain. As expected, gel permeation assays found the purified protein mainly present as dimers, the functional oligomeric state of LC8. Gel retardation assays and atomic force microscopy proved the ability of the fusion protein to interact and condense pDNA. Zeta potential measurements indicated that LC8 with DNA binding domain (LD4) has an enhanced capacity to interact and condense pDNA, generating positively charged complexes. Transfection of cultured HeLa cells confirmed the ability of the LD4 to facilitate pDNA uptake and indicate the involvement of the retrograde transport in the intracellular trafficking of pDNA: LD4 complexes. Finally, cytotoxicity studies demonstrated a very low toxicity of the fusion protein vector, indicating the potential for in vivo applications. The study presented here is part of an effort to develop new modular shuttle proteins able to take advantage of strategies used by viruses to infect mammalian cells, aiming to provide new tools for gene therapy and DNA vaccination studies. (C) 2012 Elsevier B.V. All rights reserved.

A complex chromosomal rearrangement involving chromosomes 2, 5, and X in autism spectrum disorder

Here, we describe a female patient with autism spectrum disorder and dysmorphic features that harbors a complex genetic alteration, involving a de novo balanced translocation t(2;X)(q11;q24), a 5q11 segmental trisomy and a maternally inherited isodisomy on chromosome 5. All the possibly damaging genetic effects of such alterations are discussed. In light of recent findings on ASD genetic causes, the hypothesis that all these alterations might be acting in orchestration and contributing to the phenotype is also considered. (C) 2012 Wiley Periodicals, Inc.

Identification of unannotated exons of low abundance transcripts in Drosophila melanogaster and cloning of a new serine protease gene upregulated upon injury

Abstract Background The sequencing of the D.melanogaster genome revealed an unexpected small number of genes (~ 14,000) indicating that mechanisms acting on generation of transcript diversity must have played a major role in the evolution of complex metazoans. Among the most extensively used mechanisms that accounts for this diversity is alternative splicing. It is estimated that over 40% of Drosophila protein-coding genes contain one or more alternative exons. A recent transcription map of the Drosophila embryogenesis indicates that 30% of the transcribed regions are unannotated, and that 1/3 of this is estimated as missed or alternative exons of previously characterized protein-coding genes. Therefore, the identification of the variety of expressed transcripts depends on experimental data for its final validation and is continuously being performed using different approaches. We applied the Open Reading Frame Expressed Sequence Tags (ORESTES) methodology, which is capable of generating cDNA data from the central portion of rare transcripts, in order to investigate the presence of hitherto unnanotated regions of Drosophila transcriptome. Results Bioinformatic analysis of 1,303 Drosophila ORESTES clusters identified 68 sequences derived from unannotated regions in the current Drosophila genome version (4.3). Of these, a set of 38 was analysed by polyA+ northern blot hybridization, validating 17 (50%) new exons of low abundance transcripts. For one of these ESTs, we obtained the cDNA encompassing the complete coding sequence of a new serine protease, named SP212. The SP212 gene is part of a serine protease gene cluster located in the chromosome region 88A12-B1. This cluster includes the predicted genes CG9631, CG9649 and CG31326, which were previously identified as up-regulated after immune challenges in genomic-scale microarray analysis. In agreement with the proposal that this locus is co-regulated in response to microorganisms infection, we show here that SP212 is also up-regulated upon injury. Conclusion Using the ORESTES methodology we identified 17 novel exons from low abundance Drosophila transcripts, and through a PCR approach the complete CDS of one of these transcripts was defined. Our results show that the computational identification and manual inspection are not sufficient to annotate a genome in the absence of experimentally derived data.

Gene by environment QTL mapping through multiple trait analyses in blood pressure salt-sensitivity: identification of a novel QTL in rat chromosome 5

Background The genetic mechanisms underlying interindividual blood pressure variation reflect the complex interplay of both genetic and environmental variables. The current standard statistical methods for detecting genes involved in the regulation mechanisms of complex traits are based on univariate analysis. Few studies have focused on the search for and understanding of quantitative trait loci responsible for gene × environmental interactions or multiple trait analysis. Composite interval mapping has been extended to multiple traits and may be an interesting approach to such a problem. Methods We used multiple-trait analysis for quantitative trait locus mapping of loci having different effects on systolic blood pressure with NaCl exposure. Animals studied were 188 rats, the progenies of an F2 rat intercross between the hypertensive and normotensive strain, genotyped in 179 polymorphic markers across the rat genome. To accommodate the correlational structure from measurements taken in the same animals, we applied univariate and multivariate strategies for analyzing the data. Results We detected a new quantitative train locus on a region close to marker R589 in chromosome 5 of the rat genome, not previously identified through serial analysis of individual traits. In addition, we were able to justify analytically the parametric restrictions in terms of regression coefficients responsible for the gain in precision with the adopted analytical approach. Conclusion Future work should focus on fine mapping and the identification of the causative variant responsible for this quantitative trait locus signal. The multivariable strategy might be valuable in the study of genetic determinants of interindividual variation of antihypertensive drug effectiveness.

Crystal structures of Leptospira interrogans FAD-containing ferredoxin-NADP+ reductase and its complex with NADP+

Abstract Background Ferredoxin-NADP(H) reductases (FNRs) are flavoenzymes that catalyze the electron transfer between NADP(H) and the proteins ferredoxin or flavodoxin. A number of structural features distinguish plant and bacterial FNRs, one of which is the mode of the cofactor FAD binding. Leptospira interrogans is a spirochaete parasitic bacterium capable of infecting humans and mammals in general. Leptospira interrogans FNR (LepFNR) displays low sequence identity with plant (34% with Zea mays) and bacterial (31% with Escherichia coli) FNRs. However, LepFNR contains all consensus sequences that define the plastidic class FNRs. Results The crystal structures of the FAD-containing LepFNR and the complex of the enzyme with NADP+, were solved and compared to known FNRs. The comparison reveals significant structural similarities of the enzyme with the plastidic type FNRs and differences with the bacterial enzymes. Our small angle X-ray scattering experiments show that LepFNR is a monomeric enzyme. Moreover, our biochemical data demonstrate that the LepFNR has an enzymatic activity similar to those reported for the plastidic enzymes and that is significantly different from bacterial flavoenzymes, which display lower turnover rates. Conclusion LepFNR is the first plastidic type FNR found in bacteria and, despite of its low sequence similarity with plastidic FNRs still displays high catalytic turnover rates. The typical structural and biochemical characteristics of plant FNRs unveiled for LepFNR support a notion of a putative lateral gene transfer which presumably offers Leptospira interrogans evolutionary advantages. The wealth of structural information about LepFNR provides a molecular basis for advanced drugs developments against leptospirosis.

Modularity and evolutionary constraints in a baculovirus gene regulatory network

Abstract Background The structure of regulatory networks remains an open question in our understanding of complex biological systems. Interactions during complete viral life cycles present unique opportunities to understand how host-parasite network take shape and behave. The Anticarsia gemmatalis multiple nucleopolyhedrovirus (AgMNPV) is a large double-stranded DNA virus, whose genome may encode for 152 open reading frames (ORFs). Here we present the analysis of the ordered cascade of the AgMNPV gene expression. Results We observed an earlier onset of the expression than previously reported for other baculoviruses, especially for genes involved in DNA replication. Most ORFs were expressed at higher levels in a more permissive host cell line. Genes with more than one copy in the genome had distinct expression profiles, which could indicate the acquisition of new functionalities. The transcription gene regulatory network (GRN) for 149 ORFs had a modular topology comprising five communities of highly interconnected nodes that separated key genes that are functionally related on different communities, possibly maximizing redundancy and GRN robustness by compartmentalization of important functions. Core conserved functions showed expression synchronicity, distinct GRN features and significantly less genetic diversity, consistent with evolutionary constraints imposed in key elements of biological systems. This reduced genetic diversity also had a positive correlation with the importance of the gene in our estimated GRN, supporting a relationship between phylogenetic data of baculovirus genes and network features inferred from expression data. We also observed that gene arrangement in overlapping transcripts was conserved among related baculoviruses, suggesting a principle of genome organization. Conclusions Albeit with a reduced number of nodes (149), the AgMNPV GRN had a topology and key characteristics similar to those observed in complex cellular organisms, which indicates that modularity may be a general feature of biological gene regulatory networks.

Phylogeography of the neotropical Anopheles triannulatus complex (Diptera: Culicidae) supports deep structure and complex patterns

Abstract Background The molecular phylogenetic relationships and population structure of the species of the Anopheles triannulatus complex: Anopheles triannulatus s.s., Anopheles halophylus and the putative species Anopheles triannulatus C were investigated. Methods The mitochondrial COI gene, the nuclear white gene and rDNA ITS2 of samples that include the known geographic distribution of these taxa were analyzed. Phylogenetic analyses were performed using Bayesian inference, Maximum parsimony and Maximum likelihood approaches. Results Each data set analyzed septely yielded a different topology but none provided evidence for the seption of An. halophylus and An. triannulatus C, consistent with the hypothesis that the two are undergoing incipient speciation. The phylogenetic analyses of the white gene found three main clades, whereas the statistical parsimony network detected only a single metapopulation of Anopheles triannulatus s.l. Seven COI lineages were detected by phylogenetic and network analysis. In contrast, the network, but not the phylogenetic analyses, strongly supported three ITS2 groups. Combined data analyses provided the best resolution of the trees, with two major clades, Amazonian (clade I) and trans-Andean + Amazon Delta (clade II). Clade I consists of multiple subclades: An. halophylus + An. triannulatus C; trans-Andean Venezuela; central Amazonia + central Bolivia; Atlantic coastal lowland; and Amazon delta. Clade II includes three subclades: Panama; cis-Andean Colombia; and cis-Venezuela. The Amazon delta specimens are in both clades, likely indicating local sympatry. Spatial and molecular variance analyses detected nine groups, corroborating some of subclades obtained in the combined data analysis. Conclusion Combination of the three molecular markers provided the best resolution for differentiation within An. triannulatus s.s. and An. halophylus and C. The latest two species seem to be very closely related and the analyses performed were not conclusive regarding species differentiation. Further studies including new molecular markers would be desirable to solve this species status question. Besides, results of the study indicate a trans-Andean origin for An. triannulatus s.l. The potential implications for malaria epidemiology remain to be investigated.

Behavioral meaningful opioidergic stimulation activates kappa receptor gene expression

The periaqueductal gray (PAG) has been reported to be a location for opioid regulation of pain and a potential site for behavioral selection in females. Opioid-mediated behavioral and physiological responses differ according to the activity of opioid receptor subtypes. The present study investigated the effects of the peripheral injection of the kappa-opioid receptor agonist U69593 into the dorsal subcutaneous region of animals on maternal behavior and on Oprk1 gene activity in the PAG of female rats. Female Wistar rats weighing 200-250 g at the beginning of the study were randomly divided into 2 groups for maternal behavior and gene expression experiments. On day 5, pups were removed at 7:00 am and placed in another home cage that was distant from their mother. Thirty minutes after removing the pups, the dams were treated with U69593 (0.15 mg/kg, sc) or 0.9% saline (up to 1 mL/kg) and after 30 min were evaluated in the maternal behavior test. Latencies in seconds for pup retrieval, grouping, crouching, and full maternal behavior were scored. The results showed that U69593 administration inhibited maternal behavior (P < 0.05) because a lower percentage of kappa group dams showed retrieval of first pup, retrieving all pups, grouping, crouching and displaying full maternal behavior compared to the saline group. Opioid gene expression was evaluated using real-time reverse-transcription polymerase chain reaction (RT-PCR). A single injection of U69593 increased Oprk1 PAG expression in both virgin (P < 0.05) and lactating female rats (P < 0.01), with no significant effect on Oprm1 or Oprd1 gene activity. Thus, the expression of kappa-opioid receptors in the PAG may be modulated by single opioid receptor stimulation and behavioral meaningful opioidergic transmission in the adult female might occur simultaneously to specific changes in gene expression of kappa-opioid receptor subtype. This is yet another alert for the complex role of the opioid system in female reproduction

Cryptic species in the cosmopolitan Bugulan eritina complex (Bryozoa,Cheilostomata)

Previous analyses of the mitochondrial gene cytochrome c oxidase subunit 1 (COI) and γ-proteobacterial endosymbiont diversity have suggested that the marine bryozoan Bugula neritina is a complex of three cryptic species, namely Types S, D and N. Types D and N were previously reported to have restricted distributions along California (western USA) and Delaware and Connecticut (eastern USA), respectively, whereas Type S is considered widespread in tropical, subtropical and temperate regions due to anthropogenic transport. Here, Bayesian species delimitation analysis of a data set composed of two mitochondrial (COI and large ribosomal RNA subunit [16S]) and two nuclear genes (dynein light chain roadblock type-2 protein [DYN] and voltage-dependent anion-selective channel protein [VDAC]) demonstrated that Types S, D and N correspond to three biological species. This finding was significantly supported, in spite of the combinations of priors applied for ancestral population size and root age. Furthermore, COI sequences were used to assess the introduction patterns of the cosmopolitan Type S species. Two COI haplotypes of Type S (S1a and S1d) were found occurring at a global scale. Mantel tests showed correlation between these haplotypes and local sea surface temperature tolerance. Accordingly, the distributions of Type S haplotypes may reflect intraspecific temperature tolerance variation, in addition to the role of introduction vectors. Finally, we show that the Type N may also have been introduced widely, as this species was found for the first time in Central California and north-eastern Australia.

Systematics of the Oswaldoi Complex (Anopheles, Nyssorhynchus) in South America

Abstract Background Effective malaria control relies on accurate identification of those Anopheles mosquitoes responsible for the transmission of Plasmodium parasites. Anopheles oswaldoi s.l. has been incriminated as a malaria vector in Colombia and some localities in Brazil, but not ubiquitously throughout its Neotropical range. This evidence together with variable morphological characters and genetic differences supports that An. oswaldoi s.l. compromises a species complex. The recent fully integrated redescription of An. oswaldoi s.s. provides a solid taxonomic foundation from which to molecularly determine other members of the complex. Methods DNA sequences of the Second Internal Transcribed Spacer (ITS2 - rDNA) (n = 192) and the barcoding region of the Cytochrome Oxidase I gene (COI - mtDNA) (n = 110) were generated from 255 specimens of An. oswaldoi s.l. from 33 localities: Brazil (8 localities, including the lectotype series of An. oswaldoi), Ecuador (4), Colombia (17), Trinidad and Tobago (1), and Peru (3). COI sequences were analyzed employing the Kimura-two-parameter model (K2P), Bayesian analysis (MrBayes), Mixed Yule-Coalescent model (MYC, for delimitation of clusters) and TCS genealogies. Results Separate and combined analysis of the COI and ITS2 data sets unequivocally supported four separate species: two previously determined (An. oswaldoi s.s. and An. oswaldoi B) and two newly designated species in the Oswaldoi Complex (An. oswaldoi A and An. sp. nr. konderi). The COI intra- and inter-specific genetic distances for the four taxa were non-overlapping, averaging 0.012 (0.007 to 0.020) and 0.052 (0.038 to 0.064), respectively. The concurring four clusters delineated by MrBayes and MYC, and four independent TCS networks, strongly confirmed their separate species status. In addition, An. konderi of Sallum should be regarded as unique with respect to the above. Despite initially being included as an outgroup taxon, this species falls well within the examined taxa, suggesting a combined analysis of these taxa would be most appropriate. Conclusions: Through novel data and retrospective comparison of available COI and ITS2 DNA sequences, evidence is shown to support the separate species status of An. oswaldoi s.s., An. oswaldoi A and An. oswaldoi B, and at least two species in the closely related An. konderi complex (An. sp. nr. konderi, An. konderi of Sallum). Although An. oswaldoi s.s. has never been implicated in malaria transmission, An. oswaldoi B is a confirmed vector and the new species An. oswaldoi A and An. sp. nr. konderi are circumstantially implicated, most likely acting as secondary vectors.

Exome analysis of HIV patients submitted to dendritic cells therapeutic vaccine reveals an association of CNOT1 gene with response to the treatment

INTRODUCTION: With the aim of searching genetic factors associated with the response to an immune treatment based on autologous monocyte-derived dendritic cells pulsed with autologous inactivated HIV, we performed exome analysis by screening more than 240,000 putative functional exonic variants in 18 HIV-positive Brazilian patients that underwent the immune treatment. METHODS: Exome analysis has been performed using the ILLUMINA Infinium HumanExome BeadChip. zCall algorithm allowed us to recall rare variants. Quality control and SNP-centred analysis were done with GenABEL R package. An in-house implementation of the Wang method permitted gene-centred analysis. RESULTS: CCR4-NOT transcription complex, subunit 1 (CNOT1) gene (16q21), showed the strongest association with the modification of the response to the therapeutic vaccine (p=0.00075). CNOT1 SNP rs7188697 A/G was significantly associated with DC treatment response. The presence of a G allele indicated poor response to the therapeutic vaccine (p=0.0031; OR=33.00; CI=1.74-624.66), and the SNP behaved in a dominant model (A/A vs. A/G+G/G p=0.0009; OR=107.66; 95% CI=3.85-3013.31), being the A/G genotype present only in weak/transient responders, conferring susceptibility to poor response to the immune treatment. DISCUSSION: CNOT1 is known to be involved in the control of mRNA deadenylation and mRNA decay. Moreover, CNOT1 has been recently described as being involved in the regulation of inflammatory processes mediated by tristetraprolin (TTP). The TTP-CCR4-NOT complex (CNOT1 in the CCR4-NOT complex is the binding site for TTP) has been reported as interfering with HIV replication, through post-transcriptional control. Therefore, we can hypothesize that genetic variation occurring in the CNOT1 gene could impair the TTP-CCR4-NOT complex, thus interfering with HIV replication and/or host immune response. CONCLUSIONS: Being aware that our findings are exclusive to the 18 patients studied with a need for replication, and that the genetic variant of CNOT1 gene, localized at intron 3, has no known functional effect, we propose a novel potential candidate locus for the modulation of the response to the immune treatment, and open a discussion on the necessity to consider the host genome as another potential variant to be evaluated when designing an immune therapy study

New approaches to open problems in gene expression microarray data

In the past decade, the advent of efficient genome sequencing tools and high-throughput experimental biotechnology has lead to enormous progress in the life science. Among the most important innovations is the microarray tecnology. It allows to quantify the expression for thousands of genes simultaneously by measurin the hybridization from a tissue of interest to probes on a small glass or plastic slide. The characteristics of these data include a fair amount of random noise, a predictor dimension in the thousand, and a sample noise in the dozens. One of the most exciting areas to which microarray technology has been applied is the challenge of deciphering complex disease such as cancer. In these studies, samples are taken from two or more groups of individuals with heterogeneous phenotypes, pathologies, or clinical outcomes. these samples are hybridized to microarrays in an effort to find a small number of genes which are strongly correlated with the group of individuals. Eventhough today methods to analyse the data are welle developed and close to reach a standard organization (through the effort of preposed International project like Microarray Gene Expression Data -MGED- Society [1]) it is not unfrequant to stumble in a clinician's question that do not have a compelling statistical method that could permit to answer it.The contribution of this dissertation in deciphering disease regards the development of new approaches aiming at handle open problems posed by clinicians in handle specific experimental designs. In Chapter 1 starting from a biological necessary introduction, we revise the microarray tecnologies and all the important steps that involve an experiment from the production of the array, to the quality controls ending with preprocessing steps that will be used into the data analysis in the rest of the dissertation. While in Chapter 2 a critical review of standard analysis methods are provided stressing most of problems that In Chapter 3 is introduced a method to adress the issue of unbalanced design of miacroarray experiments. In microarray experiments, experimental design is a crucial starting-point for obtaining reasonable results. In a two-class problem, an equal or similar number of samples it should be collected between the two classes. However in some cases, e.g. rare pathologies, the approach to be taken is less evident. We propose to address this issue by applying a modified version of SAM [2]. MultiSAM consists in a reiterated application of a SAM analysis, comparing the less populated class (LPC) with 1,000 random samplings of the same size from the more populated class (MPC) A list of the differentially expressed genes is generated for each SAM application. After 1,000 reiterations, each single probe given a "score" ranging from 0 to 1,000 based on its recurrence in the 1,000 lists as differentially expressed. The performance of MultiSAM was compared to the performance of SAM and LIMMA [3] over two simulated data sets via beta and exponential distribution. The results of all three algorithms over low- noise data sets seems acceptable However, on a real unbalanced two-channel data set reagardin Chronic Lymphocitic Leukemia, LIMMA finds no significant probe, SAM finds 23 significantly changed probes but cannot separate the two classes, while MultiSAM finds 122 probes with score >300 and separates the data into two clusters by hierarchical clustering. We also report extra-assay validation in terms of differentially expressed genes Although standard algorithms perform well over low-noise simulated data sets, multi-SAM seems to be the only one able to reveal subtle differences in gene expression profiles on real unbalanced data. In Chapter 4 a method to adress similarities evaluation in a three-class prblem by means of Relevance Vector Machine [4] is described. In fact, looking at microarray data in a prognostic and diagnostic clinical framework, not only differences could have a crucial role. In some cases similarities can give useful and, sometimes even more, important information. The goal, given three classes, could be to establish, with a certain level of confidence, if the third one is similar to the first or the second one. In this work we show that Relevance Vector Machine (RVM) [2] could be a possible solutions to the limitation of standard supervised classification. In fact, RVM offers many advantages compared, for example, with his well-known precursor (Support Vector Machine - SVM [3]). Among these advantages, the estimate of posterior probability of class membership represents a key feature to address the similarity issue. This is a highly important, but often overlooked, option of any practical pattern recognition system. We focused on Tumor-Grade-three-class problem, so we have 67 samples of grade I (G1), 54 samples of grade 3 (G3) and 100 samples of grade 2 (G2). The goal is to find a model able to separate G1 from G3, then evaluate the third class G2 as test-set to obtain the probability for samples of G2 to be member of class G1 or class G3. The analysis showed that breast cancer samples of grade II have a molecular profile more similar to breast cancer samples of grade I. Looking at the literature this result have been guessed, but no measure of significance was gived before.