Cancer Research Meeting Calendar  Jordan
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

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Chuang, E. Y.
Right arrow Articles by Little, J. B.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Chuang, E. Y.
Right arrow Articles by Little, J. B.
[Cancer Research 66, 3428-3433, April 1, 2006]
© 2006 American Association for Cancer Research


Molecular Biology, Pathobiology, and Genetics

Abnormal Gene Expression Profiles in Unaffected Parents of Patients with Hereditary-Type Retinoblastoma

Eric Y. Chuang1,7, Xi Chen1, Mong-Hsun Tsai1, Hailing Yan2, Chuan-Yuan Li3, James B. Mitchell1, Hatsumi Nagasawa4, Paul F. Wilson4, Yuanlin Peng4, Markus M. Fitzek5, Joel S. Bedford4 and John B. Little6

1 Radiation Biology and Oncology Branches, Center for Cancer Research, National Cancer Institute, NIH, Bethesda, Maryland; 2 Science Applications International Corporation-Frederick, Inc., Frederick, Maryland; 3 Department of Radiation Oncology, Duke University Medical Center, Durham, North Carolina; 4 Department of Environmental and Radiological Health Sciences, Colorado State University, Fort Collins, Colorado; 5 Department of Radiation Oncology, Tufts-New England Medical Center; 6 Department of Genetics and Complex Diseases, Harvard School of Public Health, Boston, Massachusetts; and 7 Department of Electrical Engineering, National Taiwan University, Taipei, Taiwan

Request for reprints: Joel S. Bedford, Department of Environmental and Radiological Health Sciences, Colorado State University, Fort Collins, CO 80523. Phone: 970-491-7492; Fax: 970-491-7742; E-mail: joel.bedford{at}colostate.edu.


    Abstract
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 References
 
The hereditary form of retinoblastoma (Rb) is associated with a germ line mutation in one RB allele and is characterized by the occurrence of multiple, bilateral Rb tumors and a predisposition to the development of second cancers. In an earlier study, we observed an unexpected hypersensitivity to ionizing radiation in skin fibroblasts derived from unaffected parents of children with hereditary Rb. In at least four of these five families, there was no family history of Rb, indicating a new germ line mutation. We hypothesize that the increased parental cell sensitivity to radiation may reflect the presence of an as yet unrecognized genetic abnormality occurring in one or both parents of children with Rb. In the present study, we use DNA microarray technology to determine whether differences in gene expression profiles occurred in the unaffected parents of patients with hereditary Rb relative to normal individuals. Microarray analyses were validated by quantitative reverse transcription-PCR measurements. A distinct difference was observed in the patterns of gene expression between unaffected Rb parents and normal controls. By use of the prediction analysis for microarrays and principal component analysis methodologies, significant differences between the two groups were identified when as few as nine genes were analyzed. Further study of this phenomenon may offer a new insight into the genetic mechanisms of Rb and perhaps more broadly in cancer biology. (Cancer Res 2006; 66(7): 3428-33)


    Introduction
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 References
 
The hereditary form of retinoblastoma (Rb) is associated with a germ line mutation in one of the RB genes and is characterized by the occurrence of multiple, bilateral Rb tumors and a predisposition to the development of second cancers. These second tumors often occur within the radiation field in patients treated by radiation therapy (1). Interestingly, however, the family history is often negative for Rb, suggesting that a new mutation has occurred in the parental germ line. In a high percentage of cases, such new mutations arise in the paternally derived allele (2, 3). However, parental mosaicism of a mutation arising during early embryonic development of the zygote may also occur (4). Among Rb patients with high-penetrance, new germ line mutations, no known genetic abnormalities have as yet been identified in the somatic cells of their parents, although a statistically significant increase in the incidence of cancer has been reported (57).

In an earlier study, we observed an unexpected hypersensitivity to ionizing radiation in skin fibroblasts derived from unaffected parents of children with hereditary Rb (8). In at least four of these five families, there was no family history of Rb, indicating a new germ line mutation. Enhanced radiosensitivity was manifested by both an enhanced radiation-induced G1 arrest and an increased sensitivity to cell killing as measured by a clonogenic survival assay. We hypothesize that the increased parental cell sensitivity to radiation may reflect the presence of an as yet unrecognized genetic abnormality occurring in one or both parents of children with Rb. Interestingly, the radiosensitivity of two of the five probands studied fell within the reference range, consistent with earlier reports showing variability in the response of cells from patients with hereditary Rb to ionizing radiation (9, 10).

