MYC regulation of a "poor-prognosis" metastatic cancer cell state.

Gene expression signatures are used in the clinic as prognostic tools to determine the risk of individual patients with localized breast tumors developing distant metastasis. We lack a clear understanding, however, of whether these correlative biomarkers link to a common biological network that regulates metastasis. We find that the c-MYC oncoprotein coordinately regulates the expression of 13 different "poor-outcome" cancer signatures. In addition, functional inactivation of MYC in human breast cancer cells specifically inhibits distant metastasis in vivo and invasive behavior in vitro of these cells. These results suggest that MYC oncogene activity (as marked by "poor-prognosis" signature expression) may be necessary for the translocation of poor-outcome human breast tumors to distant sites.

Selecting control genes for RT-QPCR using public microarray data.

Background: Gene expression analysis has emerged as a major biological research area, with real-time quantitative reverse transcription PCR (RT-QPCR) being one of the most accurate and widely used techniques for expression profiling of selected genes. In order to obtain results that are comparable across assays, a stable normalization strategy is required. In general, the normalization of PCR measurements between different samples uses one to several control genes (e. g. housekeeping genes), from which a baseline reference level is constructed. Thus, the choice of the control genes is of utmost importance, yet there is not a generally accepted standard technique for screening a large number of candidates and identifying the best ones. Results: We propose a novel approach for scoring and ranking candidate genes for their suitability as control genes. Our approach relies on publicly available microarray data and allows the combination of multiple data sets originating from different platforms and/or representing different pathologies. The use of microarray data allows the screening of tens of thousands of genes, producing very comprehensive lists of candidates. We also provide two lists of candidate control genes: one which is breast cancer-specific and one with more general applicability. Two genes from the breast cancer list which had not been previously used as control genes are identified and validated by RT-QPCR. Open source R functions are available at http://www.isrec.isb-sib.ch/similar to vpopovic/research/ Conclusion: We proposed a new method for identifying candidate control genes for RT-QPCR which was able to rank thousands of genes according to some predefined suitability criteria and we applied it to the case of breast cancer. We also empirically showed that translating the results from microarray to PCR platform was achievable.

Chromosomal gene movements reflect the recent origin and biology of therian sex chromosomes

Mammalian sex chromosomes stem from ancestral autosomes and have substantially differentiated. It was shown that X-linked genes have generated duplicate intronless gene copies (retrogenes) on autosomes due to this differentiation. However, the precise driving forces for this out-of-X gene "movement" and its evolutionary onset are not known. Based on expression analyses of male germ-cell populations, we here substantiate and extend the hypothesis that autosomal retrogenes functionally compensate for the silencing of their X-linked housekeeping parental genes during, but also after, male meiotic sex chromosome inactivation (MSCI). Thus, sexually antagonistic forces have not played a major role for the selective fixation of X-derived gene copies in mammals. Our dating analyses reveal that although retrogenes were produced ever since the common mammalian ancestor, selectively driven retrogene export from the X only started later, on the placental mammal (eutherian) and marsupial (metatherian) lineages, respectively. Together, these observations suggest that chromosome-wide MSCI emerged close to the eutherian-marsupial split approximately 180 million years ago. Given that MSCI probably reflects the spread of the recombination barrier between the X and Y, crucial for their differentiation, our data imply that these chromosomes became more widely differentiated only late in the therian ancestor, well after the divergence of the monotreme lineage. Thus, our study also provides strong independent support for the recent notion that our sex chromosomes emerged, not in the common ancestor of all mammals, but rather in the therian ancestor, and therefore are much younger than previously thought

Characterization of root apoplastic barriers in Arabidopsis thaliana

