Relation of gene expression phenotype to immunoglobulin mutation genotype in B cell chronic lymphocytic leukemia.

The most common human leukemia is B cell chronic lymphocytic leukemia (CLL), a malignancy of mature B cells with a characteristic clinical presentation but a variable clinical course. The rearranged immunoglobulin (Ig) genes of CLL cells may be either germ-line in sequence or somatically mutated. Lack of Ig mutations defined a distinctly worse prognostic group of CLL patients raising the possibility that CLL comprises two distinct diseases. Using genomic-scale gene expression profiling, we show that CLL is characterized by a common gene expression "signature," irrespective of Ig mutational status, suggesting that CLL cases share a common mechanism of transformation and/or cell of origin. Nonetheless, the expression of hundreds of other genes correlated with the Ig mutational status, including many genes that are modulated in expression during mitogenic B cell receptor signaling. These genes were used to build a CLL subtype predictor that may help in the clinical classification of patients with this disease.

A second iron-regulatory system in yeast independent of Aft1p.

Iron homeostasis in the yeast Saccharomyces cerevisiae is regulated at the transcriptional level by Aft1p, which activates the expression of its target genes in response to low-iron conditions. The yeast genome contains a paralog of AFT1, which has been designated AFT2. To establish whether AFT1 and AFT2 have overlapping functions, a mutant containing a double aft1Deltaaft2Delta deletion was generated. Growth assays established that the single aft2Delta strain exhibited no iron-dependent phenotype. However, the double-mutant aft1Deltaaft2Delta strain was more sensitive to low-iron growth conditions than the single-mutant aft1Delta strain. A mutant allele of AFT2 (AFT2-1(up)), or overexpression of the wild-type AFT2 gene, led to partial complementation of the respiratory-deficient phenotype of the aft1Delta strain. The AFT2-1(up) allele also increased the uptake of (59)Fe in an aft1Delta strain. DNA microarrays were used to identify genes regulated by AFT2. Some of the AFT2-regulated genes are known to be regulated by Aft1p; however, AFT2-1(up)-dependent activation was independent of Aft1p. The kinetics of induction of two genes activated by the AFT2-1(up) allele are consistent with Aft2p acting as a direct transcriptional factor. Truncated forms of Aft1p and Aft2p bound to a DNA duplex containing the Aft1p binding site in vitro. The wild-type allele of AFT2 activated transcription in response to growth under low-iron conditions. Together, these data suggest that yeast has a second regulatory pathway for the iron regulon, with AFT1 and AFT2 playing partially redundant roles.

Comprehensive mutational analysis of the Moloney murine leukemia virus envelope protein.

The envelope (Env) protein of Moloney murine leukemia virus is the primary mediator of viral entry. We constructed a large pool of insertion mutations in the env gene and analyzed the fitness of each mutant in completing two critical steps in the virus life cycle: (i) the expression and delivery of the Env protein to the cell surface during virion assembly and (ii) the infectivity of virions displaying the mutant proteins. The majority of the mutants were poorly expressed at the producer cell surface, suggesting folding defects due to the presence of the inserted residues. The mutants with residual infectivity had insertions either in the amino-terminal signal sequence region, two disulfide-bonded loops in the receptor binding domain, discrete regions of the carboxy-terminal region of the surface subunit (SU), or the cytoplasmic tail. Insertions that allowed the mutants to reach the cell surface but not to mediate detectable infection were located within the amino-terminal sequence of the mature Env, within the SU carboxy-terminal region, near putative receptor binding residues, and throughout the fusion peptide. Independent analysis of select mutants in this group allowed more precise identification of the defect in Env function. Mapping of mutant phenotypes to a structural model of the receptor-binding domain provides insights into the protein's functional organization. The high-resolution functional map reported here will be valuable for the engineering of the Env protein for a variety of uses, including gene therapy.

Diversity of gene expression in adenocarcinoma of the lung.

