Isolation and characterisation of recombination events ... - Nature

7 downloads 0 Views 155KB Size Report
Acknowledgements. This work was supported by Yorkshire Cancer Research, the. Leukaemia Research Fund, and the Special Trustees of the. Leeds General ...
Leukemia (1999) 13, 1100–1107  1999 Stockton Press All rights reserved 0887-6924/99 $12.00 http://www.stockton-press.co.uk/leu

Isolation and characterisation of recombination events involving immunoglobulin heavy chain switch regions in multiple myeloma using long distance vectorette PCR (LDV-PCR) J Proffitt*, J Fenton*, G Pratt*, Z Yates and G Morgan Institute of Pathology, Algernon Firth Building, University of Leeds, Leeds, LS2 9JT, UK

Immunoglobulin class switching occurs as a result of recombination between pairs of switch region sequences located 5⬘ to each constant heavy chain gene except C␦. In the B cell neoplasm multiple myeloma, tumour cells have generally undergone class switching and often contain oncogenic sequences translocated into switch regions, presumably as a result of aberrant switch recombination. We have developed a method (LDV-PCR) which combines long distance PCR with one-sided vectorette PCR that is capable of detecting and isolating both normal and aberrant switch recombination breakpoints from multiple myeloma cell lines and primary multiple myeloma tumour material. Using LDV-PCR we have directly cloned the translocation breakpoints present in two multiple myeloma cell lines and isolated a normal productive switch recombination event from a primary tumour. Furthermore, we have isolated a novel translocation t(14;22)(q32;q12) from a primary tumour sample and have demonstrated that internal deletions within switch regions can occur in multiple myeloma cells. Compared to a Southern blotting approach, LDV-PCR is simpler and more rapid to perform, allows the simultaneous detection and isolation of recombination events, and can also be applied to amounts of DNA which are too low to permit the conventional cloning of recombination breakpoints. Keywords: multiple myeloma; switch regions; translocation; Ig heavy chain; long vectorette PCR

Introduction B lymphocytes express antigen-binding immunoglobulins (Igs) as complexes containing heavy and light chain components. Ig heavy chains consist of a variable region, produced by the recombination of multiple gene segments, which is involved in antigen recognition, and a constant region which defines the class and biological function of the molecule.1 In humans, the Ig heavy chain gene locus is located on chromosome 14q32 and contains nine functional constant region genes, encoding five isotype classes (IgM, IgD, IgG1–4, IgA1–2 and IgE), located centromeric to the variable encoding VDJ region (Figure 1a). Newly developed B cells initially express IgM (IgD is generally co-expressed with IgM as a result of mRNA differential splicing). However, after contact with antigen, they are capable of switching class from IgM/IgD to either IgG1–4, IgA1– 2 or IgE, via a process which occurs in the germinal centres of lymphoid follicles.2–4 This switch of heavy chain isotype is achieved by a DNA recombination event involving switch (S) regions, segments of DNA several kilobases in size which lie 5⬘ to each constant region gene except C␦ and consist of multiple G/C-rich pentameric repeats. Legitimate recombination between an upstream S region (generally S␮) and a downstream S region results in deletion of the intervening DNA,

Correspondence: G Morgan; Fax: 0113 2333404 *Joint first authors Received 7 January 1999; accepted 13 April 1999

normally in the form of a circle, and the repositioning of an unchanged VDJ region in front of a new constant region.5–7 Chromosomal translocations involving the IgH locus at 14q32 are commonly seen in a variety of B cell neoplasms.8 It is believed that such rearrangements cause overexpression of oncogenes as a result of their coupling to strong enhancer/locus control region (LCR) sequences present within the IgH region. Examples include translocation of the c-myc gene by t(8;14)(q24;q32) in Burkitt’s lymphoma9 and translocation of the BCL 1/cyclin D1 locus by t(11;14)(q13;q32) in mantle cell lymphoma.10 Multiple myeloma (MM) is a B cell neoplasm, the clinical features of which result from the accumulation of a single clone of Ig producing plasma cells in the bone marrow.11 Despite the fact that cytogenetic analyses in MM are difficult to perform, probably due to the low proliferation rate of malignant plasma cells, 14q32 translocations have been identified in 20–60% of cytogenetically abnormal cases.12–14 Furthermore, a more sensitive FISH assay of MM tumours has revealed 14q32 translocation in up to 74% of cases.15 These translocations appear to involve a promiscuous array of partner chromosomes, although to date, the actual identity of many of these partners remains to be defined. In contrast to many B cell tumours where translocations involving 14q32 generally occur within the VDJ region, rearrangements at 14q32 in MM appear to cluster within switch region sequences, probably as a result of aberrant (illegitimate) switch recombination.16,17 This is consistent with MM being a post germinal centre B cell neoplasm in which the cell of origin has passed through the stage of development when class switch recombination is active. Recently, several groups have developed genomic Southern blot assays capable of detecting illegitimate translocations into switch region sequences in MM.16–20 Several of these identified breakpoints have been cloned mainly from MM cell lines but also from a smaller number of primary tumour samples. Analysis of these breakpoints has identified several putative oncogenic sequences that are translocated and concomitantly overexpressed including BCL1/cyclin D1,16,17 FGFR3,20,21 and cmaf.22 However, it is clear that a significant number of translocation breakpoints remain to be identified. We have developed a vectorette-based long distance PCR strategy (LDV-PCR: see Methods) using primers which bind 3⬘ to S region sequences to successfully isolate legitimate and illegitimate switch recombination breakpoints from both MM cell lines and primary MM tumour material. Due to the orientation of the IgH locus with respect to the centromere of chromosome 14, primers flanking the 3⬘ sides of S regions should detect the 14q+ component of any illegitimate recombination event (Figure 1b and 1c). This assay is more simple and rapid to perform than Southern blotting and allows the simultaneous detection and isolation of recombination breakpoint regions. Additionally, the technique can be applied to levels of DNA which are too low to permit Southern blotting and the sub-

