Linkage on chromosome 10 of several murine retroviral integration loci associated with leukaemia

Peter Haviernik1, Stephen M. Festinb,1, Rene Opavskyc,1, Richard P. Koller1, Nighean I. Barr2, James C. Neil2 and Linda Wolff1

Leukemogenesis Section, Laboratory of Cellular Oncology, National Cancer Institute, NIH, Bethesda, MD 20892-4255, USA1
Molecular Oncology Laboratory, University of Glasgow, Glasgow G61 1QH, UK2

Author for correspondence: Linda Wolff. Fax +1 301 594 3996. e-mail lwolff{at}helix.nih.gov


   Abstract
Top
Abstract
Introduction
Methods
Results
Discussion
References
 
Mml loci have been identified as provirus integration sites among a subset of monocytic tumours induced by murine leukaemia virus (MuLV) infection of BALB/c and DBA/2 mice. These myeloid leukaemias contain a retrovirus integrated on chromosome 10 in proximity to the c-myb locus; however, c-myb expression was not altered. Detailed physical mapping enabled placement of the retroviral integration sites ~25 kb (Mml1), ~51 kb (Mml2), and ~70 kb (Mml3) upstream of the c-myb locus. Furthermore, the Fti1 (fit-1) locus, a common integration site in feline leukaemia virus-induced T cell lymphomas, was mapped upstream of Mml3. Sequence analysis of Mml1, Mml2 and Mml3 loci (39·6, 16·4 and 5·9 kb, respectively) in conjunction with the BLAST (basic local alignment search tool) homology searches against the expressed sequence tag (EST) database and the use of gene/exon prediction programs revealed potential coding sequences that were not confirmed by Northern analysis or RT–PCR. The sequences between c-myb and Fti1, which were shown to include two potential scaffold/matrix attachment regions (S/MARs), are most likely regulatory in nature. An extended search for transcribed sequences far upstream of Mml3 revealed five genes, four of which were expressed in multiple tissues in mice. These genes could not be linked to tumour formation by the virus but their homologous sequences were found on human chromosome 6, thus allowing extension of the syntenic region on mouse chromosome 10 to approximately 250 kb.


   Introduction
Top
Abstract
Introduction
Methods
Results
Discussion
References
 
During leukaemogenesis, haematopoietic cell growth and differentiation is altered by changes in expression or function of proto-oncogenes or tumour suppressors. One approach that has been utilized to identify genes that are involved in the oncogenic process is retroviral insertional mutagenesis (Jonkers & Berns, 1996 ; Kung et al., 1991 ). We have employed this approach to detect the genetic alterations involved in monocytic leukaemia through the development of a unique in vivo model. Intravenous inoculation of murine leukaemia virus (MuLV), either Moloney MuLV or amphotropic virus 4070A, in conjunction with pristane, results in insertional mutagenesis of haematopoietic progenitor cells and the ultimate outgrowth of leukaemic cells in the peritoneal cavity (Nason-Burchenal & Wolff, 1993 ; Wolff et al., 1991 ; Wolff et al., 1988 ). These tumours, called MML (murine leukaemia virus-induced myeloid leukaemia), are transplantable and can be cultivated in vitro as well. In the absence of virus, pristane-treated mice never develop these neoplasms.

Using this model, several loci of retrovirus integration have been identified. The most frequent target is c-myb (Mukhopadhyaya & Wolff, 1992 ; Nazarov & Wolff, 1995 ; Shen-Ong & Wolff, 1987 ; Wolff et al., 1991 ). Its locus is rearranged and its expression at the RNA level is constitutive due to the activity of the retrovirus promoter and/or enhancer. However, many of the leukaemias induced by the 4070A virus express normal c-myb and have integrated proviruses outside the c-myb locus. In some of these tumours insertion sites were found to be clustered on chromosome 10, approximately 25 kb upstream of c-myb in a locus named Mml1 (Koller et al., 1996 ).

Recently, we identified two other sites of integration, Mml2 and Mml3, in MML. Here we have reported their position 50–70 kb upstream of c-myb, as well as a physical map linking Mml1, Mml2, Mml3 and Fti1, and a comparison of this region to an analogous region on human chromosome 6. Interestingly, this chromosomal region seems to be a preferred area for retroviral integrations that contribute to leukaemogenesis since Fti1, a common integration site in feline leukaemia virus (FeLV)-induced T cell lymphomas, maps to this region (Barr et al., 1999 ). Another integration site, Ahi1, in Abelson MuLV-induced lymphomas is located downstream of c-myb (Jiang et al., 1994 ). For this study, we set out to determine if, in a 250 kb region upstream of c-myb, there are any genes that might be altered in expression or function thus contributing to the neoplastic disease. We concluded that, although there are several genes more than 100 kb upstream of c-myb, none appeared to be related to the oncogenic process. The region of approximately 100 kb, surrounding Mml1, 2, 3 and Fti1, is devoid of genes and presumably contains sequences involved in chromosomal structure or regulation of gene expression at a distance.


   Methods