In vascular plants, the best-known feature of a differentiated endodermal cell is the "Casparian Strip" (CS). This structure refers to a highly localized cell wall impregnation in the transversal and anticlinal walls of the cell, which surrounds the cell like a belt/ring and is tightly coordinated with respect to neighboring cells. Analogous to tight junctions in animal epithelia, CS in plants act as a diffusion barrier that controls the movement of water and ions from soil into the stele. Since its first description by Robert Caspary in 1865 there have been many attempts to identify the chemical nature of the cell wall deposition in CS. Suberin, lignin, or both have been claimed to be the important components of CS in a series of different species. However, the exact chemical composition of CS has remained enigmatic. This controversy was due to the confusion and lack of knowledge regarding the precise measurement of three developmental stages of the endodermis. The CS represent only the primary stage of endodermal differentiation, which is followed by the deposition of suberin lamellae all around the cellular surface of endodermal cells (secondary developmental stage). Therefore, chemical analysis of whole roots, or even of isolated endodermal tissues, will always find both of the polymers present. It was crucial to clarify this point because this will guide our efforts to understand which cell wall biosynthetic component becomes localized in order to form the CS. The main aim of my work was to find out the major components of (early) CS, as well as their spatial and temporal development, physiological roles and relationship to barrier formation. Employing the knowledge and tools that have been accumulated over the last few years in the model plant Arabidopsis thaliana, various histological and chemical assays were used in this study. A particular feature of my work was to completely degrade, or inhibit formation of lignin and suberin biopolymers by biochemical, classical genetic and molecular approaches and to investigate its effect on CS formation and the establishment of a functional diffusion barrier. Strikingly, interference with monolignol biosynthesis abrogates CS formation and delays the formation of function diffusion barrier. In contrast, transgenic plants devoid of any detectable suberin still develop a functional CS. The combination of all these assays clearly demonstrates that the early CS polymer is made from monolignol (lignin monomers) and is composed of lignin. By contrast, suberin is formed much later as a secondary wall during development of endodermis. These early CS are functionally sufficient to block extracellular diffusion and suberin does not play important role in the establishment of early endodermal diffusion barrier. Moreover, suberin biosynthetic machinery is not present at the time of CS formation. Our study finally concludes the long-standing debate about the chemical nature of CS and opens the door to a new approach in lignin research, specifically for the identification of the components of the CS biosynthetic pathway that mediates the localized deposition of cell walls. I also made some efforts to understand the patterning and differentiation of endodermal passage cells in young roots. In the literature, passage cells are defined as a non- suberized xylem pole associated endodermal cells. Since these cells only contain the CS but not the suberin lamellae, it has been assumed that these cells may offer a continued low-resistance pathway for water and minerals into the stele. Thus far, no genes have been found to be expressed specifically in passage cells. In order to understand the patterning, differentiation, and physiological role of passage it would be crucial to identify some genes that are exclusively expressed in these cells. In order to identify such genes, I first generated fluorescent marker lines of stele-expressed transporters that have been reported to be expressed in the passage cells. My aim was to first highlight the passage cells in a non-specific way. In order to find passage cell specific genes I then adapted a two-component system based on previously published methods for gene expression profiling of individual cell types. This approach will allow us to target only the passage cells and then to study gene expression specifically in this cell type. Taken together, this preparatory work will provide an entry point to understand the formation and role of endodermal passage cells. - Chez les plantes vasculaires, la caractéristique la plus commune des cellules différentiées de l'endoderme est la présence de cadres de Caspary. Cette structure correspond à une imprégnation localisée des parties transversales et anticlinales de la paroi cellulaire. Cela donne naissance, autour de la cellule, à un anneau/cadre qui est coordonné par rapport aux cellules voisines. De manière analogue aux jonctions serrées des épithéliums chez les animaux, les cadres de Caspary agissent chez les plantes comme barrière de diffusion, contrôlant le mouvement de l'eau et des ions à travers la racine entre le sol et la stèle. Depuis leur première description par Robert Caspary en 1865, beaucoup de tentatives ont eu pour but de définir la nature chimique de ces cadres de Caspary. Après l'étude de différentes espèces végétales, à la fois la subérine, la lignine ou les deux ont été revendiquées comme étant des composants importants de ces cadres. Malgré tout, leur nature chimique exacte est restée longtemps énigmatique. Cette controverse provient de la confusion et du manque de connaissance concernant la détermination précise des trois stades de développement de l'endoderme. Les cadres de Caspary représentent uniquement le stade primaire de différentiation de l'endoderme. Celui-ci est suivi par le second stade de différentiation, la déposition de lamelles de subérine tout autour de la cellule endodermal. De ce fait, l'analyse chimique de racines entières ou de cellules d'endoderme isolées ne permet pas de séparer les stades de différentiation primaire et secondaire et aboutit donc à la présence des deux polymères. Il est également crucial de clarifier ce point dans le but de connaître quelle machinerie cellulaire localisée à la paroi cellulaire permet l'élaboration des cadres de Caspary. En utilisant les connaissances et les outils accumulés récemment grâce à la plante modèle Arabidopsis thaliana, divers techniques histologiques et chimiques ont été utilisées dans cette étude. Un point particulier de mon travail a été de dégrader ou d'inhiber complètement la formation de lignine ou de subérine en utilisant des approches de génétique classique ou moléculaire. Le but étant d'observer l'effet de l'absence d'un de ces deux polymères sur la formation des cadres de Caspary et l'établissement d'une barrière de diffusion fonctionnelle. De manière frappante, le fait d'interférer avec la voie de biosynthèse de monolignol (monomères de lignine) abolit la formation des cadres de Caspary et retarde l'élaboration d'une barrière de diffusion fonctionnelle. Par contre, des plantes transgéniques dépourvues d'une quantité détectable de subérine sont quant à elles toujours capables de développer des cadres de Caspary fonctionnels. Mises en commun, ces expériences démontrent que le polymère formant les cadres de Caspary dans la partie jeune de la racine est fait de monolignol, et que de ce fait il s'agit de lignine. La subérine, quant à elle, est formée bien plus tard durant le développement de l'endoderme, de plus il s'agit d'une modification de la paroi secondaire. Ces cadres de Caspary précoces faits de lignine suffisent donc à bloquer la diffusion extracellulaire, contrairement à la subérine. De plus, la machinerie de biosynthèse de la subérine n'est pas encore présente au moment de la formation des cadres de Caspary. Notre étude permet donc de mettre un terme au long débat concernant la nature chimique des cadres de Caspary. De plus, elle ouvre la porte à de nouvelles approches dans la recherche sur la lignine, plus particulièrement pour identifier des composants permettant la déposition localisée de ce polymère dans la paroi cellulaire. J'ai aussi fais des efforts pour mettre en évidence la formation ainsi que le rôle des cellules de passage dans les jeunes racines. Dans la littérature, les cellules de passage sont définies comme de la cellule endodermal faisant face aux pôles xylèmes et dont la paroi n'est pas subérisée. Du fait que ces cellules contiennent uniquement des cadres de Caspary et pas de lamelle de subérine, il a été supposé qu'elles ne devraient offrir que peu de résistance au passage de l'eau et des nutriments entre le sol et la stèle. Le rôle de ces cellules de passage est toujours loin d'être clair, de plus aucun gène s'exprimant spécifiquement dans ces cellules n'a été découvert à ce jour. De manière à identifier de tels gènes, j'ai tout d'abord généré des marqueurs fluorescents pour des transporteurs exprimés dans la stèle mais dont l'expression avait également été signalée dans l'endoderme, uniquement dans les cellules de passage. J'ai ensuite développé un système à deux composants basé sur des méthodes déjà publiées, visant principalement à étudier le profil d'expression génique dans un type cellulaire donné. En recoupant les gènes exprimés spécifiquement dans l'endoderme à ceux exprimés dans la stèle et les cellules de passage, il nous sera possible d'identifier le transriptome spécifique de ces cellules. Pris dans leur ensemble, ces résultats devraient donner un bon point d'entrée dans la définition et la compréhension des cellules de passage.

