Transcriptional Response in a Sepsis Mouse Model Reflects Transcriptional Response in Sepsis Patients.

Mortality due to sepsis remains unacceptably high, especially for septic shock patients. Murine models have been used to better understand pathophysiology mechanisms. However, the mouse model is still under debate. Herein we investigated the transcriptional response of mice injected with lipopolysaccharide (LPS) and compared it to either human cells stimulated in vitro with LPS or to the blood cells of septic patients. We identified a molecular signature composed of 2331 genes with an FDR median of 0%. This molecular signature is highly enriched in regulated genes in peritoneal macrophages stimulated with LPS. There is significant enrichment in several inflammatory signaling pathways, and in disease terms, such as pneumonia, sepsis, systemic inflammatory response syndrome, severe sepsis, an inflammatory disorder, immune suppression, and septic shock. A significant overlap between the genes upregulated in mouse and human cells stimulated with LPS has been demonstrated. Finally, genes upregulated in mouse cells stimulated with LPS are enriched in genes upregulated in human cells stimulated in vitro and in septic patients, who are at high risk of death. Our results support the hypothesis of common molecular and cellular mechanisms between mouse and human sepsis.

Prenatal maternal vitamin D deficiency sex-dependently programs adipose tissue metabolism and energy homeostasis in offspring.

In utero environment is crucial to ensure normal development of the fetus and to program metabolic health throughout the life. Beside macronutrients, the role of micronutrients, including vitamin D, begins to be explore. The aim of this study was to decipher the impact of maternal vitamin D deficiency (VDD), in normal and high-fat (HF) diet context, on adipose tissue metabolism and energy homeostasis in offspring, considering sex-specific responses. Body weight, energy expenditure, and spontaneous activity was differential impacted in juvenile male and female offspring born from VDD mice. In adulthood, a HF diet combined with maternal VDD disrupted glucose homeostasis and adiposity in male offspring but not in females. Such phenotypes were associated to different transcriptomic profiles in adipose tissue, which could be related to differential modulation of plasma 17ß-estradiol concentrations. Thus, maternal VDD sex-dependently modulated metabolic fate of the offspring, especially when associated with HF diet in adulthood.

Divergent transcription is associated with promoters of transcriptional regulators.

BACKGROUND: Divergent transcription is a wide-spread phenomenon in mammals. For instance, short bidirectional transcripts are a hallmark of active promoters, while longer transcripts can be detected antisense from active genes in conditions where the RNA degradation machinery is inhibited. Moreover, many described long non-coding RNAs (lncRNAs) are transcribed antisense from coding gene promoters. However, the general significance of divergent lncRNA/mRNA gene pair transcription is still poorly understood. Here, we used strand-specific RNA-seq with high sequencing depth to thoroughly identify antisense transcripts from coding gene promoters in primary mouse tissues. RESULTS: We found that a substantial fraction of coding-gene promoters sustain divergent transcription of long non-coding RNA (lncRNA)/mRNA gene pairs. Strikingly, upstream antisense transcription is significantly associated with genes related to transcriptional regulation and development. Their promoters share several characteristics with those of transcriptional developmental genes, including very large CpG islands, high degree of conservation and epigenetic regulation in ES cells. In-depth analysis revealed a unique GC skew profile at these promoter regions, while the associated coding genes were found to have large first exons, two genomic features that might enforce bidirectional transcription. Finally, genes associated with antisense transcription harbor specific H3K79me2 epigenetic marking and RNA polymerase II enrichment profiles linked to an intensified rate of early transcriptional elongation. CONCLUSIONS: We concluded that promoters of a class of transcription regulators are characterized by a specialized transcriptional control mechanism, which is directly coupled to relaxed bidirectional transcription.

In Vivo Model for Isolating Epithelial Cells of the Anorectal Transition Zone.

Different epithelia line the body and organs and form a continuous lining of cells. The junction of two different types of epithelia represents a special region called transition zone (TZ). TZ are small areas found in numerous places in the body such as between the esophagus and the stomach, in the cervix, in the eye, and between the anal canal and the rectum. These zones are associated with diverse pathologies such as cancers; however, the cellular and molecular mechanisms involved in tumor progression are poorly investigated. We recently characterized the role of anorectal TZ cells during homeostasis and after injury using an in vivo (lineage tracing) approach. To follow TZ cells, we previously developed a mouse model of lineage tracing using cytokeratin 17 (Krt17) as a promoter and GFP as a reporter. Krt17 is expressed by TZ but also by anal glands located below the TZ in the stroma that can interfere with TZ cell population isolation and analysis afterward. In this chapter, we provide a new dissection method to remove specifically anal glands without affecting anorectal TZ cells. This protocol allows the specific dissection and isolation of anal canal, TZ, and rectum epithelia.