The global gene expression profiles for 67 human lung tumors representing 56 patients were examined by using 24,000-element cDNA microarrays. Subdivision of the tumors based on gene expression patterns faithfully recapitulated morphological classification of the tumors into squamous, large cell, small cell, and adenocarcinoma. The gene expression patterns made possible the subclassification of adenocarcinoma into subgroups that correlated with the degree of tumor differentiation as well as patient survival. Gene expression analysis thus promises to extend and refine standard pathologic analysis.

Three cell wall mannoproteins facilitate the uptake of iron in Saccharomyces cerevisiae.

Analysis of iron-regulated gene expression in Saccharomyces cerevisiae using cDNA microarrays has identified three putative cell wall proteins that are directly regulated by Aft1p, the major iron-dependent transcription factor in yeast. FIT1, FIT2, and FIT3 (for facilitator of iron transport) were more highly expressed in strains grown in low concentrations of iron and in strains in which AFT1-1(up), a constitutively active allele of AFT1, was expressed. Northern blot analysis confirmed that FIT1, FIT2, and FIT3 mRNA transcript levels were increased 60-230-fold in response to iron deprivation in an Aft1p-dependent manner. Fit1p was localized exclusively to the cell wall by indirect immunofluorescence. Deletion of the FIT genes, individually or in combination, resulted in diminished uptake of iron bound to the siderophores ferrioxamine B and ferrichrome, without diminishing the uptake of ferric iron salts, or the siderophores triacetylfusarinine C and enterobactin. FIT-deletion strains exhibited increased expression of Aft1p target genes as measured by a FET3-lacZ reporter gene or by Arn1p Western blotting, indicating that cells respond to the absence of FIT genes by up-regulating systems of iron uptake. Aft1p activation in FIT-deleted strains occurred when either ferrichrome or ferric salts were used as sources of iron during growth, suggesting that the FIT genes enhance uptake of iron from both sources. Enzymatic digestion of the cell wall resulted in the release of significant amounts of iron from cells, and the relative quantity of iron released was reduced in FIT-deletion strains. Fit1p, Fit2p, and Fit3p may function by increasing the amount of iron associated with the cell wall and periplasmic space.

Microarrays in primary breast cancer--lessons from chemotherapy studies.

Current development in molecular techniques has extended the opportunities to explore genetic alterations in malignant tissue. There is a need to improve prognostication and, in particular, to understand the mechanisms of treatment resistance in different tumours. Gene analyses by microarrays allow concomitant analyses of several genes in concert, providing new opportunities for tumour classification and understanding of key biological disturbances. This paper outlines our continuing studies exploring prognostic and, we hope, predictive factors in breast cancer therapy.

Gene expression patterns of breast carcinomas distinguish tumor subclasses with clinical implications.

The purpose of this study was to classify breast carcinomas based on variations in gene expression patterns derived from cDNA microarrays and to correlate tumor characteristics to clinical outcome. A total of 85 cDNA microarray experiments representing 78 cancers, three fibroadenomas, and four normal breast tissues were analyzed by hierarchical clustering. As reported previously, the cancers could be classified into a basal epithelial-like group, an ERBB2-overexpressing group and a normal breast-like group based on variations in gene expression. A novel finding was that the previously characterized luminal epithelial/estrogen receptor-positive group could be divided into at least two subgroups, each with a distinctive expression profile. These subtypes proved to be reasonably robust by clustering using two different gene sets: first, a set of 456 cDNA clones previously selected to reflect intrinsic properties of the tumors and, second, a gene set that highly correlated with patient outcome. Survival analyses on a subcohort of patients with locally advanced breast cancer uniformly treated in a prospective study showed significantly different outcomes for the patients belonging to the various groups, including a poor prognosis for the basal-like subtype and a significant difference in outcome for the two estrogen receptor-positive groups.

RERG is a novel ras-related, estrogen-regulated and growth-inhibitory gene in breast cancer.