Top
Abstract
Introduction
Methods
Results
Discussion
References
 
{blacksquare} Cell cultures.
The murine monoblastic cell line M1 (Liebermann & Hoffman-Liebermann, 1989 ) was maintained in RPMI 1640 medium supplemented with 10% heat-inactivated horse serum (Life Technologies). The leukaemic cell lines were cultured in Dulbecco’s modified Eagle’s medium with 10% (v/v) foetal bovine serum (Life Technologies). For cell differentiation experiments, cells were seeded at a density of 1–2x105 cells/ml in medium containing IL-6. IL-6 stocks were prepared as described previously (Schmidt et al., 2000 ).

{blacksquare} Primers and probes.
The amplimers used for PCR amplification of the probes are listed in Table 1. Other probes used in this study were a c-myb genomic BglII fragment [B probe in Mukhopadhyaya & Wolff (1992) ], an Ahi1 0·8 kb genomic PstI–HindIII fragment [Ahi-D probe in Poirier et al. (1988) ] and an 4070A env (RK107) 1·5 kb EcoRI–NheI fragment.


View this table:
[in this window]
[in a new window]
 
Table 1. Amplimers used for PCR amplification of probes

 
{blacksquare} RNA isolation, Northern analysis and RT–PCR.
Total RNA was prepared using TRIzol Reagent (Gibco–BRL) according to the manufacturer’s instructions. Samples containing 10 µg total RNA were resolved on a 1·2% agarose gel containing 0·25 M formaldehyde in a MOPS buffer system. After transfer onto nylon membrane (Nytran SuPerCharge, Schleicher & Schuell) and UV cross-linking the blot was hybridized with random-priming labelled probes (Ready-To-Go DNA Labelling Beads, Amersham Pharmacia Biotech).

The RT–PCR analysis was performed using the Titan One Tube RT–PCR Kit (Roche) according to the manufacturer’s instructions.

{blacksquare} Mml2, Mml3 cloning and library screening.
The Mml2 and Mml3 integration sites were cloned from the 30-2-9 and the 30-2-7 tumours, respectively. In the case of 30-2-9, genomic DNA was purified from primary tumour tissue and for 30-2-7 from an established cell line. Genomic DNA was digested with EcoRI and size-selected on a 0·8% agarose gel. Fragments from a range of about 3–5 kb were purified and ligated to the EcoRI cut arms of the {lambda}ZAP Express vector. These {lambda} subgenomic libraries were screened with a 4070A env probe (RK107). An insert from the positive Mml2 recombinant phage ({lambda}21.1) was recloned into pBluescriptII SK(-) and sequenced. The Mml3 recombinant phage DNA ({lambda}26·2) was sequenced directly. The genomic sequences flanking the provirus were used to design amplimers for RK114 (Mml2) or MML-3 (Mml3) PCR probes.

The mouse genomic BAC library RPCI-22 (female spleen 129S6/SvEvTac, about 212000 clones, of 154 kb average insert size, and approximately 10·9-fold genomic coverage, Children’s Hospital Oakland Research Institute, Oakland, CA, USA) (Osoegawa et al., 2000 ) was screened by Southern hybridization with either an Mml1 (RK29) or an Mml2 (RK114) probe.