Defining Anorectal Transition Zone Heterogeneity Using Single-Cell RNA Sequencing.

Special regions called transition zones (TZs) are found at numerous places in the body. TZs represent the junction between two different types of epithelia and are located between the esophagus and the stomach, in the cervix, in the eye, and between the anal canal and the rectum. TZ is a heterogeneous population, and the detailed characterization of its populations requires an analysis at the single-cell level. In this chapter, we provide a protocol to do single-cell RNA sequencing primary analysis of anal canal, TZ, and rectum epithelia.

The Contribution of Multiplexing Single Cell RNA Sequencing in Acute Myeloid Leukemia.

Decades ago, the treatment for acute myeloid leukemia relied on cytarabine and anthracycline. However, advancements in medical research have introduced targeted therapies, initially employing monoclonal antibodies such as ant-CD52 and anti-CD123, and subsequently utilizing specific inhibitors that target molecular mutations like anti-IDH1, IDH2, or FLT3. The challenge lies in determining the role of these therapeutic options, considering the inherent tumor heterogeneity associated with leukemia diagnosis and the clonal drift that this type of tumor can undergo. Targeted drugs necessitate an examination of various therapeutic targets at the individual cell level rather than assessing the entire population. It is crucial to differentiate between the prognostic value and therapeutic potential of a specific molecular target, depending on whether it is found in a terminally differentiated cell with limited proliferative potential or a stem cell with robust capabilities for both proliferation and self-renewal. However, this cell-by-cell analysis is accompanied by several challenges. Firstly, the scientific aspect poses difficulties in comparing different single cell analysis experiments despite efforts to standardize the results through various techniques. Secondly, there are practical obstacles as each individual cell experiment incurs significant financial costs and consumes a substantial amount of time. A viable solution lies in the ability to process multiple samples simultaneously, which is a distinctive feature of the cell hashing technique. In this study, we demonstrate the applicability of the cell hashing technique for analyzing acute myeloid leukemia cells. By comparing it to standard single cell analysis, we establish a strong correlation in various parameters such as quality control, gene expression, and the analysis of leukemic blast markers in patients. Consequently, this technique holds the potential to become an integral part of the biological assessment of acute myeloid leukemia, contributing to the personalized and optimized management of the disease, particularly in the context of employing targeted therapies.

Single-Cell Transcriptome Analysis of Acute Myeloid Leukemia Cells Using Methanol Fixation and Cryopreservation.

INTRODUCTION: The application of single-cell RNA sequencing has greatly improved our understanding of various cellular and molecular mechanisms involved in physiological and pathophysiological processes. However, obtaining living cells for this technique can be difficult under certain conditions. To solve this problem, the methanol fixation method appeared as a promising alternative for routine clinical use. MATERIALS AND METHODS: In this study, we selected two AML samples that had been fixed in methanol for 12-18 months. Once the cells were rehydrated, these samples were subjected to single-cell RNA sequencing. We then compared the results obtained from these samples with those obtained from the same samples cryopreserved in DMSO. RESULTS: We used a previously validated methanol fixation protocol to perform scRNA-seq on DMSO cryopreserved cells and cells fixed in methanol for more than one year. Preliminary results show that methanol fixation induces some genetic and transcriptional modification compared with DMSO cryopreservation but remains a valuable method for single-cell analysis of primary human leukemia cells. CONCLUSIONS: The initial findings from this study highlight certain resemblances in methanol fixation over a 12-month period and cryopreservation with DMSO, along with associated transcriptional level modifications. However, we observed genetic degradation in the fixation condition when extending beyond one year. Despite certain study limitations, it is evident that short-term methanol fixation can be effectively used for leukemia blast samples. Its ease of implementation holds the potential to simplify the integration of this technique into routine clinical practice.

Genome-wide analysis of familial dysautonomia and kinetin target genes with patient olfactory ecto-mesenchymal stem cells.