In the present investigation, we used microarray technology to determine whether differences in gene expression profiles occurred in skin fibroblast cell strains isolated from unaffected parents of patients with hereditary Rb compared with cells from normal individuals. Although initial experiments were carried out using pooled RNA from all family members, including probands, the more detailed analysis focused on cells from individual parents analyzed separately. Extensive cell culture collections from RB families are not generally available; thus, the availability of cells from the parents of RB patients was an unusual feature that allowed this study to be carried out. A distinct difference was observed in the patterns of gene expression between the fibroblasts from unaffected Rb parents and normal controls. By use of the prediction analysis for microarrays (PAM) and principal component analysis (PCA) methodology (11, 12), significant differences in expression profiles between the parental fibroblasts and normal controls were identified when as few as nine genes were analyzed.


    Materials and Methods
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 References
 
Cell culture. Clinical information and radiosensitivities concerning the skin fibroblast cell strains from the parents and probands from five Rb families, along with the range of radiosensitivities from numerous low passage normal human fibroblasts was reported previously (8). Briefly summarizing the radiosensitivities as measured by the dose required to reduce cell survival to 10% (D10) from this report, the cells derived from parents of Rb probands yielded D10 estimates ranging from about 1.0 to 2.3 Gy, whereas most of the cell strains derived from apparently normal individuals ranged from about 2.6 to 3.1 Gy. Interestingly, D10 values in the reference range were observed for cells from probands from families II and IV (2.7 and 2.8 Gy, respectively), whereas D10 values for cells from the probands for families I, III, and V (2.1, 2.0, and 1.3 Gy, respectively) were in the hypersensitive range seen for unaffected parents (8). The cell strains from Rb families were all established at the Harvard School of Public Health and have been designated according to the five families of origin as MF-4F, MF-2M, and MF-3R from family I; MF-6F and MF-7R from family II; MF-11F, MF-12M, and MF-10R from family III; MF-15F, MF-13M, and MF-14R from family IV; and MF-18F, MF-16M, and MF-17R from family V, where in each family designation the final code letter F, M, or R, respectively, refer to the father, mother, or proband (8). The low passage normal human fibroblast cultures used in this study were from the National Institute of General Medical Sciences Human Cell Repository at the Corriell Institute for Medical Research (CIMR) in Camden, NJ (8). For the present study, 12 low-passage fibroblast cell strains from apparently normal individuals also from the CIMR were used as controls. These strains are designated GM00041, GM00969, GM02673, GM03440, GM04501, GM04503, GM04505, GM05756, GM08333, GM08447, GM08429, and AG01521 and were age and sex matched to those from the Rb family members and studied at generally similar passage levels. Exponentially growing cells were maintained in Eagle's MEM (Life Technologies, Grand Island, NY) supplemented with 15% heat-inactivated fetal bovine serum, penicillin (50 units/mL), and streptomycin (50 µg/mL). Cell cultures grown to 80% confluence in P-100 plastic dishes were subcultured into larger P-150 plastic dishes. When these cultures reached ~70% confluence, the medium was renewed, and RNA was isolated 24 hours later when the cells had resumed exponential growth.

RNA extraction and purification. Cells were washed twice with ice cold PBS; 1 mL of Trizol reagent (Invitrogen, Carlsbad, CA) was added to cell pellets, which were vortexed and incubated for 5 minutes at 4°C. The mixtures were transferred into 1.5-mL centrifuge tubes and 200 µL of chloroform per mL of Trizol was added to partition the phases. The mixtures were iced for 5 minutes and then centrifuged (10,000 rpm, 4°C) for 10 minutes. The aqueous layers were then transferred into new 1.5-mL centrifuge tubes, and a 0.5 volume of 2-propanol was added. Samples were pelleted by centrifugation (14,000 rpm) for 30 minutes at 4°C and washed twice with 1 mL of 70% ethanol. Then, 100 µg of extracted RNA was further purified using Qiagen Rneasy Mini kit (Valencia, CA) according to the manufacturer's instructions.

