|
|
|||||||
Malignant Lymphomas |
1 Molecular Pathology Programme, Centro Nacional de Investigaciones Oncológicas, Madrid
2 Pathology Department, Hospital Ramon y Cajal, Madrid
3 Genetics and Pathology Departments, Hospital Virgen de la Salud, Toledo
4 Pathology Department, MD Anderson International Spain, Madrid
5 Hematology Department, Hospital Doce de Octubre, Madrid
6 Pathology Department, Hospital Universitario de Getafe, Madrid, Spain
Correspondence: Miguel A. Piris, Centro Nacional de Investigaciones Oncológicas, Melchor Fernández Almagro, 3, Madrid 28029, Spain. E-mail:mapiris{at}cnio.es
| ABSTRACT |
|---|
|
|
|---|
Design and Methods: With the aim of identifying surrogate markers of somatic hypermutation, a combined investigation of IgVH mutational status and expression profiles of 93 samples from patients with small B-cell lymphoma was performed.
Results: The analysis identified an somatic hypermutation signature of genes involved in the regulation of gene transcription, DNA repair and replication, and chromosome maintenance. Eight of these genes were subjected to protein analysis using tissue microarrays, for a set of 118 cases. We found a clear link between RAD51C and CDK7 protein expression and somatic hypermutation status, in that positive expression of either marker was significantly associated with a mutated status (p<0.003). We also found that positive expression of TFDP1 and POLA was significantly associated with ongoing somatic hypermutation (p<0.001). To assess the potential clinical applicability of these somatic hypermutation markers, we studied a series of cases of mantle cell lymphoma included in a tissue microarray. The expression of RCC1 and CDK7, separately and together, was found to be significantly associated with longer overall survival.
Conclusions: An somatic hypermutation signature has been identified for different types of small B-cell lymphoma. This has a potential mechanistic and diagnostic value.
Key words: lymphomas, somatic hypermutation, profiling.
| Introduction |
|---|
|
|
|---|
In humans, SHM of Ig genes occurs at rates of 10–3 to 10–4 mutations per base pair per generation, which is up to six orders of magnitude greater than the spontaneous mutation rate of most other genes.2,3 The majority of changes are single base substitutions (most commonly transitions) in the Ig variable region, starting 150–200 base pairs downstream of the promoter and continuing about 1.5 kb downstream.4,5 Although mutations occur throughout the rearranged variable regions, RGYW and WRCY (R: A or G; Y: T or C; W: A or T) motifs are preferential targets (hotspots) for SHM. Higher-order structures or differences in local sequences may also play a role in targeting SHM to specific bases, as not all potential hotspots in the same region are affected by SHM.6,7 Likewise, the vast majority of genes are not affected by SHM although similar hotspots may be present. For example the Ig constant region gene, which does not exhibit SHM, is found only a few kilobases from the variable region genes.5
Activation-induced cytidine deaminase (AID) is expressed by activated B cells in the germinal centers of peripheral lymphoid organs, and is critical to the process of SHM. However, the mechanisms of SHM are not fully understood. For example, the molecules and processes responsible for DNA unwinding, which allow AID access to the DNA template, and the signaling pathways involved in the process are unknown. High AID expression is insufficient to explain SHM in chronic lymphocytic leukemia (CLL), since strong AID mRNA expression is associated with unmutated IgVH gene status.8
Besides the biological relevance of SHM, there is a broad base of evidence to suggest that IgVH SHM is a clinically significant phenomenon in various types of B-cell lymphoma, as demonstrated by studies reporting that cases of hypermutated CLL9 and splenic marginal zone lymphoma (SMZL) have a better prognosis,10 while cases of hypermutated mantle cell lymphoma (MCL) display specific clinicopathological features (leukemic course, longer survival).11–14
There were two aims of the current study. First, we wished to identify a set of markers that could be analyzed by immunohistochemical assays, since this could have a potential clinical value for analyzing SHM status. Secondly, we wanted to understand the mechanisms and markers of SHM more thoroughly, since genes associated with high SHM status could help to elucidate the mechanisms of SHM and the co-factors involved with AID. These mechanisms and markers are currently insufficiently characterized, except, to some extent, in the case of AID, whose role was mentioned above, and ZAP70, which is associated with unmutated IgVH genes in CLL.15 However, this finding is confined to CLL: no association between ZAP70 and SHM has been observed in cases of MCL.16
We used expression profiling in 93 samples from patients with B-cell lymphoma to identify genes that could serve as markers of high SHM levels, a characteristic that has prognostic significance. The study focuses on small B-cell lymphomas in which the presence of Ig somatic hypermutation has been reported to be a biological variable of clinical significance, such as CLL, MCL and SMZL, as noted above.9–14 Furthermore, expression profiling was used in patients samples to identify genes that may play a role in the process of SHM.
Finally, the association of several markers with overall survival was assessed.
| Design and Methods |
|---|
|
|
|---|
Analysis of IgVH somatic hypermutation and ongoing mutations in patients samples
DNA was extracted from frozen tissue blocks using proteinase-K and purified by phenol-chloroform extraction. Rearranged IgVH genes were amplified by using a semi-nested polymerase chain reaction (PCR) method, as described previously.10,13,17 The first round of the PCR was performed using a mixture of six framework region 1 (FR1) VH family-specific primers and two consensus primers for the JH gene. The second round of PCR was performed in six separate reactions with one of the six VH FR1 primers and JH internal primers. Direct sequencing was performed in an ABI PRISM 310 or 3700 Genetic Analyzer (Applied Biosystems, Foster City, CA, USA) from both strands, using the same primers as in the amplification. Mutations were identified by comparison with the germ-line sequence (Ig BLAST and V BASE sequence directory). Cases with a high mutational load were defined as those in which the homology between the amplified sequence and the original sequence was less than 98%, while those cases with less than 2% variation with respect to the original sequence were considered to have a low mutational status. In all cases, analysis of IgVH SHM was performed in duplicate.
Ongoing mutation was analyzed in a subset of 31 patients with more than 2% SHM, comprising 14 CLL, 4 MCL and 13 SMZL cases. The VH fragments were amplified as outlined above and the PCR products were purified and cloned into the pCR2.1-TOPO vector (Invitrogen, CA, USA). Several colonies were sequenced from each case. The sequences obtained were compared with that of the wild type and the number of stable and ongoing mutations was determined. A case was considered to have ongoing mutation if at least two of its colonies had a mutation that was not present in the other colonies. For evaluation of intraclonal heterogeneity, a mutation was considered to be confirmed if it was observed more than once in VH gene molecular clones from the same tumor specimen. Only confirmed mutations were considered as evidence of intraclonal heterogeneity. Unconfirmed mutations, a substitution mutation observed in only one of the VH gene molecular clones from the same tumor specimen, were disregarded because they could have been caused by a Taq polymerase error.
RNA isolation, cDNA microarray target preparation and hybridization
Total RNA was extracted from frozen tumor samples using the Trizol reagent (Invitrogen) and RNA was purified and treated with RNase-free DNase I using the RNeasy kit (Qiagen Inc., Valencia, CA, USA). Next, 1–5 µg of target RNA was amplified using T-7 in vitro transcription18 and 2.5 µg of amplified RNA were directly labeled with cyanine 5-conjugated dUTP or cyanine 3-conjugated dUTP (Amersham, Uppsala, Sweden). The reference sample used was 2.5 µg of amplified RNA from the Universal Human Reference RNA (Stratagene, La Jolla, CA, USA).
Microarray studies were carried out using the CNIO OncoChip and labeling and hybridizations were performed as previously described.18–20 The cloned sequences of all the genes included in the OncoChip and the reproducibility of the expression data (measured by quantitative PCR) of multiple genes have been verified.18,20–23 Scanning and image analysis were performed using a Scanarray 5000 XL (GSI Lumonics, Kanata, Ontario, Canada) and GenePix Pro Software (Axon Instruments Inc., Union City, CA, USA), respectively.
Data analysis and normalization of microarray data
Raw microarray data were processed as previously described.19,20,22 Data from CLL and SMZL cases were normalized against the average expression of each gene from five reactive lymph node and five normal spleen samples, respectively. Normalization was carried out only for those genes for which at least 50% of the data were available in the control samples. All other genes were excluded from analysis. Data from MCL samples were normalized against data from tonsillar mantle cells purified using magnetic beads.24 Raw and normalized data of the genes studied are presented in the Online Supplementary Appendix.
Statistical analyses
To identify the genes of importance in distinguishing cases with a high or low IgVH SHM load, Welchs t-statistic, which does not assume equal variances,25 was calculated, and false discovery rates were obtained from tests of 100,000 permutations. Genes were considered significant if they were differentially expressed between patients with high and low SHM burdens to a false discovery rate < 0.20.
The biological functions of genes were assigned using the Gene Ontology (www.geneontology.org) and Genecards26,27 (http://bioinformatics.weizmann.ac.il/cards) databases. Gene pathways were analyzed with the help of Gene Set Enrichment Analysis (GSEA) version 2.0.1. The gene-set database included Biocarta pathways,28 clusters of functionally related, coregulated genes identified by our unsupervised clustering and molecular signatures defined for lymphoma in the Staudt molecular signature database (http://lymphochip.nih.gov/signaturedb).29
Fishers exact test was used to investigate associations between SHM marker expression and SHM status and with ongoing mutation status in all tumor types. The statistical significance of relationships between Ki67 and SHM markers was evaluated using Pearsons
2 test. Differences were considered statistically significant for p values <0.05. Overall survival was analyzed in patients with different levels of expression of SHM markers using the Mantel-Cox test. Univariate receiver operation curves (ROC) were analyzed to identify the specificity of marker genes. All statistical analyses were done using SPSS version 13.0 (SPSS Inc., Chicago, IL, USA).
Immunohistochemistry
Antibodies against CDK7, DP1 (TFDP1), HMGB2, POLA, PRIM1, RAD51C, RCC1 (CHC1) and RFC4 were used to determine protein expression employing a previously described protocol.24,30 Suppliers and dilutions of antibodies for immunohistochemistry were as follows: Ki67 at 1:50 (MIB1 antibody from DAKO, Glostrup, Denmark), CDK7 at 1:225, HMGB2 at 1:10, RCC1 (CHC1) at 1:150, RFC4 at 1:150, RAD51C at 1:30 and POLA at 1:100 (Santa Cruz Biotechnology, Santa Cruz, CA, USA) and DP1 (TFDP1) at 1:60, and DNA primase (PRIM1) at 1:35 (NeoMarkers, Lab Vision, Fremont, CA, USA). To analyze SHM status and marker association with ongoing mutation, samples were evaluated semi-quantitatively by a panel of pathologists. Given that the 60 cases of MCL used for survival analysis were included in a tissue microarray amenable to automatic quantification, the expression of the markers studied was ranked as strong, weak or negative, as determined by the ARIOL semi-automated computerized training system (http://www.aicorp.com/products/02path.htm). The cellular location of each marker was determined using information from published studies and defined to the ARIOL training system. On the basis of this, the system was trained by an expert team of technicians and pathologists to classify the staining on tissue microarray slides as strong, weak or negative for each marker selection. The automated selection was then verified to ensure the robustness of results. Each cellular location was counted as a single point and, depending on its intensity, was assigned to the strong, weak or negative category. The operation is similar to flow cytometry in tissue sections, in which quantitative measurements (e.g., number of cells in a core from a patients sample) and qualitative measurements (e.g., intensity of staining) are made. Cores with low quality staining were considered not evaluable.
| Results |
|---|
|
|
|---|
2% or > 2% SHM, respectively; false discovery rate <0.20). The top 58 genes found to be upregulated in cases with a high mutational status are involved in cell cycling and DNA replication (PRIM1, RFC4, HOXB7, CCNA2, AIM2, BUB3, PCNA, CDK2, JUNB, MCM3, CHEK1, RFC3, POLA, WEE1), DNA repair (CDK7, RAD51C, POLS, HMG2, RAD54B, FRAP1, PRKDC), chromosome condensation (CHC1) and transcription regulation (TFDP1, POLA, MED6, TCEA2, E2F5). As a validation step, we found the same genes in the top ranks using GSEA (Table 1).
|
|
|
|
|
|
|
| Discussion |
|---|
|
|
|---|
Conceptually, SHM must consist of at least two main steps: (i) unwinding of DNA during replication or transcription to allow access to the SHM machinery; (ii) mutation of the DNA by a polymerase during replication or by a DNA repair enzyme. A large proportion of the genes identified as SHM markers or as being involved in the SHM mechanism (Tables 1 and 2, respectively) play a role in DNA repair, replication, transcription and cell cycling. Most of these genes had not previously been associated with SHM and so represent novel findings. The composition of the SHM signature identified here seems to confirm that SHM is indeed associated with a characteristic set of changes in the cell machinery in charge of cell cycling, transcription and DNA repair and replication, when considering different B-cell lymphoma types together. This contrasts with the findings when considering only CLL, in which SHM has been mainly found to be associated with changes in B-cell receptor signaling genes, such as ZAP70 or others.31–33 Although theoretically, malignant transformation could have taken place after silencing of the hypermutation process in a post-germinal center B cell, these findings seem to show that in fact hypermutated tumors, when considered as a whole, preserve a distinctive signature. An additional finding supporting this observation is the increased expression of some of the genes composing this SHM signature (TFDP1 and POLA) in cases with ongoing SHM, in which the SHM machinery is still active. Nevertheless, the limited number of cases analyzed for ongoing SHM impedes generalization regarding the value of the genes associated with ongoing SHM in conditions such as CLL, in which ongoing SHM is more infrequent.
DNA replication appears to play a crucial role in SHM given that a large number of critical genes involved in the process are overexpressed in cases with high SHM status. These genes include PRIM1, POLA, RFC4, CDK7, TFDP1 and RCC1. POLA is a replication polymerase in a complex with DNA primase (PRIM1). It is expressed in hypermutating cells but not in resting B cells,34 as our study of ongoing mutation confirms (Table 3B). PRIM1, which is responsible for synthesizing the small RNA primers during discontinuous DNA replication, is significantly associated with SHM in microarray analysis and tends to show higher protein expression in cases with > 2% SHM (Table 2). CHC1 (RCC1) is involved in regulating the onset of chromosome condensation in the S phase. Replication factor C (RFC4) is a DNA polymerase accessory protein whose function is to elongate primed DNA templates through the action of DNA polymerase. It acts as a clamp loader to enable the polymerase to bind to DNA.
CDK7 and RAD51C showed significantly greater expression in cases with higher SHM, and both genes were significantly associated with SHM in validation studies using protein expression. Furthermore nodal marginal zone lymphomas with a high SHM burden also showed increased CDK7 mRNA expression (data not shown). Both genes can activate p53, which leads to cell-cycle arrest during which DNA repair can progress. CDK7, a cyclin-dependent kinase, is a mediator of cell-cycle progression through activation by binding to cyclins. Additionally, CDK7 is a component of the transcription factor TFIIH, which helps control transcription by RNA polymerase II and possesses DNA repair and helicase activities. RAD51C is involved in the homologous recombination repair pathway of double-stranded DNA breaks that arise during DNA replication or are induced by DNA-damaging agents, and in meiotic recombination. Its overexpression may contribute to SHM by repairing the DNA after mutations have been induced, perhaps during DNA replication, and by allowing the cells to progress through the cell cycle. Finally, in terms of DNA damage repair, microarray analysis also revealed overexpression of PRKDC in IgVH mutated cases. PRKDC is a serine/threonine-protein kinase involved in DNA non-homologous end joining, which is required for double-strand-break repair and V(D)J recombination. PRKDC is also involved in the modulation of transcription. However, the role of PRKDC is unclear, as SHM is known to be able to proceed essentially unaffected by deficient DNA-PK activity.35
Interestingly, these findings show that some specific components of protein complexes involved in DNA repair and transcription control, such as CDK7 and RAD51C, appear to be selectively increased in association with SHM. The observations require additional confirmation, but do raise an interesting hypothesis to be explored in functional assays.
B-cell translocation gene 1 (BTG1) is a negative regulator of cell proliferation whose expression was increased in mutated cases. The t(8;12)(q24;q22) BTG1/MYC translocation has been identified in CLL. TFDP1 has been associated in this study with ongoing mutation. It is a component of the E2F/DP transcription factor complex and a regulator of a number of genes whose products are involved in cell-cycle regulation and DNA replication. Its expression can be induced by p53 and it is also responsible for indirect activation of p53 through its role in the cell cycle. HMG2 may also be involved in the final ligation step in DNA end-joining processes of DNA double-strand-break repair, and V(D)J recombination. p53-dependent DNA-repair may also play a role in the upregulation of a variety of transcription factors that are associated with high SHM status in patients samples.
This SHM signature does not include AID, an enzyme required for SHM but largely absent in multiple lymphoproliferative processes with increased SHM, as shown by Pasqualucci and co-workers.36
Finally, a potential clinical application of these findings was explored, since some of the markers identified here as being associated with SHM are significantly associated with improved overall survival in patients, with MCL, underlining the fact that SHM is a prognostic factor in the major types of low-grade B-cell lymphoma, including MCL.14 Taking into consideration all these findings, studies can be developed to examine the role of the specific genes detected here in the generation of IgVH SHM in small B-cell lymphomas, and to analyze their potential prognostic utility in larger series of standardized patients.
| Acknowledgments |
|---|
| Footnotes |
|---|
LT: experimental procedures, contribution to and/or analysis of microarray data, contribution of essential new reagents, manuscript preparation; MA: experimental procedures, contribution to and/or analysis of microarray data, contribution of essential new reagents, manuscript preparation; MG-C: immunohistochemistry evaluation; RV: data analysis and interpretation, contribution to and/or analysis of microarray data; PA: experimental procedures, data analysis and interpretation; AS-A: experimental procedures; JFG: immunohistochemistry evaluation; AR, FIC, NM, ER-B, MM: contribution to and/or analysis of microarray data; MAP: immunohistochemistry evaluation, contribution to and/or analysis of microarray data, contribution of essential new reagents, manuscript preparation.
The authors reported no potential conflicts of interest.
Received for publication March 3, 2008. Revision received April 17, 2008. Accepted for publication April 23, 2008.
| References |
|---|
|
|
|---|
: changes induced in IFN-
sensitive and resistant cutaneous T cell lymphoma. J Interf Cyt Res 2004;24:185-95.[CrossRef][ISI][Medline]
B targets BCL2 and BIRC5/survivin characterizes small B-cell and aggressive B-cell lymphomas, respectively. J Pathol 2005;206:123-34.[CrossRef][ISI][Medline]
| ||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | TABLE OF CONTENTS | ARCHIVE | SUBSCRIPTIONS |