Familial dysautonomia (FD) is a rare inherited neurodegenerative disorder. The most common mutation is a c.2204+6T>C transition in the 5' splice site (5'ss) of IKBKAP intron 20, which causes a tissue-specific skipping of exon 20, resulting in lower synthesis of IKAP/hELP1 protein. To better understand the specificity of neuron loss in FD, we modeled the molecular mechanisms of IKBKAP mRNA splicing by studying human olfactory ecto-mesenchymal stem cells (hOE-MSCs) derived from FD patient nasal biopsies. We explored how the modulation of IKBKAP mRNA alternative splicing impacts the transcriptome at the genome-wide level. We found that the FD transcriptional signature was highly associated with biological functions related to the development of the nervous system. In addition, we identified target genes of kinetin, a plant cytokinin that corrects IKBKAP mRNA splicing and increases the expression of IKAP/hELP1. We identified this compound as a putative regulator of splicing factors and added new evidence for a sequence-specific correction of splicing. In conclusion, hOE-MSCs isolated from FD patients represent a promising avenue for modeling the altered genetic expression of FD, demonstrating a methodology that can be applied to a host of other genetic disorders to test the therapeutic potential of candidate molecules.

Population transcriptomics: insights from Drosophila simulans, Drosophila sechellia and their hybrids.

Sequence differentiation has been widely studied between populations and species, whereas interest in expression divergence is relatively recent. Using microarrays, we compared four geographically distinct populations of Drosophila simulans and a population of Drosophila sechellia, and interspecific hybrids. We observed few differences between populations, suggesting a slight population structure in D. simulans. This structure was observed in direct population comparisons, as well as in interspecific comparisons (hybrids vs. parents, D. sechellia vs. D. simulans). Expression variance is higher in the French and Zimbabwean populations than in the populations from the ancestral range of D. simulans (Kenya and Seychelles). This suggests a large scale phenomenon of decanalization following the invasion of a new environment. Comparing D. simulans and D. sechellia, we revealed 304 consistently differentially expressed genes, with striking overrepresentation of genes of the cytochrome P450 family, which could be related to their role in detoxification as well as in hormone regulation. We also revealed differences in genes involved in Juvenile hormone and Dopamine differentiation. We finally observed very few differentially expressed genes between hybrids and parental populations, with an overrepresentation of X-linked genes.

An evaluation of the role of gene expression in the prediction and diagnosis of ventilator-associated pneumonia.

BACKGROUND: The SepsiChip project explored transcriptional modulation associated with ventilator-associated pneumonia (VAP) in patients admitted to the intensive care unit for trauma. Genome-wide expression analysis may help to identify potential diagnostic markers for diseases. The current study examined the changes in blood transcriptome during VAP. METHODS: The authors prospectively included 165 trauma patients, and 41 developed VAP. Whole blood samples were collected at admission and at VAP. To predict VAP, the admission samples were compared by microarray in patients who did or did not develop VAP. To identify diagnosis markers, paired samples of 35 patients who developed VAP were analyzed. Using NanoString (Seattle, WA), the results were confirmed in the patients who developed VAP. Trauma patients who did not develop VAP served as controls to eliminate a time effect. RESULTS: The injury severity scores of the patients who did or did not develop VAP were 36 and 29, respectively. No predictive biomarker was identified. For patients who developed VAP, a transcriptional signature was identified between the two sampling times. However, this signature was a generalized pattern related to trauma, independent of the infectious process. Genes involved in the proinflammatory response were down-regulated in the patients who developed VAP, but this difference was not statistically significant. CONCLUSIONS: In contrast to clinical assessment, transcriptional analysis of whole blood samples cannot predict or diagnose VAP in trauma patients. Differentiating infection from inflammation seems challenging.

The contribution of single-cell analysis of acute leukemia in the therapeutic strategy.

