|
||
Address correspondence to Wendy F. Davidson, Greenebaum Cancer Center and Dept. of Microbiology and Immunology, University of Maryland School of Medicine, 15601 Crabbs Branch Way, Rockville, MD 20855. Phone: (301) 517-0324; Fax: (301) 517-0344; email: wdavi002{at}umaryland.edu
| Abstract |
|---|
|
|
|---|
Key Words: autoimmunity B cell lymphoma
Abbreviations used in this paper: ANA, antinuclear antibody; BCR, B cell receptor; BL, Burkitt lymphoma; BLL, Burkitt-like lymphoma; CBL, centroblastic; CDR, complementarity-determining region; CL, cardiolipin; dsDNA, double-stranded DNA; FBL, follicular B cell lymphoma; FWR, framework region; GC, germinal center; GL, germline; IBL, immunoblastic; MALT, mucosal-associated lymphoid tissue; MZL, splenic marginal zone lymphoma; PC, plasma cell; PCA, principal component analysis; PCT, plasmacytoma; PL, plasmacytoid lymphoma; PPC, phosphorylcholine; R mutation, replacement mutation; RF, rheumatoid factor; S mutation, silent mutation; SBL, small B cell lymphoma; SJL, SJL-ß2M/ lymphoma; ssDNA, single-stranded DNA.
| Introduction |
|---|
|
|
|---|
gene groups (3). In spite of this biased repertoire, the autoimmune responses to nuclear antigen by lpr/gld and other autoimmune mice resemble normal, secondary antibody responses to foreign antigen, and affinity maturation of the autoimmune ANA responses appears to be the result of antigen drive (4, 5). Thus, within the B cell compartment of lpr and gld mice, there is evidence for loss of tolerance to and chronic stimulation by DNAhistone complexes.
GL mutations in Fas/Fasl also greatly increase the risk of B cell neoplasia in humans and mice. B cell lymphomas in autoimmune lymphoproliferative syndrome patients are predominantly follicular and phenotypically diverse (6). In mice, lymphomas develop between 6 and 12 mo of age, and tumor incidence is strain related. By 1 yr of age,
30% of C3H-gld and 60% of BALB-gld mice have monoclonal outgrowths of B cells in the spleen and LN that metastasize to nonlymphoid organs (2). These tumors can be passaged to purity in immunodeficient scid mice and are lethal for their hosts (2). The scid-passaged lpr and gld B cell lymphomas are predominantly plasmacytoid in morphology, are
dull+ CD23 CD21 Mac-1+, have undergone Ig isotype switching, and spontaneously secrete Ig (2), suggesting that they derive from antigen-experienced, possibly germinal center (GC)-selected B cells.
Lymphomagenesis is a multistep process involving both genetic and environmental factors. There is mounting evidence that chronic antigen stimulation might be an important environmental factor. In support of this, there are reports that within the Ig VH and VL genes of many B cell tumors and clonal progeny of nontransformed cells there is a counterselection for replacement mutations (R mutations) in the framework regions (FWRs) of the B cell receptor (BCR) and generally higher ratios of R versus silent mutations (S mutations) in the complementarity-determining regions (CDRs; references 711). The persistence of these mutational patterns in the face of ongoing somatic mutation suggests that there is strong selective pressure to preserve the structural integrity of the BCR. Moreover, analyses of intraclonal variation in human follicular lymphomas revealed genealogic relations between subclones, indicative of antigen-driven clonal evolution (12, 13). Other compelling evidence that chronic antigen stimulation might be a factor in B cell lymphomagenesis comes from studies of mucosal-associated lymphoid tissue (MALT) B cell lymphomas. These tumors are thought to derive from marginal zone B cells and develop at sites of chronic inflammation caused by autoimmune disease or infection (1419). As examples, patients with Sjögren's syndrome or Hashimoto's thyroiditis have an increased risk of developing MALT lymphomas of the salivary and thyroid glands, respectively (1517). Notably, in both patient populations, MALT lymphomas have been reported that produce rheumatoid factor (RF) and exhibit nonrandom use of IgVH and VL genes, suggesting that autoreactive B cells chronically stimulated by nominal self-antigen might be selectively targeted for transformation (14, 17, 18). Of interest, C57BL/6-lpr mice deficient in all T lineage cells also develop B lymphomas that produce RF (20).
The goals of this study were twofold. First, we sought to determine whether B cell transformation in gld mice is a random event or is skewed toward chronically stimulated, autoreactive B cells. Second, we undertook gene expression profiling using oligonucleotide-based microarrays to determine whether the gld plasmacytoid tumors exhibit a transcriptional profile that might be useful for diagnosis and relevant to pathogenesis.
| Materials and Methods |
|---|
|
|
|---|
Tumors.
Early passage number (P1 or P2) plasmacytoid lymphomas (PLs), including 13 BALB-gld and 6 C3H-gld, were transplanted i.p. into 46-wk-old scid mice. Sera, spleens, and mesenteric LNs were harvested from tumor-bearing mice. Single cell suspensions were made for FACS analysis and isolation of DNA and RNA. Southern hybridization analysis for IgH and IgL gene rearrangements was used to establish clonality (2). The proportion of tumor cells in spleen and LN samples was determined by FACS analysis (BD Biosciences). The majority of the tumor cell preparations contained >90% CD19+ B cells with high forward and side scatter. Primary tumor cells were obtained from four mice with advanced tumors and were enriched by MACS sorting (Miltenyi Biotec) of CD19+ cells. These cells were predominantly CD23 CD21 Mac-1+ with high forward and side scatter. Non-gld lymphomas used in array analyses were mostly from NFS.V+ congenic mice (21), but included plasmacytoma (PCT) cell lines provided by M. Potter (National Cancer Institute, Bethesda, MD) and primary lymphomas of SJL/J-ß2M/ mice.
Cultures.
Sorted CD19+ CD23 tumor cells were cultured at 107/ml in supplemented RPMI medium for 48 h in the presence or absence of 25 µg/ml LPS (Sigma-Aldrich).
ELISA for Ig Isotypes, ANA, Single-Stranded DNA (ssDNA), and RF.
Ig isotypes and ANA, anti-ssDNA, and RF titers in sera from tumor-bearing mice and tumor tissue culture supernatants were determined using a mouse Ig isotyping ELISA kit (BD Biosciences) and ANA, ssDNA, and RF ELISA kits (Alpha Diagnostic International). Antibodies to dsDNA were detected by ELISA as described previously (2). Nuclear staining was detected using a kit (Antibodies Inc.).
Cloning and Characterization of Ig Genes.
Total RNA was isolated from tumor cells using TRIzol reagent (Invitrogen). First-strand cDNA was generated using oligo (dT) or random primers and SuperScript (Invitrogen). PCR amplification of tumor VH and VL genes was performed using a set of degenerate primers and conditions listed in Table S1 (see Online Supplemental Material section below). Selected PCR products were cloned using a TOPO TA cloning kit (Invitrogen). Multiple independent plasmid isolates for VH and VL were sequenced for each tumor. Ig gene sequences were searched for their closest GL sequence match against a DNA and protein database (http://www.ncbi.nlm.nih.gov/BLAST), a database for Ig genes (http://www.ncbi.nlm.nih.gov/igblast), and the Celera database. The probability that an excess of R mutations in CDR1 and CDR2 regions or a scarcity of R mutations in FWRs occurred by chance was calculated using a multinomial distribution model and a published program (http://www.stat.stanford.edu/immunoglobin; reference 22).
Oligonucleotide Microarrays.
Microarray chips, printed by the Microarray Research Facility at NIAID (http://madb.niaid.nih.gov), comprised
6,700 mouse gene targets represented by 70-mer oligonucleotides purchased from Compugen. Total tumor cell RNA was extracted by the TRIzol protocol and further purified using the RNeasy column-based protocol (QIAGEN). A reference sample was created by extracting total RNA from 10 well-characterized hematopoietic cell lines and making a pool comprising equal amounts of RNA from each line (Table S2; see Online Supplemental Material section below). The labeling of cDNA, chip hybridization, and chip scanning were performed as described previously (23). The microarray data accession number is GSE1908.
Microarray Data Analysis.
The microarray dataset was organized and analyzed using SAS software (SAS Institute). From a total of 6,800 genes, 5,770 genes had <30% of missing values across all 134 samples. For these genes, the missing values were replaced with the overall mean of each individual gene, and this no-missing-value dataset was used in the subsequent principal component analysis (PCA). Student's t test and the nonparametric Kruskall-Wallis test were used to compare the expression levels of genes among different lymphoma groups. Hierarchical clustering (average linkage method) of selected significant genes against lymphoma samples was performed using a program available in the NIAID microarray database application (http://madb.niaid.nih.gov).
Quantitative PCR.
Primers and conditions for quantitative PCR are given in Table S3 (see Online Supplemental Material section below).
Online Supplemental Material.
The supplemental materials include the following: (a) primers and conditions for PCR amplification of VH and VL genes (Table S1); (b) cell lines used to generate the standard for microarrays (Table S2); (c) primers for quantitative PCR (Table S3); (d) translated amino acid sequences of VH and VL for gld tumors (Fig. S1); (e) VH-DH-JH rearrangement of tumor VH CDR3 regions (Fig. S2); (f) hierarchical clustering of 415 genes that discriminate gld PL from other mouse lymphomas (Fig. S3); (g) the same data as shown in Fig. S3 with all genes identified (Fig. S4); (h) hierarchical clustering of gene expression data for 99 genes that distinguish gld PLs from other PC neoplasms (Fig. S5); and (i) the same data as shown in Fig. S5 with all genes identified (Fig. S6). Tables S1S3 and Figs. S1S6 are available at http://www.jem.org/cgi/content/full/jem.20041575/DC1.
| Results |
|---|
|
|
|---|
dull CD23 CD21 CD5dull or CD5, a phenotype remarkably similar to that of the nontransformed B cells that selectively accumulate in old lpr and gld mice (2, 24, 25). The Ig isotype produced by the tumors was determined by testing sera from tumor-bearing scid mice and supernatants from cultures of sorted tumor cells. All 16 sera tested contained significantly elevated levels of Ig, and isotype switching had occurred in 11 (Table I). One tumor, 424, produced polyreactive antibodies of different isotypes that were reactive in ELISA assays for dsDNA, phosphorylcholine (PPC), and cardiolipin (CL). Tumor 600 consistently secreted both IgG2a and IgG3 in several scid passages (Table I). Ig was not detected in sera from control scid mice (not depicted). The sera also were assayed for ANA, RF, and antibodies to dsDNA, ssDNA, PPC, and CL. 9 of 16 (56%) sera tested positive for ANA by ELISA and for nuclear staining detected by immunohistochemistry (Table I). Within this group, Ig from tumor 426 also reacted strongly with ssDNA and was weakly positive for RF. Only antibodies from 221, 424, and 2685 reacted with dsDNA. antibodies from two ANA tumors, 217 and 421, also reacted with ssDNA. Although negative for ANA activity by ELISA, antibodies from 425 showed weak reactivity with nuclei, weak reactivity with ssDNA, and had high titers of RF. Tumor 8850 produced high titers of RF and anti-ssDNA antibodies, and tumors 205 and 208 produced RF only. Cross-reactivity between RF and ssDNA has been reported previously (26). With the exception of 424, no tumors produced antibodies reactive with PPC or CL. In all, 15 of 17 Ig-secreting tumors produced RF and/or antibodies reactive with nuclear antigens. 4 of 17 tumors produced antibodies that appeared to be polyreactive. Results obtained with tumor culture supernatants conformed to those from sera of scid mice bearing the same tumor (not depicted).
|
2, P < 0.05) from that reported for normal BALB/c B cells (51%, 30 of 59 sequences; reference 29). A similar trend was observed in the C3H-gld tumors, where 66.7% (four of six) used J558 genes.
|
Evidence for Nonrandom Use of DH and JH Genes in PLs.
Analysis of 18 tumor CDR3 DH sequences revealed evidence of DD fusions and significant skewing in DH family use (Table II and Fig. S2). Fusion of two D segments, rare in normal B cells, was observed in six tumors. Altogether, 15 of 18 tumors (83.3%) used a D segment from the DSP2 family, whereas 4 tumors used D genes from the DFL family. In contrast, in a previous study, 50% to 55% of BALB/c follicular B cells used DSP family members (29). Notably, many tumor VH CDR3 sequences contained one or more R and N residues, a frequent feature of ANA (3). In many cases, the R and N residues were generated by N region nucleotides (Fig. S2). One motif, [(W)LR(R)n], found recurrently in anti-DNA antibodies (3, 5, 30), was present in the VH CDR3 region of tumors 426, 593, 600, and 2685, a finding consistent with the production of ANA and anti-DNA antibodies by these tumors (Table I).
Biased usage of JH segments was also observed in the tumors, particularly in the BALB-gld group (Table II and Fig. S2). Overall, the JH distribution was 5.6% JH1, 33.3% JH2, 11.1% JH3, and 50% JH4. By comparison, the JH family usage in normal BALB/c B2 cells was reported to be 9% JH1, 25% JH2, 32% JH3, and 34% JH4 (29). The tumor and B2 cell JH distributions differ significantly (
2, P < 0.0001). The JH4 bias was even more prominent (57%) if only the 14 tumors with demonstrated autoreactivity were considered (Tables I and II).
Biased Use of VL Gene Families and J
Genes in PLs.
In normal B cells, the pattern of VL family usage generally corresponds to the estimated family size (31). Two exceptions are the V
1 and V
4 families, which each account for
13% of the expressed VL genes (31, 32). Analyses of PL VL sequences revealed a pattern different from that in normal B cells (Table II). Although 10 V
families were represented among the PLs, the following three families were overrepresented: V
4/5 (15.8%), V
21 (15.8%), and V
12/13 (10.5%). Although the sample size is small, preferential subfamily use within these families also was observed. Tumors 355 and 426 used V
12-44, and tumors 208 and 2685 used 21-7 (Table II). Notably, both of these V
subfamilies are represented frequently in anti-DNA antibodies (3, 5).
Similar to the VH genes, the VL genes also showed a trend toward increased rearrangement to distal J genes (Table II). JK1, JK2, and JK5 each were used by 5 of 18 tumors (27.8%), and JK4 was used by 3 of 18 tumors (16.7%). In normal mouse B cells, the frequency of rearrangement to JK loci is JK1 > JK2 > JK5 >> JK4, with
35% rearranging to JK4 (13%) or JK5 (22%; reference 31). The frequency of JK4 or JK5 usage by the tumors (45%) differed significantly from that of normal B cells (
2, P = 0.036).
Evidence for Antigen Selection in PLs.
Patterns of DNA mutation in Ig genes can provide valuable insights into the maturational history of the tumor precursors and the importance of BCR integrity, specificity, and affinity in tumorigenesis (711). To investigate the frequency and nature of the mutations in tumor Ig genes, the distribution of R versus S mutations in the FWR, CDR1, and CDR2 of V genes was analyzed (22). Only data for VL genes is shown because most of the mouse VL GL genes have been sequenced and the GL assignments for the tumors can be made with greater confidence. In this analysis, a significantly higher frequency of R mutations in CDRs than expected by chance (PCDR < 0.05) is considered evidence for positive, antigen-mediated selection of cells with the mutated Ig. Conversely, a frequency of R mutations in FWRs that is statistically lower than expected by chance (PFWR < 0.05) implies that there is pressure to maintain the structural integrity of Ig, possibly because antigen is providing important growth/survival signals to the tumor cells (7).
As shown in Table III, two tumors had GL VL sequences and the rest had 110 R mutations. CDR hypermutation was observed in 4 of 18 tumors and FWR conservation was observed in 4 of 18 tumors. Tumors 142 and 425 exhibited both VL CDR hypermutation and VL FWR preservation. Notably, six VL sequences had R mutations in CDR1 and/or CDR2 that gave rise to arginine (R) or asparagine (N) residues important for DNA binding (3). Overall, 6 of 18 tumors had mutated CDR and/or conserved FWR in VL sequences, suggesting that at least one third of the B cells undergoing transformation were antigen selected. All six tumors in this group were Ig isotype switched and those producing antibodies (five out of six) were reactive with nuclear antigens or RF. To test for evidence of ongoing antigen-mediated selection, sequences from up to 15 clones from each tumor were examined for intraclonal variation. This was observed only in tumor 593 and at low frequency. 2 of 15 clones had a single amino acid change at the same position in VH CDR1 (not depicted).
|
To determine how PLs relate molecularly to these lymphoma types and to PCT, we compared their gene expression profiles using a 6,680 gene chip. We also studied a series of primary PC-related neoplasms from SJL-ß2M/ mice that have not been studied previously at the molecular level and whose place in the scheme of terminal B cell differentiation is not known. These tumors resemble those that arise in aging SJL mice, but they develop at an accelerated rate (unpublished data).
PCA of the microarray data revealed that the PL tumors were remarkably similar to each other in terms of gene expression. In the majority of cases, they also were readily distinguished from lymphoma classes that originate from cells that have not yet entered the terminal progression program of PCs (Fig. 1). Further analyses of these data identified 415 genes whose expression levels by t test distinguished PLs from BLLs, BLs, SBLs, CBLs, FBLs, MZLs, and IBLs (Figs. S3 and S4). Most of these genes were expressed at significantly lower levels in the PL than the other lymphoma types, consistent with the well-established down-regulation of many genes that occurs as B cells mature toward secretory stages of development. Fewer than 20 genes were selectively expressed at high levels in PL. One gene of particular interest in this group was Apobec1, which encodes an RNA editing protein that also may function as a DNA mutator (34, 35).
|
|
To further examine the developmental relationships among the three sets of PC neoplasms, we interrogated a series of 30 genes modulated in expression late in B cell differentiation (36, 37). The matrix generated by hierarchical clustering of both the genes and the arrays is shown in Fig. 3 A. All PCTs but one fell on one branch of the tree generated by the array samples, whereas PL and SJL neoplasms again occupied the second branch. 13 of the 21 PL (61%) and 5 of the 10 SJL tumors were tightly grouped, with the remainder being interspersed outside these immediate alignments. The genes also fell into several groups. In one major group, 10 genes generally were expressed at the highest levels by PL, at somewhat lower levels by many of the SJL tumors, and at very low levels by the PCT. Among these were the genes encoding CD45, CD22, CD19, CCR7, ICSBP1, and FCGR2B (Fig. 3 A). A smaller group of genes, including Irf4, Sdc1, and Xrcc5, were expressed at highest levels in PCT and at similarly lower levels in PL and SJL tumors (Fig. 3 A).
|
Primary and scid-passaged PLs Exhibit Similar Patterns of Gene Expression.
To ensure that the remarkably similar gene expression profiles observed in the scid-passaged PLs were not a reflection of tumor selection in the scid environment, we compared gene expression in scid-passaged BALB-gld PLs, isolated primary BALB-gld PL cells, and BALB/c splenic B cells by quantitative PCR. The genes selected for this comparison were chosen from those that distinguish the scid-passaged PLs from other B cell lineage tumors and were either significantly increased (Fig. 4 A) or decreased (Fig. 4 B) in expression relative to purified BALB/c splenic B cells. Notably, these genes included candidates of potential importance in tumorigenesis (Apobec-1 [34, 35], Tssc4 [40], Btg2 [41], Bin1 [42], Prkcb [43], and Mad1l1 [44]), lymphoma survival (BAFF [4547] and IL-10 [48, 49]), and autoimmunity (IL-10 [50, 51] and Ifi202/203 [52]). Il6st, which encodes the gp130 component of the IL-6 receptor (53), also was significantly decreased in expression in primary and passaged tumors (Fig. 4 B). For some genes, such as Apobec1, there were quantitative differences between levels expressed by primary and passaged tumors, but in general, the distinguishing patterns of gene expression for PLs were in evidence before transplantation (Fig. 4, A and B).
|
| Discussion |
|---|
|
|
|---|
The biased pattern of VH and VL gene use observed in the PL is reminiscent of the distinctive repertoires of autoreactive B cells specific for nuclear antigen. Like the PL, anti-dsDNA antibodies and ANA from MRL-lpr, MRL-+/+, (NZBxNZW)F1, and (SWRxNZW)F1 mice also disproportionately used the J558 and 7183 families and the J558.18, J558.33, J558.40, J558.47, and J558.52 subfamilies (3, 5). These VH genes also tended to pair preferentially with a restricted set of VL genes including members of the V
8, V
2, V
4/5, V
9/10, V
21, and V
23 groups (3, 5). These similarities and the autospecificities of the PL strongly imply that they derive from the peripheral autoreactive B cell pool characteristic of lpr/gld mice.
DD fusions and preferential use of downstream JH and JL genes are common features of anti-DNA and antinucleosome antibodies in humans and mice with systemic lupus erythematosus (3, 5459). The preferential use of upstream D genes and JH4 is suggestive of secondary D-JH rearrangements (54, 55, 59), whereas preferential use of downstream JL genes is indicative of secondary rearrangement of IgL loci or receptor editing (5658). The high proportion of PLs using JH4 (55%) and the increased incidence of rearrangement of VL to J
4 or J
5 strongly imply that the majority of PLs derive from autoreactive B cells that have undergone repeated D-JH and/or V-JL rearrangements. Secondary VH rearrangements using D-JH4 also may occur in PLs by a mechanism proposed previously by Klonowski and Monestier (55) and Monestier and Zouali (59). Although chimeric VH chains indicative of secondary rearrangement were not observed in the PL, it is possible that they preferentially used the embedded heptamer at the 3' end of FWR3 that results in virtually seamless rearrangements (6062).
Why predominantly Ig isotypeswitched, autoreactive B cells with the hallmarks of receptor revision appear to be selectively targeted for transformation in gld mice is intriguing. One explanation is that transformation occurs randomly in a B cell population that is remarkably skewed toward receptor-edited autoreactive B cells. Although hybridomas derived from lpr splenic B cells are enriched for autoreactive specificities (35), the actual representation of autoreactive B cells within the unmanipulated B cell pool in tumor-free gld mice is unknown. Another possible reason for selected transformation of autoreactive B cells is that repeated activation of Aid during somatic hypermutation or class switch recombination increases the risk of transforming genetic mutations (63, 64). BCR specificity also might be a factor in the transformation process. Autoreactive B cells in lpr or gld mice in all likelihood are chronically activated and have an extended lifespan. These two factors are conducive to the accumulation of transforming genetic defects and the clonal expansion of mutated cells. Indeed, there is compelling evidence that chronic stimulation of B cells by infectious agents such as hepatitis C virus or Helicobacter pylori can culminate in lymphoma development (6567). Notably, these lymphomas often produce RF and other auto antibodies (11, 14, 1719, 68). This might be explained by molecular mimicry or, alternatively, RF-producing tumors may develop from RF B cells driven to clonally expand by immune complexes (14, 17, 18, 67). With the caveat that the assignments of GL genes may not be absolutely correct, the Ig VH and VL regions of most PLs exhibited multiple nucleotide substitutions in FWRs and CDRs. This observation, together with Ig isotype switching in most tumors and a gene expression profile consistent with early PC differentiation, strongly implies that PLs derive from antigen-experienced B cells that have encountered a GC-like environment. Somatic mutation may occur in follicular GCs or extrafollicularly at the T zone red pulp border as reported for RF B cells (69). ANA-producing B cells and immune complexbinding RF B cells might be particularly prone to chronic stimulation because they can present multiple nucleosomal epitopes to T helper cells and thereby generate polyclonal T cell help (70).
The relative importance of antigen drive in lymphomagenesis and lymphoma survival can be implied from analysis of R and S mutations in the predicted amino acid sequences of IgH and IgL FWRs and CDRs (22). Using this type of analysis on IgL genes, 33% of PLs showed evidence of selection for FWR sequence conservation or diversification of CDR sequence, suggesting selection by antigen. The proportion of tumors with the hallmarks of antigen selection increased to 61% if the R/S mutation analysis of Ig VH genes also was included (not depicted). The lack of evidence for intraclonal variation among individual PLs may indicate a lack of pressure to enhance binding affinity, antigen independence after transformation, or a tumor microenvironment not conducive to further mutation. Although chronic stimulation by auto antigen might be an etiologic factor in lymphomagenesis in gld mice, the importance of continued antigen stimulation for proliferation and survival of PLs is unclear. The reduced Ig expression on PLs may preclude antigen-driven proliferation but still facilitate the delivery of survival signals.
Earlier histologic and phenotypic analyses of PL tumors propagated in immunodeficient mice indicated that the majority of the tumors resembled activated B cells in phenotype and had a plasmacytoid morphology, although a proportion of cells differentiated into typical PCs (2). This phenotype suggested that most of the tumor cells were arrested at an early stage in PC differentiation. Gene expression profiling and quantitative PCR established that the PLs indeed are less differentiated than PCs or PCTs and readily distinguishable from less mature B cell lineage tumors. The similarities in gene expression patterns between the scid-passaged PLs and primary PLs provides assurance that the transplanted PLs are not undergoing appreciable changes in gene expression as a result of selection or a change from the gld to the scid microenvironment. Given the rarity of murine immature PC tumors, the PLs will be valuable tools for studying gene expression during PC differentiation and possibly for developing diagnostic transcriptional profiles.
Microarray analyses also led to the identification of genes differentially expressed in PLs, including some that might contribute to the transformation process. Genes that were significantly up-regulated relative to BALB/c follicular B cells included Apobec1, Prkcb, BAFF (Tnfst 13b), Ifi202, and IL-10. Apobec1 is of particular interest because it was selectively overexpressed in PLs and levels were particularly high in primary tumors. APOBEC1, a cytidine deaminase that constitutes the catalytic component of an RNA editing complex, has homology to AID, an activation-induced cytidine deaminase required for somatic hypermutation, and class switch recombination in B cells (34, 35). Constitutive expression of Aid in transgenic mice leads to hypermutation of TCR and Myc genes and development of T cell lymphomas (71). Although the function of APOBEC1 in B cells is unknown, it may have oncogenic potential, as overexpression in liver leads to hepatocellular carcinomas (72). High expression of Prkcb, although not unique to PLs, might be significant in view of the role of protein kinase C ß in survival signaling via NF-
B and the finding that inhibitors of protein kinase C ß induce cell death in lymphomas with dysregulated expression of the gene (73). Both BAFF and IL-10 can be important autocrine growth and survival factors for autoreactive B cells and B cell lineage lymphomas (4551). IL-10 also may protect lymphomas from immune attack through its immunomodulatory effects on T cells and macrophages (50). Overexpression of Ifi202 in PLs relative to normal BALB/c B cells warrants further exploration considering the association of this gene with susceptibility to systemic lupus erythematosus and the survival of autoreactive B cells in NZB mice (52). Ifi202 expression might be regulated by IFN or IL-6 (74); however, the consistently low level of Il6st transcripts in PLs suggests that Ifi202 regulation in PLs might be IL-6 independent. In agreement with the down-regulation of IL6st, IL-6 treatment of PL cell lines failed to induce phosphorylation of relevant JAKs or STATs (unpublished data). The block in IL-6 signaling might be a factor in the arrest of PLs at an early stage in PC development.
Genes of interest that were selectively down-regulated in PLs and may have a role in tumorigenesis included Bin1, Mad1l1, and Btg2. Although defective expression of these genes has been described in a variety of nonlymphoid cancers, this is the first report of their down-regulation in B lymphomas (41, 42, 44). Bin1, which was greatly reduced in expression in both primary and transplanted PLs, encodes a Myc-interacting protein with tumor suppressor properties (42, 75). BIN1 frequently is greatly reduced or absent in malignant cells, and this deficit correlates with resistance to cell death (42, 75). Defective stimulation of Myc-induced apoptosis might be an important factor in progression to malignancy in PLs.
In summary, we have demonstrated that the PLs that develop spontaneously in gld mice derive predominantly from antigen-driven autoreactive B cells, are arrested at an early stage of PC differentiation, and exhibit a transcriptional profile that distinguishes them from other murine B cell neoplasms and also might be relevant to pathogenesis. This model will be invaluable for further investigation of the relationship between autoimmunity and B lymphomagenesis and for elucidating changes in gene expression relevant to early events in PC differentiation.
| Acknowledgments |
|---|
This work was supported by NIH grant CA82872-02 to W.F. Davidson.
The authors have no conflicting financial interests.
Submitted: 5 August 2004
Accepted: 25 October 2004
| References |
|---|
|
|
|---|
1 Cohen, P.L., and R.A. Eisenberg. 1991. Lpr and gld: single gene models of systemic autoimmunity and lymphoproliferative disease. Annu. Rev. Immunol. 9:243269.[CrossRef][Medline]
2 Davidson, W.F., T. Giese, and T.N. Fredrickson. 1998. Spontaneous development of plasmacytoid tumors in mice with defective FasFas ligand interactions. J. Exp. Med. 187:18251838.
3 Radic, M.Z., and M.G. Weigert. 1994. Genetic and structural evidence for antigen selection of anti-DNA antibodies. Annu. Rev. Immunol. 12:487520.[CrossRef][Medline]
4 Shlomchik, M.J., A.H. Aucoin, D.S. Pisetsky, and M.G. Weigert. 1987. Structure and function of anti-DNA autoantibodies derived from a single autoimmune mouse. Proc. Natl. Acad. Sci. USA. 84:91509154.
5 Tillman, D.M., N.T. Jou, R.J. Hill, and T.N. Marion. 1992. Both IgM and IgG anti-DNA antibodies are the products of clonally selective B cell stimulation in (NZB x NZW)F1 mice. J. Exp. Med. 176:761779.
6 Straus, S.E., E.S. Jaffe, J.M. Puck, J.K. Dale, K.B. Elkon, A. Rosen-Wolff, A.M. Peters, M.C. Sneller, C.W. Hallahan, J. Wang, et al. 2001. The development of lymphomas in families with autoimmune lymphoproliferative syndrome with germline Fas mutations and defective lymphocyte apoptosis. Blood. 98:194200.
7 Lossos, I.S., C.Y. Okada, R. Tibshirani, R. Warnke, J.M. Vose, T.C. Greiner, and R. Levy. 2000. Molecular analysis of immunoglobulin genes in diffuse large B-cell lymphomas. Blood. 95:17971803.
8 Lossos, I.S., A.A. Alizadeh, M.B. Eisen, W.C. Chan, P.O. Brown, D. Botstein, L.M. Staudt, and R. Levy. 2000. Ongoing immunoglobulin somatic mutation in germinal center B cell-like but not in activated B cell-like diffuse large cell lymphomas. Proc. Natl. Acad. Sci. USA. 97:1020910213.
9 Kuppers, R., K. Rajewsky, and M.L. Hansmann. 1997. Diffuse large cell lymphomas are derived from mature B cells carrying V region genes with a high load of somatic mutation and evidence of selection for antibody expression. Eur. J. Immunol. 27:13981405.[Medline]
10 Vescio, R.A., J. Cao, C.H. Hong, J.C. Lee, C.H. Wu, M. Der Danielian, V. Wu, R. Newman, M.Z. Radic, and M.G. Weigert. 1994. Genetic and structural evidence for antigen selection of anti-DNA antibodies. Annu. Rev. Immunol. 12:487520.[CrossRef][Medline]
11 De Re, V., S. De Vita, A. Marzotto, M. Rupolo, A. Gloghini, B. Pivetta, D. Gasparotto, A. Carbone, and M. Boiocchi. 2000. Sequence analysis of the immunoglobulin antigen receptor of hepatitis C virus-associated non-Hodgkin lymphomas suggests that the malignant cells are derived from the rheumatoid factor-producing cells that occur mainly in type II cryoglobulinemia. Blood. 96:35783584.
12 Bahler, D.W., and R. Levy. 1992. Clonal evolution of a follicular lymphoma: evidence for antigen selection. Proc. Natl. Acad. Sci. USA. 89:67706774.
13 Matolcsy, A., E.J. Schattner, D.M. Knowles, and P. Casali. 1999. Clonal evolution of B cells in transformation from low- to high-grade lymphoma. Eur. J. Immunol. 29:12531264.[CrossRef][Medline]
14 Mariette, X. 2001. Lymphomas complicating Sjögren's syndrome and hepatitis C virus infection may share a common pathogenesis: chronic stimulation of rheumatoid factor B cells. Ann. Rheum. Dis. 60:10071010.
15 Hyjek, E., W.J. Smith, and P.G. Isaacson. 1988. Primary B-cell lymphoma of salivary glands and its relationship to myoepithelial sialadenitis. Hum. Pathol. 19:766776.[Medline]
16 Hyjek, E., and P.G. Isaacson. 1988. Primary B cell lymphoma of the thyroid and its relationship to Hashimoto's thyroiditis. Hum. Pathol. 19:13151326.[Medline]
17 De Re, V., S. De Vita, D. Gasparotto, A. Marzotto, A. Carbone, G. Ferraccioli, and M. Boiocchi. 2002. Salivary gland B cell lymphoproliferative disorders in Sjögren's syndrome present a restricted use of antigen receptor gene segments similar to those used by hepatitis C virus-associated non-Hodgkins's lymphomas. Eur. J. Immunol. 32:903910.[CrossRef][Medline]
18 Martin, T., C. Weber, H. Levallois, N. Labouret, A. Soley, S. Koenig, A.S. Korganow, and J.L. Pasquali. 2000. Salivary gland lymphomas in patients with Sjögren's syndrome may frequently develop from rheumatoid factor B cells. Arthritis Rheum. 43:908916.[CrossRef][Medline]
19 Eng, M.A., U. Kallemuchikkal, and P.D. Gorevic. 2000. Hepatitis C virus, autoimmunity and lymphoproliferation. Mt. Sinai J. Med. 67:120132.[Medline]
20 Peng, S.L., M.E. Robert, A.C. Hayday, and J. Craft. 1996. A tumor-suppressor function for Fas (CD95) revealed in T celldeficient mice. J. Exp. Med. 184:11491154.
21 Hartley, J.W., S.K. Chattopadhyay, M.R. Lander, L. Taddesse-Heath, Z. Naghashfar, H.C. Morse III, and T.N. Fredrickson. 2000. Accelerated appearance of multiple B cell lymphoma types in NFS/N mice congenic for ecotropic murine leukemia viruses. Lab. Invest. 80:159169.[Medline]
22 Lossos, I.S., R. Tibshirani, B. Narasimhan, and R. Levy. 2000. The inference of antigen selection on Ig genes. J. Immunol. 165:51225126.
23 Shin, M.S., T.N. Fredrickson, J.W. Hartley, T. Suzuki, K. Agaki, and H.C. Morse III. 2004. High throughput retroviral tagging for identification of genes involved in initiation and progression of mouse splenic marginal zone lymphomas. Cancer Res. 64:44194427.
24 Reap, E.A., M.L. Piecyk, A. Oliver, E.S. Sobel, T. Waldschmidt, P.L. Cohen, and R.A. Eisenberg. 1996. Phenotypic abnormalities of splenic and bone marrow B cells in lpr and gld mice. Clin. Immunol. Immunopathol. 78:2129.[CrossRef][Medline]
25 Mandik-Nayak, L., S.J. Seo, C. Sokol, K.M. Potts, A. Bui, and J. Erikson. 1999. MRL-lpr/lpr mice exhibit a defect in maintaining developmental arrest and follicular exclusion of antidouble-stranded DNA B cells. J. Exp. Med. 189:17991814.
26 Theofilopoulos, A., R.S. Balderas, L. Hang, and F.J. Dixon. 1983. Monoclonal IgM rheumatoid factors derived from MRL/Mp-lpr/lpr mice. J. Exp. Med. 158:901919.
27 Haines, B.B., C.V. Angeles, A.P. Parmelee, P.A. McLean, and P.H. Brodeur. 2001. Germline diversity of the expressed BALB/c VhJ558 gene family. Mol. Immunol. 38:918.[CrossRef][Medline]
28 Brodeur, P.H., and R. Riblet. 1984. The immunoglobulin heavy chain variable region (Igh-V) locus in the mouse. I. One hundred Igh-V genes comprise seven families of homologous genes. Eur. J. Immunol. 14:922930.[Medline]
29 Kantor, A.B., C.E. Merrill, L.A. Herzenberg, and J.L. Hillson. 1997. An unbiased analysis of V(H)-D-J(H) sequences from B-1a, B-1b, and conventional B cells. J. Immunol. 158:11751186.[Abstract]
30 Krishnan, M.R., N.T. Jou, and T.N. Marion. 1996. Correlation between the amino acid position of arginine in VH-CDR3 and specificity for native DNA among autoimmune antibodies. J. Immunol. 157:24302439.[Abstract]
31 Kalled, S.L., and P.H. Brodeur. 1991. Utilization of V
families and V
exons. Implications for the available B cell repertoire. J. Immunol. 147:31943200.[Abstract]
32 Whitcomb, E.A., and P.H. Brodeur. 1998. Rearrangement and selection in the developing V
repertoire of the mouse: an analysis of the usage of two V
gene segments. J. Immunol. 160:49044913.
33 Morse, H.C., T. McCarty, C.-F. Qi, E.A. Torrey, Z. Naghashfar, S.K. Chattopadhyay, T.N. Fredrickson, and J.W. Hartley. 2004. B lymphoid neoplasms of mice: characteristics of naturally occurring and engineered diseases and relationships to human disorders. Adv. Immunol. 81:97121.
34 Nakamuta, M., K. Oka, J. Krushka, K. Kobayashi, M. Yamamoto, W.H. Li, and L. Chan. 1995. Alternative mRNA splicing and differential promoter utilization determine tissue-specific expression of the apolipoprotein B mRNA-editing protein (Apobec1) gene in mice. Structure and evolution of Apobec1 and related nucleoside/nucleotide deaminases. J. Biol. Chem. 270:1304213056.
35 Harris, R.S., S.K. Petersen-Mahrt, and M.S. Neuberger. 2002. RNA editing enzyme APOBEC1 and some of its homologs can act as DNA mutators. Mol. Cell. 10:12471253.[CrossRef][Medline]
36 Tarte, K., F. Zhan, J. De Vos, B. Klein, and J. Shaughnessy. 2003. Gene expression profiling of plasma cells and plasmablasts: toward a better understanding of the late stages of B-cell differentiation. Blood. 102:592600.
37 Underhill, G.H., D. George, E.G. Bremer, and G.S. Kansas. 2003. Gene expression profiling reveals a highly specialized genetic program of plasma cells. Blood. 101:40134021.
38 Hromas, R., S.J. Collins, D. Hickstein, W. Raskind, L.L. Deaven, P. O'Hara, F.S. Hagen, and K. Kaushansky. 1991. A retinoic acid-responsive human zinc finger gene, MZF-1, preferentially expressed in myeloid cells. J. Biol. Chem. 266:1418314187.
39 Li, L., and S.N. Cohen. 1996. Tsg101: a novel tumor susceptibility gene isolated by controlled homozygous functional knockout of allelic loci in mammalian cells. Cell. 85:319329.[CrossRef][Medline]
40 Day, C.D., N.J. Smilinich, G.V. Fitzpatrick, P.J. deJong, T.B. Shows, and M.J. Higgins. 1999. The imprinted domain in mouse distal Chromosome 7: reagents for mutagenesis and sequencing. Mamm. Genome. 10:182185.[CrossRef][Medline]
41 Guehenneux, F., L. Duret, M.B. Callanan, R. Bouhas, R. Hayette, C. Berthet, C. Samarut, R. Rimokh, A.M. Birot, Q. Wang, et al. 1997. Cloning of the mouse BTG3 gene and definition of a new gene family (the BTG family) involved in the negative control of the cell cycle. Leukemia. 11:370375.[CrossRef][Medline]
42 Ge, K., J. Duhadaway, D. Sakamuro, R. Wechsler-Reya, C. Reynolds, and G.C. Prendergast. 2000. Losses of the tumor suppressor BIN1 in breast carcinoma are frequent and reflect deficits in programmed cell death capacity. Int. J. Cancer. 85:376383.[CrossRef][Medline]
43 Su, T.T., B. Guo, Y. Kawakami, K. Sommer, K. Chae, L.A. Humphries, R.M. Kato, S. Kang, L. Patrone, R. Wall, et al. 2002. PKC-beta controls I kappa B kinase lipid raft recruitment and activation in response to BCR signaling. Nat. Immunol. 3:780786.[Medline]
44 Tsukasaki, K., C.W. Miller, E. Greenspun, S. Eshaghian, H. Kawabata, T. Fujimoto, M. Tomonaga, C. Sawyers, J.W. Said, and H.P. Koeffler. 2001. Mutations in the mitotic check point gene, MAD1L1, in human cancers. Oncogene. 20:33013305.[CrossRef][Medline]
45 Moreaux, J., E. Legouffe, E. Jourdan, P. Quittet, T. Reme, C. Lugagne, P. Moine, J. Rossi, B. Klein, and K. Tarte. 2004. BAFF and APRIL protect myeloma cells from apoptosis induced by interleukin 6 deprivation and dexamethasone. Blood. 103:31483157.
46 He, B., A. Chadburn, E. Jou, E.J. Schattner, D.M. Knowles, and A. Cerutti. 2004. Lymphoma B cells evade apoptosis through the TNF family members BAFF/BLyS and APRIL. J. Immunol. 172:32683279.
47 Novak, A.J., J.R. Darce, B.K. Arendt, B. Harder, K. Henderson, W. Kindsvogel, J.A. Gross, P.R. Greipp, and D.F. Jelinek. 2004. Expression of BCMA, TACI, and BAFF-R in multiple myeloma: a mechanism for growth and survival. Blood. 103:689694.
48 Masood, R., Y. Zhang, M.W. Bond, D.T. Scadden, T. Moudgil, R.E. Law, M.H. Kaplan, B. Jung, B.M. Espina, and Y. Lunardi-Iskandar. 1995. Interleukin-10 is an autocrine growth factor for acquired immunodeficiency syndrome-related B-cell lymphoma. Blood. 85:34233430.
49 Lu, Z.Y., X.G. Zhang, C. Rodriguez, J. Wijdenes, Z.J. Gu, B. Morel-Fournier, J.L. Harousseau, R. Bataille, J.F. Rossi, and B. Klein. 1995. Interleukin-10 is a proliferation factor but not a differentiation factor for human myeloma cells. Blood. 85:25212527.
50 Groux, H., and F. Cottrez. 2003. The complex role of interleukin-10 in autoimmunity. J. Autoimmun. 20:281285.[CrossRef][Medline]
51 Lopatin, U., X. Yao, R.K. Williams, J.J. Bleesing, J.K. Dale, D. Wong, J. Teruya-Feldstein, S. Fritz, M.R. Morrow, I. Fuss, et al. 2001. Increases in circulating and lymphoid tissue interleukin-10 in autoimmune lymphoproliferative syndrome are associated with disease expression. Blood. 97:31613170.
52 Rozzo, S.J., J.D. Allard, D. Choubey, T.J. Vyse, S. Izui, G. Peltz, and B.L. Kotzin. 2001. Evidence for an interferon-inducible gene, Ifi202, in the susceptibility to systemic lupus. Immunity. 15:435443.[CrossRef][Medline]
53 Saito, M., K. Yoshida, M. Hibi, T. Taga, and T. Kishimoto. 1992. Molecular cloning of a murine IL-6 receptor-associated signal transducer, gp130, and its regulated expression in vivo. J. Immunol. 148:40664071.[Abstract]
54 Kompfner, E., P. Oliveira, A. Montalbano, and A.J. Feeney. 2001. Unusual germline DSP2 gene accounts for all apparent V-D-D-J rearrangements in newborn, but not adult, MRL mice. J. Immunol. 167:69336938.
55 Klonowski, K.D., and M. Monestier. 2000. Heavy chain revision in MRL mice: a potential mechanism for the development of autoreactive B cell precursors. J. Immunol. 165:44874493.
56 Radic, M.Z., J. Erikson, S. Litwin, and M. Weigert. 1993. B lymphocytes may escape tolerance by revising their antigen receptors. J. Exp. Med. 177:11651173.
57 Brard, F., M. Shannon, E.L. Prak, S. Litwin, and M. Weigert. 1999. Somatic mutation and light chain rearrangement generate autoimmunity in antisingle-stranded DNA transgenic MRL/lpr mice. J. Exp. Med. 190:691704.
58 Li, H., Y. Jiang, E.L. Prak, M. Radic, and M. Weigert. 2001. Editors and editing of anti-DNA receptors. Immunity. 15:947957.[CrossRef][Medline]
59 Monestier, M., and M. Zouali. 2002. Receptor revision and systemic lupus. Scand. J. Immunol. 55:425431.[CrossRef][Medline]
60 Chen, C., Z. Nagy, E.L. Prak, and M. Weigert. 1995. Immunoglobulin heavy chain gene replacement: a mechanism of receptor editing. Immunity. 3:747755.[CrossRef][Medline]
61 Wilson, P.C., K. Wilson, Y.J. Liu, J. Banchereau, V. Pascual, and J.D. Capra. 2000. Receptor revision of immunoglobulin heavy chain variable region genes in normal human B lymphocytes. J. Exp. Med. 191:18811894.
62 Fanning, L., F.E. Bertrand, C. Steinberg, and G.E. Wu. 1998. Molecular mechanisms involved in receptor editing at the Ig heavy chain locus. Int. Immunol. 10:241246.
63 Smit, L.A., R.J. Bende, J. Aten, J.E. Guikema, W.M. Aarts, and C.J. van Noesel. 2003. Expression of activation-induced cytidine deaminase is confined to B-cell non-Hodgkin's lymphomas of germinal-center phenotype. Cancer Res. 63:38943898.
64 Greeve, J., A. Philipsen, K. Krause, W. Klapper, K. Heidorn, B.E. Castle, J. Janda, K.B. Marcu, and R. Parwaresch. 2003. Expression of activation-induced cytidine deaminase in human B-cell non-Hodgkin lymphomas. Blood. 101:35743580.
65 Weng, W.K., and S. Levy. 2003. Hepatitis C virus (HCV) and lymphomagenesis. Leuk. Lymphoma. 44:11131120.[CrossRef][Medline]
66 Isaacson, P.G. 1999. Mucosa-associated lymphoid tissue lymphoma. Semin. Hematol. 36:139147.[Medline]
67 Crump, M., M. Gospodarowicz, and F.A. Shepherd. 1999. Lymphoma of the gastrointestinal tract. Semin. Oncol. 26:324337.[Medline]
68 Dammacco, F., D. Sansonno, C. Piccoli, V. Racanelli, F.P. D'Amore, and G. Lauletta. 2000. The lymphoid system in hepatitis C virus infection: autoimmunity, mixed cryoglobulinemia, and overt B-cell malignancy. Semin. Liver Dis. 20:143157.[CrossRef][Medline]
69 Shlomchik, M.J., C.W. Euler, S.C. Christensen, and J. William. 2003. Activation of rheumatoid factor (RF) B cells and somatic hypermutation outside of germinal centers in autoimmune-prone MRL/lpr mice. Ann. NY Acad. Sci. 987:3850.[Medline]
70 Doyle, H.A., J. Yan, B. Liang, and M.J. Mamula. 2001. Lupus autoantigens: their origins, forms, and presentation. Immunol. Res. 24:131147.[CrossRef][Medline]
71 Okazaki, I.M., H. Hiai, N. Kakazu, S. Yamada, M. Muramatsu, K. Kinoshita, and T. Honjo. 2003. Constitutive expression of AID leads to tumorigenesis. J. Exp. Med. 197:11731181.
72 Yamanaka, S., M.E. Balestra, L.D. Ferrell, J. Fan, K.S. Arnold, S. Taylor, J.M. Taylor, and T.L. Innerarity. 1995. Apolipoprotein B mRNA-editing protein induces hepatocellular carcinoma and dysplasia in transgenic animals. Proc. Natl. Acad. Sci. USA. 92:84838487.
73 Su, T.T., B. Guo, Y. Kawakami, K. Sommer, K. Chae, L.A. Humphries, R.M. Kato, S. Kang, L. Patrone, R. Wall, et al. 2002. PKC-ß controls I
B kinase lipid raft recruitment and activation in response to BCR signaling. Nat. Immunol. 3:780786.[Medline]
74 Pramanik, R., T.N. Jorgensen, H. Xin, B.L. Kotzin, and D. Choubey. 2004. Interleukin-6 induces expression of Ifi202, an interferon-inducible candidate gene for lupus susceptibility. J. Biol. Chem. 279:1612116127.
75 Elliott, K., D. Sakamuro, A. Basu, W. Du, W. Wunner, P. Staller, S. Gaubatz, H. Zhang, E. Prochownik, M. Eilers, and G.C. Prendergast. 1999. Bin1 functionally interacts with Myc and inhibits cell proliferation via multiple mechanisms. Oncogene. 18:35643573.[CrossRef][Medline]
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| TABLE OF CONTENTS |
|