Error models for microarray intensities

We derive the additive-multiplicative error model for microarray intensities, and describe two applications. For the detection of differentially expressed genes, we obtain a statistic whose variance is approximately independent of the mean intensity. For the post hoc calibration (normalization) of data with respect to experimental factors, we describe a method for parameter estimation.

GENE SET ENRICHMENT ANALYSIS MADE SIMPLE

Among the many applications of microarray technology, one of the most popular is the identification of genes that are differentially expressed in two conditions. A common statistical approach is to quantify the interest of each gene with a p-value, adjust these p-values for multiple comparisons, chose an appropriate cut-off, and create a list of candidate genes. This approach has been criticized for ignoring biological knowledge regarding how genes work together. Recently a series of methods, that do incorporate biological knowledge, have been proposed. However, many of these methods seem overly complicated. Furthermore, the most popular method, Gene Set Enrichment Analysis (GSEA), is based on a statistical test known for its lack of sensitivity. In this paper we compare the performance of a simple alternative to GSEA.We find that this simple solution clearly outperforms GSEA.We demonstrate this with eight different microarray datasets.

Short communication: duplication in the 5'-flanking region of the beta-lactoglobulin gene is linked to the BLG A allele

beta-Lactoglobulin (beta-LG) is the major whey protein in the milk of cows and other ruminants. It is well established that the predominant genetic variants beta-LG A and B are differentially expressed. Extensive investigation of the genetic variation in the promoter region of the BLG gene revealed the existence of specific haplotypes associated with the A and B variants. However, the genetic basis for the differentially expressed BLG A and B alleles is still elusive. In this study additional genetic variation further upstream in the 5'-flanking region of the BLG gene was identified, including 6 single nucleotide substitutions, a single nucleotide deletion, and a 7-bp duplication. Comparison of DNA sequences showed that the investigated 5'-flanking region is highly conserved between ruminants, and the duplication g.-1885_-1879dupCTCTCGC and the substitution g.-1888A>G are only found in the BLG A and D alleles in cattle. The cytosine at position g.-1957 and the thymines at positions g.-2008 and g.-2049 are only found in BLG B alleles of cattle. It is suggested that the described genetic variability contributes to the differential allelic expression of the BLG gene.

Systemic analysis of microarray data sets towards identifying genes regulating root development

Nitrogen and water are essential for plant growth and development. In this study, we designed experiments to produce gene expression data of poplar roots under nitrogen starvation and water deprivation conditions. We found low concentration of nitrogen led first to increased root elongation followed by lateral root proliferation and eventually increased root biomass. To identify genes regulating root growth and development under nitrogen starvation and water deprivation, we designed a series of data analysis procedures, through which, we have successfully identified biologically important genes. Differentially Expressed Genes (DEGs) analysis identified the genes that are differentially expressed under nitrogen starvation or drought. Protein domain enrichment analysis identified enriched themes (in same domains) that are highly interactive during the treatment. Gene Ontology (GO) enrichment analysis allowed us to identify biological process changed during nitrogen starvation. Based on the above analyses, we examined the local Gene Regulatory Network (GRN) and identified a number of transcription factors. After testing, one of them is a high hierarchically ranked transcription factor that affects root growth under nitrogen starvation. It is very tedious and time-consuming to analyze gene expression data. To avoid doing analysis manually, we attempt to automate a computational pipeline that now can be used for identification of DEGs and protein domain analysis in a single run. It is implemented in scripts of Perl and R.

DISSECTION OF STRESS RESPONSE NETWORKS REGULATING MULTIPLE STRESSES IN RICE