Comparative transcriptomics of rice reveals an ancient pattern of response to microbial colonization.

Glomalean fungi induce and colonize symbiotic tissue called arbuscular mycorrhiza on the roots of most land plants. Other fungi also colonize plants but cause disease not symbiosis. Whole-transcriptome analysis using a custom-designed Affymetrix Gene-Chip and confirmation with real-time RT-PCR revealed 224 genes affected during arbuscular mycorrhizal symbiosis. We compared these transcription profiles with those from rice roots that were colonized by pathogens (Magnaporthe grisea and Fusarium moniliforme). Over 40% of genes showed differential regulation caused by both the symbiotic and at least one of the pathogenic interactions. A set of genes was similarly expressed in all three associations, revealing a conserved response to fungal colonization. The responses that were shared between pathogen and symbiont infection may play a role in compatibility. Likewise, the responses that are different may cause disease. Some of the genes that respond to mycorrhizal colonization may be involved in the uptake of phosphate. Indeed, phosphate addition mimicked the effect of mycorrhiza on 8% of the tested genes. We found that 34% of the mycorrhiza-associated rice genes were also associated with mycorrhiza in dicots, revealing a conserved pattern of response between the two angiosperm classes.

STAT3 links IL-22 signaling in intestinal epithelial cells to mucosal wound healing.