Microarray fabrication (oligo arrays). The oligonucleotides used in this study are 60-mer with a 5'-C6 amino modification for flexible and universal covalent binding of the oligo to special substrates. Two oligo sets were kindly provided by Compugen (Jamesburg, NJ), including the Human Oligo Library set and the Human Cancer/Apoptosis Common Oligo Mini-Library set. The Human Oligo Library set contained 18,861 oligonucleotides, which represents 18,664 LEADS clusters plus 197 controls (one oligo has been designed for each gene). The Human Cancer/Apoptosis Common Oligo Mini-Library set contained 276 oligos, of which three oligos were designed for one gene. The total number of oligos used in the study was 19,137. The gene list and detailed oligo information can be found on http://www.labonweb.com/chips.

By use of an OmniGrid arrayer (GeneMachines, San Carlos, CA), 25 pL of 50 µmol/L oligonucleotide solutions in 3x SSC were spotted onto poly-L-lysine–coated slides (National Cancer Institute ROSP Human Oligo 19K Array) according to Eisen and Brown (13). Details of this array can be found on http://nciarray.nci.nih.gov. Post-processing was achieved by UV cross-linking for 600 (x100) mJ inside a Stratalinker (Stratagene, La Jolla, CA), boiling in water for 2 minutes, and dehydrating in 95% ethanol before spin drying. Alternatively, after UV cross-linking, printed slides were blocked for 15 minutes with succinic anhydride before boiling and then spin-dried.

Microarray fabrication (cDNA arrays). Microarrays contained 16,200 human cDNA clones and were prepared from two different clone sets, including the Incyte UniGEM2 set (Fremont, CA) and Research Genetics Named Genes set (Huntsville, AL). These cDNA clones are enriched for known genes. All 16K cDNAs were spotted onto poly-L-lysine–coated slides (National Cancer Institute ROSP 16K Human Array) according to Eisen and Brown (13) using an OmniGrid arrayer (GeneMachines).

RNA preparation and labeling. The labeling method was as described by Chuang et al. (14); 20 and 40 µg of total RNA were labeled with cyanine 3-dUTP (Cy3) and cyanine 5-dUTP (Cy5; Perkin-Elmer, Boston, MA), respectively, using SUPER-SCRIPT II and Oligo (dT) 12-18 from Invitrogen. The cDNA probes were compared with a reference probe that was generated from a universal human reference RNA (Stratagene), which consisted of RNAs isolated from 10 cell lines. This reference probe was used for both oligo and cDNA arrays.

Probe cleanup and array hybridization. The labeled probes were first cleaned by using Micro Bio-Spin 6 Chromatography Columns (Bio-Rad, Hercules, CA) and then were further purified using Microcon-30 spin columns (Millipore, Bedford, MA) with three washes of 400 µL TE. The final elution volume was to 17 µL.

The arrays were prehybridized with prehybridization buffer (5x SSC, 0.1% SDS, 1% bovine serum albumin) at 42°C for at least 1 hour. Slides were washed in deionized water followed by 2-propanol. Cy5- and Cy3-labeled samples were mixed with 1 µL each of human Cot1-DNA (10 µg/µL; Invitrogen), poly(A) (8-10 µg/µL; Amersham Pharmacia Biotech, Piscataway, NJ), and yeast tRNA (4 µg/µL; Ambion, Austin, TX) for hybridization. The mixed samples were denatured at 100°C for 1 minute, and after the addition of 20 µL of 2x hybridization buffer (50% formamide, 10x SSC, 0.2% SDS), the entire sample (40 µL) was loaded onto the slides. Hybridization was carried out in a humid chamber overnight in a 42°C water bath.

Washing, scanning, and image analysis. Following hybridization, the slides were washed in first wash buffer (2x SSC, 0.1% SDS), second wash buffer (1x SSC, 0.1% SDS), and third wash buffer (0.2x SSC) each for 4 minutes, then in 0.05x SSC for 1 minute. The slides were subsequently dried by brief centrifugation at 650 rpm at room temperature and scanned on a GenePix 4000A scanner (Axon Instruments, Inc., Foster City, CA) at 10-µm resolution. The Cy5- and Cy3-labeled samples were scanned at 635 and 532 nm, respectively. The resulting TIFF images were analyzed by GenePix Pro software (Axon Instruments). The ratios of the sample intensity to the reference red (Cy5)/green (Cy3) intensity for all targets were determined, and ratio normalization was done to normalize the center of ratio distribution to 1.0.