Important food crops like rice are constantly exposed to various stresses that can have devastating effect on their survival and productivity. Being sessile, these highly evolved organisms have developed elaborate molecular machineries to sense a mixture of stress signals and elicit a precise response to minimize the damage. However, recent discoveries revealed that the interplay of these stress regulatory and signaling molecules is highly complex and remains largely unknown. In this work, we conducted large scale analysis of differential gene expression using advanced computational methods to dissect regulation of stress response which is at the heart of all molecular changes leading to the observed phenotypic susceptibility. One of the most important stress conditions in terms of loss of productivity is drought. We performed genomic and proteomic analysis of epigenetic and miRNA mechanisms in regulation of drought responsive genes in rice and found subsets of genes with striking properties. Overexpressed genesets included higher number of epigenetic marks, miRNA targets and transcription factors which regulate drought tolerance. On the other hand, underexpressed genesets were poor in above features but were rich in number of metabolic genes with multiple co-expression partners contributing majorly towards drought resistance. Identification and characterization of the patterns exhibited by differentially expressed genes hold key to uncover the synergistic and antagonistic components of the cross talk between stress response mechanisms. We performed meta-analysis on drought and bacterial stresses in rice and Arabidopsis, and identified hundreds of shared genes. We found high level of conservation of gene expression between these stresses. Weighted co-expression network analysis detected two tight clusters of genes made up of master transcription factors and signaling genes showing strikingly opposite expression status. To comprehensively identify the shared stress responsive genes between multiple abiotic and biotic stresses in rice, we performed meta-analyses of microarray studies from seven different abiotic and six biotic stresses separately and found more than thirteen hundred shared stress responsive genes. Various machine learning techniques utilizing these genes classified the stresses into two major classes' namely abiotic and biotic stresses and multiple classes of individual stresses with high accuracy and identified the top genes showing distinct patterns of expression. Functional enrichment and co-expression network analysis revealed the different roles of plant hormones, transcription factors in conserved and non-conserved genesets in regulation of stress response.

Expression and hypoxic regulation of the endothelin system in endocrine cells of human and rat pancreatic islets

CONTEXT: The success of pancreatic islet transplantation depends largely on the capacity of the islet graft to survive the initial phase immediately after transplantation until revascularization is completed. Endothelin-1 (ET-1) is a strong vasoconstrictor which has been involved in solid organ graft failure but is also known to be a potent mitogenic/anti-apoptotic factor which could also potentially enhance the survival of the transplanted islets. OBJECTIVE: Characterization of the endothelin system with regard to a potential endothelin agonist/antagonist treatment. DESIGN: Regulated expression of the endothelin system in human and rat pancreatic islets and beta-cell lines was assessed by means of immunohistochemistry, competition binding studies, western blot, RT-PCR, real-time PCR and transplant studies. RESULTS: ET-1, ETA- and ETB-receptor immunoreactivity was identified in the endocrine cells of human and rat pancreatic islets. The corresponding mRNA was detectable in rat beta-cell lines and isolated rat and human pancreatic islets. Competition binding studies on rat islets revealed binding sites for both receptor types. ET-1 stimulated the phosphorylation of mitogen-activated protein kinase, which was prevented by ETA- and ETB-receptor antagonists. After exposure to hypoxia equal to post-transplant environment oxygen tension, mRNA levels of ET-1 and ETB-receptor of human islets were robustly induced whereas ETA-receptor mRNA did not show significant changes. Immunostaining signals for ET-1 and ETA-receptor of transplanted rat islets were markedly decreased when compared to native pancreatic sections. CONCLUSIONS: In pancreatic islets, ET-1 and its receptors are differentially expressed by hypoxia and after transplantation. Our results provide the biological basis for the study of the potential use of endothelin agonists/antagonists to improve islet transplantation outcome.

Integrative analysis of RUNX1 downstream pathways and target genes