After decades during which the treatment of acute myeloblastic leukemia was limited to variations around a skeleton of cytarabine/anthracycline, targeted therapies appeared. These therapies, first based on monoclonal antibodies, also rely on specific inhibitors of various molecular abnormalities. A significant but modest prognosis improvement has been observed thanks to these new treatments that are limited by a high rate of relapse, due to the intrinsic chemo and immune-resistance of leukemia stem cell, together with the acquisition of these resistances by clonal evolution. Relapses are also influenced by the equilibrium between the pro or anti-tumor signals from the bone marrow stromal microenvironment and immune effectors. What should be the place of the targeted therapeutic options in light of the tumor heterogeneity inherent to leukemia and the clonal drift of which this type of tumor is capable? Novel approaches by single cell analysis and next generation sequencing precisely define clonal heterogeneity and evolution, leading to a personalized and time variable adapted treatment. Indeed, the evolution of leukemia, either spontaneous or under therapy selection pressure, is a very complex phenomenon. The model of linear evolution is to be forgotten because single cell analysis of samples at diagnosis and at relapse show that tumor escape to therapy occurs from ancestral as well as terminal clones. The determination by the single cell technique of the trajectories of the different tumor sub-populations allows the identification of clones that accumulate factors of resistance to chemo/immunotherapy ("pan-resistant clones"), making possible to choose the combinatorial agents most likely to eradicate these cells. In addition, the single cell technique identifies the nature of each cell and can analyze, on the same sample, both the tumor cells and their environment. It is thus possible to evaluate the populations of immune effectors (T-lymphocytes, natural killer cells) for the leukemia stress-induced alteration of their functions. Finally, the single cells techniques are an invaluable tool for evaluation of the measurable residual disease since not only able to quantify but also to determine the most appropriate treatment according to the sensitivity profile to immuno-chemotherapy of remaining leukemic cells.

Response of Downy Oak (Quercus pubescens Willd.) to Climate Change: Transcriptome Assembly, Differential Gene Analysis and Targeted Metabolomics.

Global change scenarios in the Mediterranean basin predict a precipitation reduction within the coming hundred years. Therefore, increased drought will affect forests both in terms of adaptive ecology and ecosystemic services. However, how vegetation might adapt to drought is poorly understood. In this report, four years of climate change was simulated by excluding 35% of precipitation above a downy oak forest. RNASeq data allowed us to assemble a genome-guided transcriptome. This led to the identification of differentially expressed features, which was supported by the characterization of target metabolites using a metabolomics approach. We provided 2.5 Tb of RNASeq data and the assembly of the first genome guided transcriptome of Quercus pubescens. Up to 5724 differentially expressed transcripts were obtained; 42 involved in plant response to drought. Transcript set enrichment analysis showed that drought induces an increase in oxidative pressure that is mitigated by the upregulation of ubiquitin-like protein protease, ferrochelatase, oxaloacetate decarboxylase and oxo-acid-lyase activities. Furthermore, the downregulation of auxin biosynthesis and transport, carbohydrate storage metabolism were observed as well as the concomitant accumulation of metabolites, such as oxalic acid, malate and isocitrate. Our data suggest that early metabolic changes in the resistance of Q. pubescens to drought involve a tricarboxylic acid (TCA) cycle shunt through the glyoxylate pathway, galactose metabolism by reducing carbohydrate storage and increased proteolytic activity.

Effect of landiolol on sex-related transcriptomic changes in the myocardium during sepsis.

OBJECTIVES: The aims of this study are to better understand phenotypic differences between male and female rats during sepsis, to characterise the contribution of the beta1-adrenergic blocker landiolol to septic cardiomyopathy and to determine why landiolol induces divergent effects in males and females. METHODS: The myocardial transcriptional profiles in male and female Wistar rats were assessed after the induction of sepsis by cecal ligation and puncture and addition of landiolol. RESULTS: Our results showed major differences in the biological processes activated during sepsis in male and female rats. In particular, a significant decrease in processes related to cell organisation, contractile function, ionic transport and phosphoinositide-3-kinase/AKT (PI3K/AKT) signalling was observed only in males. The transcript of ATPase sarcoplasmic/endoplasmic reticulum Ca2+ transporting 3 (SERCA3) was sex-differently regulated. In males, landiolol reversed several signalling pathways dysregulated during sepsis. The expression level of genes encoding tubulin alpha 8 (TUBA8) and myosin heavy chain 7B (MYH7) contractile proteins, phosphatase 2 catalytic subunit alpha (PPP2CA), G protein-coupled receptor kinase 5 (GRK5) and A-kinase anchoring protein 6 (AKAP6) returned to their basal levels. In contrast, in females, landiolol had limited effects. CONCLUSION: In males, landiolol reversed the expression of many genes that were deregulated in sepsis. Conversely, sepsis-induced deregulation of gene expression was less pronounced in females than in males, and was maintained in the landiolol-treated females. These findings highlight important sex-related differences and confirm previous observations on the important benefit of landiolol intake on cardiac function in male rats.

