| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Tumor Biology |
Institute for Pathology, University of Basel, CH-4004 Basel, Switzerland [F. J., J. R., P. S., G. S., M. J. M., H. M.]; Deutsches Krebsforschungszentrum, Abt. Theoretische Bioinformatik, D-69120 Heidelberg, Germany [R. D.]; Department of Electrical Engineering and Computer Science, University of California at Berkeley, Berkeley, California 94720 [C. H. P.]; Computational Biology Branch, National Center for Biotechnology Information, NIH, Bethesda, Maryland 20894 [A. A. S.], and Cancer Genetics Branch, National Human Genome Research Institute, NIH, Bethesda, Maryland 20892 [O. P. K.]
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
The study of colorectal cancer by Vogelstein et al. (5) suggested that the progression of colorectal cancer can be described by a chain of four genetic events: (a) mutation or deletion of the adenopolyposis coli, gene on chromosome 5q; (b) mutation of the K-ras gene on chromosome 12p; (c) deletion of the deleted in colon cancer gene on chromosome 18q; and (d) mutation or deletion of the p53 gene on chromosome 17p. When the first of these events occurs, the chance of the second event occurring increases, and when the second event occurs the chances of the third event increases, and so on. These events are irreversible, in that once an event occurs, it is never reversed in the future. In mathematical terms, this provides a path model for tumor progression. This model suggests that the cell starts as normal and proceeds along a path with four vertices representing different genetic changes. The presence of all four events appears to be an indicator of colorectal cancer because adenomas do not show all four alterations (5) . Whereas the path model suggests a most likely order of occurrence, colorectal cancer is really associated with the accumulation of the genetic changes on the path. In fact, the genetic changes will not always occur in the order of the path, but the path defines a preferred order.
Models for tumor progression pathways would be of obvious value to
define gene loci relevant for the early diagnosis or treatment of
cancer. However, attempts to find similar path models for other types
of cancer have not been successful. The main reason for this is that
cancer is genetically heterogeneous, even in tumors that are considered
to be clinically or pathologically homogeneous by all current tests
(6)
. Previous molecular studies of RCC have
demonstrated a high degree of genetic heterogeneity
(7, 8, 9)
. Therefore, it is important to explore models that
are not dependent on a single linear progression. The natural way of
generalizing path models to capture heterogeneity leads to tree models
(Fig. 1)
. Previously, two classes of tree models for oncogenesis, called
branching and distance-based trees, have been proposed by Desper
et al. (10
, 11)
. In these models, it is
believed that genetic events do not occur in a random fashion. Once an
event occurs, it increases the probability of other events occurring.
However, the connection between one event and the next will be specific
and directly causal in some instances, whereas in other cases, the
later events occur seemingly at random because of the basic genetic
instability in a tumor cell.
|
|
Many CGH studies have been based on the analysis of only a limited number of tumors (n < 50). Often genetic changes in tumors appear to be random and may be caused by the inherent genetic instability of cancer cells. Large data sets are needed to try to clarify the difference between the nonrandom and the random genetic changes. In this retrospective study, CGH was applied to screen for genomic alterations in a large set of formalin-fixed, paraffin-embedded tissue of CRCCs, the most frequent subtype of renal tumors. One major aim of this paper was therefore to build a model of the oncogenesis process in CRCC from CGH data using a tree-like model (10) .
| MATERIALS AND METHODS |
|---|
|
|
|---|
Comparative Genomic Hybridization
DNA Preparation.
Tumor DNA was extracted from formalin-fixed and paraffin-embedded
primary tumors. Specimens were trimmed to ensure a minimum of 75%
tumor cells in the sample. Tissue preparation and DNA extraction was as
described (9
, 20
, 21)
. One µg of tumor DNA was
nick-translated by using a commercial kit (BioNick kit; Life
Technologies, Gaithersburg, MD) and Spectrum Green- dUTPs (Vysis,
Inc., Downers Grove, IL) for direct labeling of tumor DNA. Spectrum
Red-labeled normal reference DNA (Vysis) was used for cohybridization.
CGH and Digital Image Analysis.
The hybridization mixture consisted of 200 ng Spectrum Green-labeled
tumor DNA, 200 ng Spectrum Red-labeled normal reference DNA, and 20
µg of Cot-1 DNA (Life Technologies) dissolved in 10 µl
hybridization buffer [50% formamide, 10% dextran sulfate,
2 x SSC (pH 7.0)]. Hybridization, image acquisition,
image analysis, and control experiments were exactly as described
(20
, 21)
. At least four observations/autosome and two
observations/sex chromosome were included in each analysis. Each CGH
experiment included a tumor cell line (Spectrum Green MPE-600, Vysis)
with known aberrations (positive control) and a hybridization of two
differentially labeled sex-mismatched normal DNAs to each other
(negative control). Sex-mismatched normal controls were also used to
test the ability of each metaphase batch to allow for a linear
relationship between fluorescence intensities and DNA sequence copy
numbers. Metaphases were used only if the color ratio of sex-mismatched
normal DNAs was
0.66 at the X chromosome (14)
.
The thresholds used for definition of DNA sequence copy number gains and losses were based on results of CGH analyses of formalin-fixed normal tissues (14 , 20 , 22 , 23) . A gain of DNA sequences was recorded at chromosomal regions where the hybridization resulted in a tumor:normal ratio >1.20. Overrepresentations were considered amplifications when the fluorescence ratio values exceeded 1.5 in a subregion of a chromosome arm. A loss of DNA sequences was recorded where the tumor:normal ratio was <0.80. To define an aberration, additionally it was required that the first SD was above (gain) or below (deletion) 1.00. Because in negative control hybridizations the mean green:red ratio occasionally exceeded the fixed 1.2 cutoff level in normal tissues at 1p, 16p, 19, and 22, these guanine-cytosine-rich regions, known to produce false positive results by CGH, were excluded from all analyses (14) . CGH profiles were only evaluated if the quality of the CGH experiment was appropriate (14) .
| Statistics |
|---|
|
|
|---|
| Oncogenetic Trees |
|---|
|
|
|---|
It is neither practical nor useful to make a model of all 72 possible events. Most of these events occur only in a few tumors and are probably truly random genetic changes because of the inherent instability of tumor cells or because of false positives attributable to imperfections in the CGH analysis. Therefore, it is preferable to choose a small set of events that occur sufficiently frequently to appear nonrandom. We implemented the well-known method of Brodeur et al. (26) to select a set of events that appear to occur nonrandomly in this data set of 116 RCC tumors. The method of Brodeur et al. uses random simulation to derive a distribution of event lists for tumors under the null hypothesis that all events occur randomly. To generate a sample from the null distribution, one needs a prior probability for each event to occur. We derived a prior distribution by assuming that gains and losses are equally likely, and the probability of a gain/loss on an arm is proportional to the length of that arm; we used one of the tables of chromosome arm lengths provided in Ref. 27 . We generated 10,000 replicates of sets of 116 tumors sampled from the null distribution. For each replicate, a score is computed relating each event to its prior probability, and the maximum score, which represents a measure of how often the most frequent event occurs in random data, is recorded. Scores are also computed for the 72 events in the real data, and an event is considered nonrandom if its score in the real data exceeds the 95th percentile of the 10,000 maximum scores from the random replicates. The method of Brodeur et al. (26) selected 12 events (-3p, -4p, -4q, -6q, -8p, -9p, -13q, -18q, -Xp, +17p, +17q, and +Xp) as nonrandom events in the set of 116 RCC tumors.
Branching Trees.
Recently, we proposed a method for reconstructing a tree model of
oncogenesis (10)
based on a computer science technique
called a maximum weight branching (28)
. The tree models
derived by this method are therefore called "branching
trees." In a branching tree, there is one node called a root,
and every other node is one of the events. An edge
i > j represents a hypothetical
cause-and-effect relationship meaning that the occurrence of event
i makes the occurrence of j more likely.
Branching trees are natural generalizations of Vogelsteins path model
(5)
because there can be multiple edges coming out of each
node that intuitively represent the heterogeneous possibilities for how
oncogenesis can progress. The choice of which edges to include in the
branching tree is based on a weight function that takes into account
how often each event occurs and how often each pair of events occur
together in the same tumor. If there are n events and 1
root, the tree will always have n edges. We choose the tree
for which the sum of the weights on the n edges is
maximized. Although the number of possible trees grows
exponentially as a function of n, it is possible to
choose the maximum weight tree quickly using an algorithm
developed in the field of combinatorial optimization
(28, 29, 30, 31)
.
The central theoretical result of Desper et al. (10) is that under plausible assumptions about the stochastic process of oncogenesis, and given enough samples, the maximum weight branching method provably converges to the correct tree model. Although in practice we generally do not have enough tumor samples for the convergence result to apply, we still would expect that the maximum weight branching tree has most of the edges correct. The technical assumptions needed to prove that the convergence results were inspired by some work on similar phylogenetic trees, now called Cavender-Farris trees because of the seminal papers (32 , 33) .
Distance-based Trees.
A different class of tree models that uses existing phylogenetic tools
more directly was described previously by Desper et al.
(11)
. The trees constructed by this approach have all of
the events at leaves, whereas the internal nodes are hidden, unnamed
events, much like a phylogenetic tree has the existing species as
leaves, and the hypothetical common ancestors as internal nodes (Fig. 1)
. Another significant difference is that, in these trees, each edge
has a length, and we draw the trees so that horizontal distance along
an edge is proportional to length. The trees are constructed by first
defining an n x n distance
matrix that describes for each pair of events whether they tend to
occur together or not. The second step is to use existing phylogenetic
methods to find the phylogenetic tree that best fits the distance
matrix. For this reason the trees constructed by this method are called
"distance-based trees." The problem of fitting a tree to a matrix
of distances has been extensively studied in phylogenetics
(34, 35, 36, 37)
. The distance between two events i and
j is the sum of the distances on the edges between them in
the tree.
Typically, distance-based methods produce unrooted trees, or root the resulting phylogeny by an outgroup. We added an extra leaf, to serve as the root of the tree. In our model, each edge e = (i,j) had an associated probability pe, representing the conditional probability that event j occurs, given that event i occurs. To use a distance method, there must be an underlying additive metric. As probabilities are multiplicative in this model when directed away from the root, we used the logarithm to transform the edge probabilities to an additive metric, such that de = -log pe for each edge e.
Under this transformation, the distance from the root to any other event i is -log (pi), where pi is the probability of observing event i. Letting pij be the probability of observing both events i and j, the distance from i to j is defined as -2 log pij + log pi + log pj.
Although distance-based trees are not a direct generalization of the Vogelstein path model, they may provide a statistically robust answer to the basic questions:
(a) Which events are early? Those near the root. The distance to the root is the negative logarithm of the probability, so a shorter distance means a higher probability, i.e., an earlier event.
(b) Which events mark subclasses of tumors? Those which cluster together in subtrees.
The relative advantages and disadvantages of the branching trees and distance-based trees are described at length in the "Discussion" section. In practice, we find it useful to construct both trees and to compare the predictions they make about oncogenesis.
For the distance-based tree presented here, first we used the method of Desper et al. (11) to construct a distance matrix. Then we used the Fitch and Neighbor programs from the PHYLIP package (version 3.5c; Ref. 38 ) to fit the distance matrix to phylogenetic tree topologies with distances on each edge. The edge lengths were converted back to underlying probabilities, and we found that the tree inferred by the Fitch program better fit the probabilities of event occurrences in the original data.
| RESULTS |
|---|
|
|
|---|
The number of aberrations was not associated with histological grade or pT stage. To screen for specific chromosomal alterations that might be linked to tumor progression, the association between the most frequent aberrations, histological grade, and pT stage was analyzed. Chromosome 6q loss was detected in 14 of 48 pT1/2 tumors (29%) and in 31 of 67 pT3/4 RCCs (46%), but this trend did not reach significance (P = 0.06). Interestingly, gains of chromosomes 5q, 17p, and 17q were associated with good tumor differentiation. Chromosome 5q gains were detected in 10 of 46 grade 1/2 tumors (22%), and in 6 of 70 grade 3/4 tumors (9%; P = 0.05). Chromosome 17p gains were seen in 67% of grade 1/2 tumors, but in only 38% of grade 3/4 tumors (P = 0.008). Gains of chromosome 17q were detected in 49% of grade 1/2 tumors and in 29% of grade 3/4 tumors (P = 0.03). Chromosome 9p, 13q, and 18q losses were not associated with tumor grade or stage (contingency table analysis).
Construction of Oncogenetic Trees for Renal Carcinoma Progression
Using the CGH data of 116 CRCC, a branching tree and a
distance-based tree were constructed. The branching tree shown in Fig. 2
is based on a maximum branching in a weighted graph. A branching
tree has a root vertex representing a normal cell, whereas the other
vertices represent observed CNAs of interest. An edge from CNA
i to CNA j indicates that the occurrence of
i increases the probability for the occurrence of
j in branching trees. The branching tree is generated purely
from considering individual probabilities and pairwise joint
probabilities. A strongly correlated pair is usually presented in the
tree with an edge from one to the other. The placement of leaves
adjacent to the root, such as -3p and -8p, indicates that these two
events are not more strongly correlated with some events than with
other events. It leaves unclear whether -3p and -8p are important
early events or not. In contrast, -4q, which is next to the root and
the root of a large subtree, is clearly predicted to be an important
early event. The edge lengths in the branching tree are irrelevant.
|
|
| DISCUSSION |
|---|
|
|
|---|
Mathematical models for tumor progression are of interest to characterize pathways whereby oncogenesis proceeds. For example, the path model is a well-established hypothesis for tumor progression of colorectal cancer (5) . However, analysis of the complex CGH data suggests that a path model is not useful for other types of cancer (41) . It has been assumed that the genetic events take place in a more tree-like than path-like pattern (10) . The analogy between evolving tumors and evolving species had been observed by Buetow et al. (42) , who used phylogenetic methods to genetically classify a set of liver tumor samples. We recently developed a tree model to analyze genetic progression of solid tumors, and showed that an algorithm based on maximum-weight branching in a graph correctly infers the topology of the tree. This group of tree models has been called branching trees (10) . Using tree-fitting algorithms developed by researchers in phylogeny, we constructed a different type of tree model called distance-based trees, in which events are leaves of the tree (11) .
Remarkably, both tree models shared many relevant properties and allowed us to derive new hypotheses for CRCC progression: (a) there are important subtrees emanating from a common root; (b) a subtree with -13q, -9p, and -18q; (c) a subtree with -6q, +17q, +17p; (d) a close relationship between -4q and -4p; and (e) -3p and -4q are important early events near the root.
The subtrees in our models predict that there may be at least two subclasses of CRCC loosely associated with -4q: one subclass in which +17q is tightly linked with and precedes +17p in a side branch closely related to -6q; and another subclass that is characterized by the events -9p, -13q, and -18q.
The event -4q is closer to the first subclass, but may represent a third subclass because it appears in its own subtree, or it may be closer to independent as indicated by the long branch leading to it. Losses of 4q were observed in 58% of RCCs by LOH studies (43 , 44) but in 1050% by CGHs (3 , 16 , 21 , 45) . Interestingly, losses on 4q were detected in 4677% of colon, cervix, and hepatocellular carcinoma, suggesting the presence of tumor suppressor genes on this chromosomal arm (41 , 46 , 47) .
Interestingly, this CGH analysis revealed a significant association between 13q losses and poor prognosis. Previous studies examining LOH have shown that 13q allelic imbalances can occur in CRCC. The rates of 13q LOH ranged from 4 to 57% (3 , 45 , 48, 49, 50, 51, 52) . Our finding that chromosome 13q loss is associated with poor patient outcome is consistent with a previous study of Lai et al. (52) . This group found that LOH of the retinoblastoma gene, located on chromosome 13q14, was detectable in 4 of 7 informative RCCs (57%) and was associated with high tumor grade, high stage, and poor prognosis. This raises the possibilities that inactivation of one or more gene(s) on 13q might be relevant for CRCC progression.
In this CGH analysis of 116 CRCC, we confirmed previous studies suggesting that 9p losses are related to short recurrence-free survival in pT3 CRCC (3 , 53) . In this study, 9p losses were only weakly related to poor tumor-specific survival in pT3 CRCCs by CGH. Although this trend did not reach significance, a microsatellite analysis from our laboratory demonstrated a clear relationship between LOH on 9p and poor prognosis in papillary RCC (54) and in CRCC.4 According to our preliminary mapping analyses, a gene on 9p13-14 might be relevant for both types of RCC.
Losses of the short arm of chromosome 8 have been recently considered as important events in RCC progression because 8p losses were associated with higher tumor grade in nonpapillary CRCC (55) . LOH on 8p was detected in 1833% of RCC by microsatellite analysis (44 , 55) . The extremely long branch to -8p in the distance-based tree suggests that this event is more likely to be a late effect than an early cause. This finding also corresponds to previous biological observations. We have recently studied CRCC metastases and their primary tumors by CGH (9) . Loss of 8p was a frequent finding in metastases but was quite rare in corresponding primary tumors. This finding is consistent with the observation derived from the distance-based tree model that 8p loss is a late event in primary CRCC. If one assumes that 8p losses are late events according to the present tree analysis, and our previous study found them associated with metastases, it is tempting to speculate that the presence of 8p loss in the primary tumors would be predictive of metastatic progression. However, cell clones with 8p loss may represent minor cell populations in the primary tumor and are therefore not detectable by CGH because CGH allows only detection of clonal copy number aberrations, which are present in >5075% of the tumor cells (13 , 14) . Therefore, one would need other techniques (such as FISH) to detect these individual, possibly metastasis-prone cells.
Both tree models have important advantages and disadvantages. The distance-based method considers all pairwise correlations simultaneously (11) . If an event is very weakly correlated with a large number of events, it will be pulled to the subtree containing those events. Such a phenomenon explains the fact that -3p is pulled toward the center of the distance-based tree. The relatively long length of the edge from the root to the first central node suggests a general phenomenon that the occurrence of any one event makes any other event more likely. Comparing how we use the branching and distance-based tree models to make predictions illustrates some of the strengths of the distance-based method. For both types of trees, early events should be close to the root. The distance-based approach quantifies this precisely, and allows us to infer an order of events by ranking according to distance to the root.
In the branching tree approach (10) , distances are measured less precisely by numbers of edges. A CNA is adjacent to the root in the branching, but has no subbranches (such as -8p in the renal cancer branching tree). When a CNA is a branch adjacent to the root in the branching tree (such as -3p and -8p) it is unclear whether this is because it is an early event that precedes the others but is not correlated especially well with any of them, or whether it is a relatively late event that is also not well correlated with a subclass of tumors. In this case, the distance-based tree clarifies that -3p is of the first type, and -8p is of the second type.
Both methods tend to cluster CNAs that occur together in subtrees, but the distance-based method is more robust in that it preserves information about the co-occurrence of all pairs of CNAs, whereas the branching tree only shows the best-correlated pairs. One advantage of the branching tree is that edges between CNAs lead to direct predictions of cause-and-effect relationships. In practice, we constructed both types of trees. We looked for similarities in which events were near the root to predict the early CNAs. We looked for similar clustering in subtrees to predict which events tend to occur together and may mark genetically homogeneous subsets. The two tree construction methods are encoded in free software.5 The branching tree method is entirely contained in the software we distribute. For the distance-based method, our software goes as far as producing a distance matrix suitable for input to PHYLIP (38) .
Since the publication in 1988 by Vogelstein et al. (5) of the path model for colorectal cancer, it has been hoped that similar models could be found for other types of cancer; but this search has been largely fruitless. Because the CGH technique allowed a genome-wide view of genomic alterations in solid tumors, it has been hoped that CGH will clarify which genetic changes are nonrandom; this has had mixed results (16) . In both cases it appears that a central problem is that many types of solid tumors, such as CRCC, are very heterogeneous and appear to have complex patterns of genetic changes. Furthermore, it cannot be excluded that small deletions <10 Mb are missed by CGH (12) . However, mathematical modeling based on the primary CGH data are desirable to propose models more general than path models, to elucidate which genetic changes are most worthy of additional study, and to suggest some good hypotheses for additional investigation. To get useful results from mathematical modeling, it is necessary to collect large tumor data sets such as the one described herein, so that small numbers of truly random genetic events do not obscure the basic patterns.
In summary, we proposed in this study tree models for CRCC oncogenesis that lead to very similar predictions, although they are constructed by rather different mathematical methods. Two interesting hypotheses that emerge from our study are: (a) -4q is an important early event in CRCC worthy of additional study; and (b) there are two subclasses of CRCCs, one of which is marked by -9p, -13q, and -18q, and the second of which is marked by -6q, +17q, and +17p. The results of the present study and previous analyses (3) indicate that chromosome 13q and 9p are chromosomal loci, which may harbor important tumor suppressor genes for RCC progression. More generally, we described a study illustrating the utility of two techniques to construct tree models of oncogenesis that may be applicable to other large CGH data sets.
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
1 Supported by the Swiss National Science
Foundation (97.3100-050752), Krebsforschung Schweiz (KFS 367-9-1996),
Basler Krebsliga (3/97) and the National Science Foundation Grant
CCR-9626361 (to C. H. P.) ![]()
2 To whom requests for reprints should be
addressed, at Institute for Pathology, University of Basel,
Schönbeinstrasse 40, CH-4003 Basel, Switzerland. Phone: 061-265
28 90; Fax: 061-265 31 94; E-mail: hmoch{at}uhbs.ch ![]()
3 RCC, renal cell carcinoma; LOH, loss of
heterozygosity; FISH, fluorescence in situ
hybridization; CGH, comparative genomic hybridization; CRCC,
conventional (clear cell) renal carcinoma; CNA, copy number
aberration. ![]()
4 P. Schraml, K. Struckmann, R. Bednar, W. Fu, T.
Gasser, K. Wilber, J. Kononen, G. Sauter, M. J. Mihatsch, and Moch, H.
CDKN2A mutation analysis, protein expression and deletion mapping of
chromosome 9p in conventional (clear-cell) renal carcinomas: Evidence
for a new tumor suppressor gene proximal to CDKN2A, submitted for
publication. ![]()
5 E-mail address: desper@ncbi.nlm.nih.gov. or
schaffer@helix.nih.gov. ![]()
Received 1/27/00. Accepted 9/12/00.
| REFERENCES |
|---|
|
|
|---|
-fetoprotein production. Gastroenterology, 110: 184-192, 1996.[Medline]
This article has been cited by other articles:
![]() |
M. I. Zhou, R. L. Foy, V. C. Chitalia, J. Zhao, M. V. Panchenko, H. Wang, and H. T. Cohen Jade-1, a candidate renal tumor suppressor that promotes apoptosis PNAS, August 2, 2005; 102(31): 11035 - 11040. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Rahnenfuhrer, N. Beerenwinkel, W. A. Schulz, C. Hartmann, A. von Deimling, B. Wullich, and T. Lengauer Estimating cancer survival and clinical outcome based on genetic tumor progression scores Bioinformatics, May 15, 2005; 21(10): 2438 - 2446. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Struckmann, P. Schraml, R. Simon, K. Elmenhorst, M. Mirlacher, J. Kononen, and H. Moch Impaired Expression of the Cell Cycle Regulator BTG2 Is Common in Clear Cell Renal Cell Carcinoma Cancer Res., March 1, 2004; 64(5): 1632 - 1638. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Rennstam, M. Ahlstedt-Soini, B. Baldetorp, P.-O. Bendahl, A. Borg, R. Karhu, M. Tanner, M. Tirkkonen, and J. Isola Patterns of Chromosomal Imbalances Defines Subgroups of Breast Cancer with Distinct Clinical Features and Prognosis. A Study of 305 Tumors by Comparative Genomic Hybridization Cancer Res., December 15, 2003; 63(24): 8861 - 8868. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. A. Gatenby and E. T. Gawlinski The Glycolytic Phenotype in Carcinogenesis and Tumor Invasion: Insights through Mathematical Models Cancer Res., July 15, 2003; 63(14): 3847 - 3854. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. A. Gatenby and B. R. Frieden Application of Information Theory and Extreme Physical Information to Carcinogenesis Cancer Res., July 1, 2002; 62(13): 3675 - 3684. [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 |