Signal transducer and activator of transcription (STAT) 3 is a pleiotropic transcription factor with important functions in cytokine signaling in a variety of tissues. However, the role of STAT3 in the intestinal epithelium is not well understood. We demonstrate that development of colonic inflammation is associated with the induction of STAT3 activity in intestinal epithelial cells (IECs). Studies in genetically engineered mice showed that epithelial STAT3 activation in dextran sodium sulfate colitis is dependent on interleukin (IL)-22 rather than IL-6. IL-22 was secreted by colonic CD11c(+) cells in response to Toll-like receptor stimulation. Conditional knockout mice with an IEC-specific deletion of STAT3 activity were highly susceptible to experimental colitis, indicating that epithelial STAT3 regulates gut homeostasis. STAT3(IEC-KO) mice, upon induction of colitis, showed a striking defect of epithelial restitution. Gene chip analysis indicated that STAT3 regulates the cellular stress response, apoptosis, and pathways associated with wound healing in IECs. Consistently, both IL-22 and epithelial STAT3 were found to be important in wound-healing experiments in vivo. In summary, our data suggest that intestinal epithelial STAT3 activation regulates immune homeostasis in the gut by promoting IL-22-dependent mucosal wound healing.

Gene finding in the chicken genome.

BACKGROUND: Despite the continuous production of genome sequence for a number of organisms, reliable, comprehensive, and cost effective gene prediction remains problematic. This is particularly true for genomes for which there is not a large collection of known gene sequences, such as the recently published chicken genome. We used the chicken sequence to test comparative and homology-based gene-finding methods followed by experimental validation as an effective genome annotation method. RESULTS: We performed experimental evaluation by RT-PCR of three different computational gene finders, Ensembl, SGP2 and TWINSCAN, applied to the chicken genome. A Venn diagram was computed and each component of it was evaluated. The results showed that de novo comparative methods can identify up to about 700 chicken genes with no previous evidence of expression, and can correctly extend about 40% of homology-based predictions at the 5' end. CONCLUSIONS: De novo comparative gene prediction followed by experimental verification is effective at enhancing the annotation of the newly sequenced genomes provided by standard homology-based methods.

Cellular origin of T-cell lymphomas.

In this issue of Blood, Iqbal et al, having compiled gene expression profiles from >300 peripheral T-cell lymphomas, expand previous findings on the diagnostic value of molecular signatures that correlate with different histological types of T-cell lymphomas. They report the discovery of 2 molecular subgroups of peripheral T-cell lymphomas, not otherwise specified (PTCL, NOS), characterized by high expression of either GATA-binding protein 3 (GATA-3) or t-box 21 (TBX21) transcription factors and corresponding target genes, with the GATA3 subgroup being associated with distinctly worse prognosis. In an independent study, Wang et al(2) also show that GATA3 expression in a subset of PTCL, NOS identifies a subgroup of patients with inferior survival.

Epigenetic features of human mesenchymal stem cells determine their permissiveness for induction of relevant transcriptional changes by SYT-SSX1.

BACKGROUND: A characteristic SYT-SSX fusion gene resulting from the chromosomal translocation t(X;18)(p11;q11) is detectable in almost all synovial sarcomas, a malignant soft tissue tumor widely believed to originate from as yet unidentified pluripotent stem cells. The resulting fusion protein has no DNA binding motifs but possesses protein-protein interaction domains that are believed to mediate association with chromatin remodeling complexes. Despite recent advances in the identification of molecules that interact with SYT-SSX and with the corresponding wild type SYT and SSX proteins, the mechanisms whereby the SYT-SSX might contribute to neoplastic transformation remain unclear. Epigenetic deregulation has been suggested to be one possible mechanism. METHODOLOGY/PRINCIPAL FINDINGS: We addressed the effect of SYT/SSX expression on the transcriptome of four independent isolates of primary human bone marrow mesenchymal stem cells (hMSC). We observed transcriptional changes similar to the gene expression signature of synovial sarcoma, principally involving genes whose regulation is linked to epigenetic factors, including imprinted genes, genes with transcription start sites within a CpG island and chromatin related genes. Single population analysis revealed hMSC isolate-specific transcriptional changes involving genes that are important for biological functions of stem cells as well as genes that are considered to be molecular markers of synovial sarcoma including IGF2, EPHRINS, and BCL2. Methylation status analysis of sequences at the H19/IGF2 imprinted locus indicated that distinct epigenetic features characterize hMSC populations and condition the transcriptional effects of SYT-SSX expression. CONCLUSIONS/SIGNIFICANCE: Our observations suggest that epigenetic features may define the cellular microenvironment in which SYT-SSX displays its functional effects.