A critical regulator of Bcl2 revealed by systematic transcript discovery of lncRNAs associated with T-cell differentiation.

Normal T-cell differentiation requires a complex regulatory network which supports a series of maturation steps, including lineage commitment, T-cell receptor (TCR) gene rearrangement, and thymic positive and negative selection. However, the underlying molecular mechanisms are difficult to assess due to limited T-cell models. Here we explore the use of the pro-T-cell line P5424 to study early T-cell differentiation. Stimulation of P5424 cells by the calcium ionophore ionomycin together with PMA resulted in gene regulation of T-cell differentiation and activation markers, partially mimicking the CD4-CD8- double negative (DN) to double positive (DP) transition and some aspects of subsequent T-cell maturation and activation. Global analysis of gene expression, along with kinetic experiments, revealed a significant association between the dynamic expression of coding genes and neighbor lncRNAs including many newly-discovered transcripts, thus suggesting potential co-regulation. CRISPR/Cas9-mediated genetic deletion of Robnr, an inducible lncRNA located downstream of the anti-apoptotic gene Bcl2, demonstrated a critical role of the Robnr locus in the induction of Bcl2. Thus, the pro-T-cell line P5424 is a powerful model system to characterize regulatory networks involved in early T-cell differentiation and maturation.

The nature of the electron acceptor (MnIV/NO3) triggers the differential expression of genes associated with stress and ammonium limitation responses in Shewanella algae C6G3.

Shewanella algae C6G3 can dissimilatively reduce nitrate into ammonium and manganese oxide (MnIV) into MnII. It has the unusual ability to anaerobically produce nitrite from ammonium in the presence of MnIV. To gain insight into their metabolic capabilities, global mRNA expression patterns were investigated by RNA-seq and qRT-PCR in cells growing with lactate and ammonium as carbon and nitrogen sources, and with either MnIV or nitrate as electron acceptors. Genes exhibiting higher expression levels in the presence of MnIV belonged to functional categories of carbohydrate, coenzyme, lipid metabolisms and inorganic ion transport. The comparative transcriptomic pattern between MnIV and NO3 revealed that the strain presented an ammonium limitation status with MnIV, despite the presence of a non-limiting concentration of ammonium under both culture conditions. In addition, in the presence of MnIV, ntrB/nrtC regulators, ammonium channel, nitrogen regulatory protein P-II, glutamine synthetase and asparagine synthetase glutamine-dependent genes were over-represented. Under the nitrate condition, the expression of genes involved in the synthesis of several amino acids was increased. Finally, the expression level of genes associated with the general stress response was also amplified in both conditions and among them, katE, a putative catalase/peroxidase present on several Shewanella genomes, was highly expressed with a median value relatively higher in the MnIV condition.

A transcriptomic signature predicting septic outcome in patients undergoing autologous stem cell transplantation.

Autologous hematopoietic stem cell transplantation is the standard treatment for multiple myeloma and relapsed or refractory lymphomas. After autologous hematopoietic stem cell transplantation, hematologic reconstitution and infectious complications are the two most critical issues. Although many patients develop infectious complications after therapeutic intensification, it remains impossible to predict infection for each individual. The goal of this work was to determine and identify a predictive transcriptomic signature of systemic inflammatory response syndrome and/or sepsis in patients receiving autologous hematopoietic stem cell transplantation. High-throughput transcriptomic and bioinformatics analysis were performed to analyze gene expression modulation in peripheral blood mononuclear cells in 21 patients undergoing autologous hematopoietic stem cell transplantation for hematological malignancies (lymphoma or multiple myeloma). Transcriptomic analysis of peripheral blood mononuclear cells samples collected just after conditioning regimen identified an 11-gene signature (CHAT, CNN3, ANKRD42, LOC100505725, EDAR, GPAT2, ENST00000390425, MTRM8, C6orf192, LOC10289230, and XLOC-005643) that was able to early predict (at least 2-7 days before its occurrence) the development of systemic inflammatory response syndrome or sepsis. The possibility of systemic inflammatory response syndrome or sepsis occurrence early prediction (2-7 days before occurrence) opens up new therapeutic strategies based on preemptive antibiotic and/or antifungal prophylaxis adapted to the specific risk profile of each patient.