Switch recombination in multiple myeloma J Proffitt et al

1101

Figure 1 Diagram of the IgH locus showing the location of primers used in LDV-PCR. (a) Map of the IgH locus following VDJ recombination but prior to switch recombination. In B cells, functional heavy chain is expressed from just one of the two copies of the IgH locus present (the productive allele). (b) Diagram indicating the location of nested primer pairs m1:m2, g1:g2 and a1:a2 relative to switch repeat regions ␮, ␥1– 4 and ␣1–2, respectively. Due to strong sequence homology primer pair g1:g2 binds 3⬘ to S␥1–4, and primer pair a1:a2 binds 3⬘ to both S␣1 and S␣2. (c) Schematic diagram detailing the generation of possible LDV-PCR products for primer pair g1:g2. Germline DNA or recombined DNA is digested with restriction enzyme R and vectorette units V are ligated. Two rounds of nested PCR using primers directed towards V in combination with g1 and g2 generate LDV-PCR products containing the DNA between g2 and the closest 5⬘ recognition site for R.

sequent generation of genomic libraries for conventional cloning.

Materials and methods

Cell culture Human MM cell lines JJN3 and U266B1 were grown (37°C, 5% Co2) in RPMI 1640 medium containing glutamax-1 (GIBCO BRL, Paisley, UK) and supplemented with 10% heatinactivated fetal bovine serum, 100 U/ml penicillin, and 100 ␮g/ml streptomycin. Cells were maintained in culture at densities between 0.5 and 1.5 × 106 cells/ml.

Preparation of genomic DNA High molecular weight genomic DNA was prepared from MM cell lines and MM patient samples by proteinase K digestion, phenol–chloroform extraction, and ethanol precipitation as described by Sambrook et al.23 Peripheral blood samples and bone marrow aspirates were collected from MM patients in EDTA.

Construction of vectorette libraries Standard PCR requires sequence information from both ends of the fragment to be amplified. Vectorette PCR is a method which allows DNA fragments to be amplified when the sequence of only one end is known.24 Such one-sided PCR allows the amplification of any uncharacterised sequence adjacent to a known region. The technique works by ligating short synthetic double-stranded oligonucleotide linkers (termed vectorette units, these are identical in core sequence but are made with ends compatible for different restriction enzymes) to DNA which has been digested with a particular restriction enzyme. Two rounds of nested PCR can then be performed using primers directed towards these vectorette units in combination with known sequence-specific primers. This generates a PCR product containing the DNA between the inner specific primer and the closest recognition site for the enzyme used for target DNA digestion. Specifically, we digested 1 ␮g of target genomic DNA separately with 20 U of BamHI, ClaI, EcoRI, HindIII, PvuII and RsaI in a 50 ␮l reaction volume as recommended by the enzyme manufacturer (GIBCO BRL). Ligation of the appropriate vectorette unit (Genosys Biotechnologies, Pampisford, UK) to the digested DNA was then achieved by the addition of 3 pmol vectorette units, 1 ␮l 100 mm ATP, 1 ␮l 100 mm DTT, and 1 U T4 DNA ligase (GIBCO BRL), and incubation at 20°C for 60 min followed by 37°C for 30 min. This temperature cyc-

Switch recombination in multiple myeloma J Proffitt et al

1102

ling was repeated three times in total. Restriction enzyme binding sites are reformed when target DNA fragments ligate to each other but not when they ligate to the appropriate vectorette unit, and therefore this temperature cycling increases the proportion of target DNA correctly linked to vectorette units. After incubation, 200 ␮l of sterile water was added to the reaction and the resulting vectorette library stored in aliquots at −20°C.