Real-time reverse transcription-PCR assay. After an initial expression analysis, seven clones were selected for validation of cDNA microarray data by use of a reverse transcription-PCR (RT-PCR) assay. In addition, one clone (GADPH) was selected as a control gene. Quantitative RT-PCR was done by taking 1 mg of total RNA, which was reverse transcribed with Superscript II enzyme (Invitrogen). Real-time PCR was done with the AB1 Prism 7900HT sequence detection system with reagents purchased from PE Applied Biosystems (Foster City, CA). Primer sequences were designed using Primer Express software (PE Applied Biosystems). Amplicons were 50 to 150 bp long.

Data analysis. Data were corrected by median centering. Clustering was done by use of Eisen's Cluster and TreeView. For unsupervised clustering, genes were chosen based on the PAM result and were from the entire sample set. PCA was done to visualize the similarity of the log-transformed expression ratios of genes. PCA reduces the dimensionality of the data (expression profile of 19,000 genes) into three dimensions by summarizing the most important (i.e., defining) parts while simultaneously filtering out noise.

Supervised class prediction analyses were done by applying PAM software hosted on National Cancer Institute's mAdb server (http://nciarray.nci.nih.gov/). The method of the nearest shrunken centroids identifies a subgroup of genes that best characterizes a predefined class (12). The prediction error was calculated by means of 10-fold cross-validation. The model is fitted on 90% of the samples, and the class of the remaining 10% is then predicted. This process is repeated 10 times, with each part playing the role of the test samples and the error of all 10 parts added together to compute the overall error. The error within the validation set reflects the number of samples wrongfully predicted to be in this set.


    Results
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 References
 
As previously reported (8), all of the cell strains derived from the unaffected parents in the Rb families, as well as strains from three of the five probands (families 1, 3, and 5), showed a moderate degree of hypersensitivity to {gamma}-radiation. All of the probands were classified as having hereditary-type Rb, based on the occurrence of multiple bilateral tumors. With the exception of Family 3, family histories were negative for Rb or other known cancer-susceptibility disorders. Thus, at least four of the cases represented new germ line mutations; the specific mutations were previously described (8). The father in family 3 had one eye enucleated at the age of 2 years for a single, unilateral Rb tumor.

In a pilot experiment, RNA samples extracted from all 14 cell strains derived from the Rb family members were pooled and compared directly with pooled RNA from 11 cell strains from apparently normal individuals by microarray analysis. A subset from the microarray analysis was then further analyzed using quantitative RT-PCR. Fifteen genes were identified through cDNA arrays, and 22 through oligo arrays that showed >2-fold increased or 2-fold decreased expression in the pooled RNA samples from the Rb families compared with the normal controls (data not shown). Among these genes, six showed significant changes in expression in both cDNA and oligo arrays. These six genes (i.e., IGFBP5, IGFBP7, ALDH1A3, JAG1, COL4A2, and COLA5A3) were selected for further analysis by quantitative real-time RT-PCR. One additional clone (GAPDH) was employed as an unrelated housekeeping gene. The results from the cDNA microarray and quantitative RT-PCR analyses are shown in Fig. 1 . As can be seen in Fig. 1, there was excellent agreement between differences in expression for these genes as measured by the cDNA microarray and RT-PCR analyses. Similar agreement was shown for the oligo microarrays (data not shown). The excellent agreement between both cDNA and oligo microarray data and RT-PCR analysis thus validates the techniques for measuring expression levels in the present study. It is noteworthy that the expression level of the GAPDH housekeeping gene was not sufficiently higher than that necessary to reach the criterion of a 2-fold change considered to be significantly different from normal background fluctuations.


Figure 1
View larger version (9K):
[in this window]
[in a new window]
 