Using microarray analysis, we identified a unique ras superfamily gene, termed RERG (ras-related and estrogen-regulated growth inhibitor), whose expression was decreased or lost in a significant percentage of primary human breast tumors that show a poor clinical prognosis. Importantly, high RERG expression correlated with expression of a set of genes that define a breast tumor subtype that is estrogen receptor-positive and associated with a slow rate of tumor cell proliferation and a favorable prognosis for these cancer patients. RERG mRNA expression was induced rapidly in MCF-7 cells stimulated by beta-estradiol and repressed by tamoxifen treatment. Like Ras, RERG protein exhibited intrinsic GDP/GTP binding and GTP hydrolysis activity. Unlike Ras proteins, RERG lacks a known recognition signal for COOH-terminal prenylation and was localized primarily in the cytoplasm. Expression of RERG protein in MCF-7 breast carcinoma cells resulted in a significant inhibition of both anchorage-dependent and anchorage-independent growth in vitro and inhibited tumor formation in nude mice. These features of RERG are strikingly different from most Ras superfamily GTP-binding pro-teins and suggest that the loss of RERG expression may contribute to breast tumorigenesis.

Use of cDNA microarrays to analyze dioxin-induced changes in human liver gene expression.

One mechanism by which cells adapt to environmental changes is by altering gene expression. Here, we have used cDNA microarrays to identify genes whose expression is altered by exposure to the environmental contaminant 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). The goal of our study was to enhance our understanding of toxicity mediated through the pathway by which TCDD stimulates gene expression. To model this toxicity response, we exposed human hepatoma (HepG2) cells to TCDD (10 nM for 18 h) and analyzed mRNA by two-color fluorescent hybridization to cDNA sequences immobilized on glass microscope slides (2.5 x 7.5 cm) covering a surface area of 2.25 cm(2). We analyzed approximately one-third of the genes expressed in HepG2 cells and found that TCDD up- or down-regulates 112 genes two-fold or more. Most changes are relatively subtle (two- to four-fold). We verified the regulation of protooncogene cot, XMP, and human enhancer of filamentation-1 (HEF1), genes involved in cellular proliferation, as well as metallothionein, plasminogen activator inhibitor (PAI1), and HM74, genes involved in cellular signaling and regeneration. To characterize the response in more detail, we performed time-course, dose-dependence studies, and cycloheximide experiments. We observed direct and indirect responses to TCDD implying that adaptation to TCDD (and other related environmental stimuli) is substantially more complex than we previously realized.

Regulation of CSF1 promoter by the SWI/SNF-like BAF complex.

The mammalian BAF complex regulates gene expression by modifying chromatin structure. In this report, we identify 80 genes activated and 2 genes repressed by the BAF complex in SW-13 cells. We find that prior binding of NFI/CTF to the NFI/CTF binding site in CSF1 promoter is required for the recruitment of the BAF complex and the BAF-dependent activation of the promoter. Furthermore, the activation of the CSF1 promoter requires Z-DNA-forming sequences that are converted to Z-DNA structure upon activation by the BAF complex. The BAF complex facilitates Z-DNA formation in a nucleosomal template in vitro. We propose a model in which the BAF complex promotes Z-DNA formation which, in turn, stabilizes the open chromatin structure at the CSF1 promoter.

Haa1, a protein homologous to the copper-regulated transcription factor Ace1, is a novel transcriptional activator.

The Saccharomyces cerevisiae genome contains a predicted gene, YPR008w, homologous to the gene encoding the copper-activated transcription factor Ace1. The product of the YPR008w gene, designated Haa1, regulates the transcription of a set of yeast genes, many of which encode membrane proteins. Two main target genes of Haa1 are the multidrug resistance gene YGR138c and the YRO2 homolog to the plasma membrane Hsp30. Haa1 is localized to the nucleus. Haa1-induced expression of YGR138c and YRO2 appears to be direct. Induction of HAA1 using a GAL1/HAA1 fusion gene resulted in rapid galactose-induced expression of both HAA1 and target genes. Although Haa1 has a sequence very similar to the Cu-activated DNA binding domain of Ace1, expression of Haa1 target genes was found to be independent of the copper status of cells. Haa1 does not exhibit metalloregulation in cells incubated with a range of transition metal salts. Haa1 does not exhibit any cross-talk with Ace1. Overexpression of Haa1 does not compensate for cells lacking a functional Ace1. The lack of metalloregulation of Haa1 despite the strong sequence similarity to the copper regulatory domain of Ace1 is discussed.