BACKGROUND: The RUNX1 transcription factor gene is frequently mutated in sporadic myeloid and lymphoid leukemia through translocation, point mutation or amplification. It is also responsible for a familial platelet disorder with predisposition to acute myeloid leukemia (FPD-AML). The disruption of the largely unknown biological pathways controlled by RUNX1 is likely to be responsible for the development of leukemia. We have used multiple microarray platforms and bioinformatic techniques to help identify these biological pathways to aid in the understanding of why RUNX1 mutations lead to leukemia. RESULTS: Here we report genes regulated either directly or indirectly by RUNX1 based on the study of gene expression profiles generated from 3 different human and mouse platforms. The platforms used were global gene expression profiling of: 1) cell lines with RUNX1 mutations from FPD-AML patients, 2) over-expression of RUNX1 and CBFbeta, and 3) Runx1 knockout mouse embryos using either cDNA or Affymetrix microarrays. We observe that our datasets (lists of differentially expressed genes) significantly correlate with published microarray data from sporadic AML patients with mutations in either RUNX1 or its cofactor, CBFbeta. A number of biological processes were identified among the differentially expressed genes and functional assays suggest that heterozygous RUNX1 point mutations in patients with FPD-AML impair cell proliferation, microtubule dynamics and possibly genetic stability. In addition, analysis of the regulatory regions of the differentially expressed genes has for the first time systematically identified numerous potential novel RUNX1 target genes. CONCLUSION: This work is the first large-scale study attempting to identify the genetic networks regulated by RUNX1, a master regulator in the development of the hematopoietic system and leukemia. The biological pathways and target genes controlled by RUNX1 will have considerable importance in disease progression in both familial and sporadic leukemia as well as therapeutic implications.

Novel functional profiling approach combining reverse phase protein microarrays and human 3-D ex vivo tissue cultures: expression of apoptosis-related proteins in human colon cancer

Cancer is caused by a complex pattern of molecular perturbations. To understand the biology of cancer, it is thus important to look at the activation state of key proteins and signaling networks. The limited amount of available sample material from patients and the complexity of protein expression patterns make the use of traditional protein analysis methods particularly difficult. In addition, the only approach that is currently available for performing functional studies is the use of serial biopsies, which is limited by ethical constraints and patient acceptance. The goal of this work was to establish a 3-D ex vivo culture technique in combination with reverse-phase protein microarrays (RPPM) as a novel experimental tool for use in cancer research. The RPPM platform allows the parallel profiling of large numbers of protein analytes to determine their relative abundance and activation level. Cancer tissue and the respective corresponding normal tissue controls from patients with colorectal cancer were cultured ex vivo. At various time points, the cultured samples were processed into lysates and analyzed on RPPM to assess the expression of carcinoembryonic antigen (CEA) and 24 proteins involved in the regulation of apoptosis. The methodology displayed good robustness and low system noise. As a proof of concept, CEA expression was significantly higher in tumor compared with normal tissue (p<0.0001). The caspase 9 expression signal was lower in tumor tissue than in normal tissue (p<0.001). Cleaved Caspase 8 (p=0.014), Bad (p=0.007), Bim (p=0.007), p73 (p=0.005), PARP (p<0.001), and cleaved PARP (p=0.007) were differentially expressed in normal liver and normal colon tissue. We demonstrate here the feasibility of using RPPM technology with 3-D ex vivo cultured samples. This approach is useful for investigating complex patterns of protein expression and modification over time. It should allow functional proteomics in patient samples with various applications such as pharmacodynamic analyses in drug development.

Mammary epithelial-specific knockout of the ephrin-B2 gene leads to precocious epithelial cell death at lactation

The family of Eph receptor tyrosine kinases and their membrane bound ligands, the ephrins, are involved in a wide variety of morphogenic processes during embryonic development and adult tissue homeostasis. Receptor-ligand interaction requires direct cell-cell contact and results in forward and reverse signaling originating from the receptor and ligand, respectively. We have previously shown that EphB4 and ephrinB2 are differentially expressed during the development of the adult mammary parenchyma. Overexpression of EphB4 in the mammary epithelium of transgenic mice leads to perturbations in mammary epithelial morphology, motility and growth. To investigate the role of ephrinB2 signaling in mammary gland biology, we have established transgenic mice exhibiting conditional ephrinB2 knockout in the mammary epithelium. In homozygote double transgenic CreLox mice, specific knockout of ephrinB2 occurred in the mammary epithelium during the first pregnancy-lactating period. Abolishing ephrinB2 function led to severe interference with the architecture and functioning of the mammary gland at lactation. The morphology of the transgenic lactating glands resembled that of involuting controls, with decreased epithelial cell number and collapsed lobulo-alveolar structures. Accordingly, massive epithelial cell death and expression of involution-specific genes were observed. Interestingly, in parallel to cell death, significant cell proliferation was apparent, suggestive of tissue regeneration.