Figure 1. Validation of cDNA microarray data by RT-PCR. Expression ratios between the pooled RB samples relative to the pooled normal samples were compared using cDNA microarrays (filled columns) or quantitative RT-PCR (open columns) for each of seven selected genes. A ratio >1 indicates upregulation in the pooled RB samples, whereas a ratio below –1 indicates down-regulation. GAPDH was used as an unrelated housekeeping gene whose relative expression level, as expected, was not significantly different from normal variations in expression levels, in keeping with the 2-fold change criterion usually applied to judge significant changes.

 
Because the focus of this study is on possible differences in gene expression for cells from a normal control population relative to the parents of Rb patients and not the patients themselves, a comparison of gene expression in cells from individuals rather than the pooled averages was considered likely to be more informative. Thus, in the following experiments, RNA samples from the nine cell strains derived from the parents in the Rb families and the 12 apparently normal cell strains were analyzed individually by use of oligo microarrays. To compare gene expression profiles between the RB and normal RNA samples based on global gene expression data, the PCA method was employed. PCA reduces the dimensionality of the 19,000 genes in this study to only very few dimensions, by summarizing the most defining parts while simultaneously filtering out noise components. With the analysis of all 19,000 genes, the difference between the RB parental and normal samples was immediately evident; the RB and normal samples formed two distinct groups on the right and left side of the first PCA dimension (data not shown). This clear separation of samples, which implied the existence of an inherent difference between them in their gene expression profiles, led to our next step, in which we identified a much smaller set of genes to be used for future grouping.

We applied the method of PAM, in which samples were divided into training and testing sets instead of training sets alone. This was designed to avoid bias prediction and to estimate error rate. Among the 19,000 genes examined, nine sets of genes were identified, all with the lowest misclassification error of 0.048. The number of genes in those sets ranged from 9 to 289 (Table 1 ). Genes in the smaller sets are wholly included in the larger sets. These sets were then examined to determine the smallest number of genes that can separate clearly the two groups of samples.


View this table:
[in this window]
[in a new window]
 
Table 1. PAM training/cross-validation result

 
PCA analysis of the above nine sets of genes showed that the one with 42 genes most effectively separated the Rb cell strains from the normal strains. These results are shown in Fig. 2A . As can be seen, the nine samples from the Rb parents (green spheres) clustered closely together and were distinctly different from those of the 12 normal controls (red spheres). The set with only nine genes also showed clear separation between the RB parents and normal controls, with the exception of a single control strain (GM00969), which seems closer to the cluster of RB samples than normal controls (Fig. 2B).


Figure 2
View larger version (40K):
[in this window]
[in a new window]
 
Figure 2. PCA mapping of the subset of 42 genes (A) or 9 genes (B) shown in Table 1 and hierarchical clustering of RNA expression observed for RB parents and normal individuals using Pearson correlation and complete linkage (C). Expression of the 42 genes is shown on a logarithmic scale. Green, for the nine retinoblastoma family parents; red, for the 12 normal control individuals. The first three dimensions of both PCA maps covered the majority of the defining part of the expression data: 91.9% for 9 genes and 80.5% for 42 genes.

 
The 42-gene set was further analyzed by use of the hierarchical cluster mapping method (Fig. 2C). Samples can be grouped into two distinct clusters, one RB parent cluster and one normal cluster. Each cluster displayed a pattern that seems to be different from the other, and these patterns were consistent among most of the samples within each group. Interestingly, normal control strain GM00969, while being clustered into the normal group, displayed some patterns that were more similar to the RB samples than to the remainder of the control group. The examination of the array cluster indicated that this individual had the least linkage to the rest of the control group and was positioned on the border of the cluster of RB samples. The same individual seemed to be outside of the core cluster by PCA analysis (Fig. 2A and B).


    Discussion
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 References
 
Microarray experiments result in a large number of data points that require extensive data mining and analysis to retrieve useful information for resolving biological questions. Therefore, the development of robust analytic tools for microarray data is very important. PCA has been widely employed as an analytic tool for this purpose. For example, Huang et al. (15) used PCA to distinguish different hepatotoxic agents and toxicity end points based upon gene expression data. PCA was also used to identify a set of genes that may play a role in tumor development (16). The analysis of global expression patterns by this approach in the present study yielded a clear indication that skin fibroblasts isolated from parents of patients with Rb have an inherently different gene expression profile than cell strains derived from apparently normal individuals. The subsequent identification of 42 genes or 9 genes provides a simpler method to differentiate the two groups and makes such separation more obvious. Identification of these small subsets of genes from global gene expression analysis may be used in the future as new biomarkers, including perhaps increased rates of germ line mutation at the Rb or other loci, possible changes in systems involved in genome stability, or increased sensitivity to radiation or other genotoxic agents.