PCR primers and long distance PCR amplification of vectorette libraries (LDV-PCR) Primers were designed to bind to the 3⬘ flanking sequences of S␮,25 S␥1–4,26 and S␣1–2.27,28 Primers m1 (5⬘-TGTCAGGGCTGAGTGGCGTTC-3⬘) and m2 (5⬘-CGTTCTGAGTGCCCTCACTACTTG-3⬘) anneal 3⬘ of S␮, primers g1 (5⬘TCCCAGTGTCCTGCATTACTTC-3⬘) and g2 (5⬘-TTCTGGAGGCTCAGTCACCAC-3⬘) anneal 3⬘ of S␥1–4, and primers a1 (5⬘-AACCCAGCACAGAGAGGCCT-3⬘) and a2 (5⬘-AGAGAGGCCTGGTGACAGCC-3⬘) anneal 3⬘ of S␣1–2. Primers m1, g1 and a1 were used separately in combination with a vectorette unit binding primer (Genosys) to amplify vectorette libraries. Each PCR reaction contained 200 ␮m each dNTP, 60 mm Tris-SO4 (pH 9.1), 18 mm (NH4,2SO4 1.4 mm MgSO4, 12.5 pmoles each primer, 2 ␮l vectorette library, and 1 ␮l Elongase enzyme mix (a mixture of thermostable DNA polymerases optimised for long distance PCR amplification, GIBCO BRL) in a 50 ␮l volume under mineral oil. After an initial incubation at 95°C for 3 min, each reaction was subjected to 30 cycles consisting of 95°C for 45 s, 60°C for 45 s, and 68°C for 5 min. After the final cycle each reaction was incubated at 68°C for 10 min. A second round of nested PCR was performed using primers m2, g2, and a2 in combination with a vectorette unit binding nested primer (Genosys). First round PCR reactions were diluted 1/1000 and 2 ␮l of this dilution was used as a template in PCR reactions performed exactly as detailed above except that 35 cycles were performed.

Cloning and DNA sequence analysis of LDV-PCR products LDV-PCR amplification products were purified from 1% agarose gels using an Agarose Gel DNA Extraction Kit (Boehringer, Lewes, UK) according to the manufacturer’s protocol. Purified fragments were blunt-end cloned into the plasmid pT7Blue using a Perfectly Blunt Cloning Kit (Novagen, Nottingham, UK). Fluorescent DNA sequencing of purified PCR fragments (60–90 ng) using internally hybridising primers, or sequencing of plasmid subclones (300–400 ng) using pUC/M13 universal and reverse sequencing primers was performed using a Dye Terminator Cycle Sequencing Kit (Perkin Elmer, Warrington, UK) and an ABI 373A DNA Sequencer according to the manufacturer’s instructions. DNA database homology searches were performed using the BLAST program (NCBI internet site).

Standard DNA amplifications A 7.5 kb fragment specific for the ␤-globin gene was amplified using the following primers: Beg1 (5⬘-CTGCTGAAAGAGATGCGGTGG-3⬘) and Beg2 (5⬘-GCACTGGCTTAG-

GAGTTGGACT-3⬘), and 35 cycles of the following conditions: denaturation at 95°C for 45 s, annealing at 58°C for 45 s, and extension at 68°C for 8 min. Each 50 ␮l PCR reaction contained 200 ␮m dNTPs, 60 mm Tris-SO4 (pH 9.1), 18 mm (NH4)2SO4 2 mm MgSO4, 12.5 pmoles each primer, 100 ng genomic DNA, and 1 ␮l Elongase. Amplification of genomic DNA sample (100 ng/reaction) or a human monochromosomal somatic cell hybrid DNA panel (40 ng/reaction)29 was performed using primers gmp0–4 as described above but using 30 cycles of the following conditions: 95°C for 30 s, 56°C for 30 s, and 68°C for 30 s. Results

Development of long distance vectorette-PCR (LDVPCR) using primers located 3⬘ to IgH switch regions We have designed a strategy which combines ligation of vectorette linker units to restriction endonuclease-digested DNA fragments with long distance PCR (LDV-PCR) to amplify and clone recombination events occurring within IgH locus switch region sequences (Figure 1a). Due to strong sequence homology around the regions encoding all four ␥ isotypic subclasses and the two ␣ subclasses, seven out of the eight classical switch regions can be screened by LDV-PCR for rearrangements using just three pairs of nested primers (m1:m2, g1:g2 and a1:a2 which bind specifically 3⬘ to S␮, S␥1–4 and S␣1–2 regions, respectively) in combination with a common pair of nested vectorette primers (Figure 1b and c). When an appropriate restriction enzyme recognition site is within range, LDV-PCR will amplify a band of predictable size from germline DNA, and in informative cases, a band of different size from DNA template that has been subject to recombination. Sequencing of non-germline bands allows legitimate recombination and illegitimate recombination events to be differentiated (Figure 1c). PvuII vectorette libraries were constructed using genomic DNA isolated from myeloma cell lines, myeloma patient blood and bone marrow, and normal blood as a germline control. Nested vectorette-PCR amplification of these libraries with primer paris m1:m2, g1:g2 and a1:a2 produced single bands of the expected size, 250 bp, 750 bp and 610 bp, respectively, indicating that these primers function correctly in vectorette PCR (data not shown). Similarly the predicted sizes of PCR products following nested vectorette-PCR amplification of genomic germline DNA for all six vectorette libraires and the three primer pairs have been calculated (Table 1). Additionally, DNA isolated from all sources is a suitable template for long distance PCR as amplification with primers (Beg1 and Beg2) which bind within the ␤-globin gene reproducibly generated a 7.5 kb product (data not shown). Long distance PCR and vectorette PCR can be successfully combined as demonstrated by the generation of specific amplification products between approximately 2 kb and 5 kb in size using the 3⬘S␮, 3⬘S␥ and 3⬘S␣ specific primer pairs and vectorette libraries prepared from MM cell line or patient genomic DNA (Figure 2 and Figure 3, lanes 2, 4 and 6).

LDV-PCR amplification of switch region translocations from multiple myeloma cell lines Two MM cell lines, JJN3 and U266B1, have been previously shown to contain translocations of chromosomes 16 and 11

Switch recombination in multiple myeloma J Proffitt et al

Table 1

Switch region

S␮ S␥1 S␥2 S␥3 S␥4 S␣1–2

1103

LDV-PCR amplification of genomic (control) DNA

Vectorette library HindIII

BamHI

EcoRI

RsaI

ClaI

PvuII

4960 4580 3240 3810 3020 NA

NA 3230 NA 2460 1720 NA

NA NA NA NA NA NA

3280 2380 1650 1340 1470 860

NA NA NA NA NA NA

250 750 750 750 740 610

The predicted size of PCR products after two rounds of LDV-PCR using the six vectorette libraries and primer pairs m1:m2 (␮), g1:g2 (␥) and al:a2 (␣) are shown to the nearest 10 bp. When an appropriate restriction enzyme recognition site is not within range (⬎5 kb) there is no value shown (NA).

Figure 3 Isolation of translocations into switch region sequences for two MM cell lines. Vectorette libraries were constructed using DNA isolated from normal blood (CONTROL), U266B1 MM cell line (U266), and JJN3 MM cell line (JJN3) using the restriction enzymes indicated. CONTROL and U266 libraries were subject to amplification using primer pair a1:a2 (lanes 1 and 2). CONTROL and JJN3 libraries were amplified using primer pair m1:m2 (lanes 3–6). The 5.0 kb JJN3/ClaI product (lane 6) was digested with RsaI (lane 7). Figure 2 LDV-PCR amplification of genomic DNA. HindIII, BamHI and RsaI vectorette libraries were made from genomic DNA isolated from MM patient bone marrow (BM), and MM patient blood (PB). Nested LDV-PCR was performed on these libraries using primer pair g1:g2. PCR products were analysed as described above.

into IgH switch regions ␮ and ␣1 respectively using a genomic Southern blot assay.17 To demonstrate that these same translocation breakpoints could be detected and directly isolated using LDV-PCR, vectorette libraries (BamHI, ClaI, EcoRI, HindIII, PvuII and RsaI) were prepared using genomic DNA isolated from JJN3 cells, U266B1 cells and normal blood, and used in LDV-PCR with primer pairs m1:m2, g1:g2 and a1:a2. For cell line JJN3, amplification of DNA using primer pair m1:m2 from the ClaI and RsaI vectorette libraries reproducibly gave rise to fragments 5.0 kb and 3.8 kb in size, respectively, which were not produced from the corresponding germline vectorette libraries (Figure 3). Digestion of the ClaI product with RsaI produced a 3.8 kb product of exactly the same size as the RsaI fragment strongly suggesting that PCR is consistently initiating from the same position (Figure 3, lane 7). The JJN3-specific 3.8 kb (RsaI) product was cloned into the vector pT7Blue and subjected to fluorescent DNA sequencing. Sequence from the vectorette end of the fragment matched exactly with that reported as the JJN3 t(14;16)(q32;q23) breakpoint by Bergsagel et al17 (data not shown). Amplification of the U266B1 EcoRI vectorette library with primer pair a1:a2 produced a 3.5 kb band which was not seen with the germline EcoRI vectorette library (Figure 3). Direct sequencing of this 3.5 kb fragment with a primer designed to

hybridise upstream to the reported U266B1 breakpoint sequence confirmed that this fragment contained 11q13 sequences recombined to S␣1 (data not shown).

Cloning non-germline switch region DNA fragments from multiple myeloma patient bone marrow samples In MM patients a single clone of malignant plasma cells accumulates in the bone marrow (generally 10% of cells or greater). These cells are present at much lower levels in the peripheral blood. Performing LDV-PCR using the 3⬘S␮, 3⬘S␥ and 3⬘S␣ specific primer pairs on vectorette libraries constructed from genomic DNA isolated from MM patient bone marrow and matched blood as a control should allow detection and isolation of switch recombination events (both legitimate and illegitimate) which have occurred in the myeloma clone against a background of germline bands which will be common to both blood and bone marrow. Vectorette libraries (BamHI, ClaI, EcoRI, HindIII, PvuII, RsaI) were constructed using DNA isolated from bone marrow and peripheral blood samples obtained from four patients with multiple myeloma. Patients 1, 2 and 3 were shown to be IgG producing and patient 4 to be IgA producing. All four cases contained greater than 20% malignant plasma cells within the bone marrow as estimated by immunophenotypic analysis (data not shown). Using LDV-PCR bone marrow specific bands could be repeatably obtained for all four patients (Figure 4). For patient 1, amplification of HindIII libraries with primer

Switch recombination in multiple myeloma J Proffitt et al

1104

Figure 4 LDV-PCR analysis of DNA isolated from peripheral blood and bone marrow of multiple myeloma patients. Vectorette libraries were constructed with genomic DNA isolated from blood (PB) and bone marrow (BM) samples obtained from four patients with multiple myeloma. With patients 1, 2 and 3, LDV-PCR amplification of HindIII vectorette libraries with primer pair g1:g2 produced fragments specific to the bone marrow 1.9 kb, 1.8 kb and 2.7 kb in size, respectively (marked with arrows), in addition to the expected 3.3 kb germline band present in both blood and bone marrow. For patient 4, amplification of (PB) and (BM) HindIII libraries with primer pair m1:m2 produced a 2.3 kb bone marrowspecific fragment (marked with an arrow).

pair g1:g2 produced a 1.9 kb fragment specific to the bone marrow, in addition to an expected 3.3 kb germline band present in both blood and bone marrow (Figure 4). The 1.9 kb bone marrow-specific fragment was cloned and partially sequenced. DNA sequence obtained from the vectorette end of the fragment matched exactly with the reported germline sequence of the S␮ flanking region starting at a HindIII site located approximately 900 bp upstream from the start of the S␮ pentameric repeats (data not shown). Sequencing the other end of the product revealed a perfect match with known germline S␥1 3⬘ flanking sequence initiating at the g2 primer binding site. These data imply that the 1.9 kb HindIII-vectorette product represents a myeloma clone-specific legitimate switch from ␮ to ␥1. This is consistent with patient 1 being an IgG producing myeloma. Non-germline, bone marrow-specific fragments 2.3 kb and 2.7 kb in size were also generated for patients 3 and 4 respectively with HindIII vectorette libraries (Figure 4). Partial sequence analysis of the 2.3 kb patient 4-specific product identified 5⬘ S␮ sequence initiating from the same 5⬘ HindIII site used by the 1.9 kb patient 1-specific product at the the vectorette end of the fragment and 3⬘ S␮ sequence starting from the m2 primer at the opposite end. Additionally, sequencing the 2.7 kb patient 3-specific product characterised 5⬘ S␥2 sequence starting from a HindIII site at the vectorette end of the fragment and 3⬘ S␥2 sequence from the primer g2 at the other. In germline, primers m2 and g2 are separated from the closest 5⬘ HindIII sites by approximately 5.0 kb and 3.3 kb of intervening S␮ and S␥2 DNA, respectively. Therefore, it would appear that internal deletions of about 2.7 kb and 0.6 kb in S␮ and S␥2 have occurred within the myeloma clone present in patients 3 and 4, either by the looping out and deletion of DNA within a single switch region, or by trans-recombination involving the same switch region from both the productive and non-productive IgH alleles. Internal switch region deletions have not been reported in MM cells to date. However, a deletion of several kilobases within S␮ has been observed in a B cell line derived from BALB/c mice.30

Analysis of a novel t(14;22) breakpoint isolated from a primary multiple myeloma tumour sample LDV-PCR using primers g1:g2 on HindIII vectorette libraries made from patient 2 blood and bone marrow resulted in a

bone marrow-specific 1.8 kb product, almost identical in size to the 1.9 kb product produced using the same primers for patient 1 (Figure 4). However, sequencing of this 1.8 kb fragment from the vectorette end identified 84 bp of sequence which demonstrated no significant homology to any published chromosome 14 sequences. This 84 bp region was fused to a sequence which matches exactly with the reported germline 5⬘ S␮ flanking sequence starting from a point approximately 300 bp 5⬘ to the start of the S␮ pentameric repeat region (Figure 5). Analysis of sequence data generated from the opposite end of this 1.8 kb product identified a 100% match (except for a small region containing deletions of 1 bp and 7 bp: see Discussion) to germline S␥2 3⬘ flanking sequence starting at the g2 primer (Figure 5). This implies that this 1.8 kb fragment has the structure 5⬘-84 bp unknown sequence (Ch ?)-␮ switch sequence (5⬘S␮)-␥2 switch sequence (3⬘ S␥2)-3⬘ and has been generated by at least two separate recombination events (see Figure 5). Two PCR primers (gmp1 and gmp2) were designed to bind to this unknown 84 bp sequence, and PCR was performed on a human monochromosomal somatic cell hybrid DNA panel.29 A single band of the expected size was obtained from human positive control DNA and from the hybrid containing human chromosome 22, but not from mouse and hamster negative control DNA or from any of the other hybrids containing chromosomes 1–21, X and Y (data not shown). Consistent with this result, a DNA database homology search revealed that the 5⬘ 84 bp sequence obtained from the 1.8 kb patient 2 product matched 100% to a chromosome 22 PAC clone (DJ130H16) which has recently been sequenced by the St Louis Genome Sequencing Center. The PAC DJ130H16 sequence has been mapped to 22q12 (Dr J Collins, Sanger Centre, personal communication), which indicates that the patient 2 translocation can be further defined as t(14;22)(q32;q12). Interestingly, analysis of PAC DJ130H16 sequence reveals that the breakpoint in patient 2 occurs within a gene, specifically between exons 2 and 3 of a 12 exon gene which shows similarity to squid retinal-binding protein.31 Furthermore, a single round of PCR with a primer (gmp0) designed to bind to PAC clone DJ130H16 sequence upstream from that contained within the 1.8 kb patient 2specific product and either of two primers (gmp3 and gmp4) which hybridise to the 5⬘ S␮ flanking region (Figure 5) produced bands of 155 bp and 288 bp when used with patient 2 bone marrow DNA, but was negative with patient 2 blood or normal germline DNA (data not shown). These results strongly suggest that the 1.8 kb HindIII-vectorette product represents a

Switch recombination in multiple myeloma J Proffitt et al

Figure 5 Structure and partial sequence of the 1.8 kb patient 2specific product. The structure of the 1.8 kb fragment is shown as a schematic diagram. The amount of DNA sequence obtained from each end of the fragment is indicated by horizontal lines (a) and (b). The exact location of the 5⬘S␮-3⬘S␥2 breakpoint has not yet been determined. (a) Sequence obtained from the vectorette end of the 1.8 kb patient-specific fragment is shown (lower case, accession number AJ010127) aligned with sequence from chromosome 22 PAC clone DJ130H16 (upper case, accession number AC004997). The breakpoint from 100% match with DJ130H16 sequence to 100% match with germline 5⬘S␮ sequence25 is marked by vertical arrows (there is 2 bp homology between DJ130H16 and 5⬘S␮ sequences at this point). PCR primers gmp0–4 are marked in bold and underlined, sequence corresponding to the ligated vectorette unit is shown in bold, and the HindIII recognition site is overlined. (b) Sequence from the 3⬘ end of the 1.8 kb patient 2-specific fragment is shown in lower case. This sequence matches exactly with that reported for germline 3⬘S␥226 except for two deletions of 1 bp and 7 bp (upper case). Sequence corresponding to PCR primer g2 is marked in bold and underlined.

breakpoint specific for the myeloma clone in patient 2 for a chromosome 14–22 translocation. To date, translocation of chromosome 22 sequences into the IgH locus has not been reported in any myeloma cell lines or primary myeloma tumours. Discussion We have applied LDV-PCR to isolate legitimate and illegitimate switch recombination breakpoints from both MM cell

lines and primary MM tumour material. A number of other B cell neoplasms such as sporadic Burkitt’s lymphoma and HIVassociated non-Hodgkin’s lymphoma can present with translocations in Ig switch region sequences.32 Therefore, it should also be possible to apply LDV-PCR to detect and isolate breakpoints in these cases. It should be noted that although translocation breakpoints in the IgH locus can be identified, the actual deregulated oncogene is often located many kb away from the breakpoint so no part of the gene need be contained on the vectorette sequence isolated. Our LDV-PCR data predict the presence of a novel t(14;22)(q32;q12) breakpoint in patient 2. To date, such a translocation has not been detected cytogenetically in MM cases. However, cytogenetics fails to detect a number of 14q+ translocations which have been identified at the molecular level as occurring repetitively in MM,20–22 including translocation of the FGFR3 gene by t(4;14)(p16;q32) and translocation of the c-maf gene via t(14;16)(q32;q23). It is possible that the extreme telomeric position of the IgH gene could explain the failure of classical cytogenetic techniques to detect certain translocations involving this locus in MM, and that this could be the case for t(14;22)(q32;q12) as well. A recent comparative genomic hybridisation analysis of chromosomal gains and losses in MM identified gains of 22q occurring in approximately 10% of cases.33 This is consistent with the presence of one or more oncogenes important in the pathogenesis of MM on the long arm of chromosome 22 that could be activated by translocation to 14q32. A recent study by Stec et al34 has shown that, in addition to causing FGFR3 overexpression, a number of t(4;14)(p16;q32) translocations in MM also result in the merging of a candidate gene for Wolf–Hirschhorn syndrome (WHSC1) with the IgH locus and the possible generation of WHSC1-IgH fusion transcripts. Since the breakpoint in patient 2 occurs within a gene we are currently investigating whether t(14;22)(q32;q12) results in the expression of hybrid or truncated transcripts involving this gene and the IgH locus. To our knowledge, t(14;22)(q32;q12) represents the first translocation in MM to be isolated initially from a primary tumour rather than a MM cell line. Analyses of translocations in MM cell lines are likely to contain a bias towards those favouring the formation of cells capable of growth in vitro. For example, t(11;14)(q13;q32) involving the BCL1/cyclin D1 oncogene has been detected in more than 25% of MM cell lines analysed to date,17,18 but a recent study by Fonseca et al35 observed t(11;14)(q13;q32) in only 13 of 263 (4.9%) primary MM cases, and of these only three patients had the translocation at the time of initial presentation. The ability to rapidly isolate breakpoints from primary tumour material will help to determine the full range and prevalence of translocations in MM. The t(14;22) breakpoint in patient 2 is located approximately 300 bp upstream from the start of the S␮ repeats. Other studies have previously shown recombination sites to be distributed throughout the S␮ region and its flanking regions.36,37 Our finding is consistent with their conclusion that recombination does not necessarily have to take place with the pentameric repeat region iself. Downstream of the t(14;22) breakpoint there is a joining of 5⬘S␮ to 3⬘S␥2 sequences indicating that two recombination events have occurred in patient 2, presumably on the non-productive IgH allele as myeloma cells from patient 2 still express monoclonal IgG. It is not possible to determine the order in which these two recombination events occurred. However, work on murine plasmacytoma cell lines (an animal model for MM) has established that c-

1105

Switch recombination in multiple myeloma J Proffitt et al

1106

myc often translocates first into S␮, and then ‘remodels’ into the structure c-myc-S␮-S␣ via a second recombination event.38 This results in a closer juxtaposition of c-myc to the 3⬘ IgH LCR, and therefore it has been postulated that this remodeling results in enhanced expression of c-myc. It is possible that a similar selection event has happened in patient 2, which would have important implications for the pathogenesis of MM. Sequence analysis of the vectorette end of the 1.8 kb patient 2-specific product revealed perfect matches with reported germline chromosome 22 (PAC clone DJ130H16) and chromosome 14 (5⬘S␮) sequences. However, sequencing from the opposite end 3⬘S␥2) of the product identified two deletions of 1 bp and 7 bp compared to the germline sequence reported for 3⬘S␥2.26 The presence of these deletions were confirmed by the direct sequencing of patient 2-specific products generated from several independently performed PCR reactions and therefore are unlikely to be the result of errors during PCR amplification. Mutations including deletions and substitutions have been found in a number of recombined switch regions and it has been suggested that error-prone DNA replication is involved in the process of switch recombination.30,36,39 The molecular mechanism of immunoglobulin switch recombination is not clear and contradictory evidence supporting different models has been produced. For example, Li et al39 have analysed extrachromosomal constructs containing two switch regions which have been transfected into B cell lines and observed mutations on either but not on both sides of the switch junctions produced. In contrast, analysis of S␮/S␥3 recombination breakpoints isolated by PCR from mitogen-stimulated murine B cells identified mutations on both sides of the crossover point.30 Both these studies analyse recombination breakpoints which have been generated in vitro. LDV-PCR allows the generation of PCR products containing switch recombination breakpoints that have been generated in vivo (eg the 1.9 kb patient 1-specific product which represents a S␮/S␥1 recombination). Analysis of this product and the isolation of other legitimate switches that have occurred in germinal centres could help to distinguish between the different models developed to explain how switch recombination occurs. In conclusion, we have demonstrated the ability of LDVPCR to rapidly isolate and characterise both legitimate and illegitimate breakpoints from switch recombination events. This should assist both in the elucidation of the mechanism of normal switch recombination, and in the identification of novel oncogenic sequences which are translocated to the heavy chain locus in MM and possibly other B cell neoplasms. Acknowledgements This work was supported by Yorkshire Cancer Research, the Leukaemia Research Fund, and the Special Trustees of the Leeds General Infirmary. We also thank the UK HGMP Resource Centre for providing the human somatic cell hybrid DNA panel, Fraser Lewis for the rapid synthesis of oligonucleotide primers, and Gina Scott for help with fluorescent DNA sequencing. References 1 Schatz DG, Oettinger MA, Schlissel MS. V(D)J recombination: molecular biology and regulation. Ann Rev Immunol 1992; 10: 359–383.