Changes in the transcriptional profile of transporters in the intestine along the anterior-posterior and crypt-villus axes.

BACKGROUND: The purpose of this work was to characterize the expression of drug and nutrient carriers along the anterior-posterior and crypt-villus axes of the intestinal epithelium and to study the validity of utilizing whole gut tissue rather than purified epithelial cells to examine regional variations in gene expression. RESULTS: We have characterized the mRNA expression profiles of 76 % of all currently known transporters along the anterior-posterior axis of the gut. This is the first study to describe the expression profiles of the majority of all known transporters in the intestine. The expression profiles of transporters, as defined according to the Gene Ontology consortium, were measured in whole tissue of the murine duodenum, jejunum, ileum and colon using high-density microarrays. For nine transporters (Abca1, Abcc1, Abcc3, Abcg8, Slc10a2, Slc28a2, Slc2a1, Slc34a2 and Slc5a8), the mRNA profiles were further measured by RT-PCR in laser micro-dissected crypt and villus epithelial cells corresponding to the aforementioned intestinal regions. With respect to differentially regulated transporters, the colon had a distinct expression profile from small intestinal segments. The majority (59 % for p cutoff < or = 0.05) of transporter mRNA levels were constant across the intestinal sections studied. For the transporter subclass "carrier activity", which contains the majority of known carriers for biologically active compounds, a significant change (p < or = 0.05) along the anterior-posterior axis was observed. CONCLUSION: All nine transporters examined in laser-dissected material demonstrated good replication of the region-specific profiles revealed by microarray. Furthermore, we suggest that the distribution characteristics of Slc5a8 along the intestinal tract render it a suitable candidate carrier for monocarboxylate drugs in the posterior portion of the intestine. Our findings also predict that there is a significant difference in the absorption of carrier-mediated compounds in the different intestinal segments. The most pronounced differences can be expected between the adjoining segments ileum and colon, but the differences between the other adjoining segments are not negligible. Finally, for the examined genes, profiles measured in whole intestinal tissue extracts are representative of epithelial cell-only gene expression.

Regulation of gammadelta versus alphabeta T lymphocyte differentiation by the transcription factor SOX13.

alphabeta and gammadelta T cells originate from a common, multipotential precursor population in the thymus, but the molecular mechanisms regulating this lineage-fate decision are unknown. We have identified Sox13 as a gammadelta-specific gene in the immune system. Using Sox13 transgenic mice, we showed that this transcription factor promotes gammadelta T cell development while opposing alphabeta T cell differentiation. Conversely, mice deficient in Sox13 expression exhibited impaired development of gammadelta T cells but not alphabeta T cells. One mechanism of SOX13 function is the inhibition of signaling by the developmentally important Wnt/T cell factor (TCF) pathway. Our data thus reveal a dominant pathway regulating the developmental fate of these two lineages of T lymphocytes.

Histological and histochemical evaluation of human oral mucosa constructs developed by tissue engineering

Reconstruction of large oral mucosa defects is often challenging, since the shortage of healthy oral mucosa to replace the excised tissues is very common. In this context, tissue engineering techniques may provide a source of autologous tissues available for transplant in these patients. In this work, we developed a new model of artificial oral mucosa generated by tissue engineering using a fibrin-agarose scaffold. For that purpose, we generated primary cultures of human oral mucosa fibroblasts and keratinocytes from small biopsies of normal oral mucosa using enzymatic treatments. Then we determined the viability of the cultured cells by electron probe quantitative X-ray microanalysis, and we demonstrated that most of the cells in the primary cultures were alive and had high K/Na ratios. Once cell viability was determined, we used the cultured fibroblasts and keratinocytes to develop an artificial oral mucosa construct by using a fibrin-agarose extracellular matrix and a sequential culture technique using porous culture inserts. Histological analysis of the artificial tissues showed high similarities with normal oral mucosa controls. The epithelium of the oral substitutes had several layers, with desmosomes and apical microvilli and microplicae. Both the controls and the oral mucosa substitutes showed high suprabasal expression of cytokeratin 13 and low expression of cytokeratin 10. All these results suggest that our model of oral mucosa using fibrin-agarose scaffolds show several similarities with native human oral mucosa.