Transcriptional Profiling of Midguts Prepared from Trypanosoma/T. congolense-Positive Glossina palpalis palpalis Collected from Two Distinct Cameroonian Foci: Coordinated Signatures of the Midguts' Remodeling As T. congolense-Supportive Niches.

Our previous transcriptomic analysis of Glossina palpalis gambiensis experimentally infected or not with Trypanosoma brucei gambiense aimed to detect differentially expressed genes (DEGs) associated with infection. Specifically, we selected candidate genes governing tsetse fly vector competence that could be used in the context of an anti-vector strategy, to control human and/or animal trypanosomiasis. The present study aimed to verify whether gene expression in field tsetse flies (G. p. palpalis) is modified in response to natural infection by trypanosomes (T. congolense), as reported when insectary-raised flies (G. p. gambiensis) are experimentally infected with T. b. gambiense. This was achieved using the RNA-seq approach, which identified 524 DEGs in infected vs. non-infected tsetse flies, including 285 downregulated genes and 239 upregulated genes (identified using DESeq2). Several of these genes were highly differentially expressed, with log2 fold change values in the vicinity of either +40 or -40. Downregulated genes were primarily involved in transcription/translation processes, whereas encoded upregulated genes governed amino acid and nucleotide biosynthesis pathways. The BioCyc metabolic pathways associated with infection also revealed that downregulated genes were mainly involved in fly immunity processes. Importantly, our study demonstrates that data on the molecular cross-talk between the host and the parasite (as well as the always present fly microbiome) recorded from an experimental biological model has a counterpart in field flies, which in turn validates the use of experimental host/parasite couples.

Omics analysis of mouse brain models of human diseases.

The identification of common gene/protein profiles related to brain alterations, if they exist, may indicate the convergence of the pathogenic mechanisms driving brain disorders. Six genetically engineered mouse lines modelling neurodegenerative diseases and neuropsychiatric disorders were considered. Omics approaches, including transcriptomic and proteomic methods, were used. The gene/protein lists were used for inter-disease comparisons and further functional and network investigations. When the inter-disease comparison was performed using the gene symbol identifiers, the number of genes/proteins involved in multiple diseases decreased rapidly. Thus, no genes/proteins were shared by all 6 mouse models. Only one gene/protein (Gfap) was shared among 4 disorders, providing strong evidence that a common molecular signature does not exist among brain diseases. The inter-disease comparison of functional processes showed the involvement of a few major biological processes indicating that brain diseases of diverse aetiologies might utilize common biological pathways in the nervous system, without necessarily involving similar molecules.

Gene profiling reveals specific oncogenic mechanisms and signaling pathways in oncocytic and papillary thyroid carcinoma.

The oncogenic pathways in mitochondrial-rich thyroid carcinomas are not clearly understood. To investigate the possible implication of mitochondrial abundance in the genesis of thyroid tumors, we have explored the gene expression profile of six oncocytic carcinomas and six mitochondrial-rich papillary carcinomas using cDNA-microarray technology. A supervised approach allowed us to identify 83 genes differentially expressed in the two types of carcinoma. These genes were classified according to their ontologic profiles. Three genes, NOS3, alpha-actinin-2 and alpha-catenin, suspected of playing a role in tumor genesis, were explored by quantitative RT-PCR analysis and immunohistochemistry. Of the 59 genes overexpressed in papillary carcinomas, 51% were involved in cell communication. Of the 24 genes overexpressed in oncocytic carcinomas, 84% were involved in mitochondrial and cellular metabolism. Our results suggest that mitochondrial respiratory chain complexes III and IV play a significant role in the regulation of reactive oxygen species production by oncocytic tumors.

The chemosensitivity to therapy of childhood early B acute lymphoblastic leukemia could be determined by the combined expression of CD34, SPI-B and BCR genes.

We have identified genes differentially expressed in childhood early B acute lymphoblastic leukemia at diagnosis, according to chemosensitivity. Chemosensitive (M1) and chemoresistant (M3) patients present <5% and >25% of residual leukemic blasts at 21 days of treatment, respectively. The expression profiles of 4205 genes for 32 patients included in the FRALLE93 protocol have been determined using microarray. From differential analysis, CD34, SPI-B and BCR distinguished M1 from M3 patients using microarray and RT-PCR data. Linear discriminant analysis (LDA) and cross-validation show that the combined expression of these three genes classify and predict correctly around 90% and 80% of patients, respectively.