2 Harriman W, Volk H, Defranoux N, Wabl M. Immunoglobulin class switch recombination. Ann Rev Immunol 1993; 11: 361– 384. 3 Wabl M, Steinberg C. Affinity maturation and class switching. Curr Opin Immunol 1996; 8: 89–92. 4 Stavnezer J. Immunoglobulin class switching. Curr Opin Immunol 1996; 8: 199–205. 5 vonSchwedler U, Jack H, Wabl M. Circular DNA is a product of the immunoglobulin class switch rearrangement. Nature 1990; 345: 452–456. 6 Matsuoka M, Yoshida K, Maeda T, Usuda S, Sakano H. Switch circular DNA formed in cytokine-treated mouse splenocytes: evidence for intramolecular DNA deletion in immunoglobulin class switching. Cell 1990; 62: 135–142. 7 Iwasato T, Shimizu A, Honjo T, Yamagishi H. Circular DNA is excised by immunoglobulin class switch recombination. Cell 1990; 62: 143–149. 8 Korsmeyer SJ. Chromosomal translocations in lymphoid malignancies reveal novel protooncogenes. Ann Rev Immunol 1992; 10: 785–807. 9 Dalla-Favera R, Bregni M, Erickson D, Patterson D, Gallo RC, Croce CM. Human c-myc oncogene is located on the region of chromosome 8 that is translocated in Burkitt lymphoma cells. Proc Natl Acad Sci USA 1982; 79: 7824–7827. 10 Williams ME, Swerdlow SH, Rosenberg CL, Arnold A. Chromosome 11 translocation breakpoints at the PRAD1/CyclinD1 gene locus in centrocytic lymphoma. Leukemia 1993; 7: 241–245. 11 Hallek M, Bergsagel PL, Anderson KC. Multiple myeloma: increasing evidence for a multistep transformation process. Blood 1998; 91: 3–21. 12 Nishida K, Taniwaki M, Misawa S, Abe T. Nonrandom rearrangement of chromosome 14 at band q32.33 in human lymphoid malignancies with mature B cell phenotype. Cancer Res 1989; 49: 1275–1281. 13 Taniwaki M, Nishida K, Takashima T, Nakagawa H, Fujii H, Tamaki T, Shimazaki C, Horiike S, Misawa S, Abe T, Kashima K. Nonrandom chromosomal rearrangements of 14q32.3 and 19p13.3 and preferential deletion of 1p in 21 patients with multiple myeloma and plasma cell leukemia. Blood 1994; 84: 2283–2290. 14 Calasanz MJ, Cigudosa JC, Odero MD, Ferreira C, Ardanaz MT, Fraile A, Carrasco JL, Sole F, Cuesta B, Gullon A. Cytogenetic analysis of 280 patients with multiple myeloma and related disorders: primary breakpoints and clinical correlations. Genes Chromosomes Cancer 1997; 18: 84–93. 15 Nishida K, Tamura A, Nakazawa N, Ueda Y, Abe T, Matsuda F, Kashima K, Taniwaki M. The Ig heavy chain gene is frequently involved in chromosomal translocations in multiple myeloma and plasma cell leukemia as detected by in situ hybridization. Blood 1997; 90: 526–534. 16 Chesi M, Bergsagel PL, Brents LA, Smith CM, Gerhard DS, Kuehl WM. Dysregulation of cyclin D1 by translocation into an IgH gamma switch region in two multiple myeloma cell lines. Blood 1996; 88: 674–681. 17 Bergsagel PL, Chesi M, Nardini E, Brents LA, Kirby SL, Kuehl WM. Promiscuous translocations into immunoglobulin heavy chain switch regions in multiple myeloma. Proc Natl Acad Sci USA 1996; 93: 13931–13936. 18 Bergsagel PL, Nardini E, Brents L, Chesi M, Kuehl WM. IgH translocations in multiple myeloma: a nearly universal event that rarely involves c-myc. Curr Topics Microbiol Immunol 1997; 224: 283–287. 19 Iida S, Rao PH, Butler M, Corradini P, Boccadoro M, Klein B, Chaganti RSK, Dalla-Favera R. Deregulation of MUM1/IRF4 by chromosomal translocation in multiple myeloma. Nat Genet 1997; 17: 226–230. 20 Richelda R, Ronchetti D, Baldini L, Cro L, Viggiano L, Marzella R, Rocchi M, Otsuki T, Lombardi L, Maiolo AT, Neri A. A novel chromosomal translocation t(4;14)(p16.3;q32) in multiple myeloma involves the fibroblast growth factor receptor 3 gene. Blood 1997; 90: 4062–4070. 21 Chesi M, Nardini E, Brents LA, Schrock E, Ried T, Kuehl WM, Bergsagel PL. Frequent translocation t(4;14)(p16.3;q32.3) in multiple myeloma is associated with increased expression and activat-