Although the PCA analyses showed that there were distinct differences in gene expression profiles between the Rb parents and normal samples, one of the normal samples (GM00969) showed a genetic profile approaching that of the Rb parents. This raises the question as to whether other cell strains may show moderate hyper-radiosensitivity and may share certain genetic characteristics with the Rb parents. A corollary question is the extent to which the changes we have observed in specific gene expression may be related to the radiosensitivity phenotype rather than the Rb susceptibility phenotype. Interestingly, control strain GM00969 that displayed some gene expression patterns similar to the RB samples was also moderately radiosensitive; however, two other control strains showed similar hypersensitivity, whereas their expression patterns clearly clustered with the RB parents (data not shown). We are currently further investigating this question, including specifically an examination of gene expression profiles in the Rb probands, two of which were reported by Fitzek et al. to be normal (not hypersensitive) with respect to radiation sensitivity (8).

Among the 42 genes that were shown by PCA analysis to be the best set of genes for distinguishing normal control fibroblasts from RB family member–derived fibroblasts, NM23A, MCM5, HOXB2, HOXD10, HOXC10, polo-like kinase 1, and E2F1 are especially interesting because each of these seven genes has been shown to be involved in cancer development. All are associated with various aspects of cellular growth. These include cancer metastasis (NM23A; ref. 17), DNA replication (MCM5; refs. 18, 19), cellular differentiation (HOXB2, HOXD10, and HOXC10; refs. 20, 21), mitosis (polo-like kinase 1; ref. 22), and G1-S phase transition (E2F1; ref. 23). In all these cases, expression seemed to be decreased in cells from RB family members.

The exact biological significance of this gene expression pattern is not clear. The reduced expression of NM23A in the RB parents seems plausible because it had been implicated as a tumor suppressor gene. The decreased expression of the other genes is, however, puzzling and counterintuitive as most of these genes have been associated with cellular proliferation and cancer development (18, 2426). One might expect the expression of these genes to be increased in RB-associated fibroblasts, as the RB gene negatively regulates genes that are involved in cellular proliferation (e.g., E2F1), but in the present study, we are dealing with cells from unaffected parents who are not RB heterozygotes and are presumably normal for both copies of the RB gene. Another interesting observation is the absence of known DNA repair genes from the 42-gene list. Such genes might be among one group of logical candidates for parental hypersusceptibility to the development of new germ line RB mutations or hypersensitivity to radiation. This finding suggests that the moderate radiosensitivity observed in cells from the RB families is not directly related to significant alterations in the expression of a known DNA repair gene.

DNA microarray analysis done on clinical specimens provides information directly applicable to cancer diagnosis as well as disease outcome (2732). Gene expression profiling of leukemia patients was employed to distinguished acute myeloid leukemia (AML) from acute lymphoblastic leukemia (ALL); in addition, it showed that a subset of leukemia with morphologically apparent ALL had a similar gene expression profile and therapy response to AML (28). Two classes of breast tumors with or without BRCA gene mutations can be classified by different gene expression profiles (29). In view of the different origins and diversity of genetic backgrounds for the normal and RB cell strains we used in the present study, it is worth noting that despite the different origins and genetic backgrounds for the cells used in the BRCA study, the different expression profiles were, nevertheless, clear (29). In another study, molecular profiles in esophageal squamous cell carcinoma were found to differ in familial cases compared with sporadic cases (32). A few studies have related tumor gene expression profiles to treatment responses and outcomes after surgery, radiation, or cytotoxic agents (3336). A total of 92 genes selected as a discriminatory set predicted all responders after docetaxel treatment in 24 breast cancers (33). Genes related to tumor chemosensitivity were reported in patients with esophageal cancer (36). Fukuda et al. reported that genes related to apoptosis (BIRC2 and CASP6), inflammation (COX-2), DNA metabolism (CD73), cell growth (PLAU), cell adhesion (CDH1 and CDH3), transcription (MLL3), and cell cycle (CDK6) were expressed differently in radioresistant esophageal cancer cases compared with radiosensitive cases (35).

The present studies showing that there are significant differences in gene expression in somatic cells from unaffected parents of Rb patients compared with normal controls suggest that these parents may harbor some as yet unidentified genetic abnormality. Further study of this phenomenon may offer a new insight into the genetic mechanisms of Rb and perhaps more broadly in cancer biology.


    Acknowledgments
 