Elevated albumin in retinas of monkeys with experimental glaucoma.

PURPOSE: To establish the identity of a prominent protein, approximately 70 kDa, that is markedly increased in the retina of monkeys with experimental glaucoma compared with the fellow control retina, the relationship to glaucoma severity, and its localization in the retina. METHODS: Retinal extracts were subjected to 2-D gel electrophoresis to identify differentially expressed proteins. Purified peptides from the abundant 70 kDa protein were analyzed and identified by liquid chromatography/mass spectrometry/mass spectrometry (LC/MS/MS) separation, and collision-induced dissociation sequencing. Protein identity was performed on MASCOT (Matrix Science, Boston, MA) and confirmed by Western blot. The relationship between the increase in this protein and glaucoma severity was investigated by regression analyses. Protein localization in retina was evaluated by immunohistochemistry with confocal imaging. RESULTS: The abundant protein was identified as Macaca mulatta serum albumin precursor (67 kDa) from eight non-overlapping proteolytic fragments, and the identity was confirmed by Western blot. The average increase in retinal albumin content was 2.3 fold (P = 0.015). In glaucoma eyes, albumin was localized to some neurons of the inner nuclear layer, in the inner plexiform layer, and along the vitreal surface, but it was only found in blood vessels in control retinas. CONCLUSIONS: Albumin is the abundant protein found in the glaucomatous monkey retinas. The increased albumin is primarily localized to the inner retina where oxidative damage associated with experimental glaucoma is known to be prominent. Since albumin is a major antioxidant, the increase of albumin in the retinas of eyes with experimental glaucoma may serve to protect the retina against oxidative damage.

LDOC1, A NOVEL BIOMARKER OF PROGNOSIS IN CHRONIC LYMPHOCYTIC LEUKEMIA

In chronic lymphocytic leukemia (CLL), one of the best predictors of outcome is the somatic mutation status of the immunoglobulin heavy chain variable region (IGHV) genes. Patients whose CLL cells have unmutated IGHV genes have a median survival of 8 years; those with mutated IGHV genes have a median survival of 25 years. To identify new prognostic biomarkers and molecular targets for therapy in untreated CLL patients, we reanalyzed the raw data from four published gene expression profiling microarray studies. Of 88 candidate biomarkers associated with IGHV somatic mutation status, we identified LDOC1 (Leucine Zipper, Down-regulated in Cancer 1), as one of the most significantly differentially expressed genes that distinguished mutated from unmutated CLL cases. LDOC1 is a putative transcription factor of unknown function in B-cell development and CLL pathophysiology. Using a highly sensitive quantitative RT-PCR (QRT-PCR) assay, we confirmed that LDOC1 mRNA was dramatically down-regulated in mutated compared to unmutated CLL cases. Expression of LDOC1 mRNA was also vii strongly associated with other markers of poor prognosis, including ZAP70 protein and cytogenetic abnormalities of poor prognosis (deletions of chromosomes 6q21, 11q23, and 17p13.1, and trisomy 12). CLL cases positive for LDOC1 mRNA had significantly shorter overall survival than negative cases. Moreover, in a multivariate model, LDOC1 mRNA expression predicted overall survival better than IGHV mutation status or ZAP70 protein, among the best markers of prognosis in CLL. We also discovered LDOC1S, a new LDOC1 splice variant. Using isoform-specific QRT-PCR assays that we developed, we found that both isoforms were expressed in normal B cells (naïve > memory), unmutated CLL cells, and in B-cell non-Hodgkin lymphomas with unmutated IGHV genes. To investigate pathways in which LDOC1 is involved, we knocked down LDOC1 in HeLa cells and performed global gene expression profiling. GFI1 (Growth Factor-Independent 1) emerged as a significantly up-regulated gene in both HeLa cells and CLL cells that expressed high levels of LDOC1. GFI1 oncoprotein is implicated in hematopoietic stem cell maintenance, lymphocyte development, and lymphomagenesis. Our findings indicate that LDOC1 mRNA is an excellent biomarker of overall survival in CLL, and may contribute to B-cell differentiation and malignant transformation.