Identification of tumor-associated antigens by large-scale analysis of genes expressed in human colorectal cancer.

Despite the high prevalence of colon cancer in the world and the great interest in targeted anti-cancer therapy, only few tumor-specific gene products have been identified that could serve as targets for the immunological treatment of colorectal cancers. The aim of our study was therefore to identify frequently expressed colon cancer-specific antigens. We performed a large-scale analysis of genes expressed in normal colon and colon cancer tissues isolated from colorectal cancer patients using massively parallel signal sequencing (MPSS). Candidates were additionally subjected to experimental evaluation by semi-quantitative RT-PCR on a cohort of colorectal cancer patients. From a pool of more than 6000 genes identified unambiguously in the analysis, we found 2124 genes that were selectively expressed in colon cancer tissue and 147 genes that were differentially expressed to a significant degree between normal and cancer cells. Differential expression of many genes was confirmed by RT-PCR on a cohort of patients. Despite the fact that deregulated genes were involved in many different cellular pathways, we found that genes expressed in the extracellular space were significantly over-represented in colorectal cancer. Strikingly, we identified a transcript from a chromosome X-linked member of the human endogenous retrovirus (HERV) H family that was frequently and selectively expressed in colon cancer but not in normal tissues. Our data suggest that this sequence should be considered as a target of immunological interventions against colorectal cancer.

Modeling stochasticity and robustness in gene regulatory networks.

MOTIVATION: Understanding gene regulation in biological processes and modeling the robustness of underlying regulatory networks is an important problem that is currently being addressed by computational systems biologists. Lately, there has been a renewed interest in Boolean modeling techniques for gene regulatory networks (GRNs). However, due to their deterministic nature, it is often difficult to identify whether these modeling approaches are robust to the addition of stochastic noise that is widespread in gene regulatory processes. Stochasticity in Boolean models of GRNs has been addressed relatively sparingly in the past, mainly by flipping the expression of genes between different expression levels with a predefined probability. This stochasticity in nodes (SIN) model leads to over representation of noise in GRNs and hence non-correspondence with biological observations. RESULTS: In this article, we introduce the stochasticity in functions (SIF) model for simulating stochasticity in Boolean models of GRNs. By providing biological motivation behind the use of the SIF model and applying it to the T-helper and T-cell activation networks, we show that the SIF model provides more biologically robust results than the existing SIN model of stochasticity in GRNs. AVAILABILITY: Algorithms are made available under our Boolean modeling toolbox, GenYsis. The software binaries can be downloaded from http://si2.epfl.ch/ approximately garg/genysis.html.

Breast cancer prognostic classification in the molecular era: the role of histological grade.

Breast cancer is a heterogeneous disease with varied morphological appearances, molecular features, behavior, and response to therapy. Current routine clinical management of breast cancer relies on the availability of robust clinical and pathological prognostic and predictive factors to support clinical and patient decision making in which potentially suitable treatment options are increasingly available. One of the best-established prognostic factors in breast cancer is histological grade, which represents the morphological assessment of tumor biological characteristics and has been shown to be able to generate important information related to the clinical behavior of breast cancers. Genome-wide microarray-based expression profiling studies have unraveled several characteristics of breast cancer biology and have provided further evidence that the biological features captured by histological grade are important in determining tumor behavior. Also, expression profiling studies have generated clinically useful data that have significantly improved our understanding of the biology of breast cancer, and these studies are undergoing evaluation as improved prognostic and predictive tools in clinical practice. Clinical acceptance of these molecular assays will require them to be more than expensive surrogates of established traditional factors such as histological grade. It is essential that they provide additional prognostic or predictive information above and beyond that offered by current parameters. Here, we present an analysis of the validity of histological grade as a prognostic factor and a consensus view on the significance of histological grade and its role in breast cancer classification and staging systems in this era of emerging clinical use of molecular classifiers.