A genomic region containing the wild-type Mml2 locus was cloned using a BAC DNA positive clone from RPCI-22 library (#401-G9) by digestion with BamHI. Size-selected and purified fragments of 15–20 kb were ligated into the {lambda}DASHII vector. The subgenomic {lambda} library was screened using an Mml2 PCR probe (RK114) and the insert from a positive phage ({lambda}II-4a) was subcloned into a plasmid. Similarly, the Mml3 locus was cloned from the same BAC (#401-G9) after BamHI digestion and size selection. Fragments of 4–8 kb were cloned directly into the pBluescriptIISK(-) vector and the library was screened with an Mml3 PCR probe (MML-3) to obtain a positive clone (pSK-MML3).

{blacksquare} Plasmid, {lambda}, BAC and genomic DNA isolation and Southern analysis.
Plasmid DNA was prepared with the Wizard Plus SV Miniprep System (Promega), {lambda} DNA was purified from lysate using the Lambda Maxi Kit (Qiagen), BAC DNA was isolated with the Plasmid Maxi Kit (Qiagen) and genomic DNA was prepared using the Wizard Genomic DNA Purification Kit (Promega). Ten µg of DNA were digested with various restriction endonucleases, separated on a 0·8% agarose gel in a TAE system, transferred onto nylon membrane and hybridized as described for Northern analysis.

{blacksquare} Pulsed field gel electrophoresis (PFGE).
PFGE was performed on the CHEF-DR II System (Bio-Rad) in 1% agarose in 0·5x TBE buffer at 6 V/cm for 20 h with initial switch time 50 s and final switch time 90 s.

{blacksquare} DNA sequencing and bioinformatics.
DNA was sequenced using the BigDye Terminator Cycle Sequencing Kit (Applied Biosystems) and the samples were run on the ABI PRIZM 377 sequencer at the DNA Sequencing MiniCore Facility (NCI). Assembly of obtained sequence was performed by the Sequencher program (Gene Codes) and homology searches were done by the BLAST (basic local alignment search tool) program at the NCBI (Altschul et al., 1997 ). For masking the repetitive elements the RepeatMasker program was used (A. F. A. Smit and P. Green at http://ftp.genome.washington.edu/cgi-bin/RepeatMasker). Sequences were analysed by the RUMMAGE sequence annotation system (http://gen100.imb-jena.de/rummage). The search for matrix attachment regions (MAR) was carried out using the MAR-Wiz program (Singh et al., 1997 ) (www.futuresoft.org/MAR-Wiz/).


   Results
Top
Abstract
Introduction
Methods
Results
Discussion
References
 
Cloning and mapping of novel integration sites Mml2 and Mml3 on chromosome 10 between c-myb and Fti1
In 4070A retrovirus-induced promonocytic leukaemias in DBA/2 mice, c-myb is often the target of retrovirus insertional mutagenesis (Wolff et al., 1991 ). However, many tumours have retroviral integration sites outside the immediate c-myb locus. The cloning of another novel integration site, Mml1, rearranged due to proviral integration in ten tumours, was reported (Koller et al., 1996 ) and we later identified three additional tumours with alterations at this site (data not shown). Recently, proviral integration sites other than c-myb or Mml1 were identified in MML by Southern blot analysis using a virus env-specific probe. Two of these, 30-2-7 and 30-2-9, contain a single provirus integrated in their genome (Fig. 1). Genomic DNA comprising proviral and flanking sequences was cloned and probes were designed that were specific for the new cellular proviral integration sites, Mml2 and Mml3. Rearrangements at these loci have not been detected in any other tumours so far, but subsequent mapping near Mml1 and c-myb indicated that the effects of all of these integrations might be similar.



View larger version (38K):
[in this window]
[in a new window]
 
Fig. 1. Southern analysis of proviral integration in new MML. Ten µg of genomic DNA from a cell line derived from tumour 30-2-7 and from primary tumour 30-2-9, as well as from BALB/c normal liver was digested with EcoRI and separated on a 0·8% agarose gel. After blotting, the membrane was hybridized with a 4070A env probe. Arrows on the left indicate the positions of marker fragments, those on the right mark the genomic fragments containing provirus.

 
Screening of a mouse BAC genomic library, RPCI-22, using an Mml1 probe, RK29 yielded several positive clones (C6, G9, K9, D17, K14, L15). Interestingly, they were found to hybridize with probes Mml2 and Mml3 (Fig. 2). In addition, clones A7 and J8 were identified as being positive for Mml2 and Mml3, but not Mml1 or c-myb (Fig. 2). Fti1 had previously been shown to be in the same chromosomal region as Mml1, although the distance between the integration sites was not determined (Barr et al., 1999 ). A Fti1 probe hybridized to all of the BAC clones shown in Fig. 2. We therefore concluded that all of the proviral integration sites Mml1, 2, 3 and Fti1 are located in the same chromosomal region. Using endonuclease digestion with several rare cutting enzymes, PFGE, blotting, and hybridization with integration site-specific probes, we constructed a physical map of this region as shown in Fig. 2. According to our analysis, Mml1, Mml2, Mml3 and Fti1 are situated approximately 25, 51, 70 and 81–90 kb, respectively, upstream of the c-myb gene.



View larger version (10K):
[in this window]
[in a new window]
 
Fig. 2. Map of chromosome 10 segment. The physical map was built using the PFGE of G9 BAC DNA digested with the restriction endonucleases and hybridization with different probes. The BAC contig map was composed from hybridization experiments of individual BAC clones with various probes. The genetic map was compiled from two of the above analyses. Fti1, a common retrovirus integration site in FeLV-induced T cell lymphomas (Tsujimoto et al., 1993 ); Mml1, a common retrovirus integration site in 13 MuLV-induced myeloid leukaemias (Koller et al., 1996 ); Mml2, Mml3, retrovirus integration sites each found in one MuLV-induced myeloid leukaemia.

 
Effect of Mml1 integration on c-myb expression
Since Mml1 is located near the c-myb locus and c-myb is activated in a large majority of monocytic tumours (Wolff, 1996 ), it is possible that proviral insertion in Mml1 has a positive influence on c-myb expression at the mRNA level. To our surprise, in some of the tumours (30A-2-6 and 30C-13) c-myb was not expressed, while in others (30-3-14, 30-4-18, 30C-19) c-myb was expressed (Fig. 3A, B). In the instances where c-myb was expressed, we wondered whether its expression could be modulated in response to normal signalling. It is known that in myeloid cells, c-myb is expressed during proliferation and down-regulated during differentiation (Coll et al., 1983 ; Duprey & Boettiger, 1985 ; Liebermann & Hoffman-Liebermann, 1989 ; Sheiness & Gardinier, 1984 ; Wolff et al., 1996 ). If proviral insertion resulted in constitutive expression, it would be expected to be refractory to differentiation-associated down-regulation. We analysed a c-myb-positive Mml1 cell line 30-4-18 following treatment with IL-6, an inducer of differentiation. Down-regulation of c-myb occurred in these cells in response to the differentiation signal as it did in the M1 monoblastic cell line and a phenotypically similar tumour with virus integration in the 3' end of the c-myb locus (45-16). (Fig. 3B). In this case, deregulated transcription of c-myb could not be an explanation for the sustained proliferation of leukaemia 30-4-18.



View larger version (34K):
[in this window]
[in a new window]
 
Fig. 3. Expression analysis of c-myb. (A) Expression profile of c-myb in the leukaemic cells with the retrovirus integration in the Mml1 locus. The samples of total RNA from various leukaemias with proviral insertion in the Mml1 were analysed by Northern blotting and hybridization with a c-myb cDNA probe. M1 cells were used as a positive control. (B) Down-regulation of c-myb. Northern blot analysis of c-myb expression in total RNA from leukaemic cells with a proviral integration in either the c-myb (45-16) or the Mml1 (30-4-18) loci without (-) or after (+) 24 h treatment with IL-6. M1 cells were used as a control of c-myb down-regulation that can be observed even with a truncated version of this gene (45-15). The GAPDH expression served as an internal normalizing control.

 
Sequence, informatics and expression analysis
Because a correlation between integration in Mml1 and deregulated c-myb expression was not established, we decided to search for genes that could be potential oncogenes or tumour suppressor(s) in the general vicinity of Mml1, Mml2, Mml3 and Fti1. Probes obtained in close proximity to the integration sites (RK29, ML2D-D, MML-3, mu fit-1) did not detect any gene expression in normal tissues including spleen, thymus, bone marrow, heart, liver, kidney or in tumour-derived cell lines, with or without integrations in these loci (data not shown). Due to an inability to detect any transcribed sequences close to the integration sites, we decided to sequence the areas surrounding the proviral integrations and search for genes using bioinformatic sequence analysis tools. Because of the possibility that retrovirus-induced transcriptional deregulation can occur from a long distance, large regions were sequenced (Jonkers & Berns, 1996 ).

To identify potential genes in the vicinity around c-myb, Mml1, Mml2, Mml3 and Fti1, 39·6 kb in the Mml1 region was sequenced including between Mml1 and the c-myb, and sequences around the Mml2 and Mml3 (16·4 kb and 5·9 kb, respectively). To sequence the region between Mml1 and c-myb we used lambda clones subcloned from a P1 clone (#5362, Genome Systems). For sequencing around Mml2 and Mml3, plasmid clones derived from a G9 BAC (RPCI-22, #401-G9) were used. After sequencing and prior to comparative sequence analysis, repetitive elements were filtered using the RepeatMasker program (A. F. A. Smit and P. Green, unpublished data, University of Washington Genome Center, http://ftp.genome.washington.edu/cgi-bin/RepeatMasker). The extent of masking varied. It was 27·4, 28 and 38·9% for Mml1, Mml2 and Mml3, respectively (Table 2).


View this table:
[in this window]
[in a new window]
 
Table 2. Sequencing statistics

 
Using the BLAST algorithm (Altschul et al., 1997 ) to search the expressed sequence tag (EST) database, four ESTs with coding capacity sequence were identified near the integration sites. Two ESTs were found in Mml1 and two in Mml2 sequences as shown in Table 3 and Fig. 4(A, B). In a second approach, gene/exon prediction programs were used to identify potential exons. Sequences were analysed by the RUMMAGE sequence annotation system (Taudien et al., 2000 ) (http://gen100.imb-jena.de/rummage), which includes several exon prediction programs (GenScan, Grail, MZEF and Xpound). The GenScan program predicts a complex gene structure comprising several exons (Burge & Karlin, 1997 ). The other programs search for individual exons of different reliability based on their probability score (Thomas & Skolnick, 1994 ; Xu et al., 1994 ; Zhang, 1997 ). No exon was predicted by all programs, so we analysed only those identified by at least two programs. Two were located in the Mml1 region and two in Mml2 (Fig. 4A, B). Next, RNAs from embryos, several normal tissues and tumours were analysed for transcription by Northern blot analysis using PCR-derived probes covering the EST regions or by RT–PCR with amplimers derived from predicted exon sequences. However, as in the case of probes that were proximal to integration sites, we again did not detect any transcripts. In support of this finding, we did not find any CpG islands (clusters over 400 bp with at least 55% GC content) in the region, which would have suggested the presence of a gene(s). Although no transcripts were detected, this result does not exclude the possibility that these regions are associated with genes that may be expressed at early stages of tumour formation.


View this table:
[in this window]
[in a new window]
 
Table 3. EST homologues in Mml

 


View larger version (8K):
[in this window]
[in a new window]
 
Fig. 4. Features of the Mml1 (A)and Mml2 (B) locus sequences. Several features including predicted exons, EST homologues and predicted S/MARs are depicted in relationship to virus integration sites. Virus integration sites are shown by vertical lines and the orientation of provirus is indicated by an arrow.

 
Since retroviruses can activate genes from a long distance, the search was extended upstream of Mml3. Shotgun cloning (Deininger, 1983 ) and low redundancy sequencing (Bouck et al., 1998 ) of BAC clone J8 (RPCI-22, #471-J8) resulted in approximately 130 kb of sequence (Table 2). Within this sequence, regions of homology to several ESTs were found. The identified EST sequences were used to probe for expression by Northern blot analysis on mouse tissues and tumour cell lines and our characterization is summarized in Table 4. In addition, their presence is noted on the map in Fig. 2. Although probes to two of the ESTs (A2, A239) detected message in some non-haematopoietic tissues of the mouse, they did not detect message in the spleen or thymus or tumours with integrations in Mml1, 2 or 3 or c-myb. Probes to the other EST detected transcripts in all of the above types of tumours as well as in normal tissues; however, there was no alteration in size or abundance of RNA when comparing normal versus neoplastic cells. Therefore, we concluded that none of these ESTs appears to be involved in tumour development. One EST homologue, A115, also shows homology to the HBS1L/eRFS gene (Wallrapp et al., 1998 ). The protein encoded by this gene possesses GTP-binding activity and is related in primary structure to proteins involved in translational machinery. Human HBS1L spans 90 kb and is located approximately 130 kb upstream of c-myb on chromosome 6. Our genomic clone (A115) covers approximately 160 bp of mRNA sequence, which, from our comparison of the human mRNA (acc.# XM 051825) and genomic sequences (acc.# NT 027047, contig containing c-myb and HBS1L), appears to be homologous to the second exon of the human gene.


View this table:
[in this window]
[in a new window]
 
Table 4. EST homologues in J8 BAC

 
We compared our sequences using BLAST to a recently published draft of the human genome sequence (http://www.ncbi.nlm.nih.gov/genome/seq/) and identified regions homologous to Mml on human chromosome 6 upstream of human c-myb. This was not a surprise because this part of mouse chromosome 10 is considered to be syntenic with an analogous region on human chromosome 6. We were interested in the extent of homology based on collinearity of individual sequences. According to our data, the order of loci on mouse chromosome 10 is c-mybMml1–Mml2–Mml3–Fti1–A239–A182–A115/A2. We used the c-myb gene, which was localized on human chromosome 6 (acc.#. NT 025741) to position 2801411–2839158, as a reference. The order of the loci on the human chromosome was determined to be c-mybMml1–Mml2–Mml3–A115–A274 (A274 was not mapped in the mouse BAC using Southern analysis), thus expanding the syntenic region to almost 250 kb.

As indicated above, our search for exons in close proximity to Mml1, Mml2 and Mml3 did not lead to direct identification of any genes influenced by proviral integration at these loci. The lack of detectable coding regions in the vicinity might suggest that these regions harbour regulatory or structural elements. Chromatin-organizing structures, such as scaffold/matrix attachment regions (S/MARs) are found in AT-rich areas. Although the overall AT content for the 39·6 kb Mml1 region, as shown in Fig. 4A, is 56·5%, we found several areas that have a relatively high AT content, between 61 and 65%. The MAR-Wiz program (Singh et al., 1997 ) (www.futuresoft.org/MAR-Wiz/) predicted the S/MAR around the position of 18·3 kb which is next to the provirus insertion site in one Mml1 tumour – V46–34 at a position around 17·9 kb. Since S/MARs often contain specific recognition sequences defined by a MAR/SAR recognition signature (MRS), we searched for these elements. The signature is made up of two degenerate sequences, an 8 bp sequence (AATAAYAA) and a 16 bp sequence (AWWRTAANNWWGNNNC), found within 200 bp of each other (van Drunen et al., 1999 ). Although we could identify several 8 bp sequences throughout the Mml1 sequence, we could not identify the 16 bp sequences. The Mml2 sequence of 16·4 kb (Fig. 4B) with AT content of 59·2% also contains several AT-rich areas of 66–72% AT. In one of them the MAR-Wiz predicted an S/MAR around 11·3 kb, but we did not find any complete MRS. We conclude from this search for S/MARs that there is some evidence for S/MARs. Although the data are only partially predictive, they warrant further experimentation to determine directly whether or not S/MARs are present.


   Discussion
Top
Abstract
Introduction
Methods
Results
Discussion
References
 
Our analysis of a set of monocytic leukaemias established a target region for retroviral integration that spans approximately 60 kb and is located between the Fti1 and c-myb loci. Although we found several genes transcribed from a region more than 80 kb upstream of c-myb, there was no evidence within this 80 kb area for the presence of transcribed sequences. The lack of genes in a region of this size is not unusual. Data obtained from sequencing of the human genome show that exons comprise only about 1–1·5% of the genome sequence and 20% of the genome sequences are contained in islands, larger than 500 kb, without a gene. However, gene density varies between 43·5 kb per gene for the gene-rich chromosome 19 and 200 kb per gene for two gene-poor chromosomes, 13 and Y (Venter et al., 2001 ). With this in mind it would not be unusual to find an 80 kb plus segment of DNA that is devoid of genes.

Repeated integration of proviruses into this chromosome 10 region, in association with tumorigenesis, suggests an important function for these sequences, despite the absence of transcribed genes in this area, that might be altered in expression due to the provirus. The proviruses at these loci are like many others that appear to be involved in neoplastic disease but have not yet been assigned a function (Jonkers & Berns, 1996 ). The fact that the region is devoid of genes, and the MAR-Wiz program predicted two potential S/MARs, might suggest that at least part of the area contains structure-related regulatory sequences. Matrix attachment elements organize chromatin topologically, allowing spatial and temporal gene regulation and the nuclear matrix is implicated in the control of development and differentiation (Nepveu, 2001 ; Stein et al., 1995 ). Interestingly, alterations in nuclear architecture have been associated with the neoplastic condition (Deppert et al., 2000 ); there are transforming genes that have the ability to bind to nuclear matrix including SV40 Large T antigen and mutant p53, the latter of which was found to bind directly to the S/MAR (Deppert, 2000 ). In addition, S/MARs are AT rich and some proteins, including the mixed lineage leukaemia (MLL) protein, contain an AT hook domain, which binds to these sequences (Broeker et al., 1996 ; Caslini et al., 2000 ). Moreover, MLL fusion proteins have been implicated in myeloid leukaemia in man (Dimartino & Cleary, 1999 ). It is conceivable, therefore, that the virus could perturb the binding of DNA to matrix through the physical disruption or could introduce regulatory proteins such as those that bind to the LTR. Further studies will be required to determine the function of the sequences predicted to be S/MARs by the MAR-Wiz program.

One must still consider the possibility that regulatory regions associated with these integration sites could function to influence c-myb expression. Integrations have now been found in a wide region (approximately 150 kb) both upstream (Fti1, Mml3, Mml2, Mml1) of c-myb, as well as downstream (Ahi1). Until recently, downstream integrations were only found in lymphoid tumours. However, Blaydes et al. (2001) also discovered provirus insertions in this region in a series of BXH-2 myeloid leukaemias.

In the present study we considered the possibility that integrations affected Myb transcription by looking to see if, in established tumour cell lines, provirus integration in Mml1 caused increased expression of c-myb (Fig. 3A, B) and found that some cell lines did not express it. However, our inability to find increased expression in the leukaemias may not rule out the possibility that integration into this area influenced c-myb expression during an early phase of the disease prior to differentiation. Our previous data suggest that mature MMLs (with integrations in c-myb) develop from pre-leukaemic cells in the bone marrow or spleen that are immature (Nason-Burchenal & Wolff, 1993 ). If pre-leukaemic cells involving Mml1, 2, 3 are also immature cells, their mechanism of endogenous c-myb gene regulation could be very different from that of end-stage tumour cells, which are mature monocytic cells. Integrations that affect c-myb could influence early stage progression, for example by providing increased proliferation potential or an anti-apoptotic role, since these are assigned functions of c-Myb. If this is the case, one would have to predict that, at later stages of progression, these functions are provided by additional oncogenic event(s). The proximity of this proto-oncogene to the Mml integration sites continues to raise the question of whether these integration sites affect c-myb.

Overall, the data presented here indicate that the 80 kb region upstream of c-myb, which contains three separate retroviral integration sites associated with myeloid tumours, may be devoid of genes, but could possess regulatory function. Since other integration sites associated with lymphoid and myeloid tumours have also been found in this chromosomal region, it will be important to determine how sequences in this area can influence the neoplastic state.


   Acknowledgments
 
We thank Paul Jolicoeur for kindly providing the Ahi-1 probe. We are also grateful to Lawrence Kirschner for providing the PFGE apparatus for the physical mapping. We would like to thank Mark Miller for analysis of the sequencing reactions.


   Footnotes
 
b Present address: Department of Biology, Hamilton College, Clinton, NY 13323, USA.

c Present address: Division of Human Cancer Genetics, Ohio State University, Columbus, OH 43210, USA.


   References
Top
Abstract
Introduction
Methods
Results
Discussion
References
 
Altschul, S. F., Madden, T. L., Schaffer, A. A., Zhang, J., Zhang, Z., Miller, W. & Lipman, D. J. (1997). Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Research 25, 3389-3402.[Abstract/Free Full Text]

Barr, N. I., Stewart, M., Tsatsanis, C., Fulton, R., Hu, M., Tsujimoto, H. & Neil, J. C. (1999). The fit-1 common integration locus in human and mouse is closely linked to MYB. Mammalian Genome 10, 556-559.[Medline]

Blaydes, S. M., Kogan, S. C., Truong, B. T., Gilbert, D. J., Jenkins, N. A., Copeland, N. G., Largaespada, D. A. & Brannan, C. I. (2001). Retroviral integration at the Epi1 locus cooperates with Nf1 gene loss in the progression to acute myeloid leukemia. Journal of Virology 75, 9427-9434.[Abstract/Free Full Text]

Bouck, J., Miller, W., Gorrell, J. H., Muzny, D. & Gibbs, R. A. (1998). Analysis of the quality and utility of random shotgun sequencing at low redundancies. Genome Research 8, 1074-1084.[Abstract/Free Full Text]

Broeker, P. L., Harden, A., Rowley, J. D. & Zeleznik-Le, N. (1996). The mixed lineage leukemia (MLL) protein involved in 11q23 translocations contains a domain that binds cruciform DNA and scaffold attachment region (SAR) DNA. Current Topics in Microbiology and Immunology 211, 259-268.[Medline]

Burge, C. & Karlin, S. (1997). Prediction of complete gene structures in human genomic DNA. Journal of Molecular Biology 268, 78-94.[Medline]

Caslini, C., Alarcon, A. S., Hess, J. L., Tanaka, R., Murti, K. G. & Biondi, A. (2000). The amino terminus targets the mixed lineage leukemia (MLL) protein to the nucleolus, nuclear matrix and mitotic chromosomal scaffolds. Leukemia 14, 1898-1908.[Medline]

Coll, J., Saule, S., Martin, P., Raes, M. B., Lagrou, C., Graf, T., Beug, H., Simon, I. E. & Stehelin, D. (1983). The cellular oncogenes c-myc, c-myb and c-erb are transcribed in defined types of avian hematopoietic cells. Experimental Cell Research 149, 151-162.[Medline]

Deininger, P. L. (1983). Random subcloning of sonicated DNA: application to shotgun DNA sequence analysis. Analytical Biochemistry 129, 216-223.[Medline]

Deppert, W. (2000). The nuclear matrix as a target for viral and cellular oncogenes. Critical Reviews in Eukaryotic Gene Expression 10, 45-61.[Medline]

Deppert, W., Gohler, T., Koga, H. & Kim, E. (2000). Mutant p53: ‘gain of function’ through perturbation of nuclear structure and function? Journal of Cellular Biochemistry S35, 115-122.

Dimartino, J. F. & Cleary, M. L. (1999). Mll rearrangements in haematological malignancies: lessons from clinical and biological studies. British Journal of Haematology 106, 614-626.[Medline]

Duprey, S. P. & Boettiger, D. (1985). Developmental regulation of c-myb in normal myeloid progenitor cells. Proceedings of the National Academy of Sciences, USA 82, 6937-6941.[Abstract]

Jiang, X., Villeneuve, L., Turmel, C., Kozak, C. A. & Jolicoeur, P. (1994). The Myb and Ahi-1 genes are physically very closely linked on mouse chromosome 10. Mammalian Genome 5, 142-148.[Medline]

Jonkers, J. & Berns, A. (1996). Retroviral insertional mutagenesis as a strategy to identify cancer genes. Biochimica et Biophysica Acta 1287, 29-57.[Medline]

Koller, R., Krall, M., Mock, B., Bies, J., Nazarov, V. & Wolff, L. (1996). Mml1, a new common integration site in murine leukemia virus-induced promonocytic leukemias maps to mouse chromosome 10. Virology 224, 224-234.[Medline]

Kung, H.-J., Boerkoel, C. & Carter, T. H. (1991). Retroviral mutagenesis of cellular oncogenes: a review with insights into the mechanisms of insertional activation. Current Topics in Microbiology and Immunology 171, 1-25.[Medline]

Liebermann, D. A. & Hoffman-Liebermann, B. (1989). Proto-oncogene expression and dissection of the myeloid growth to differentiation developmental cascade. Oncogene 4, 583-592.[Medline]

Mukhopadhyaya, R. & Wolff, L. (1992). New sites of proviral integration associated with murine promonocytic leukemias and evidence for alternate modes of c-myb activation. Journal of Virology 66, 6035-6044.[Abstract]

Nason-Burchenal, K. & Wolff, L. (1993). Activation of c-myb is an early bone-marrow event in a murine model for acute promonocytic leukemia. Proceedings of the National Academy of Sciences, USA 90, 1619-1623.[Abstract]

Nazarov, V. & Wolff, L. (1995). Novel integration sites at the distal 3' end of the c-myb locus in retrovirus-induced promonocytic leukemias. Journal of Virology 69, 3885-3888.[Abstract]

Nepveu, A. (2001). Role of the multifunctional CDP/Cut/Cux homeodomain transcription factor in regulating differentiation, cell growth and development. Gene 270, 1-15.[Medline]

Osoegawa, K., Tateno, M., Woon, P. Y., Frengen, E., Mammoser, A. G., Catanese, J. J., Hayashizaki, Y. & de Jong, P. J. (2000). Bacterial artificial chromosome libraries for mouse sequencing and functional analysis. Genome Research 10, 116-128.[Abstract/Free Full Text]

Poirier, Y., Kozak, C. & Jolicoeur, P. (1988). Identification of a common helper provirus integration site in Abelson murine leukemia virus-induced lymphoma DNA. Journal of Virology 62, 3985-3992.[Medline]

Schmidt, M., Nazarov, V., Stevens, L., Watson, R. & Wolff, L. (2000). Regulation of the resident chromosomal copy of c-myc by c-Myb is involved in myeloid leukemogenesis. Molecular and Cellular Biology 20, 1970-1981.[Abstract/Free Full Text]

Sheiness, D. & Gardinier, M. (1984). Expression of a proto-oncogene (proto-myb) in hemopoietic tissues of mice. Molecular and Cellular Biology 4, 1206-1212.[Medline]

Shen-Ong, G. L. & Wolff, L. (1987). Moloney murine leukemia virus-induced myeloid tumors in adult BALB/c mice: requirement of c-myb activation but lack of v-abl involvement. Journal of Virology 61, 3721-3725.[Medline]

Singh, G. B., Kramer, J. A. & Krawetz, S. A. (1997). Mathematical model to predict regions of chromatin attachment to the nuclear matrix. Nucleic Acids Research 25, 1419-1425.[Abstract/Free Full Text]

Stein, G. S., van Wijnen, A. J., Stein, J., Lian, J. B. & Montecino, M. (1995). Contributions of nuclear architecture to transcriptional control. International Review of Cytology 162A, 251-278.

Taudien, S., Rump, A., Platzer, M., Drescher, B., Schattevoy, R., Gloeckner, G., Dette, M., Baumgart, C., Weber, J., Menzel, U. & Rosenthal, A. (2000). RUMMAGE – a high-throughput sequence annotation system. Trends in Genetics 16, 519-520.[Medline]

Thomas, A. & Skolnick, M. H. (1994). A probabilistic model for detecting coding regions in DNA sequences. IMA Journal of Mathematics Applied in Medicine and Biology 11, 149-160.[Medline]

Tsujimoto, H., Fulton, R., Nishigaki, K., Matsumoto, Y., Hasegawa, A., Tsujimoto, A., Cevario, S., O’Brien, S. J., Terry, A., Onions, D. & Neil, J. C. (1993). A common proviral integration region, fit-1, in T-cell tumors induced by myc-containing feline leukemia viruses. Virology 196, 845-848.[Medline]

van Drunen, C. M., Sewalt, R. G., Oosterling, R. W., Weisbeek, P. J., Smeekens, S. C. & van Driel, R. (1999). A bipartite sequence element associated with matrix/scaffold attachment regions. Nucleic Acids Research 27, 2924-2930.[Abstract/Free Full Text]

Venter, J. C., Adams, M. D., Myers, E. W., Li, P. W., Mural, R. J. and others (2001). The sequence of the human genome. Science 291, 1304–1351.[Abstract/Free Full Text]

Wallrapp, C., Verrier, S. B., Zhouravleva, G., Philippe, H., Philippe, M., Gress, T. M. & Jean-Jean, O. (1998). The product of the mammalian orthologue of the Saccharomyces cerevisiae HBS1 gene is phylogenetically related to eukaryotic release factor 3 (eRF3) but does not carry eRF3-like activity. FEBS Letters 440, 387-392.[Medline]

Wolff, L. (1996). Myb-induced transformation. Critical Reviews in Oncogenesis 7, 245-260.[Medline]

Wolff, L., Mushinski, J. F., Shen-Ong, G. L. & Morse III, H. C. (1988). A chronic inflammatory response. Its role in supporting the development of c-myb and c-myc related promonocytic and monocytic tumors in BALB/c mice. Journal of Immunology 141, 681-689.[Abstract/Free Full Text]

Wolff, L., Koller, R. & Davidson, W. (1991). Acute myeloid leukemia induction by amphotropic murine retrovirus (4070A): clonal integrations involve c-myb in some but not all leukemias. Journal of Virology 65, 3607-3616.[Medline]

Wolff, L., Koller, R., Bies, J., Nazarov, V., Hoffman, B., Amanullah, A., Krall, M. & Mock, B. (1996). Retroviral insertional mutagenesis in murine promonocytic leukemias: c-myb and Mml1. Current Topics in Microbiology and Immunology 211, 191-199.[Medline]

Xu, Y., Mural, R., Shah, M. & Uberbacher, E. (1994). Recognizing exons in genomic sequence using GRAIL II. Genetic Engineering 16, 241-253.[Medline]

Zhang, M. Q. (1997). Identification of protein coding regions in the human genome by quadratic discriminant analysis. Proceedings of the National Academy of Sciences, USA 94, 565-568.[Abstract/Free Full Text]

Received 15 October 2001; accepted 7 December 2001.