DYSREGULATION OF MEOX2 FOLLOWING WT1 MUTATION IN KIDNEY DEVELOPMENT AND WILMS TUMORIGENESIS

Wilms tumor (WT) is a childhood tumor of the kidney and a productive model for understanding the role of genetic alteration and interactions in tumorigenesis. The Wilms tumor gene 1 (WT1) is a transcriptional factor and one of the few genes known to have genetic alterations in WT and has been shown be inactivated in 20% of WTs. However, the mechanisms of how WT1 mutations lead to Wilms tumorigenesis and its influence on downstream genes are unknown. Since it has been established that WT1 is a transcriptional regulator, it has been hypothesized that the loss of WT1 leads to the dysregulation of downstream genes, in turn result in the formation of WTs. To identify the dysregulated downstream genes following WT1 mutations, an Affymetrix GeneChip Human Genome Array was previously conducted to assess the differentially expressed genes in the WT1-wildtype human and WT1-mutant human WTs. Approximately 700 genes were identified as being significantly dysregulated. These genes were further prioritized based on their statistical significance, fold change, chromosomal region, spatial pattern of gene expression and known or putative cellular functions. Mesenchyme homeobox 2 (MEOX2) was one of the most significantly upregulated genes in WT1-mutant WT. MEOX2 is known to play a role in cell proliferation, apoptosis, and differentiation. In addition to its biological roles, it is expressed during early kidney development in the condensed mesenchyme similar to WT1. Furthermore, the use of the Match® web-based tool from the BIOBASE Biological Data base identified a significant predicted WT1 binding site within the first intron of MEOX2. The similarity in spatial gene expression in the developing kidney and the significant predicted WT1 binding site found in the first intron of MEOX2 lead to the development of my hypothesis that MEOX2 is upregulated via a WT1-dependent manner. Here as a part of my master’s work, I have validated the Affymetrix GeneChip Human Genome Array data using an independent set of Wilms tumors. MEOX2 remained upregulated in the mutant WT1 Wilms tumor by 41-fold. Wt1 and Meox2 gene expression were assessed in murine newborn kidney; both Wt1 and Meox2 were expressed in the condensed, undifferentiated metanephric mesenchyme. I have shown that the in vivo ablation of Wt1 during embryonic development at embryonic day (E) 13.5 resulted in the slight increase of Meox2 gene expression by two fold. In order to functionally demonstrate the effect of the loss of Wt1 on Meox2 gene expression in undifferentiated metanephric mesenchyme, I have generated a kidney mesenchymal cell line to genetically ablate Wt1 in vitro by adenoviral infection. The ablation of Wt1 in the kidney mesenchymal cell line resulted in the upregulation of Meox2 by 61-fold. Moreover, the upregulation of Meox2 resulted in the significant induction of p21 and Itgb5. In addition to the dysregulation of these genes the ablation of Wt1 in the kidney mesenchymal cells resulted in decrease in cell growth and loss of cellular adherence. However, it is uncertain whether the upregulation of Meox2 caused this particular cellular phenotype. Overall, I have demonstrated that the upregulation of Meox2 is Wt1-dependent during early kidney development.

