| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Advances in Brief |
1 University Childrens Hospital, Division of Oncology, Zurich, Switzerland; 2 Seminar for Statistics, Swiss Federal Institute of Technology, Zurich, Switzerland; and 3 Olga Hospital, Department of Pediatric Oncology/Hematology, Stuttgart, Germany
| ABSTRACT |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
Whereas no consistent and unique genetic alterations have been identified in embryonal rhabdomyosarcoma, recent molecular findings allow additional subdivision of the alveolar rhabdomyosarcoma group. A great majority of the alveolar rhabdomyosarcoma tumors are associated with unique chromosomal translocations (fusion-positive alveolar rhabdomyosarcoma). Of these,
75% are t(2;13)(q35;q14), and
25% are t(1;13)(p36;q14), leading to the fusion of either the transcription factor PAX3 in the former case or the related transcription factor PAX7 in the latter case to another transcription factor named FKHR (FOXO1A; reviewed in Ref. 3
). The resulting fusion products are composed of the NH2-terminal DNA-binding domain of PAX3 or PAX7 and the COOH-terminal transactivation domain of FKHR. These fusion proteins are more potent transcriptional activators than the wild-type PAX proteins, suggesting that PAX-specific target genes may be deregulated and thereby involved in the oncogenic transformation of the fusion-positive alveolar rhabdomyosarcoma cells.
The remaining 20% of alveolar rhabdomyosarcoma are translocation-negative (fusion-negative alveolar rhabdomyosarcoma) and form a more heterogenous, unexplored group, of which the unambiguous classification and discrimination from embryonal rhabdomyosarcoma based on classical methods such as histology still remains challenging due to the lack of immunohistochemical or ultrastructural markers. Furthermore, this group may also include cases with cryptic PAX3(7) fusions or alternative fusions with other members of the forkhead family of transcription factors (4) that cannot be identified by classical PCR diagnostics. However, clinical outcome of this subtype is as unfavorable as the one from fusion-positive alveolar rhabdomyosarcoma (5) and needs increased treatment intensity compared with embryonal rhabdomyosarcoma. Therefore, novel insight into the underlying genetics of this and the other rhabdomyosarcoma subgroups is clearly needed.
Because gene expression profiling has been shown to reliably classify a range of tumors (6, 7, 8) , we hypothesized that it should allow us to improve the molecular characterization of rhabdomyosarcoma, thereby facilitating future diagnostic classification. Hence, we performed a microarray study on a cohort of rhabdomyosarcoma biopsy samples (n = 29) and show that embryonal rhabdomyosarcoma as well as fusion-positive alveolar rhabdomyosarcoma and fusion-negative alveolar rhabdomyosarcoma have transcriptional signatures that allow identification of these distinct subtypes. Interestingly, the fusion-positive alveolar rhabdomyosarcoma-specific signature led to the identification of a novel variant translocation, which would not have been detected by specific PCR screening.
| Materials and Methods |
|---|
|
|
|---|
RNA Purification and Target cRNA Preparation
To prepare total RNA, 50100 mg of frozen tissue sections were directly homogenized in Trizol (Invitrogen, Groningen, the Netherlands). RNA precipitates were resuspended in diethyl pyrocarbonate-treated water and stored at 80°C. Only RNAs with rRNA ratios > 1.2 as measured by Bioanalyzer (Agilent, Palo Alto, CA) were used.
RNA amplification and labeling were done according to the GeneChip eukaryotic small sample target labeling technical note (Affymetrix, Santa Clara, CA). Briefly, 250 ng of total RNA was used for first and second strand cDNA synthesis with oligo-dT primers containing a T7 RNA polymerase promotor site (Microsynth, Balgach, Switzerland) using the SuperScript cDNA synthesis customer kit (Invitrogen). cRNA was prepared by in vitro transcription using the BioArray High Yield RNA Transcript Labeling kit (Enzo, Farmingdale, NY). Of the resulting cRNA, 250 ng was used for a second round of first-strand cDNA synthesis with random hexamer primers (Promega), followed by second strand cDNA synthesis with the T7-oligo-dT primers. Labeled target cRNA was then prepared in a second round of in vitro transcription using biotin-labeled CTP and UTP (Enzo).
Array Hybridization and Scanning
Fifteen µg of labeled cRNA was fragmented at 94°C for 35 min in a buffer containing 40 mM Tris-acetate (pH 8.1), 100 mM potassium acetate, and 30 mM magnesium acetate. Fragmented cRNA was mixed with a mixture containing 0.1 M 4-morpholinepropanesulfonic acid (pH 6.56.7), 1 M NaCl, 20 mM EDTA, 0.01% Tween-20, 0.5 mg/ml acetylated bovine serum albumin, 0.1 mg/ml herring sperm DNA, 1x eukaryotic hybridization control mix, and 50 pM B2 control oligo (Affymetrix). This mixture was heated to 95°C for 5 min, followed by incubation at 45°C for 5 min and centrifugation for 5 min at maximum speed. The probe array was incubated at 45°C for 16 h at constant rotation (60 rpm). Subsequent washings and stainings were done in an Affymetrix GeneChip Fluidics Station using the protocol EukGen W2V4. Stained arrays were scanned in a Gene Array Scanner (Agilent), and expression values for each gene were calculated by Affymetrix GeneChip software.
Statistical Analysis
Clustering.
Unsupervised two-way hierarchical clustering was done with the software D-chip (9)
. Expression values were base 2 logarithmized and normalized. Before clustering, genes were filtered using a variation coefficient (SD/mean) of at least 0.058 yielding 4168 nonredundant genes, which were used for an average-linkage hierarchical clustering based on Pearson correlation coefficients. Other variation coefficients gave similar results.
Gene Selection.
Supervised testing for genes with differential expression was done for all three of the groups simultaneously by computing between-group to within-group sums of squares, based on logarithmized and normalized data. Raw P values were computed from an F-distribution with 2 and 26 degrees of freedom, and they were then adjusted by the Benjamini-Hochberg false discovery rate, according to the implementation in the statistical software bundle R (package, multitest). For dimension reduction and visualization of this gene set, a principal component analysis was conducted. In binary problems, gene selection was based on Students t test statistic. Raw P values were obtained from t-distributions with corresponding degrees of freedom before the Benjamini-Hochberg false discovery rate adjustment was performed.
Classification.
In the absence of genuine test sets, the predictive potential for the classification of all three of the groups simultaneously, as well as for the binary translocation status, was evaluated by leave-one-out cross-validation. In each cross-validation run, we identified the 10, 25, 50, 100, 200, and 500 most predictive genes according to their t- or F-statistic, and used them as input for three established methods in microarray classification: a support vector machine with radial basis kernel, the 1-nearest-neighbor rule, and diagonal linear discriminant analysis. We observed the best prediction results, however, by applying an alternative method called Pelora (10)
, a classifier that, in a forward search, constructs features (so-called supervised clusters, which are meta-genes that summarize the predictively most informative individual genes) and uses them for prediction with a penalized form of logistic regression.
3'-Race and Sequencing
3'-Race was performed with the GeneRacer kit (Invitrogen). Briefly, first-strand cDNA was reverse transcribed from 1 µg of total RNA using SuperScript II RNase H reverse transcription (Invitrogen) and the GeneRacer Oligo-dT Primer and then amplified using a PAX3-specific forward primer (5'-AACGGCAGGCCGCTGCCCAACC-3') and the GeneRacer 3' Primer, followed by a nested PCR reaction using a second PAX3-specific forward primer (5'-ATCGTGGAGATGGCCCACCACG-3') and the GeneRacer 3' nested reverse primer. In the first reaction an initial denaturation at 94°C for 5 min was followed by 5 cycles at 94°C for 30 s and at 72°C for 6 min, by 5 cycles at 94°C for 30 s, at 71°C for 30 s, and at 72°C for 6 min and by 25 cycles at 94°C for 30 s, at 70°C for 30 s, and at 72°C for 6 min, followed by a final elongation at 72°C for 10 min. In the second reaction the same scheme was applied, but annealing temperatures were 70°C, 68°C, and 66°C. PCR products were cloned into the pCR-XL-TOPO cloning vector (Invitrogen) and sequenced (Dye Terminator Sequencing kit; Applied Biosystems, Applera Europe B.V., Rotkreuz, Switzerland).
Transactivation Assay
The full-length cDNA of PAX3/NCOA1 was amplified with the expand long template PCR system (Roche Applied Science, Penzberg, Germany). The cDNA was cloned into the pcDNA3 expression vector. Similar constructs were generated for PAX3, PAX3/FKHR, and PAX3 without the transactivation domain. Expression constructs were cotransfected together with a reporter plasmid containing the luciferase gene downstream of multimerized PAX3 DNA-binding sites (6xCD19 DNA-binding sites; Ref. 11
) and a plasmid containing the lacZ gene into 293T cells. Forty-eight h after transfection, cells were lysed in reporter lysis buffer (Promega, Madison, WI). ß-Galactosidase and luciferase activities were determined with the corresponding assay systems (Promega). Luciferase activity values were normalized with the ß-galactosidase activity values. Experiments were performed at least twice, each in triplicate.
Subcellular Localization of PAX3 Proteins
Constructs were transfected into 293T cells as described above and nuclear and cytosolic extracts prepared 48 h after transfection (12)
. Extracts were analyzed on Western blots using an anti-PAX3/7 antibody (clone N-19; Santa Cruz Biotechnology, Santa Cruz, CA) and visualized with the enhanced chemiluminescence system (Amersham Biosciences, Piscataway, NJ).
| Results |
|---|
|
|
|---|
To get an overview of the variability and similarity within the transcriptomes of these 29 rhabdomyosarcoma samples, we first performed an unsupervised two-way hierarchical clustering. Genes that did not show variability in expression among the rhabdomyosarcoma samples were excluded from the clustering (for filtering criteria see "Materials and Methods"). 4168 nonredundant probe sets passed the applied filtering procedure and were used for the analysis. In the resulting hierarchical clustering, the 10 translocation-positive rhabdomyosarcoma samples are clearly separated from the 19 translocation-negative ones (Fig. 1A)
. We observed a part of several hundred genes, which are consistently up-regulated in the translocationpositive samples when compared to the translocation-negative ones (see magnification of a stretch with 271 consistently differentially expressed genes) and which produce a very dominant signal. Indeed, a leave-one-out cross validation study for estimating the misclassification risk in out-of-sample classification yielded the benchmark of 0% error-rate with a wide range of genes applying four different statistical methods including a newly developed method based on gene groups called Pelora (Ref. 10
; Fig. 1B
). This demonstrates that fusion-positive alveolar rhabdomyosarcoma are characterized by an expression signature, which separates this group with very high accuracy from other rhabdomyosarcoma. Within the translocation-negative group, the fusion-negative alveolar rhabdomyosarcoma group together with the embryonal rhabdomyosarcoma.
|
To define the gene signatures that characterize each subgroup in more detail, genes that discriminate most significantly were selected by pair-wise comparison of different tumor subgroups with different P values. The number of genes that discriminate most significantly are summarized in Fig. 2A
, and a selection of 30 genes is illustrated in Fig. 2B
. The strongest signature was found for genes discriminating translocation-positive from translocation-negative samples with 299 genes at P < 0.001. The weakest signature discriminates fusion-negative alveolar rhabdomyosarcoma from embryonal rhabdomyosarcoma samples with only 28 genes at P < 0.01. The complete gene lists for these signatures can be found in Supplemental Tables 25. Taken together, these findings extend the results from the hierarchical clustering and demonstrate that rhabdomyosarcoma subgroups can be characterized by distinct gene expression signatures.
|
To test this hypothesis, we performed a series of 3'- rapid amplification of cDNA ends experiments. Using PAX3- and PAX7-specific forward primers and universal reverse primers in two rounds of nested PCR, we obtained a PCR product of
3 kbp when using PAX3-specific forward primers. Sequence analysis of this product showed that there was an in-frame fusion transcript present consisting of PAX3 and the nuclear receptor coactivator NCOA1. The translocation break point of the fusion is located in intron 6 of PAX3 and in intron 12 of NCAO1 (splice variant 2) leading to a fusion of the first 6 exons of PAX3 to the last 9 exons of the NCOA1 (splice variant 2; Fig. 3C
). Hence, the fusion protein lacks exon 7 of PAX3, which is always present in PAX3(7)/FKHR translocation products. The resulting fusion protein is 913 AA long and is composed of the 338 NH2-terminal AA of PAX3 and the 557 COOH-terminal AA of NCOA1 (AA 843-1399) plus a unique AA at the break point. Therefore, as predicted originally, the chimeric protein includes both DNA-binding domains of PAX3, paired-box and homeodomain, and a long COOH-terminal part of NCOA1 presumably acting as transactivation domain (Fig. 3D)
.
|
|
10-fold increase in transactivation over PAX3 (Fig. 4C)
18.5 fold over PAX3-TA. Because similar transactivation efficiencies have been obtained for PAX3/FKHR using the same reporter sequences in a previous study (13)
, we conclude that PAX3/NCOA1 can act as transcription factor on PAX3-regulated genes and thereby is involved, directly or indirectly, in the generation of the translocation-specific expression profile. | Discussion |
|---|
|
|
|---|
In our analysis, strongest and most consistent differences in the transcriptional pattern were found to discriminate between translocation-positive and -negative tumors. This signature is made up of several hundred genes, which are consistently up-regulated in the translocation-positive samples. This robust gene signature suggests that at least some of the genes in this pattern may be direct or indirect targets of PAX3(7)/FKHR. Moreover, the pattern is present in all of the translocation-positive samples irrespective of the PAX gene involved in the translocation, suggesting a similar influence of the two translocation products on the transcriptome. This is in accordance with the finding that PAX3 and PAX7 bind to the same DNA sequences (11) and, therefore, might have an overlapping spectra of target genes.
In an attempt to additionally define these target genes, transient or stable ectopic expression of PAX3/FKHR was carried out in a range of cell lines (HeLa, Ruch3, 293T, and Rh1; data not shown) but failed to induce the expression of these genes. These results suggest that PAX3/FKHR expression alone is not sufficient to induce the gene expression pattern identified in fusion-positive alveolar rhabdomyosarcoma biopsy samples. Cellular context such as the differentiation state or cell culture effects may be important parameters influencing the activation of target genes by PAX3/FKHR.
Up to now, few genes have been suggested to constitute targets for the PAX3/FKHR translocation product, among them NCAM (16)
, PDGFR-
(17)
, c-met (18)
, bcl-xl (19)
, as well as MYCN and arginino succinate synthetase (20)
. It is interesting to note that from these NCAM, MYCN, and argininosuccinate synthetase were also part of our expression signature. In contrast, no consistently higher expression levels in translocation-positive samples were found for PDGFR-
, c-met, and bcl-xl.
Our signature contains several genes thus far not associated with the development of rhabdomyosarcoma. Among them are numerous transcriptional regulators with either enhanced (AP2beta and enhancer-of-split) or reduced (HoxC6, HoxC10, ZFH4, and ZFP36) expression in the translocation-positive samples. Translocation-positive samples also displayed prominent expression of two proteases (ADAM10 and elastase), which could play a role in the enhanced metastatic potential of these cells. However, from a therapeutic point of view undoubtedly the most interesting genes with differential expression were found in the group of receptor proteins. In translocation-positive samples we found up-regulation by a factor 10100-fold of cannabinoid receptor 1, FGFR2, and FGFR4 (see also Ref. 7
),
-2 adrenergic receptor, as well as the acetylcholine receptor ß3. In contrast, high expression in translocation-negative samples was found for the epidermal growth factor receptor. None of these genes are expressed in normal muscle at significant levels (data not shown). For at least two of these receptors (cannabinoid receptor 1 and epidermal growth factor receptor) experimental drugs are available suggesting novel treatment strategies for rhabdomyosarcoma. This is especially remarkable because our data derive from tumor biopsies, unlike the few other expression profiling studies for rhabdomyosarcoma that are available in the literature and that are mostly restricted to rhabdomyosarcoma cell lines. In the only other study involving tumor biopsies, rhabdomyosarcoma samples were profiled against different small, round blue cell tumors (7)
, and no attempt was made to discriminate rhabdomyosarcoma subtypes.
The significance of the identified translocation-positive signature was additionally underscored with the identification of a novel PAX3 partner gene, namely NCOA1. NCOA1 is a member of the p160 steroid receptor coactivator gene family, which is composed of the three members NCOA1 (SRC-1), NCOA2 (SRC-2), and SRC-3. These coactivators interact with ligand-bound nuclear receptors to recruit histone acetyltransferases and methyltransferases to specific enhancer/promotor regions, which facilitates chromatin remodelling, assembly of general transcription factors, and transcription of target genes. In the COOH-terminal part, present in the fusion protein, two transcriptional activation domains, TAD1 and TAD2, are located. They are responsible for interaction with general transcriptional cointegrators such as p300 and CBP (TAD1), as well as with histone methyltransferases, coactivator-associated arginine methyltransferase 1 and PMRT1 (TAD2). In our assays, PAX3/NCOA1 proved to be a transcriptional activator comparable with PAX3/FKHR additionally supporting the notion that deregulated PAX3 target genes are important for tumor development.
Because the relevant expression signature was not identified in the remaining fusion-negative alveolar rhabdomyosarcoma samples, we conclude that these samples are indeed free of cryptic translocations involving PAX3(7) as found in some alveolar rhabdomyosarcoma initially classified as fusion-negative alveolar rhabdomyosarcoma (4) and that the fusion-negative alveolar rhabdomyosarcoma might be biologically more related to embryonal rhabdomyosarcoma than to fusion-positive alveolar rhabdomyosarcoma. However, our finding of a novel PAX3 translocation in an alveolar rhabdomyosarcoma sample suggests that other PAX3(7) translocations may exist in alveolar rhabdomyosarcoma. Therefore, expression profiling might be superior to specific PCR diagnostics because it also allows the identification and recognition of variant translocations. The identification of subtype-specific gene expression signatures from in vivo rhabdomyosarcoma biopsies may allow a faster and more precise molecular classification of rhabdomyosarcoma samples in the future. Certainly, it will assist in classifying cases of the translocation-negative alveolar subtype where histopathological criteria are ambiguous. Moreover, the genes that form the expression signatures displaying differential expression do also hint at novel pathways involved in rhabdomyosarcoma development and, hence, suggest new treatment strategies that could help to improve current treatment regimens.
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
Note: Supplementary data for this article can be found at Cancer Research Online (http://cancerres.aacrjournals.org). DNA sequence reported here is available from genbank, accession no. AY633656. Original array data are available from array express (www.ebi.ac.uk/arrayexpress) under accession no. E-MEXP-121.
Requests for reprints: Beat Schäfer, Division of Oncology, University Childrens Hospital, Steinwiesstrasse 75, CH-8032 Zürich, Switzerland. E-mail: beat.schaefer{at}kispi.unizh.ch
Received 3/ 9/04. Revised 6/ 1/04. Accepted 7/ 7/04.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
S. Oesch, D. Walter, M. Wachtel, K. Pretre, M. Salazar, M. Guzman, G. Velasco, and B. W. Schafer Cannabinoid receptor 1 is a potential drug target for treatment of translocation-positive rhabdomyosarcoma Mol. Cancer Ther., July 1, 2009; 8(7): 1838 - 1845. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Nishijo, Q.-R. Chen, L. Zhang, A. T. McCleish, A. Rodriguez, M. J. Cho, S. I. Prajapati, J. A.L. Gelfond, G. B. Chisholm, J. E. Michalek, et al. Credentialing a Preclinical Mouse Model of Alveolar Rhabdomyosarcoma Cancer Res., April 1, 2009; 69(7): 2902 - 2911. [Abstract] [Full Text] [PDF] |
||||
![]() |
Bin Chang, Li Juan Pang, Yan Qi, Chun Xia Liu, Ying Cao, Hong An Li, Wen Hao Hu, Jin Fang Jiang, Wen Jie Zhang, and Feng Li PAX--FKHR Fusion Genes and AChR-{gamma} in Chinese Patients With Rhabdomyosarcoma: Diagnosis Using Formalin-Fixed Archival Tissues International Journal of Surgical Pathology, February 1, 2009; 17(1): 6 - 15. [Abstract] [PDF] |
||||
![]() |
K. T. Montone, F. G. Barr, P. J. Zhang, M. D. Feldman, and V. A. LiVolsi Embryonal and Alveolar Rhabdomyosarcoma of Parameningeal Sites in Adults: A Report of 13 Cases International Journal of Surgical Pathology, February 1, 2009; 17(1): 22 - 30. [Abstract] [PDF] |
||||
![]() |
E. Davicioni, M. J. Anderson, F. G. Finckenstein, J. C. Lynch, S. J. Qualman, H. Shimada, D. E. Schofield, J. D. Buckley, W. H. Meyer, P. H.B. Sorensen, et al. Molecular Classification of Rhabdomyosarcoma--Genotypic and Phenotypic Determinants of Diagnosis: A Report from the Children's Oncology Group Am. J. Pathol., February 1, 2009; 174(2): 550 - 564. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Amstutz, M. Wachtel, H. Troxler, P. Kleinert, M. Ebauer, T. Haneke, C. Oehler-Janne, D. Fabbro, F. K. Niggli, and B. W. Schafer Phosphorylation Regulates Transcriptional Activity of PAX3/FKHR and Reveals Novel Therapeutic Possibilities Cancer Res., May 15, 2008; 68(10): 3767 - 3776. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. G. Barr and W. H. Meyer Rhabdomyosarcoma: An Overview of Biology, Clinical Features, and the Current Children's Oncology Group Studies ASCO Educational Book, January 1, 2008; 2008(1): 465 - 470. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Merlino and C. Khanna Fishing for the origins of cancer Genes & Dev., June 1, 2007; 21(11): 1275 - 1279. [Full Text] [PDF] |
||||
![]() |
D. M. Langenau, M. D. Keefe, N. Y. Storer, J. R. Guyon, J. L. Kutok, X. Le, W. Goessling, D. S. Neuberg, L. M. Kunkel, and L. I. Zon Effects of RAS on the genesis of embryonal rhabdomyosarcoma Genes & Dev., June 1, 2007; 21(11): 1382 - 1395. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Miremadi, M. Z. Oestergaard, P. D.P. Pharoah, and C. Caldas Cancer genetics of epigenetic genes Hum. Mol. Genet., April 15, 2007; 16(R1): R28 - R49. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. Davicioni, F. Graf Finckenstein, V. Shahbazian, J. D. Buckley, T. J. Triche, and M. J. Anderson Identification of a PAX-FKHR Gene Expression Signature that Defines Molecular Classes and Determines the Prognosis of Alveolar Rhabdomyosarcomas. Cancer Res., July 15, 2006; 66(14): 6936 - 6946. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Taulli, C. Scuoppo, F. Bersani, P. Accornero, P. E. Forni, S. Miretti, A. Grinza, P. Allegra, M. Schmitt-Ney, T. Crepaldi, et al. Validation of Met as a Therapeutic Target in Alveolar and Embryonal Rhabdomyosarcoma. Cancer Res., May 1, 2006; 66(9): 4742 - 4749. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Wachtel, T. Runge, I. Leuschner, S. Stegmaier, E. Koscielniak, J. Treuner, B. Odermatt, S. Behnke, F. K. Niggli, and B. W. Schafer Subtype and Prognostic Classification of Rhabdomyosarcoma by Immunohistochemistry J. Clin. Oncol., February 10, 2006; 24(5): 816 - 822. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. P. Breitfeld and W. H. Meyer Rhabdomyosarcoma: New Windows of Opportunity Oncologist, August 1, 2005; 10(7): 518 - 527. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Cancer Research | Clinical Cancer Research |
| Cancer Epidemiology Biomarkers & Prevention | Molecular Cancer Therapeutics |
| Molecular Cancer Research | Cancer Prevention Research |
| Cancer Prevention Journals Portal | Cancer Reviews Online |
| Annual Meeting Education Book | Meeting Abstracts Online |