Switch recombination in multiple myeloma J Proffitt et al

22

23 24

25 26 27 28 29

30

31

ing mutations of fibroblast growth factor receptor 3. Nat Genet 1997; 16: 260–264. Chesi M, Bergsagel PL, Shonukan OO, Martelli ML, Brents LA, Chen T, Schrock E, Ried T, Kuehl WM. Frequent dysregulation of the c-maf proto-oncogene at 16q23 by translocation to an Ig locus in multiple myeloma. Blood 1998; 91: 4457–4463. Sambrook J, Fritsch EF, Maniatis T. Molecular Cloning: A Laboratory Manual. Cold Spring Harbor Laboratory Press, Cold Spring Harbor, NY, 1989. Riley J, Butler R, Ogilvie D, Finniear R, Powell S, Anand R, Smith JC, Markham AF. A novel rapid method for the isolation of terminal sequences from yeast artificial chromosome clones. Nucleic Acids Res 1990; 18: 2887–2890. Mills FC, Brooker JS, Camerini-Otero RD. Sequences of human immunoglobulin switch regions: implication for recombination and transcription. Nucleic Acids Res 1990; 18: 7305–7316. Mills FC, Mitchell MP, Harindranath N, Max EE. Human IgS gamma regions and their participation in sequential switching to IgE. J Immunol 1995; 155: 3021–3036. Islam K, Baskin B, Nilsson L, Hammarstrom L, Sideras P, Smith CIE. Molecular analysis of IgA deficiency: evidence for impaired switching to IgA. J Immunol 1994; 152: 1442–1452. Flanagan JG, Lefranc MP, Rabbitts TH. Mechanisms of divergence and convergence of the human immunoglobulin alpha 1 and alpha2 constant region gene sequences. Cell 1984; 36: 681–688. Kelsall DP, Rooke L, Warne D, Bonzyk M, Cullin L, Cox S, West L, Povey S, Spurr NK. Development of a panel of monochromosomal somatic cell hybrids for rapid gene mapping. Ann Hum Genet 1995; 59: 233–241. Du J, Zhu Y, Shanmugam A, Kenter AL. Analysis of immunoglobulin S gamma 3 recombination breakpoints by PCR: implications for the mechanism of isotype switching. Nucleic Acids Res 1997; 25: 3066–3073. Ozaki K, Terakita A, Ozaki M, Hara R, Hara T, Hara-Nishimura I,