Differential gene expression in multistep tumorigenesis using an {\it in vitro\/} human cell model

Using a human terato-carcinoma cell line, PA-1, the functional role of the oncogenes and tumor suppressor gene involved in the multistep process of carcinogenesis have been analyzed. The expression of AP-2 was strongly correlated with the susceptibility to ras transformation. The differential responsiveness to growth factors between stage 1 ras resistant cells and stage 2 ras susceptible cells was observed, indicating that the ability of stage 2 cells to respond to the mutated ras oncogenes in transformation correlated with the ability to be stimulated by certain growth factors. Using differential screening of cDNA libraries, a number of differentially expressed cDNA clones was isolated. One of those, clone 12, is overexpressed in ras transformed stage 3 cells. The amino acid sequence of clone 12 is almost identical to a mouse LLrep3 gene that was growth-regulated, and 78% similar to a yeast ribosomal protein S4. These results suggest that the S4 gene may be involved in regulation of growth. Clone 9 is expressed in stage 1 ras resistant cells (3.5-kb and 3.0-kb transcripts) but the expression of this clone in stage 2 ras susceptible cells and stage 3 ras-transformed cells is greatly diminished. The expression of this cDNA clone was increased to at least five fold in ras resistant cells and nontumorigenic hybrids treated with retinoic acid but not increased in retinoic acid treated ras susceptible cells, ras transformed cells and the tumorigenic segregants. Partial sequence of this clone showed no homology to the sequences in Genbank. These findings suggest that clone 9 could be a suppressor gene or the genes that are involved in the biochemical pathway of tumor suppression or neurogenic differentiation. The apparent pleiotropic effect of the loss of this suppressor gene function support Harris' proposal that tumor suppressor genes regulate differentiation. The tumor suppressor gene may act as negative regulator of tumor growth by controlling gene expression in differentiation. ^

The molecular cloning and characterization of human heparanase cDNA and the immunochemical localization of heparanase in metastatic melanomas

Heparanase, an endo-$\beta$-D-glucuronidase, has been associated with melanoma metastasis. Polyclonal antibodies directed against the murine N-terminal heparanase peptide detected a M$\sb{\rm r}\sim 97,000$ protein upon SDS-polyacrylamide gel electrophoresis of mouse melanoma and human melanoma cell lysates. In an indirect immunocytochemical study, metastatic human A375-SM and mouse B16-BL6 melanoma cells were stained with the anti-heparanase antibodies. Heparanase antigen was localized in the cytoplasm of permeabilized melanoma cells as well as at the cell surface of unpermeabilized cells. Immunohistochemical staining of frozen sections from syngeneic mouse organs containing micrometastases of B16-BL6 melanoma demonstrated heparanase localized in metastatic melanoma cells, but not in adjacent normal tissues. Similar studies using frozen sections of malignant melanomas resected from patients indicated that heparanase is localized in invading melanoma cells, but not in adjacent connective tissues.^ Monoclonal antibodies directed against murine heparanase were developed and characterized. Monoclonal antibody 10E5, an IgM, precipitated and inhibitated the enzymatic activity of heparanase. A 2.6 kb cDNA was isolated from a human melanoma $\lambda$gt11 cDNA library using the monoclonal antibody 10E5. Heparan sulfate cleavage activity was detected in the lysogen lysates from E. Coli Y1089 infected with the $\lambda$gt11 cDNA and this activity was inhibited in the presence of 10-fold excess of heparin, a potent inhibitor of heparanase. The nucleotide sequence of the cDNA was determined and insignificant homology was found with the gene sequences currently known. The cDNA hybridized to a 3.2-3.4 kb mRNA in human A375 melanoma, WI-38 fibroblast, and THP-1 leukemia cells using Northern blots.^ Heparanase expression was examined using Western and Northern blots. In comparison to human A375-P melanoma cells, the quantity of 97,000 protein recognized by the polyclonal anti-heparanase antibodies doubled in the metastatic variant A375-SM cells and the quantity of 3.2-3.4 kb mRNA doubled in A375MetMix, a metastatic variant similar to A375-SM cells. In B16 murine melanoma cell, the intensity of the 97,000 protein increased more than 2 times comparing with B16-F1 cells. The extent in the increase of the protein and the mRNA levels is comparable to the change of heparanase activity observed in those cells.^ In summary, the studies suggest that (a) the N-terminus of the heparanase molecule in mouse and human is antigenically related; (b) heparanase antigens are localized at the cell surface and in the cytoplasm of metastatic human and mouse melanoma cells; (c) heparanase antigens are localized in invasive and metastatic murine and human melanomas in vivo, but not in adjacent normal tissues; (d) heparanase molecule appeared to be differentially expressed at the transcriptional as well as at the translational level; and (e) the size of human heparanase mRNA is 3.2-3.4 kilobase. ^

