| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Advances in Brief |
Carcinogenesis Division [T. K., S. Y., E. O., N. W., T. S., T. U.] and Medical Genomics Center [T. O., M. O.], National Cancer Center Research Institute, Tokyo 104-0045, and Department of Pathology, Osaka City University Medical School, Osaka 545-8585 [K. M., S. F.], Japan
| ABSTRACT |
|---|
|
|
|---|
8000 rat genes. By hierarchical clustering analysis, the seven carcinomas were classified into two groups that exactly coincided with the DMBA-induced and the PhIP-induced groups. The correlation coefficient between the two groups was 0.63, and those between any carcinomas within each group ranged from 0.78 to 0.95. In addition, characteristic clusters of genes were also identified that highlighted distinct and common characteristics of both groups. Seventeen genes were down-regulated in the DMBA and up-regulated in the PhIP-induced groups. Thirty-three genes were regulated in the opposite manner. Our results indicated that gene expression profiles in cancers reflect their etiology and suggested a possibility that etiology of cancers could be retrospectively estimated from their expression profiles. | Introduction |
|---|
|
|
|---|
cDNA microarray technology has provided new information in many aspects of tumor biology (6) . The expression profiles in cancers were shown to reflect their histological characteristics, clinical outcomes, and the responses to treatment (7, 8, 9, 10, 11, 12) . In addition, the expression profiles were used to develop clinical biomarkers (13, 14, 15) . However, no reports are available regarding the correlation between the etiology of human cancers and their gene expression profiles, except for hepatocellular carcinomas caused by hepatitis viruses (16) . This is because the etiologic factors of human cancers are complex and very difficult to determine.
For this reason, animal models are expected to provide good resources, as animal cancers are induced by defined protocols and the carcinogenic factors involved are clear. Rat mammary carcinomas can be reproducibly induced by well-established methods using DMBA3 (17) , PhIP (18) , and NMU (19) . The induced mammary carcinomas predominantly originate from mammary ducts (17) and show smaller phenotypic variations (20) and fewer genetic alterations than the human counterpart (21 , 22) . In addition, normal mammary glands can be purely isolated from inguinal fat by the gland isolation technique.
In this study, we analyzed expression profiles of rat mammary carcinomas induced by DMBA and those induced by PhIP, a representative food-born carcinogen (23) , using a cDNA microarray, and we classified these carcinomas based on their gene expression profiles and estimated whether the profiles reflect their etiology.
| Materials and Methods |
|---|
|
|
|---|
GeneChip Protocol.
GeneChip Rat Genome U34A arrays, which contained probe sets for
8000 rat genes, were purchased from Affymetrix (Santa Clara, CA). Five µg of total RNA were used as starting material for cDNA preparation. The first-strand cDNA was synthesized with the SuperScript II reverse transcriptase (Invitrogen, Groningen, The Netherlands) and a T7-(dT)24 primer (Amersham Pharmacia Biotech, Buckinghamshire, United Kingdom). The double-strand cDNA was synthesized with Escherichia coli RNase H, E. coli DNA polymerase I, and E. coli DNA ligase (Toyobo, Tokyo, Japan). Biotin-labeled cRNA was prepared using a HighYield RNA transcript labeling kit (Affymetrix). After in vitro transcription, the unincorporated nucleotides were removed using the RNeasy Mini kit (Qiagen, Valencia, CA). Twenty µg of each labeled cRNA were fragmented, and its quality was assessed by gel electrophoresis.
The labeled cRNA was placed in a hybridization mixture containing four biotinylated hybridization controls (BioB, BioC, BioD, and Cre), as recommended by the manufacturer, and the U34A arrays were hybridized for 16 h at 40°C at constant rotation (60 rpm). After washing, the GeneChip arrays were stained with streptavidin-phycoerythrin conjugate (Molecular probes, Eugene, OR) and then scanned with a GeneArray scanner (Hewlett-Packard, Palo Alto, CA). The scanned images were processed using an Affymetrix GeneChip Analysis Suite (Version 4.0.1) to obtain the "average difference value" for each probe set. The image from each array was scaled so that the average of average difference values of all probe sets was adjusted to 1000. The values of "fold change" were obtained with normalization for all probe sets. Scaled average difference values and fold change data from each GeneChip array were exported to flat text files and used for statistical analysis. The complete set of data are available at our web site.4
Data Analysis.
Statistical comparisons were performed on the 2000 transcripts with the highest intensity from each GeneChip array. Before clustering, the fold change values of the genes were log10 transformed. The log-transformed fold change values were filtered for presence in six or more of the seven carcinomas. Average-linkage hierarchical clustering of an uncentered Pearson correlation similarity matrix was applied with the program Cluster, and the figures were generated with the program TreeView (25)
. Similarities of gene expression profiles between two given carcinomas were assessed by the Pearson correlation coefficient (25)
.
Quantitative RT-PCR.
cDNA was synthesized from total RNA with oligo (dT)1218 primer and Superscript II reverse transcriptase (Invitrogen). Quantitative RT-PCR analysis was performed using a iCycler iQ detection system (Bio-Rad Laboratories, Hercules, CA) with Sybr Green PCR Core Reagents (Applied Biosystems). The sequences of the primers and annealing temperature are listed in Table 1
. The number of molecules of a specific gene in a sample was measured by comparing its amplification with the amplifications of standard samples that contained 101-106 copies of the gene and was normalized to that of cyclophilin (26)
. Fold change of the gene in a carcinoma was calculated by dividing the normalized value of the carcinoma by that of the normal control mammary duct.
|
| Results |
|---|
|
|
|---|
Gene Expression Profiling of the Mammary Carcinomas.
The four DMBA-induced and three PhIP-induced mammary carcinomas were analyzed by GeneChip microarrays. From the 2000 genes, a set of 1564 genes was selected for cluster analysis with the criterion that the data of each gene in the set were present in six or seven carcinomas. Hierarchical cluster analysis separated these seven carcinomas into two main branches in concordance with the carcinogens used to induce them (Fig. 1)
. The correlation coefficients between two of the four DMBA-induced carcinomas ranged from 0.88 to 0.95 and those between two of the three PhIP-induced carcinomas ranged from 0.78 to 0.80. The correlation coefficient between the DMBA-induced and the PhIP-induced carcinomas was 0.63. These results showed that, although the gene expression profiles of the chemically induced rat mammary carcinomas were relatively similar, distinct gene expression profiles existed in the DMBA-induced and the PhIP-induced mammary carcinomas.
|
On the other hand, clusters C1, C2, and C3 consisted of genes up-regulated in the DMBA group while down-regulated in the PhIP group, and genes in cluster C3 were more clearly up-regulated than those in clusters C1 and C2. All of the 33 genes in cluster C3 showed significantly higher expression levels in the DMBA group than in the PhIP group (P < 0.01 by Students t test). Glutathione S-transferase IV(Gstp1, Rn.44821), cathepsin D, and MHC genes were in this cluster (Fig. 1C)
.
Cluster B consisted of genes down-regulated in both the DMBA and PhIP groups, whereas cluster D consisted of genes up-regulated in both groups. Of the 83 genes in cluster B, 37 (45%) were down-regulated with fold change values >2.5 in all seven mammary carcinomas, including tissue inhibitor of metalloproteinase-1(Timp1, Rn.25754) and transforming growth factor-ß-inducible early response gene (Tieg, Rn.2398). Of the 109 genes in cluster D, 68 (62%) were up-regulated with fold changes >2.5 in all seven mammary carcinomas, including cyclin D1 (Ccnd1, Rn.9471) and DNA methyltransferase (Dnmt1, Rn.6955).
Confirmation of Data from GeneChip Analysis by Quantitative RT-PCR Analysis.
The data obtained by the GeneChip analysis was confirmed by quantitative RT-PCR analysis of 12 genes, which were selected as those with typical expression in each cluster (Table 2)
. Only when fold changes in GeneChip analysis were between -1.3 and 1.6, changes in opposite directions were occasionally observed by quantitative RT-PCR analysis. Specific changes observed in clusters A1, B, C3, and D were reproducible in the quantitative RT-PCR analysis.
|
| Discussion |
|---|
|
|
|---|
8000 rat genes in four DMBA-induced and three PhIP-induced rat mammary carcinomas revealed that the carcinomas could be classified into two carcinogen-specific groups by nonbiased hierarchical analysis. Any two carcinomas within the DMBA-induced or PhIP-induced groups showed much higher correlation coefficients (>0.78) than the average among all of the carcinomas (Fig. 1)As for mechanisms of how expression profiles specific to DMBA- or PhIP-induced mammary carcinomas were generated, the following mechanisms can be considered. First, carcinogens are known to induce genetic alterations relatively specific to each carcinogen (28, 29, 30) . These genetic alterations would affect their downstream signaling pathways, which would be reflected in the expression profiles in the resultant tumors. DMBA is known to induce A:T to T:A transversions predominantly (31) , whereas PhIP induces G:C to T:A transversions and G:C deletions (24) , and these differences are considered to lead to different target genes. The H-ras mutation frequencies and incidences of loss of heterozygosity are reported to be different between the PhIP-induced and DMBA-induced mammary carcinomas (22 , 32, 33, 34) . Secondly, carcinogens are known to induce specific responses in normal cells (35 , 36) . DMBA has a stronger effect on serum prolactin and estradiol levels than PhIP, whereas PhIP shows stronger suppression of apoptosis in the mammary glands (37) . Even after repeated clonal selections during carcinogenesis, there is a possibility that some of the specific responses might be retained.
Genes differentially expressed between DMBA- and PhIP-induced mammary carcinomas are expected to be related with the molecular pathway specific to each group. Genes up-regulated in the DMBA-induced carcinomas and down-regulated in PhIP-induced carcinomas (cluster C3) included genes related to immune reactions such as a rat major histocompatibility gene and tumor rejection antigen 1 (38) . This suggested that stronger immune responses were induced in the DMBA-induced carcinomas. Genes up- or down-regulated commonly in DMBA- and PhIP-induced mammary carcinomas were expected to include the genes that were related to phenotypes as mammary carcinomas. Increased expression of cyclin D1 observed in both groups was in accordance with the previous report (39) .
We showed the presence of etiology-specific expression profiles in chemically induced rat mammary carcinomas. This suggested the presence of etiology-specific expression profiles in human cancers also, and this could provide clues to identifying carcinogenic agents in a certain cancer based on its gene expression profile.
| FOOTNOTES |
|---|
1 Supported by a Grant-in-Aid for the Second-term Cancer Control Strategy and by a Grant-in-Aid for the Millennium Genome Project from the Ministry of Health, Labor and Welfare. ![]()
2 To whom requests for reprints should be addressed, at Carcinogenesis Division, National Cancer Center Research Institute, 5-1-1, Tsukiji, Chuo-ku, Tokyo, 104-0045, Japan. Phone: 81-3-3547-5240; Fax: 81-3-5565-1753; E-mail: tushijim{at}ncc.go.jp ![]()
3 The abbreviations used are: DMBA, 7,12-dimethylbenz[a]anthracene; PhIP, 2-amino-1-methyl-6-phenylimidazo[4,5-b]pyridine; NMU, N-methyl-N-nitrosourea; RT-PCR, reverse transcription-PCR. ![]()
4 Internet address: http://www.ncc.go.jp/research/rat-genome/. ![]()
Received 2/15/02. Accepted 5/17/02.
| REFERENCES |
|---|
|
|
|---|
5'-GGA-3' mutation of the Apc gene in rat colon tumors induced by 2-amino-1-methyl-6-phenylimidazo[4,5-b]pyridine. Proc. Natl. Acad. Sci. USA, 92: 910-914, 1995.This article has been cited by other articles:
![]() |
L. Shan, M. Yu, and E. G. Snyderwine Global Gene Expression Profiling of Chemically Induced Rat Mammary Gland Carcinomas and Adenomas Toxicol Pathol, December 1, 2005; 33(7): 768 - 775. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Thangaraju, M. Rudelius, B. Bierie, M. Raffeld, S. Sharan, L. Hennighausen, A-M. Huang, and E. Sterneck C/EBP{delta} is a crucial regulator of pro-apoptotic gene expression during mammary gland involution Development, November 1, 2005; 132(21): 4675 - 4685. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Yamashita, K. Wakazono, T. Nomoto, Y. Tsujino, T. Kuramoto, and T. Ushijima Expression Quantitative Trait Loci Analysis of 13 Genes in the Rat Prostate Genetics, November 1, 2005; 171(3): 1231 - 1238. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. Currier, S. E. Solomon, E. G. Demicco, D. L. F. Chang, M. Farago, H. Ying, I. Dominguez, G. E. Sonenshein, R. D. Cardiff, Z.-X. J. Xiao, et al. Oncogenic Signaling Pathways Activated in DMBA-Induced Mouse Mammary Tumors Toxicol Pathol, October 1, 2005; 33(6): 726 - 737. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. M. Chan, X. Lu, F. M. Merchant, J.D. Iglehart, and P. L. Miron Gene expression profiling of NMU-induced rat mammary tumors: cross species comparison with human breast cancer Carcinogenesis, August 1, 2005; 26(8): 1343 - 1353. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Yamashita, S. Suzuki, T. Nomoto, Y. Kondo, K. Wakazono, Y. Tsujino, T. Sugimura, T. Shirai, Y. Homma, and T. Ushijima Linkage and Microarray Analyses of Susceptibility Genes in ACI/Seg Rats: A Model for Prostate Cancers in the Aged Cancer Res., April 1, 2005; 65(7): 2610 - 2616. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. Shan, M. Yu, and E. G. Snyderwine Gene expression profiling of chemically induced rat mammary gland cancer Carcinogenesis, February 1, 2005; 26(2): 503 - 509. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Nishida, S. Mine, T. Utsunomiya, H. Inoue, M. Okamoto, H. Udagawa, T. Hanai, and M. Mori Global Analysis of Altered Gene Expressions during the Process of Esophageal Squamous Cell Carcinogenesis in the Rat: A Study Combined with a Laser Microdissection and a cDNA Microarray Cancer Res., January 15, 2005; 65(2): 401 - 409. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Qiu, L. Shan, M. Yu, and E. G. Snyderwine Deregulation of the Cyclin D1/Cdk4 Retinoblastoma Pathway in Rat Mammary Gland Carcinomas Induced by the Food-derived Carcinogen 2-Amino-1- methyl-6-phenylimidazo[4,5-b]pyridine Cancer Res., September 15, 2003; 63(18): 5674 - 5678. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. El-Bayoumy, B. A. Narayanan, D. H. Desai, N. K. Narayanan, B. Pittman, S. G. Amin, J. Schwartz, and D. W. Nixon Elucidation of molecular targets of mammary cancer chemoprevention in the rat by organoselenium compounds using cDNA microarray Carcinogenesis, September 1, 2003; 24(9): 1505 - 1514. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Abe, S. Yamashita, T. Kuramoto, Y. Hirayama, T. Tsukamoto, T. Ohta, M. Tatematsu, M. Ohki, T. Takato, T. Sugimura, et al. Global expression analysis of N-methyl-N'-nitro-N-nitrosoguanidine-induced rat stomach carcinomas using oligonucleotide microarrays Carcinogenesis, May 1, 2003; 24(5): 861 - 867. [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 |