Comparative genome-scale analysis of gene expression profiles in T cell lymphoma cells during malignant progression using a complementary DNA microarray.

Using a cDNA microarray, we compared the expression of approximately 8000 genes between two unique, clonally related T cell lines derived from different stages of a progressive T cell lymphoma involving skin. A total of 180 genes was found to be differentially expressed at the RNA level by a factor of fivefold or greater. Compared with the cells from the earlier, clinically indolent stage of the lymphoma, 56 genes were up-regulated, whereas 124 genes were down-regulated in the cells from the advanced, clinically aggressive stage lymphoma. The functions of approximately 65% of these genes are currently unknown. The 22 genes with a known function that were up-regulated in the advanced lymphoma cells included several genes involved in promotion of cell proliferation and survival as well as drug resistance. The 42 functionally characterized genes that were down-regulated in the advanced lymphoma cells included negative regulators of cell activation and cell cycle, and mediators of cell adhesion, apoptosis, and genome integrity. The differential expression identified by the cDNA microarray analysis was confirmed for selected genes by reverse transcription-polymerase chain reaction and Northern blotting. The identified differences in gene expression may be related to the differences in behavior between the early and advanced stages of the T cell lymphoma and point to directions for further investigations into mechanisms of lymphoma progression.

Modulation of cellular and viral gene expression by the latency-associated nuclear antigen of Kaposi's sarcoma-associated herpesvirus.

Kaposi's sarcoma-associated herpesvirus (KSHV), also called human herpesvirus 8 (HHV-8), is the likely etiological agent of Kaposi's sarcoma and primary effusion lymphoma. Common to these malignancies is that tumor cells are latently infected with KSHV. Viral gene expression is limited to a few genes, one of which is the latency-associated nuclear antigen (LANA), the product of ORF73. Examination of the primary sequence of LANA reveals some structural features reminiscent of transcription factors, leading us to hypothesize that LANA may regulate viral and cellular transcription during latency. In reporter gene-based transient transfection assays, we found that LANA can have either positive or negative effects on gene expression. While expression of a reporter gene from several synthetic promoters was increased in the presence of LANA, expression from the human immunodeficiency virus (HIV) long terminal repeat (LTR)-and from NF-kappaB-dependent reporter genes-was reduced by LANA expression. In addition, the promoter of KSHV ORF73 itself is activated up to 5.5-fold by LANA. This autoregulation may be important in tumorigenesis, because two other genes (v-cyclin and v-FLIP) with likely roles in cell growth and survival are also controlled by this element. To identify cellular genes influenced by LANA, we employed cDNA array-based expression profiling. Six known genes (and nine expressed sequence tags) were found to be upregulated in LANA-expressing cell lines. One of these, Staf-50, is known to inhibit expression from the HIV LTR; most of the other known genes are interferon inducible, although the interferon genes themselves were not induced by LANA. These data demonstrate that LANA expression has effects on cellular and viral gene expression. We suggest that, whether direct or indirect in origin, these effects may play important roles in the pathobiology of KSHV infection.

The Stanford Microarray Database.

The Stanford Microarray Database (SMD) stores raw and normalized data from microarray experiments, and provides web interfaces for researchers to retrieve, analyze and visualize their data. The two immediate goals for SMD are to serve as a storage site for microarray data from ongoing research at Stanford University, and to facilitate the public dissemination of that data once published, or released by the researcher. Of paramount importance is the connection of microarray data with the biological data that pertains to the DNA deposited on the microarray (genes, clones etc.). SMD makes use of many public resources to connect expression information to the relevant biology, including SGD [Ball,C.A., Dolinski,K., Dwight,S.S., Harris,M.A., Issel-Tarver,L., Kasarskis,A., Scafe,C.R., Sherlock,G., Binkley,G., Jin,H. et al. (2000) Nucleic Acids Res., 28, 77-80], YPD and WormPD [Costanzo,M.C., Hogan,J.D., Cusick,M.E., Davis,B.P., Fancher,A.M., Hodges,P.E., Kondu,P., Lengieza,C., Lew-Smith,J.E., Lingner,C. et al. (2000) Nucleic Acids Res., 28, 73-76], Unigene [Wheeler,D.L., Chappey,C., Lash,A.E., Leipe,D.D., Madden,T.L., Schuler,G.D., Tatusova,T.A. and Rapp,B.A. (2000) Nucleic Acids Res., 28, 10-14], dbEST [Boguski,M.S., Lowe,T.M. and Tolstoshev,C.M. (1993) Nature Genet., 4, 332-333] and SWISS-PROT [Bairoch,A. and Apweiler,R. (2000) Nucleic Acids Res., 28, 45-48] and can be accessed at http://genome-www.stanford.edu/microarray.