Studies on the adhesion properties of the liver-metastatic murine RAW117 large-cell lymphoma cells: Roles of $\beta\sb3$ integrin

Metastasis is the major cause of death in cancer patients. Since many cancers show organ-preference of metastasis, elucidation of the underlying mechanisms of metastasis will benefit diagnosis or treatment of metastatic diseases. Adhesion mechanisms are thought to be involved in organ-preference of metastasis, because metastatic cells show organ preference in adhering to organ-derived microvascular endothelial cells. The adhesion molecules in this process remain largely unidentified. I have examined a series of murine RAW117 large-cell lymphoma cells variants selected in vivo for liver-colonizing properties ($\rm{H10{>>}L17>P}$). The highly liver-metastatic H10 cells were found to differentially express much higher levels of integrin $\alpha\rm\sb{v}\beta\sb3$ than L17 or P cells. H10 cells also adhered at higher rates to vitronectin and fibronectin than to fibrinogen, fibrin, laminin and type I collagen, and adhered at significantly higher rates to (GRGDS)$\sb4$ than to monomeric RGD-peptides. In contrast, P and L17 cells did not adhere well to the above substrates. H10 cells also spread well on vitronectin and migrated toward vitronectin concentration gradients. Pretreament of H10 cells with anti-$\beta\sb3$ monoclonal antibodies resulted in significant decreases in adhesion of H10 cells to vitronectin and immobilized (GRGDS)$\sb4$, and reduced the formation of experimental liver metastases in syngeneic Balb/c mice.^ Adhesion of RAW117 cells under hydrodynamic shear stresses was also studied because tumor cell adhesion occurs under fluid shear stresses in target organ microvessels. Similar to their properties found with static adhesion assays, H10 cells stabilized their hydrodynamic adhesion to vitronectin, fibronectin and (GRGDS)$\sb4$ much more quickly than P or L17 cells. Unlike their static adhesion properties, RAW117 cells showed differential adhesion stabilization to liver-sinusoidal endothelial cell-derived extracellular matrix ($\rm{H10{>>}L17>P}$). Although not supporting static adhesion of RAW117 cells, monomeric RGD-peptides mediated adhesion stabilization of H10 cells but not L17 or P cells. Integrin $\rm\alpha\sb{v}\beta\sb3$ was found to be involved in stabilizing H10 cell adhesion to vitronectin, (GRGDS)$\sb4$, monomeric RGD-peptide R1, and liver sinusoidal endothelial cell-derived extracellular matrix.^ This study is the first to provide evidence that integrin $\rm\alpha\sb{v}\beta\sb3$ is differentially expressed in liver-metastatic lymphoma cells and involved in differential adhesion of these cells. The results indicate that strong static adhesion and especially the unique hydrodynamic adhesion of RAW117 cells to the RGD-containing substrates correlate with liver-metastatic potentials. Thus, integrin $\rm\alpha\sb{v}\beta\sb3$ may play an important role in liver-preferential metastasis of RAW117 large-cell lymphoma cells. ^