Grant support: Intramural Research Program of the NIH, National Cancer Institute, Center for Cancer Research; Department of Energy grants DE-FG02-03ER63635 (C-Y. Li and E.Y. Chuang), DE-FG03-01ER63235 (J.S. Bedford), and DE-FG02-05ER64089 (J.B. Little); NIH grant CA81512 (C-Y. Li); NASA grant NAG2-1629 (C-Y. Li); and National Aeronautics and Space Administration grant NAG 9-1569 (J.S. Bedford).

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.

We thank Christy L. Warner for technical assistance.

Received 8/10/05. Revised 1/11/06. Accepted 1/31/06.


    References
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 References
 

  1. Wong FL, Boice JD, Jr., Abramson DH, et al. Cancer incidence after retinoblastoma. Radiation dose and sarcoma risk. JAMA 1997;278:1262–7.[Abstract/Free Full Text]
  2. Dryja TP, Morrow JF, Rapaport JM. Quantification of the paternal allele bias for new germline mutations in the retinoblastoma gene. Hum Genet 1997;100:446–9.[Medline]
  3. Kato MV, Ishizaki K, Shimizu T, et al. Parental origin of germ-line and somatic mutations in the retinoblastoma gene. Hum Genet 1994;94:31–8.[CrossRef][Medline]
  4. Sasaki MS, Kato M, Toguchida J, et al. Somatic and germinal mutations of tumor-suppressor genes in the development of cancer. J Radiat Res (Tokyo) 1991;32 Suppl 2:266–76.
  5. Olsen JH, Boice JD, Jr., Seersholm N, Bautz A, Fraumeni JF, Jr. Cancer in the parents of children with cancer. N Engl J Med 1995;333:1594–9.[Abstract/Free Full Text]
  6. Pang D, McNally R, Kelsey A, Birch JM. Cancer incidence and mortality among the parents of a population-based series of 2604 children with cancer. Cancer Epidemiol Biomarkers Prev 2003;12:538–44.[Abstract/Free Full Text]
  7. Strong LC, Herson J, Haas C, et al. Cancer mortality in relatives of retinoblastoma patients. J Natl Cancer Inst 1984;73:303–11.[Medline]
  8. Fitzek MM, Dahlberg WK, Nagasawa H, Mukai S, Munzenrider JE, Little JB. Unexpected sensitivity to radiation of fibroblasts from unaffected parents of children with hereditary retinoblastoma. Int J Cancer 2002;99:764–8.[Medline]
  9. Garcia Heras J, Coco R. X-radiation-induced chromosome breakage in retinoblastoma lymphocytes. Mutat Res 1987;178:225–33.[Medline]
  10. Little JB, Nichols WW, Troilo P, Nagasawa H, Strong LC. Radiation sensitivity of cell strains from families with genetic disorders predisposing to radiation-induced cancer. Cancer Res 1989;49:4705–14.[Abstract/Free Full Text]
  11. Jolliffe IT. Principal component analysis. New York: Springer-Verlag; 1986.
  12. Tibshirani R, Hastie T, Narasimhan B, Chu G. Diagnosis of multiple cancer types by shrunken centroids of gene expression. Proc Natl Acad Sci U S A 2002;99:6567–72.[Abstract/Free Full Text]
  13. Eisen MB, Brown PO. DNA arrays for analysis of gene expression. Methods Enzymol 1999;303:179–205.[Medline]
  14. Chuang YY, Chen Y, Chandramouli GVR, et al. Gene expression after treatment with hydrogen peroxide, menadione, or t-butyl hydroperoxide in breast cancer cells. Cancer Res 2002;62:6246–54.[Abstract/Free Full Text]
  15. Huang Q, Jin X, Gaillard ET, et al. Gene expression profiling reveals multiple toxicity endpoints induced by hepatotoxicants. Mutat Res 2004;549:147–67.[Medline]
  16. Lu J, Liu Z, Xiong M, et al. Gene expression profile changes in initiation and progression of squamous cell carcinoma of esophagus. Int J Cancer 2001;91:288–94.[CrossRef][Medline]
  17. Cropp CS, Lidereau R, Leone A, et al. NME1 protein expression and loss of heterozygosity mutations in primary human breast tumors. J Natl Cancer Inst 1994;86:1167–9.[Free Full Text]
  18. Alison MR, Hunt T, Forbes SJ. Minichromosome maintenance (MCM) proteins may be pre-cancer markers. Gut 2002;50:290–1.[Free Full Text]
  19. Stoeber K, Tlsty TD, Happerfield L, et al. DNA replication licensing and human cell proliferation. J Cell Sci 2001;114:2027–41.[Abstract/Free Full Text]
  20. Giampaolo A, Sterpetti P, Bulgarini D, et al. Key functional role and lineage-specific expression of selected HOXB genes in purified hematopoietic progenitor differentiation. Blood 1994;84:3637–47.[Abstract/Free Full Text]
  21. Kato K, O'Dowd DK, Fraser SE, Smith MA. Heterogeneous expression of multiple putative patterning genes by single cells from the chick hindbrain. Dev Biol 1997;191:259–69.[CrossRef][Medline]
  22. Cogswell JP, Brown CE, Bisi JE, Neill SD. Dominant-negative polo-like kinase 1 induces mitotic catastrophe independent of cdc25C function. Cell Growth Differ 2000;11:615–23.[Abstract/Free Full Text]
  23. Kowalik TF, DeGregori J, Schwarz JK, Nevins JR. E2F1 overexpression in quiescent fibroblasts leads to induction of cellular DNA synthesis and apoptosis. J Virol 1995;69:2491–500.[Abstract]
  24. Johnson DG, Cress WD, Jakoi L, Nevins JR. Oncogenic capacity of the E2F1 gene. Proc Natl Acad Sci U S A 1994;91:12823–7.[Abstract/Free Full Text]
  25. Liu X, Erikson RL. Polo-like kinase (Plk)1 depletion induces apoptosis in cancer cells. Proc Natl Acad Sci U S A 2003;100:5789–94.[Abstract/Free Full Text]
  26. Takai N, Hamanaka R, Yoshimatsu J, Miyakawa I. Polo-like kinases (Plks) and cancer. Oncogene 2005;24:287–91.[CrossRef][Medline]
  27. Beer DG, Kardia SL, Huang CC, et al. Gene-expression profiles predict survival of patients with lung adenocarcinoma. Nat Med 2002;8:816–24.[Medline]
  28. Golub TR, Slonim DK, Tamayo P, et al. Molecular classification of cancer: class discovery and class prediction by gene expression monitoring. Science 1999;286:531–7.[Abstract/Free Full Text]
  29. Hedenfalk I, Duggan D, Chen Y, et al. Gene-expression profiles in hereditary breast cancer. N Engl J Med 2001;344:539–48.[Abstract/Free Full Text]
  30. Rosenwald A, Wright G, Chan WC, et al. The use of molecular profiling to predict survival after chemotherapy for diffuse large-B-cell lymphoma. N Engl J Med 2002;346:1937–47.[Abstract/Free Full Text]
  31. Su H, Hu N, Shih J, et al. Gene expression analysis of esophageal squamous cell carcinoma reveals consistent molecular profiles related to a family history of upper gastrointestinal cancer. Cancer Res 2003;63:3872–6.[Abstract/Free Full Text]
  32. van de Vijver MJ, He YD, van't Veer LJ, et al. A gene-expression signature as a predictor of survival in breast cancer. N Engl J Med 2002;347:1999–2009.[Abstract/Free Full Text]
  33. Chang JC, Wooten EC, Tsimelzon A, et al. Gene expression profiling for the prediction of therapeutic response to docetaxel in patients with breast cancer. Lancet 2003;362:362–9.[CrossRef][Medline]
  34. Cleator S, Ashworth A. Molecular profiling of breast cancer: clinical implications. Br J Cancer 2004;90:1120–4.[CrossRef][Medline]
  35. Fukuda K, Sakakura C, Miyagawa K, et al. Differential gene expression profiles of radioresistant oesophageal cancer cell lines established by continuous fractionated irradiation. Br J Cancer 2004;91:1543–50.[CrossRef][Medline]
  36. Kihara C, Tsunoda T, Tanaka T, et al. Prediction of sensitivity of esophageal tumors to adjuvant chemotherapy by cDNA microarray analysis of gene-expression profiles. Cancer Res 2001;61:6474–9.[Abstract/Free Full Text]




This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Chuang, E. Y.
Right arrow Articles by Little, J. B.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Chuang, E. Y.
Right arrow Articles by Little, J. B.


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