New components of a system for phosphate accumulation and polyphosphate metabolism in Saccharomyces cerevisiae revealed by genomic expression analysis.

The PHO regulatory pathway is involved in the acquisition of phosphate (P(i)) in the yeast Saccharomyces cerevisiae. When extracellular P(i) concentrations are low, several genes are transcriptionally induced by this pathway, which includes the Pho4 transcriptional activator, the Pho80-Pho85 cyclin-CDK pair, and the Pho81 CDK inhibitor. In an attempt to identify all the components regulated by this system, a whole-genome DNA microarray analysis was employed, and 22 PHO-regulated genes were identified. The promoter regions of 21 of these genes contained at least one copy of a sequence that matched the Pho4 recognition site. Eight of these genes, PHM1-PHM8, had no previously defined function in phosphate metabolism. The amino acid sequences of PHM1 (YFL004w), PHM2 (YPL019c), PHM3 (YJL012c), and PHM4 (YER072w) are 32-56% identical. The phm3 and phm4 single mutants and the phm1 phm2 double mutant were each severely deficient in accumulation of inorganic polyphosphate (polyP) and P(i). The phenotype of the phm5 mutant suggests that PHM5 (YDR452w) is essential for normal catabolism of polyP in the yeast vacuole. Taken together, the results reveal important new features of a genetic system that plays a critical role in P(i) acquisition and polyP metabolism in yeast.

Molecular portraits of human breast tumours.

Human breast tumours are diverse in their natural history and in their responsiveness to treatments. Variation in transcriptional programs accounts for much of the biological diversity of human cells and tumours. In each cell, signal transduction and regulatory systems transduce information from the cell's identity to its environmental status, thereby controlling the level of expression of every gene in the genome. Here we have characterized variation in gene expression patterns in a set of 65 surgical specimens of human breast tumours from 42 different individuals, using complementary DNA microarrays representing 8,102 human genes. These patterns provided a distinctive molecular portrait of each tumour. Twenty of the tumours were sampled twice, before and after a 16-week course of doxorubicin chemotherapy, and two tumours were paired with a lymph node metastasis from the same patient. Gene expression patterns in two tumour samples from the same individual were almost always more similar to each other than either was to any other sample. Sets of co-expressed genes were identified for which variation in messenger RNA levels could be related to specific features of physiological variation. The tumours could be classified into subtypes distinguished by pervasive differences in their gene expression patterns.

Identification of the copper regulon in Saccharomyces cerevisiae by DNA microarrays.

In Saccharomyces cerevisiae, copper ions regulate gene expression through the two transcriptional activators, Ace1 and Mac1. Ace1 mediates copper-induced gene expression in cells exposed to stressful levels of copper salts, whereas Mac1 activates a subset of genes under copper-deficient conditions. DNA microarray hybridization experiments revealed a limited set of yeast genes differentially expressed under growth conditions of excess copper or copper deficiency. Mac1 activates the expression of six S. cerevisiae genes, including CTR1, CTR3, FRE1, FRE7, YFR055w, and YJL217w. Two of the last three newly identified Mac1 target genes have no known function; the third, YFR055w, is homologous to cystathionine gamma-lyase encoded by CYS3. Several genes that are differentially expressed in cells containing a constitutively active Mac1, designated Mac1(up1), are not direct targets of Mac1. Induction or repression of these genes is likely a secondary effect of cells because of constitutive Mac1 activity. Elevated copper levels induced the expression of the metallothioneins CUP1 and CRS5 and two genes, FET3 and FTR1, in the iron uptake system. Copper-induced FET3 and FTR1 expression arises from an indirect copper effect on cellular iron pools.