32

33

34

35

36 37 38

39

Mori H, Nishimura M. Molecular characterisation and functional expression of squid retinal-binding protein: a novel species of hydrophobic ligand-binding protein. J Biol Chem 1994; 269: 3838–3845. Dalla-Favera R. Chromosomal translocations involving the c-myc oncogene in lymphoid neoplasia. In: Kirsch IR (ed). The Causes and Consequences of Chromosomal Aberrations. CRC Press: New York, 1993, pp 313–332. Cigudosa JC, Rao PH, Calasanz MJ, Odero MD, Michaeli J, Jhanwar SC, Chaganti RSK. Characterization of nonrandom chromosomal gains and losses in multiple myeloma by comparative genomic hybridization. Blood 1998; 91: 3007–3010. Stec I, Wright TJ, van Ommen GB, de Boer PAJ, van Haeringen A, Moorman AFM, Altherr MR, den Dunnen JT. WHSC1, a 90 kb SET domain-containing gene, expressed in early development and homologous to a Drosophila dysmorphy gene maps in the Wolf– Hirschhorn syndrome critical region and is fused to IgH in t(4;14) multiple myeloma. Hum Mol Genet 1998; 7: 1071–1082. Fonseca R, Witzig TE, Gertz MA, Kyle RA, Hoyer JD, Jalal SM, Greipp PR. Multiple myeloma and the translocation t(11;14)(q13;q32): a report on 13 cases. Br J Haematol 1998; 101: 296–301. Dunnick W, Hertz GZ, Scappino L, Gritzmacher C. DNA sequences at immunoglobulin switch region recombination sites. Nucleic Acids Res 1993; 21: 365–372. Lee C, Kondo S, Honjo T. Frequent but biased class switch recombination in the S mu flanking regions. Curr Biol 1998; 8: 227–230. Janz S, Kovalchuk AL, Muller JR, Potter M. Genetic recombinations between c-myc and IgH mu as precursors for recombinations between c-myc and IgH alpha in BALB/c plasmacytomas. Curr Topics Microbiol Immunol 1997; 224: 241–250. Li J, Daniels GA, Lieber MR. Asymmetric mutation around the recombination break point of immunoglobulin class switch sequences on extrachromosomal substrates. Nucleic Acids Res 1996; 24: 2104–2111.

1107