Previous Article | Next Article ![]()
Journal of Virology, February 2006, p. 1922-1938, Vol. 80, No. 4
0022-538X/06/$08.00+0 doi:10.1128/JVI.80.4.1922-1938.2006
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Maud Szynal,1,
Yvette Cleuter,1
Makram Merimi,1
Hugues Duvillier,1
Françoise Lallemand,1
Claude Bagnis,3
Philip Griebel,4
Christos Sotiriou,1
Arsène Burny,1
Philippe Martiat,1 and
Anne Van den Broeke1*
Laboratory of Experimental Hematology, Bordet Institute, 1000 Brussels, Belgium,1 Institute of Pathological Physiology, Charles University, Prague, Czech Republic,2 Etablissement Français du Sang, 13009 Marseille, France,3 Vaccine and Infectious Disease Organization, S7H5E3 Saskatoon, Canada4
Received 22 August 2005/ Accepted 18 November 2005
|
|
|---|
|
|
|---|
The sheep has been of particular interest as a large-animal model for studying aspects of immunology and offers a number of experimental opportunities that are not available in murine systems. Furthermore, sheep develop B-cell leukemia following experimental transmission of bovine leukemia virus (BLV), a complex retrovirus that is structurally and functionally related to human T-cell leukemia virus type 1 (HTLV-1) (6, 42, 78). This virus-associated leukemia model in sheep has been extensively studied as an in vivo approach to understand the molecular basis of human leukemogenic processes.
Transactivating oncoretroviruses (BLV and HTLV-1) induce tumors after long latency by using poorly understood mechanisms that involve Tax (6, 12, 45, 75). Both TaxHTLV-1 and TaxBLV were originally identified as transcriptional transactivators of viral expression through the indirect binding to Tax-responsive elements located in the viral long terminal repeat (4, 10, 76). Both proteins mediate the transformation of rat embryonic fibroblasts in cooperation with Ha-ras, and injection of these cells induces tumors in nude mice (57, 77). In addition, TaxHTLV-1 immortalizes human T lymphocytes (17), while TaxBLV was recently shown to disrupt the homeostatic control of ovine B cells (65). Altogether, Tax proteins are thought to be the major viral contributors to the leukemogenic processes leading to either T-cell leukemia (HTLV-1) or B-cell malignancy (BLV). Because Tax expression is not required to maintain the transformed phenotype (69, 71), it is believed to act at early stages in the multistep process leading to full malignancy. TaxHTLV-1 has been extensively studied and was found to affect a variety of cellular functions in T cells essentially through its ability to transcriptionally regulate cellular gene expression and functionally inactivate proteins involved in cell cycle progression and DNA repair (14, 36). Using a medium-scale microarray, genes related to apoptosis, the cell cycle, and DNA repair; signaling factors; immune modulators; cytokines; growth factors; and adhesion molecules were identified as TaxHTLV-1 responsive (50).
We have recently provided evidence to support the oncogenic potential of TaxBLV in ovine B cells, and we suggest an important role for NF-
B-dependent pathways in TaxBLV-associated B-cell leukemia (19, 65, 70). So far, however, little is known about the interactions of TaxBLV with cellular proteins, and there is no evidence for TaxBLV-mediated transcriptional changes associated with abnormal B-cell growth. Although there is compelling evidence that TaxBLV is an essential contributor in the initial steps of oncogenesis, it is not sufficient for transformation, and cumulative changes are necessary for the development of overt leukemia. In this study, we examined gene expression changes in response to TaxBLV in an attempt to identify genes that play a role in the molecular events leading to leukemia. We have taken a cross-species ovine-to-human gene-profiling approach using human cDNA microarrays to identify TaxBLV-associated transcriptional changes in two sets of ovine B-cell populations. The rationale behind this study based on cross-species hybridization was not to generate an exhaustive list of differentially expressed genes but rather to provide initial clues for the functional roles of genes that participate in B-cell transformation.
|
|
|---|
-irradiated murine CD154 (mCD154) L cells (0.4 x 106 cells/well) (a kind gift from Troy Randall) in AIM-V medium supplemented with 10% fetal bovine serum, 1 mM sodium pyruvate, 2 mM glutamine, nonessential amino acids, kanamycin (100 µg/ml) (Invitrogen), recombinant human interleukin-2 (IL-2) (20 U/ml), IL-4 (40 U/ml), IL-7 (20 U/ml), and IL-15 (20 U/ml) (Peprotech) (19). Clone2 and Clone2LTAXSN B-cell clones were established from ovine jejunal Peyer's patch (18) and maintained in similar conditions with 2% fetal bovine serum and IL-4 (20 U/ml). All cultures were transferred every 3 to 4 days to fresh medium, cytokines, and mCD154 L cells. The supernatant from PG13-derived retrovirus producer cell lines was used to transduce the SBL (sheep B-lymphocyte) cells as previously described (69). All cultures were incubated at 37°C in a 5% CO2 humidified atmosphere. Viable cell numbers were measured by trypan blue exclusion. Apoptotic cell death was assessed by a terminal deoxynucleotidyltransferase-mediated dUTP-biotin nick end-labeling assay according to the manufacturer's instructions (Roche), with the following modifications: cells were fixed in 2% paraformaldehyde and stored at 20°C in 70% ethanol before labeling (65). RNA extraction and microarray probe preparation. RNA was extracted using TriPure (Roche) according to the manufacturer's instructions. Quality was assessed using the Agilent 2100 Bioanalyzer (Agilent Technology). Total RNA was linearly amplified according to a method described previously by Eberwine (11), with minor modifications. Briefly, RNA was reverse transcribed using a 63-nucleotide synthetic primer containing the T7 RNA polymerase binding site 5'-GGCCAGTGAATTGTAATACGACTCACTATAGGGAGGCGG(T)24-3'. Second-strand cDNA synthesis was obtained with RNase H, Escherichia coli DNA polymerase I, and E. coli DNA ligase (Invitrogen). cDNA was blunt ended with T4 DNA polymerase (Invitrogen). Double-stranded cDNA was transcribed using T7 polymerase (Ambion), yielding amplified antisense RNA which was purified using RNeasy Mini-Columns (QIAGEN). Total RNA from the Universal Human Reference (Stratagene) was amplified and used as a reference for microarray analysis. The cDNA microarray chips consisted of 10,336 total features and were manufactured at the Bordet Institute microarray facility using IMAGE clones. Three micrograms of RNA was reverse transcribed and labeled using cyanine 5-conjugated- or cyanine 3-conjugated dUTP. Hybridization was performed in 5x SSC (1x SSC is 0.15 M NaCl plus 0.015 M sodium citrate) and 25% formamide for 14 to 16 h at 42°C. Slides were washed, dried, and scanned using an Axon 4000a laser scanner. A detailed protocol for RNA amplification, cDNA probe labeling, and hybridization is available at http://nciarray.nci.nih.gov/reference/index.shtml. Genepix software (Axon Instruments) was used to analyze the raw data, which were then uploaded to a relational database maintained by the Bordet Institute.
Data analysis. Images of scanned slides were inspected for artifacts, and aberrant spots and slide regions were flagged for exclusion from analysis. Log ratios for each spot were calculated as follows. In each channel, the signal was calculated as foreground median minus background median. If the signal was less than 300 in any single channel, the signal value in that channel was set to 300. If the signal was less than 300 in both channels, the spot was flagged as unreliable and not used in any further analysis. For all nonflagged spots, a log ratio was calculated as log2 [(red signal)/(green signal)]. The log ratios were then normalized within each array by subtracting from each array the median log ratio value across the spots on the array. The criteria for identifying genes with different levels of expression among the samples were based on the variance of their normalized log ratios across all experiments. The variance of the log ratio of each gene was compared with the median of all the variance, where the median is taken over all genes. A chi-squared statistic was computed in which the theoretical variance of the log ratios under true repetition was estimated by the median of the computed variances, and all genes with significance levels at a P value of >0.01 were filtered out.
Northern blot analysis. Ten micrograms of total RNA was incubated for 10 min in a solution containing 6 M deionized glyoxal, 50% dimethyl sulfoxide, and 0.1 M phosphate buffer (pH 7.0) on ice followed by 3 min at 55°C. RNA was separated by electrophoresis through a 1% agarose gel containing 10 mM phosphate buffer, transferred onto nylon membranes (Amersham), and UV cross-linked. Hybridization was performed in a solution containing 50% deionized formamide and 5x SSC containing a 32P-labeled TAXBLV or glyceraldehyde-3-phosphate dehydrogenase (G3PDH) probe as previously described (69). Human cDNA probes (LENG5, FOSL2, CART1, methionine aminopeptidase 2 [METAP2], and butyrate response factor 2 [BRF2]) were obtained by T3/T7 PCR amplification of IMAGE clones (Invitrogen) selected according to GenBank accession number from the Invitrogen dbEST clone database. Signals were quantified by densitometric analysis using the Bio-Profil Bio-1D program, Windows Application V99.03.
Quantitative RT-PCR. cDNA was synthesized by reverse transcription (RT) of total RNA from ovine B cells for 45 min at 42°C using Moloney murine leukemia virus reverse transcriptase (Promega). Real-time quantitative PCR was performed using SYBR green dye (SYBR green PCR master mix; Applied Biosystems) and either ovine (Mcl-1) or bovine (PTK2 and c-FOS) sequence-specific primers (Mcl-1 forward, ACGGCTTTCCAAGGCATG; Mcl-1 reverse, CCATCACTCGAGACAAAGACTTGA; PTK2 forward, CCAAATGGAGCCAGTGAACCT; PTK2 reverse, AAGCACGTGGCCTGCTATG; c-FOS forward, TCAATGACCCTGAGCCCAAG; c-FOS reverse, TCAGCCTTCAGCTCCATGCT). The amplification was performed in an ABI prism 7900 sequence detector system (Applied Biosystems) using 40 cycles of a two-step PCR (15 s at 95°C and 60 s at 60°C) after an initial activation step (95°C for 10 min). Melting curves from 60°C to 99°C were assessed to evaluate specificity. Serial dilutions of purified amplicons served to generate standard melting curves. Amplification of ovine ß-actin mRNA as an endogenous control was used to standardize the amount of sample added to the reactions (ß-actin forward, CCATCCTGCGTCTGGACCT; ß-actin reverse, CCACGTTCCGTGAGGATCTT).
|
|
|---|
chain (
c)-common cytokines, provide an excellent in vitro model for investigating cellular processes resulting from TaxBLV expression. It was deemed appropriate to first use a homogeneous B-cell population for designing optimal experimental conditions and evaluating the feasibility of the ovine-to-human microarray approach. We previously demonstrated the clear supportive role of TaxBLV on the growth of Clone2, a surface immunoglobulin M (IgM) (sIgM)-positive (sIgM+) jejunal Peyer's patch-derived B-cell clone dependent on mCD154 and IL-4 (18, 65). Clone2LTAXSN was derived following retroviral vector-mediated delivery of TaxBLV into Clone2. Both Clone2LTAXSN, which contains detectable TaxBLV protein, and the control parental Clone2 were further used for gene expression profiling. In the present study, we hybridized cDNA from ovine B-cell cultures to human cDNA-spotted high-density microarrays comprising 10,336 genes/expressed sequence tags (ESTs). Total RNA was extracted from Clone2LTAXSN and Clone2 cells at three time points posttransduction (days 125, 153, and 182) and treated as described in Materials and Methods. Amplified RNA was reverse transcribed and labeled with either cyanine 5-conjugated (Clone2LTAXSN) or cyanine 3-conjugated (Clone2) dUTP. Our experimental design thus included three replicate hybridizations with three independent RNA extractions referred to as hybridizations A, B, and C. The data were first processed as described in detail in Materials and Methods. Only the features that hybridized in at least two out of the three microarray experiments were processed for further analysis.
Ovine cDNA from the B-cell clones hybridized to 6,484 (62.54%) of the human gene targets on a microarray containing 10,336 targets (Table 1). We considered genes that were up- or down-regulated by more than 1.5-fold to be differentially expressed. This cutoff was defined relatively arbitrarily, but as described below, it was independently validated using alternative techniques for measuring relative gene expression. Finally, across all experiments, we identified 263 genes that fulfilled the criterion of differential expression. A total of 211 genes (2.05%) were significantly up-regulated (median, >1.5), and 52 genes (0.50%) were down-regulated (median, <0.66) in Clone2LTAXSN compared to Clone2 cells, respectively. In the gene list presented here, we restricted the data set to a selection of 103 up-regulated (Table 2) and 23 down-regulated (Table 3) genes, according to existing ideas about how altered gene expression may have a functional contribution to the complex processes involved in leukemogenesis. However, to keep the benefits of an unbiased approach, it is important to provide access to the complete lists of altered genes. As our knowledge progresses, we or others might be able to make more sense of the comprehensive expression patterns, ultimately aiding in the understanding of networks that contribute to cancer. The complete microarray data set is available at http://www.bordet.be/servmed/martiat/avdb/data.htm. In Tables 2 and 3, genes are arranged according to the arithmetic means of the relative expression ratios of three independent experiments, experiments A, B, and C. We identified expression changes in genes from distinct groups, including proto-oncogenes; genes regulating apoptosis, DNA transcription, and repair; transcription factors; cell cycle regulators; kinases; phosphatases; and small Rho GTPases/GTPase-binding proteins. It is interesting that genes associated with various human malignancies, especially B-cell tumors, were largely represented among the differentially expressed genes. Other genes were novel or unexpected or had no established role in tumor development. We refer the reader to Discussion (see below) for further details regarding a selection of genes from these lists.
|
View this table: [in a new window] |
TABLE 1. Basic evaluation of hybridization and altered gene expression in Tax-positive ovine B-cells using a microarray comprising 10,336 human genes/ESTsa
|
|
View this table: [in a new window] |
TABLE 2. Selected list of genes up-regulated in Clone2LTAXSN compared to Clone2a
|
|
View this table: [in a new window] |
TABLE 3. Selected list of genes down-regulated in Clone2LTAXSN compared to Clone2a
|
![]() View larger version (50K): [in a new window] |
FIG. 1. Northern blot confirmation of microarray gene expression patterns. Northern blot analysis of five putative differentially expressed genes confirms their transcriptional up-regulation in TaxBLV-expressing B cells (+) compared to control B cells (). A and B correspond to RNA extracted from B-cell clones (Clone2LTAXSN/Clone2 cells) used in two independent microarray hybridizations, hybridizations A and B (Table 2), respectively; K stands for RNA extracted from PBMC-derived B cells (SBLTCE/SBLiE) used in microarray hybridization K (Table 4). RNA extracted from ovine B cells was analyzed using the 32P-labeled human cDNA probes BRF2, LENG5, METAP2, FOSL2, CART1, and TAXBLV, described in Materials and Methods. G3PDH was used as a control. The numbers below each pair of signals indicate the relative expression ratios quantitated by densitometry analysis. The numbers in parentheses show the ratios from the corresponding microarray data.
|
|
View this table: [in a new window] |
TABLE 4. qRT-PCR analysis of mRNA levels in ovine B cellsa
|
c-common cytokines IL-2, IL-4, IL-7, and IL-15 as described in detail in Materials and Methods. After a 2-week costimulation period, the cultures consisted exclusively of sIgM+ B cells, as indicated by flow cytometry phenotypic analysis (data not shown), and were referred to as SBL (sheep B lymphocytes). The Clone2 and polyclonal SBL populations probably share many biological similarities in terms of B-cell development stage. The continued expression of sIgM in both populations suggests that their development may be arrested prior to the memory B-cell stage. The SBL, however, represent a more heterogenous sIgM+ B-cell population that has retained its dependence on both mCD154 (T-cell-dependent development) and the
c-common cytokines for continued growth. Thus, the major difference is that SBL cells require a broader range of cytokines, which may reflect the diversity of this B-cell population.
For the transfer of tax cDNA into SBL cells, we designed a vector derived from pSFß (3), referred to as pSFTaxCMVGFP. In this vector, the BLV tax cDNA is under the control of a hybrid spleen focus-forming virus/murine embryonic stem cell virus promoter demonstrated to be active in a wide range of hematopoietic cells at different stages of differentiation. The enhanced version of the green fluorescent protein (GFP) is under the control of the cytomegalovirus promoter. The pSFiECMVGFP control vector lacks tax. In contrast to the pLTaxSN/pNUNL vectors used for the transduction of Clone2 (65), this type of vector offers the possibility of selecting GFP-positive (GFP+) cells by flow cytometry and monitoring GFP+ cells over time in mixed cultures consisting of both transduced and untransduced cells. The vectors were first evaluated in Clone2 cells, given that TaxBLV expression is known to enhance their viability. Transfer of pSFTaxCMVGFP, but not pSFiECMVGFP, confers a growth advantage to Clone2 cells, confirming that TaxBLV and not the vector itself is involved in the altered phenotype (Fig. 2). The BLV tax cDNA was then delivered into SBL cells using supernatant-mediated retroviral vector transduction, and these cultures were further maintained for 7 days in the presence of mCD154 and IL-2, IL-4, IL-7, and IL-15. The pSFTaxCMVGFP- and pSFiECMVGFP-transduced B cells, referred to as SBLTCE and SBLiE, respectively, were then examined for GFP expression by flow cytometry. The GFP+ cells, representing approximately 5% of the total cell population, were sorted (Fig. 3a), and TaxBLV expression was evaluated by Northern analysis of RNA extracted from the purified B-cell populations (Fig. 3b). We then examined the growth characteristics of these cells and found that the levels of apoptotic cell death were decreased in SBLTCE compared to SBLiE under normal costimulated culture conditions (Fig. 4a) as well as in the absence of
c-common cytokines both with and without mCD154, an essential factor required to support the continuous growth of B cells (Fig. 4b). In addition, we found that SBLTCE cells displayed increased viable cell numbers (Fig. 4c), providing clear evidence that TaxBLV alters the growth of peripheral blood-derived B cells.
![]() View larger version (33K): [in a new window] |
FIG. 2. TaxBLV provides a growth advantage to ovine B-cell clones. Clone2 cells were transduced with either pSFTaxCMVGFP or control pSFiECMVGFP retroviral vectors and cultured in the presence of mCD154 and c-common cytokines. The resulting B-cell populations comprising transduced (+) and untransduced () cells, referred to as Clone2iE+/ (a) and Clone2TCE+/ (b), were examined over time for GFP fluorescence using flow cytometry. Percentages of GFP+ cells are shown. Data presented are from one representative example of three independent experiments.
|
![]() View larger version (19K): [in a new window] |
FIG. 3. Flow cytometry sorting of retroviral vector-mediated transduced SBL cells results in highly purified B-cell populations that express TaxBLV. (a) SBL were derived from ovine PBMC cultured in the presence of mCD154 and c-common cytokines for 14 days. The resulting sIgM+ B cells were transduced with either the pSFTaxCMVGFP retroviral vector encoding both Tax and GFP or the control pSFiECMVGFP retroviral vector for the transfer of GFP alone and further cultured in the presence of mCD154 and c-common cytokines for 7 days before sorting of the GFP+ cell populations. The resulting GFP+ B-cell populations are referred to as SBLTCE and SBLiE, respectively. (b) Both the SBLTCE and SBLiE populations were examined for Tax expression by Northern analysis using a TAX probe. G3PDH was used as a control of RNA load.
|
![]() View larger version (11K): [in a new window] |
FIG. 4. TaxBLV protects PBMC-derived ovine SBL cells from apoptotic cell death and enhances viable B-cell numbers. SBL cells were cultured in the presence of mCD154 and c-common cytokines IL-2, IL-4, IL-7, and IL-15. Apoptotic cell death in both SBLTCE and SBLiE cells was examined by a terminal deoxynucleotidyltransferase-mediated dUTP-biotin nick end-labeling assay of normal costimulated cultures collected 72 and 96 h after cell passage (a) or after a 24-h culture period in the absence of c-common cytokines (cyt) in either the presence or the absence of mCD154 (+CD154/CD154) (b). Data presented in a and b are means ± standard deviations of values of three independent cultures. Viable SBL cell numbers were determined by trypan blue exclusion at each passage throughout a 40-day culture period (c). Data presented in c are from one representative assay of three independent experiments.
|
|
View this table: [in a new window] |
TABLE 5. Selected list of genes up-regulated in SBLTCE compared to SBLiE
|
|
View this table: [in a new window] |
TABLE 6. Selected list of genes down regulated in SBLTCE compared to SBLiEa
|
|
View this table: [in a new window] |
TABLE 7. Genes up-regulated in both Clone2LTAXSN and SBLTCE compared to the controlsa
|
|
View this table: [in a new window] |
TABLE 8. Genes down-regulated in both Clone2LTAXSN and SBLTCE compared to the controlsa
|
![]() View larger version (29K): [in a new window] |
FIG. 5. Distribution of TaxBLV-associated differentially expressed genes between two distinct ovine B-cell populations, Clone2LTAXSN and SBLTCE. Eighty genes differentially expressed in both B-cell clones and PBMC-derived B cells were visualized as the intersect of two rectangles, each representing the overall number of significantly altered genes (median, >1.5 or <0.66) as a result of TaxBLV expression in three microarray experiments for each B-cell population. Taking relative expression ratio median values, 211 genes were >1.5 and 52 genes were <0.66 (263) in hybridizations A, B, and C (Clone2LTAXSN/Clone2), and 292 genes were >1.5 and 149 genes were <0.66 (441) in hybridizations K, L, and M (SBLTCE/SBLiE). Sixty-five genes were >1.5 and 15 genes were <0.66 (80) in both A, B, and C and K, L, and M hybridizations (overlap).
|
![]() View larger version (25K): [in a new window] |
FIG. 6. Similar changes in overall gene transcription patterns as a result of TaxBLV expression in two different ovine B-cell populations. (a) Comparison of overall gene expression profiles generated from Clone2LTAXSN/Clone2 and SBLTCE/SBLiE populations. Log ratio medians were computed from three corresponding microarray experiments, hybridizations A, B, and C (Clone2LTAXSN/Clone2 cells) and K, L, and M (SBLTCE/SBLiE cells), respectively. b is similar to a except that median values were replaced by means. Dots above the upper oblique lines represent genes overexpressed by >2-fold, and dots below the lower oblique lines represent genes underexpressed by >2-fold.
|
|
|
|---|
Using these microarrays, we evaluated gene expression profiles of previously described B-cell clones (Clone2LTAXSN/Clone2) as well as PBMC-derived SBL cells. Both Clone2 and SBL are characterized by the surface expression of IgM, a consistent finding in leukemic B cells isolated from BLV-infected sheep. Those B-cell populations probably represent a stage in T-cell-dependent B-cell development that may occur within germinal centers, but the continued expression of sIgM would indicate that they are not centroblasts and are probably not memory B cells, since there is no evidence for isotype switching. The Clone2 and polyclonal SBL populations thus share many biological similarities in terms of B-cell development stage, but Clone2 developed a relative cytokine independence over time, probably reflecting the selection of one specific B cell during the cloning process. However, these B cells remain dependent upon both mCD154 stimulation and the presence of IL-4 to proliferate and remain viable, since removal of these stimuli results in apoptosis. Thus, these are not transformed B cells. The SBL represent a more heterogenous sIgM-positive B-cell population that has retained its dependence on both CD154 (T-cell-dependent development) and the spectrum of
c-common cytokines (IL-2, IL-4, IL-7, and IL-15) for continued growth. Thus, the major difference is that SBL cells require a broader range of cytokines, which may reflect the diversity of this B-cell population and may be associated with a decreased risk of identifying changes attributable to an individual cell type. Data are presented from SBL cells isolated from the blood of only one sheep because we observed very little variation between microarray data sets for separate blood B-cell cultures from different sheep (data not shown). This suggests that animal-to-animal variation has little impact on baseline gene expression patterns in B-cell populations generated by costimulation with CD154 and cytokines. These stimuli appear to impose a strong selection on the surviving B cells. Ultimately, to avoid the limitations of the present system for exploring genes involved in oncogenesis, it will be important to analyze Tax-induced gene expression changes in B cells directly ex vivo, but this will require vectors that are capable of delivering a transgene into resting cells. The biological relevance of these two B-cell populations for studying human B-cell malignancies is that they both represent T-cell-activated B-cell populations and may be most comparable to human follicular B cells. Thus, these sheep B-cell lines may be most relevant for understanding human follicular lymphomas.
We have previously shown that TaxBLV has a clear supportive role in the growth of these B cells, with an impact on B-cell proliferation, survival, and cell cycle phase distribution (65). In both B-cell populations, we identified a significant number of genes that underwent changes in response to TaxBLV (263 genes in Clone2LTAXSN and 441 genes in SBLTCE), with an intersect of approximately 25%. Although the samples for comparison were as closely matched as possible, analysis of gene expression data generated long gene lists, suggesting that TaxBLV expression in cultured B cells leads to direct or indirect changes in mRNA levels of hundreds of genes. Faced with this mass of data, the temptation was to search for genes that conform to existing ideas about how B cells acquire abnormal phenotypes. A selection of those genes is displayed in Tables 2, 3, and 5 to 8, but only a limited number are discussed here. The benefits of an unbiased approach, however, is lost if exploration is limited to our current framework of understanding. Therefore, in this study, although we do not provide clues as to which of these changes may be important, we provide the complete data set of differentially expressed genes. It is thus important to underline that we do not exclude a critical functional role for some of the genes that are not displayed in the shortened lists.
The first striking observation was the strong agreement between our results and those for genes previously reported as being TaxHTLV-1 responsive or TaxHTLV-1 interacting in HTLV-1-infected T-cell cultures and human T-cell lines that express TaxHTLV-1 (ATF4, FOS, FOSL2, TNFSF4, and BRF2), confirming the specificity of our approach. We identified two members of the bZip transcription factors (FOS and FOSL2), ATF4, that forms dimers with bZip proteins including Fos (22), and TNFSF4 or TaxHTLV1 transcriptionally activated glycoprotein 1. Furthermore, given that TaxHTLV1 is capable of modulating the expression of BRF1 through a cis-acting replication element (38), it is tempting to speculate that BRF2, scored as up-regulated in ovine B cells, may have a similar function.
There was a good agreement between the microarray results and both the Northern blot and qRT-PCR experiments, two methods designed to validate our data. Human cDNAs were used for probing the Northern blots, while ovine and bovine sequence-specific primers were utilized in qRT-PCR assays. In many cases, the exact severalfold changes were higher than the microarray ratios. This is commonly observed, and it is now widely accepted that those methods are a better reflection of the gene expression levels in each sample. Northern blot and qRT-PCR data thus provide a rationale for the 1.5-fold cutoff used in the microarray analysis.
Interestingly, we identified several genes differentially expressed in both the B-cell clone and the PBMC-derived B cells, emphasizing their contribution to the TaxBLV-associated abnormal B-cell growth. The observation, however, that a significant proportion of altered genes do not overlap between Clone2LTAXSN and SBLTCE may appear surprising. One possible explanation is that differences in expression patterns found in B-cell clones are partly linked to changes that occur downstream from the early TaxBLV-associated events. Furthermore, given the greater diversity of SBL cells, the finding that genes are differentially expressed in SBL cells while they are not in B-cell clones does not exclude a functional role of these genes in leukemogenic processes. We thus considered that the genes differentially expressed in both B-cell populations as well as those identified in a single B-cell population were of interest.
In the vast majority of genes, increases in expression levels did not exceed two- to threefold, with the exception of CART1 (ninefold), fibrinogen-like protein 2 (FGL-2) (fivefold), and fibrillin 1 (FBN-1) (fivefold). CART1 encodes a paired-like DNA-binding homeoprotein with transactivating transcriptional activity, but its target gene remains unknown (13). p300/CREB binding protein (CBP) acts as a coactivator to CART1 through direct interaction and lysine acetylation, enhancing its transactivation capacity (29). Given that genes encoding proteins with HAT activity, including p300/CBP, were up-regulated in our study (discussed below), it is tempting to speculate that CART1 overexpression is a secondary event and not a direct consequence of TaxBLV expression. Both FGL-2, a prothrombinase, and FBN-1, involved in the structure of connective tissue, were found to be overexpressed by fivefold. To our knowledge, aberrant expression of neither of these genes has been correlated with pathways known to be associated with deregulated cell growth. Keeping in mind that the interpretation of data is complicated by changes that may occur either downstream from or independently of the stimulus under consideration, our results do not allow primary and secondary changes in gene expression to be distinguished, and further investigation is necessary to draw conclusions from these data.
In examining the expression profiles, several alterations were found in genes involved in DNA transcription/repair and genes encoding or interacting with proteins possessing HAT activity (p300, p400, GPS2, EBF, CART1, and ATF4). Gps2 modulates bovine papillomavirus type 1 E2 expression through recruiting p300 and directly interacts with TaxHTLV1 (54), known to target p300/CBP, leading to transcriptional deregulation of several cellular genes (16). EBF activates transcription factors through p300/CBP (83), and PA2G4, a transcriptional repressor, acts via recruiting histone deacetylases (82). CART1, increased by ninefold and discussed above, recruits p300/CBP to a CART-specific subnuclear compartment, thereby modulating transcription from different promoters. Furthermore, we identified transcriptional alterations of members of the chromodomain helicase DNA-binding protein (CHD) and SWI/SNF chromatin-remodeling complex gene families (CHD2, CHD4, CHD6, SWI/SWF-related regulator of chromatin 2 [SMARCA2], and SMARCA5), genes encoding DNA-binding proteins (TCF4, SREBF2, TRIP-BR2, PC4, and BHLHB2), and genes involved in DNA damage/repair (DDB2, TOP1, APEX1, and RAD51). Suppression of DNA repair likely increases the incidence of genomic mutations and may contribute to B-cell transformation. A striking example supporting the impact of DNA damage/repair imbalance in TaxBLV-mediated leukemogenesis is the up-regulation of topoisomerase 1 (TOP1) in both B-cell populations. TOP1 is involved in multiple DNA transactions, including DNA relaxation (56), and is the fusion partner of the nucleoporin gene NUP98 in acute myeloid leukemia-associated rearrangements. NUP98-TOP1 blocks differentiation of hematopoietic precursors and induces a lethal leukemia related to the DNA-binding properties and the capacity to induce DNA lesions mediated by excess TOP1. TOP1 is the target of novel anticancer drugs, the camptothecins (21).
Among the cohort of genes with altered expression, we identified genes associated with NF-
B activation and p53 suppression, including MDM2 (20), SSRP1 (81), Trip-Br2 (28), and many others. Interestingly, PPM1D overexpression abrogates the p53 suppressor activity through the inactivation of p38-mitogen-activated protein kinase (MAPK), a product of MAPK14, down-regulated in our study (5). Translocated in liposarcoma (FUS/TLS) interacts with NF-
B-p65, and its fusion to different transcription factors is associated with various malignancies, including acute myeloid leukemia (44, 67). Our observations thus support previous findings that suggested that NF-
B activation and p53 inhibition contribute to the oncogenic effects of TaxHTLV1/TaxBLV (14, 65, 80). We also observed altered expression of numerous small Rho GTPases/GTPase-binding proteins. These small G proteins contribute to various regulatory processes and play a critical role in mechanisms associated with cancer, including cell survival, invasion, and angiogenesis (23, 39, 60). Members of this family were shown to interact with TaxHTLV1 (79). We identified Ras homolog gene family members Q, B, and H (ARHQ, ARHB, and ARHH), members 1 and 4 of the ARF family of GTP-binding proteins (ARL1 and ARL4), members of Ras oncogene family (RAB6A and RAB5B), GTP-binding protein cell division cycle 42 (CDC42), and other GTPase-binding proteins (ARL6IP, SEC61B, RACGAP1, BCAR3, RAB5EP, DAAM2, C8FW, RASAL1, RGR, DOCK2, and CHN1). Our findings thus suggest that TaxBLV may exert its largely pleiotropic effects partly through the deregulation of a broad spectrum of small Rho GTPase/cytoskeleton-related genes. Furthermore, a large number of well-known oncogenes and diverse genes known to be associated with human malignancies were deregulated, including EMP3, PDE3B, LRMP, LASP1, L7a, PTK2, MEIS1, BCAS2, LMO4, OS4, and VEGF. Although the majority were found in hematological disorders, especially B-cell malignancies (PDE3B [35, 37, 48], LRMP [1], LASP1 [7, 63], MEIS1 [15, 43, 55], and VEGF [61, 74]), others were identified in hepatocarcinoma, breast cancer, and colorectal and prostate cancers (EMP3 [41], BCAS2 [40, 72], L7a [34, 68, 73, 85], OS4 [64], and PTK2 [51]).
Interestingly, transcription patterns of antiapoptotic genes were changed. MCL1 is overexpressed in various B-cell malignancies and is found to be associated with drug resistance in B-lineage acute lymphoblastic leukemia, leading to poor response to treatment (27, 33, 53, 58). MCL1-transgenic mice develop chronic clonal B-cell lymphoma closely matching BLV-induced leukemia (84). High levels of METAP2 were detected in B-cell lymphoma subtypes derived from germinal center B cells, and METAP2 expression correlated with that of BCL6, a contributor to B-cell lymphomagenesis (2, 32, 46) up-regulated in response to TaxBLV. Furthermore, METAP2 represents a validated target for inhibition of neovascularization, suggesting that it may play a complex role in tumor progression. Finally, decreases as well as increases in cellular gene expression are important regulatory events. We have identified 52 and 149 down-regulated genes in Clone2LTAXSN and SBLTCE cells, respectively (Table 1), but only 15 genes were commonly repressed. Interestingly, among these genes, we identified the proto-oncogene ARHH, the inactivation of which promotes transformation in germinal center-derived diffuse large B-cell lymphomas (52). Also striking was the down-regulation of STMN1 in Clone2LTAXSN and DOC1 in SBLTCE. STMN1, a gene regulating cell division, is decreased in hepatoblastomas and neuroblastomas (30, 49), while DOC1 is critical for allowing cells to bypass senescence (62).
Altogether, our results identify for the first time a cohort of TaxBLV-associated transcriptionally regulated genes. The good agreement between our findings and previous observations in both human B-cell tumors and HTLV-1-associated T-cell malignancies validates cross-species chip approaches for further elucidation of molecular pathways associated with oncogenesis. The genes identified as either up- or down-regulated likely play important roles in the signaling network in which TaxBLV participates. Analysis of these genes, in terms of known function, allows several important conclusions. There is no single regulatory process that is solely targeted by TaxBLV. Rather, a network of interrelated pathways is deregulated, leading to the disruption of growth control and favoring B-cell survival. Interestingly, given that we identified changes previously associated with HTLV-1, a T-cell-tropic oncoretrovirus, as well as genes linked to non-virus-associated human B-cell malignancies, our results suggest the implication of broad mechanisms that are independent of both lineage and etiology. The power of microarray analysis is that it reveals the complexity of biological events. Cell transformation is a process that affects many different aspects of cell physiology. The complexity of cell signaling may make it difficult to identify the critical or limiting events that occur during transformation. TaxBLV-responsive genes might not be unique contributors to the malignant phenotype, and additional events are required for progression to full-blown leukemia. Analysis directed towards gene expression profiling of B-cell tumors isolated from leukemic animals will permit the elucidation of molecular pathways that promote leukemia progression. Although cross-species approaches do not permit a comprehensive analysis of gene expression patterns, we provide insights into cellular events impacting B-cell transformation and pinpoint molecular targets not identified using other methods in animal models.
We thank L. Vanderweerden for website development and Troy Randall for the generous gift of mCD154 L cells.
P.K. and M.S. contributed equally to this work. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»