Ongoing immunoglobulin somatic mutation in germinal center B cell-like but not in activated B cell-like diffuse large cell lymphomas.

B cell diffuse large cell lymphoma (B-DLCL) is a heterogeneous group of tumors, based on significant variations in morphology, clinical presentation, and response to treatment. Gene expression profiling has revealed two distinct tumor subtypes of B-DLCL: germinal center B cell-like DLCL and activated B cell-like DLCL. In a separate study, we determined that B-DLCL can also be subdivided into two groups based on the presence or absence of ongoing Ig gene hypermutation. Here, we evaluated the correlation between these B-DLCL subtypes established by the two different methods. Fourteen primary B-DLCL cases were studied by gene expression profiling using DNA microarrays and for the presence of ongoing mutations in their Ig heavy chain gene. All seven cases classified as germinal center B cell-like DLCL by gene expression showed the presence of ongoing mutations in the Ig genes. Five of the seven cases classified by gene expression as activated B cell-like DLCL had no ongoing somatic mutations, whereas, in the remaining two cases, a single point mutation was observed in only 2 of 15 and 21 examined molecular clones of variable heavy (V(H)) chain gene, respectively. These two cases were distantly related to the rest of the activated B cell-like DLCL tumors by gene expression. Our findings validate the concept that lymphoid malignancies are derived from cells at discrete stages of normal lymphocyte maturation and that the malignant cells retain the genetic program of those normal cells.

Genomic views of the immune system*.

Genomic-scale experimentation aims to view biological processes as a whole, yet with molecular precision. Using massively parallel DNA microarray technology, the mRNA expression of tens of thousands of genes can be measured simultaneously. Mathematical distillation of this flood of gene expression data reveals a deep molecular and biological logic underlying gene expression programs during cellular differentiation and activation. Genes that encode components of the same multi-subunit protein complex are often coordinately regulated. Coordinate regulation is also observed among genes whose products function in a common differentiation program or in the same physiological response pathway. Recent application of gene expression profiling to the immune system has shown that lymphocyte differentiation and activation are accompanied by changes of hundreds of genes in parallel. The databases of gene expression emerging from these studies of normal immune responses will be used to interpret the pathological changes in gene expression that accompany autoimmunity, immune deficiencies, and cancers of immune cells.

Desferrioxamine-mediated iron uptake in Saccharomyces cerevisiae. Evidence for two pathways of iron uptake.

In the yeast Saccharomyces cerevisiae, uptake of iron is largely regulated by the transcription factor Aft1. cDNA microarrays were used to identify new iron and AFT1-regulated genes. Four homologous genes regulated as part of the AFT1-regulon (ARN1-4) were predicted to encode members of a subfamily of the major facilitator superfamily of transporters. These genes were predicted to encode proteins with 14 membrane spanning domains and were from 26 to 53% identical at the amino acid level. ARN3 is identical to SIT1, which is reported to encode a ferrioxamine B permease. Deletion of ARN3 did not prevent yeast from using ferrioxamine B as an iron source; however, deletion of ARN3 and FET3, a component of the high affinity ferrous iron transport system, did prevent uptake of ferrioxamine-bound iron and growth on ferrioxamine as an iron source. The siderophore-mediated transport system and the high affinity ferrous iron transport system were localized to separate cellular compartments. Epitope-tagged Arn3p was expressed in intracellular vesicles that co-sediment with the endosomal protein Pep12. In contrast, Fet3p was expressed on the plasma membrane and was digested by extracellular proteases. These data indicate that S. cerevisiae has two pathways for ferrrioxamine-mediated iron uptake, one occurring at the plasma membrane and the other occurring in an intracellular compartment.