New Research In
Physical Sciences
Social Sciences
Featured Portals
Articles by Topic
Biological Sciences
Featured Portals
Articles by Topic
- Agricultural Sciences
- Anthropology
- Applied Biological Sciences
- Biochemistry
- Biophysics and Computational Biology
- Cell Biology
- Developmental Biology
- Ecology
- Environmental Sciences
- Evolution
- Genetics
- Immunology and Inflammation
- Medical Sciences
- Microbiology
- Neuroscience
- Pharmacology
- Physiology
- Plant Biology
- Population Biology
- Psychological and Cognitive Sciences
- Sustainability Science
- Systems Biology
Variation in gene expression patterns in follicular lymphoma and the response to rituximab
-
Contributed by David Botstein
Abstract
Analysis of the patterns of gene expression in follicular lymphomas from 24 patients suggested that two groups of tumors might be distinguished. All patients, whose biopsies were obtained before any treatment, were treated with rituximab, a monoclonal antibody directed against the B cell antigen, CD20. Gene expression patterns in the tumors that subsequently failed to respond to rituximab appeared more similar to those of normal lymphoid tissues than to gene expression patterns of tumors from rituximab responders. These findings suggest the possibility that the response of follicular lymphoma to rituximab treatment may be predicted from the gene expression pattern of tumors.
Follicular non-Hodgkin's lymphoma (NHL) is an indolent B cell malignancy with an annual incidence exceeding 10,000 cases in the United States. Although follicular lymphoma (FL) is frequently responsive to treatment, therapy is very rarely, if ever, curative. Rituximab, a chimeric IgG1 monoclonal antibody directed at the B cell antigen CD20, has become a mainstay of treatment for low-grade NHL; >400,000 patients worldwide have been treated with rituximab. Phase II trials of rituximab in patients with refractory or relapsed low grade or follicular NHL demonstrated a 50% response rate (1).
Despite this extensive clinical experience, the mechanism of action of rituximab remains unclear, as does the nature of resistance (2). Among the proposed mechanisms are antibody-dependent cell-mediated cytotoxicity (3), complement-mediated cytotoxicity (4), and direct cytotoxicity through modulating CD20 function (5–7). The association with resistance to rituximab treatment of a low-affinity variant of the Fc receptor (8) is suggestive of an immune mechanism, and remains the only plausible hint about the nature of resistance.
In this study, we examined whether gene expression profiling using cDNA microarrays could reveal biological diversity among follicular lymphomas and, more specifically, whether gene expression patterns in tumors might predict sensitivity to rituximab treatment.
Materials and Methods
Patient Characteristics.
Patients were included in this study based on the availability of freshly frozen lymph node biopsy material containing enough mRNA to allow cDNA microarray analysis. Only patients with samples that had been obtained before any systemic therapy were included. In all cases the pathological diagnosis was follicular non-Hodgkin's lymphoma [follicular small cleaved (grade 1), follicular mixed (grade 2), or follicular large cell (grade 3) histology]. Each patient received rituximab treatment with documentation of clinical outcome. In all cases, biopsy and pathology review were performed at Stanford University Medical Center. Rituximab treatment was administered either at Stanford University Medical Center or by an outside oncologist.
Microarray Procedures.
Freshly frozen lymph node samples were obtained from patients who underwent excisional biopsy at Stanford University Medical Center between 1984 and 1997, who subsequently received rituximab between 1994 and 2000 and whose clinical response to rituximab treatment had been recorded. Tonsil and spleen samples were similarly obtained from patients treated at Stanford University Medical Center in 2000 or 2001. Biopsy samples were stored frozen in optimal cutting temperature compound. Poly(A)+ mRNA was obtained from biopsy samples after homogenization of tissue with the FASTTRACK 2.0 kit (Invitrogen). An experimental cDNA probe incorporating Cy5 dye was generated from mRNA from malignant and normal lymphoid tissues; a common reference cDNA probe incorporating Cy3 dye was from mRNA derived from a panel of cell lines and probes were hybridized to cDNA microarrays as described (9, 10). Two types of microarrays were used. Some experiments in this study used Stanford Human arrays comprised of 38,431 DNA spots of 38,276 unique cDNA clones, representing ≈31,139 unique Unigene clusters of which 16,152 correspond to unique named genes. Some experiments were conducted with lymphochip (LC) microarrays comprised of 37,632 DNA spots with 32,876 unique cDNA clones, representing ≈17,622 Unigene clusters of which at least 10,250 are unique named genes. More detailed information regarding microarray methods, and data selection, and analysis, as well as searchable figures and microarray data files, can be found at http://genome-www.stanford.edu/rituximab.
Statistical Analysis.
Before statistical analysis, individual data points were median centered for each cDNA clone and filtered for data quality and signal at least 2-fold above the median in two or three of the samples in each data set, as described in the web supplement. Agglomerative hierarchical cluster analysis was applied to the gene axis and to the sample axis as described (11). Hierarchical cluster analysis of LC data revealed a technical artifact that resulted in samples segregating by the date of the experiment. Further investigation revealed that this artifact was likely caused by differences in the calibration of the two scanners used to read the arrays. Singular value decomposition was used to remove the pattern corresponding to this artifact before analysis (12) after missing data were estimated by using a K-nearest neighbors (KNN) impute algorithm with 12 nearest values (13). Supervised analysis taking into account known outcome to rituximab treatment was performed by using Wilcoxon rank sum test to generate a rank list of genes whose corresponding mRNA levels differ significantly in rituximab responders versus nonresponders (14).
Results
Patient Characteristics.
Tumor samples from 24 patients were analyzed in this study. No significant differences in age or treatment history were observed between responders and nonresponders (Table 1). All patients except one received at least one course of chemotherapy before receiving rituximab (range, 0–6 prior courses). One patient had received a shared anti-idiotype antibody as the sole prior treatment. For the treatment of FL, a course of rituximab typically consists of four weekly infusions of 375 mg/m2; 21 of 24 patients received this dosing regimen, including all of the patients in the nonresponder group. In the partial response/complete response (PR/CR) group, two patients received a single dose of 375 mg/m2 and one patient received eight weekly doses of 375 mg/m2 with a documented near CR after 6 weeks of treatment. The overall response rate (CR+PR) for the patients in this study was 50%, which is similar to the overall response rate of 60% demonstrated for FL in the pivotal Phase II trial (1). Five patients achieved CR in response to rituximab, seven patients experienced PR, and 12 patients had no or minimal response.
Gene Expression Patterns Identify Two Subtypes of Follicular Lymphoma.
An overview of the gene expression patterns from FL patients was generated by hierarchical cluster analysis of data from the first 16 FL patients analyzed together with samples from nonmalignant tonsil and spleen (Fig. 1A). The hierarchical cluster algorithm arranges tissue samples based on the degree of similarity in their gene expression patterns (11). In addition to tumor material obtained before treatment, we included tissue samples from later biopsies for four patients for whom such material was available. These repeat biopsies were obtained 2.5–6.5 years after the initial biopsy but before rituximab treatment. In all of the repeat biopsy cases, patients received one or more courses of chemotherapy between the initial and later biopsies and the histopathological diagnosis was the same for both biopsy samples.
To focus on genes differentially expressed in the samples, expression data from 2,037 unique genes whose expression differed 2-fold in at least three arrays were used for analysis. The dendrogram shown in Fig. 1A shows that the FL samples could be divided into two groups, with a subset of the FL samples exhibiting a gene expression pattern more similar to nonmalignant lymphoid tissues than to the other FL subtype. Hierarchical cluster analysis of the genes and tissue samples that contribute to the subclassification of FL samples are shown in Fig. 1B. The two subtypes of FL display gene expression patterns similar to or opposite those of nonmalignant spleen. It should be emphasized that the observed patterns of gene expression reflect all of the cells in the tumor, not only the tumor cells themselves. Interpretation of expression differences may include the differential presence or absence of other cell types, as has been found for breast cancer (9).
Rituximab Nonresponders Display Gene Expression Patterns Characteristic of Normal Lymphoid Tissue.
Having divided the FL samples into two groups based on gene expression patterns, we sought to determine whether these subtypes correspond to clinical differences in individual patients, in particular, in response to rituximab treatment. We found that the rituximab nonresponders (orange in Fig. 1) were disproportionately distributed between the two FL clusters. Most of the rituximab nonresponders clustered with normal tonsil and spleen tissue (P < 0.005, Fisher's exact test). These findings indicate that there is biological diversity in FL lymph nodes from different patients before rituximab treatment and that the specific gene expression patterns defining the two groups described may be useful in predicting outcome to rituximab treatment.
Two features of the dendrogram in Fig. 1 indicate that the clustering of tissue samples based on gene expression patterns is reflective of biologically relevant similarities between samples. First, pairs of normal spleen and tonsil tissue cluster very closely together. Second, FL lymph node biopsies from the same individual cluster together in three of the four cases (see black bars in Fig. 1A) despite the passage of time and treatment with systemic chemotherapy between biopsies. An exception was the two samples obtained from patient 13 which did not cluster together (Fig. 1A, black arrows) with the pretreatment biopsy clustering with the responders' portion of the dendrogram and the later sample displaying a gene expression pattern more similar to that of rituximab nonresponders and normal lymphoid tissue. The latter sample was predictive of the patient's actual outcome after rituximab treatment; the patient did not respond. This observation suggests that the gene expression phenotype of FL can change over time.
Genes with Significantly Different Expression In Rituximab Responders Versus Nonresponders.
To better understand differences in gene expression in involved lymph nodes from rituximab responders versus nonresponders, we used supervised statistical analysis to determine which genes had the most significant differences in expression between the two groups. Having observed that differences in gene expression in FL subtypes correlate with expression patterns in normal lymphoid tissues, samples from the original 16 patients and an additional 8 patient samples were analyzed on LC microarrays enriched for genes expressed in lymphoid cells and genes known or suspected to be important in the immune response or cancer (10). A list of genes whose expression differed between responders and nonresponders with a P value <0.005, as determined by Wilcoxon rank-sum test (14), was chosen. In cases where expression of a given gene was measured on both Stanford Human and LC microarrays, genes were included if the P value was less than 0.005 in either the Stanford Human or LC data set and a cDNA from the same Unigene cluster displayed a P value of <0.05 in the other dataset. Genes that were measured in only one data set were included if the P value is <0.005. The results of this analysis are presented in Tables 2 and 3. By these criteria, 71 genes had significantly higher expression in rituximab nonresponders versus responders; 53 named genes are on this list and 35 of the Unigene clusters represented were measured on both arrays. A total of 27 genes were more highly expressed by these criteria in rituximab responders; of these, 11 are named genes and 5 were measured on both arrays. Of note, CD20 expression was well measured and did not correlate with outcome to rituximab treatment.
Many of the genes with higher expression in tissue from rituximab nonresponders appear to be involved in cellular immune response and inflammation, specifically those encoding cytokine, tumor necrosis factor, and T cell receptor signaling, and complement proteins (Table 4). Given the limitations of available experimental evidence, this simplistic classification cannot account for the complexity of function or regulation of these genes; perhaps, mediators of cellular immune response, such T cells, macrophages, monocytes, and natural killer cells, may be relatively more abundant and/or more active in lymph node tissue of rituximab nonresponders. Several mRNAs for proteins involved in the complement cascade are more abundant in nonresponders; the complement cascade can mediate direct cytotoxicity and/or act to stimulate cell-mediated cytotoxicity by opsonization of target cells. We were unable to identify a compelling physiological link among the 11 named genes with significantly increased expression in rituximab responders.
Discussion
Variation in gene expression patterns of FL lymph nodes suggests that it may be useful to divide these tumors into two groups. We suppose that these subtypes are reflective of basic biological differences between FL samples from different individuals. As mentioned before, the differences could be the result, at least in part, of the differential admixture of cells other than the malignant B cells that comprise most of the tumor. Thus, although it is quite possible that the subtypes identified simply represent molecular heterogeneity in the malignant cells, the finding that one subtype displays an expression pattern more similar to nonmalignant lymphoid tissue raises the possibility that the subtypes may, in fact, result from interactions between the malignant B cells and other cell types in the host. On the basis of the genes differing most in their expression between the two subtypes of tumors, it is possible that the differences may involve either the participation or the activities of cells involved in host immune response to the tumors. In any case, on the basis of this preliminary study, the subtype (and the response to rituximab) appears to be determined by the time of initial diagnosis and distinctive features of the gene expression pattern of a given individual's lymphoma appear to be recognizably retained over an interval of years.
Gene expression profiling studies of various malignancies have identified previously unappreciated subtypes within accepted pathological diagnoses and, in some cases, the prognoses of these subtypes varied significantly (10, 15–17). Several models of rituximab resistance are consistent with our findings. If rituximab nonresponders are capable of mounting a cellular immune response to their tumor before therapy, the FL cells may have been selected to evade this immune response and this may translate into an ability to evade any antibody-dependent cell-mediated cytotoxicity (ADCC) response subsequently triggered by rituximab. Alternatively, it has been proposed that FL cells may be able to recruit a cellular microenvironment similar to that in the germinal center, which facilitates the growth of the lymphoma cells (18). Following this line of reasoning, the list of genes with higher expression in rituximab nonresponders may indicate that these tumors may more effectively induce a growth stimulating microenvironment and this more favorable microenvironment may make these cells less sensitive to killing by rituximab through ADCC, complement-mediated cytotoxicity, or induction of apoptosis.
Prognostic tests allowing the prediction of response to specific therapeutic agents would be of great utility, as they would allow physicians to choose therapy to maximize patient benefit while avoiding unnecessary delay, toxicity, and expense. Our findings suggest that it may be possible, at the time of the initial diagnosis of follicular lymphoma, to predict whether patients will later respond to rituximab treatment for relapsed or refractory disease. The data described here suggest that analysis of gene expression patterns in tumors may allow prediction of the sensitivity of tumors to particular antitumor agents and elucidate the biology underlying resistance to a given therapy. In fact, the possibility of dividing FL into molecularly distinct subtypes may confer additional useful information on differences among these tumors, such as response to therapies other than rituximab. Clearly, our findings require validation on a larger, independent patient cohort. A better understanding of the nature of resistance to rituximab treatment may allow more effective use of this powerful agent. For instance, it is possible that patients who are resistant to standard rituximab dosing may benefit from rituximab administered at a different dose or schedule, or may respond to rituximab in combination with other agents. The ability to prospectively identify patients who are relatively resistant to rituximab would facilitate clinical trials to determine the optimal treatment for such patients.
Acknowledgments
We thank members of the Brown and Levy laboratories for advice and thoughtful comments on data analysis, and Sandra Horning for thought-provoking discussions. LC microarrays were a much appreciated gift from A. Alizadeh, M. Diehn, and A. Whitney. This work was supported by National Cancer Institute Grants CA85129 (to P.O.B.) and CA77097 (to D.B.), National Institutes of Health Grants CA33399 and CA34233 (to R.L.), and the Howard Hughes Medical Institute. S.P.B. received fellowship support from Amgen and a Public Health Service training grant, and is a Howard Hughes Medical Institute physician postdoctoral fellow. O.A. is an Alfred P. Sloan and U.S. Department of Energy Postdoctoral Fellow in Computational Biology (Grant DE-FG03-99ER62836) and a National Human Genome Research Institute Individual Mentored Scientist Development Awardee in Genomic Research and Analysis (National Institutes of Health Grant K01 HG00038-01). P.O.B. is an Investigator of the Howard Hughes Medical Institute. R.L. is an American Cancer Society Clinical Research Professor.
Abbreviations
- FL,
- follicular lymphoma;
- LC,
- lymphochip
- Accepted December 24, 2002.
- Copyright © 2003, The National Academy of Sciences
References
- ↵
- ↵
- ↵
- ↵
- Golay J,
- Zaffaroni L,
- Vaccari T,
- Lazzari M,
- Borleri G M,
- Bernasconi S,
- Tedesco F,
- Rambaldi A,
- Introna M
- ↵
-
- Shan D,
- Ledbetter J A,
- Press O W
- ↵
- ↵
- Cartron G,
- Dacheux L,
- Salles G,
- Solal-Celigny P,
- Bardos P,
- Colombat P,
- Watier H
- ↵
- ↵
- ↵
- Eisen M B,
- Spellman P T,
- Brown P O,
- Botstein D
- ↵
- Alter O,
- Brown P O,
- Botstein D
- ↵
- Troyanskaya O,
- Cantor M,
- Sherlock G,
- Brown P,
- Hastie T,
- Tibshirani R,
- Botstein D,
- Altman R B
- ↵
- Troyanaskaya O G,
- Garber M E,
- Brown P O,
- Botstein D,
- Altman R B
- ↵
-
- Sorlie T,
- Perou C M,
- Tibshirani R,
- Aas T,
- Geisler S,
- Johnsen H,
- Hastie T,
- Eisen M B,
- van de Rijn M,
- Jeffrey S S,
- et al.
- ↵
- Garber M E,
- Troyanskaya O G,
- Schluens K,
- Petersen S,
- Thaesler Z,
- Pacyna-Gengelbach M,
- van de Rijn M,
- Rosen G D,
- Perou C M,
- Whyte R I,
- et al.
- ↵
- Ghia P,
- Granziero L,
- Chilosi M,
- Caligaris-Cappio F
Citation Manager Formats
More Articles of This Classification
Biological Sciences
Related Content
- No related articles found.
Cited by...
- Anti-CD20-interleukin-21 fusokine targets malignant B cells via direct apoptosis and NK-cell-dependent cytotoxicity
- Indolent lymphoma: follicular lymphoma and the microenvironment--insights from gene expression profiling
- CALGB 150905 (Alliance): Rituximab Broadens the Antilymphoma Response by Activating Unlicensed NK Cells
- The reliability of immunohistochemical analysis of the tumor microenvironment in follicular lymphoma: a validation study from the Lunenburg Lymphoma Biomarker Consortium
- Integration of Host-Related Signatures with Cancer Cell-Derived Predictors for the Optimal Management of Anticancer Chemotherapy
- A Pharmacologically Active Monoclonal Antibody against the Human Melanocortin-4 Receptor: Effectiveness After Peripheral and Central Administration
- Pathways of the Heart
- Notch Ligands Expressed by Follicular Dendritic Cells Protect Germinal Center B Cells from Apoptosis
- High Numbers of Tumor-Infiltrating Programmed Cell Death 1-Positive Regulatory Lymphocytes Are Associated With Improved Overall Survival in Follicular Lymphoma
- Paraffin-based 6-gene model predicts outcome in diffuse large B-cell lymphoma patients treated with R-CHOP
- Prognostic significance of host immune gene polymorphisms in follicular lymphoma survival
- Gene-Expression and Immunohistochemical Study of Specific T-Cell Subsets and Accessory Cell Types in the Transformation and Prognosis of Follicular Lymphoma
- Immunohistochemical Patterns of Reactive Microenvironment Are Associated With Clinicobiologic Behavior in Follicular Lymphoma Patients
- High numbers of tumor-infiltrating FOXP3-positive regulatory T cells are associated with improved overall survival in follicular lymphoma
- R-CHOP strikes again with survival benefit in follicular lymphoma
- Molecular Pathogenesis of Follicular Lymphoma: A Cross Talk of Genetic and Immunologic Factors
- Optimizing Rituximab in B-Cell Lymphoma
- Importance of Cellular Microenvironment and Circulatory Dynamics in B Cell Immunotherapy
- Unique Toxicities and Resistance Mechanisms Associated with Monoclonal Antibody Therapy
- The Role of B Cell-Mediated T Cell Costimulation in the Efficacy of the T Cell Retargeting Bispecific Antibody BIS20x3
- Bcl-2 Overexpression Leads to Increases in Suppressor of Cytokine Signaling-3 Expression in B Cells and De novo Follicular Lymphoma
- Stage I and II Follicular Non-Hodgkin's Lymphoma: Long-Term Follow-Up of No Initial Therapy
- Molecular Diagnostics
- Genomic Perspective and Cancer