| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Clinical Research |
1 Genetics Branch, Center for Cancer Research and 2 Nutritional Epidemiology Branch, Division of Cancer Epidemiology and Genetics, National Cancer Institute (NCI); 3 Office of Clinical and Regulatory Affairs, National Center for Complementary and Alternative Medicine and 4 Genetic Epidemiology Branch, Division of Cancer Epidemiology and Genetics, NIH; 5 Office of Centers, Training, and Resources, NCI, NIH, Bethesda, Maryland; 6 Cancer Institute, Chinese Academy of Medical Sciences, Beijing, China; and 7 Mayo Clinic College of Medicine, Rochester, Minnesota
Requests for reprints: Philip R. Taylor, Division of Cancer Epidemiology and Genetics, National Cancer Institute, NIH, EPS, 6120 Executive Boulevard, Room 7006, Bethesda, MD 20892-7236. Phone: 301-594-2932; Fax: 301-402-4489; E-mail: ptaylor{at}mail.nih.gov.
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
The precise molecular targets associated with regression and progression of squamous dysplastic lesions in the esophagus have not been identified. The present study aimed to address this knowledge gap. We used microarray analysis to assess changes in gene expression before and after treatment in histopathologically confirmed normal mucosa taken from 29 individuals, a subset of the 117 subjects in the selenomethionine and placebo arms of the original Linxian chemoprevention trial cohort. This subset included 11 individuals whose dysplasia regressed, 13 who were stable, and 5 whose dysplasia progressed.
| Materials and Methods |
|---|
|
|
|---|
Baseline and end-of-trial esophagogastroduodenoscopy exams were done using Pentax videoendoscopy equipment (Pentax Precision Instrument Corp., Orangeburg, NY). Following initial inspection, the esophageal mucosa was sprayed with 1.2% Lugol's iodine solution, after which dysplastic lesions appear unstained relative to the surrounding normal mucosa (9). All unstained lesions
5 mm in diameter were biopsied. In addition, two adjacent biopsy samples were obtained from endoscopically normal mucosa located at least 10 mm away from any grossly visible lesion. One of these latter biopsies was snap frozen in liquid nitrogen, and the other was fixed and processed along with the lesion biopsies. The fixed biopsy samples were immediately placed in 10% neutral-buffered formalin, embedded in paraffin, sectioned in 5-µm thicknesses, and stained with H&E. The biopsy slides from the baseline and end-of-trial endoscopy exams for each patient were paired and reviewed together, with the dates of the exams masked. The slides were independently reviewed in a blinded fashion by two gastrointestinal pathologists who were unaware of the intervention group assignments. Cases with discrepant results were adjudicated in a blinded fashion by a third gastrointestinal pathologist. The study was approved by the Institutional Review Boards of the Cancer Institute, Chinese Academy of Medical Sciences (Beijing, China) and the U.S. National Cancer Institute (Bethesda, MD).
Sample selection for expression analysis. Baseline and end-of-trial frozen biopsies from endoscopically normal mucosa were available from 66 individuals from the selenomethionine and placebo arms of the chemoprevention trial. The adjacent fixed biopsies were histologically normal at both time points in 29 individuals, and the frozen biopsies from these cases were used for the microarray analyses.
Microarray experiments. Total RNA was isolated from the snap-frozen biopsies using the Qiagen (Valencia, CA) RNeasy Mini kit and a Micro-H Omni (Marietta, GA) homogenizer. RNA quality was examined using the RNA 6000 Nano assay on the Agilent (Palo Alto, CA) 2100 Bioanalyzer. The Affymetrix (Santa Clara, CA) small sample labeling protocol vII was used as specified by the manufacturer (10). In short, two cycles of linear amplification were done. cDNA was first synthesized from 1 µg total RNA using a first-strand and a second-strand protocol, with the T7-oligo(dT) promoter primer used for the first-strand synthesis. The double-stranded cDNA was purified by ethanol precipitation, amplified with the Ambion (Austin, TX) MEGAscript T7 kit, and then purified again using the RNeasy kit. cRNA (2 µg), quantified by spectrophotometer, was then used to resynthesize cDNA with a random primer for the first-strand synthesis and the T7-oligo(dT) promoter primer for the second-strand synthesis. The cDNA was purified by ethanol precipitation and biotin labeled with the ENZO BioArray High Yield RNA Transcript labeling kit purchased from Affymetrix. The labeled cRNA was then purified using the RNeasy columns. The cRNA concentration was measured by a spectrophotometer and fragmented.
The GeneChip Human Genome U133A platform was used, which consists of 22,000 oligonucleotide probe sets representing 18,400 transcripts and variants, including 14,500 known genes. The GeneChips were prehybridized using 1x hybridization buffer [100 mmol/L MES, 1 mol/L Na+, 20 mmol/L EDTA, 0.01% Tween 20 (final concentration)] at 65°C, with 60 rpm rotation for 10 minutes. The prehybridization solution was removed and replaced with denatured hybridization solution containing 0.06 µg/µL fragmented cRNA (15 µg in total), 50 pmol/L control oligonucleotide B2, 1.5, 5, 25, and 100 pmol/L, respectively, of 20x eukaryotic hybridization controls (bioB, bioC, bioD, and cre), and 0.1 mg/mL herring sperm DNA. The hybridizations were done in an oven at 65°C, rotating at 60 rpm for 16 hours. Protocol EukGe_WS2v5_450 was used to wash and stain the GeneChips on Affymetrix Fluidics Station 450. The arrays were washed with a series of nonstringent (wash A: 6x saline-sodium phosphate-EDTA, 0.01% Tween 20) and stringent (wash B: 100 mmol/L MES, 0.1 mol/L Na+, 0.01% Tween 20) wash buffers and stained with streptavidin phycoerythrin [1x MES stain buffer, 2 mg/mL acetylated bovine serum albumin (BSA), 10 µg/mL streptavidin phycoerythrin]. The signal was then amplified with an antibody solution (1x MES stain buffer, 2 mg/mL acetylated BSA, 0.1 mg/mL normal goat IgG, 3 µg/mL biotinylated antibody) and stained a second time with streptavidin phycoerythrin.
Microarray image acquisition and data analysis. The probe arrays were scanned in an Affymetrix GeneChip Operating Software Scanner 3500 equipped with an argon-ion laser at the excitation wavelength of 488 nm. Signal intensities were calculated from the cellular intensity files (.cel) files, generated by the GeneChip Operating System, using justRMA, a library of Bioconductor (version 1.9) supported by R (version 1.8.1). We also used database for annotation, visualization, and integrated discovery (DAVID) and NetAffx, genomic databases, for gene annotation (11).
Statistical analysis. We analyzed mRNA expression levels in normal mucosa samples from individuals categorized as regressors, stable, or progressors based on the change in their worst histologic lesions over a period of 10 months. For subjects whose most advanced histologic diagnosis at baseline was mild dysplasia, regression was defined as no evidence of dysplasia, whereas progression was defined as moderate dysplasia, severe dysplasia, or invasive carcinoma. For subjects whose most advanced histologic diagnosis at baseline was moderate dysplasia, regression was defined as no evidence of dysplasia or mild dysplasia, whereas progression was defined as severe dysplasia or invasive carcinoma.
We examined whether the change in gene expression between T0 and T10 was associated with regression, stability, or progression of esophageal squamous dysplasia. A standard t test, often associated with microarray analysis, could not be used to analyze the data generated from this study because there are three different response groups instead of the usual two. These three ordinal categories dictated the use of ordinal regression (12). The selenomethionine treatment effects were adjusted for by placing this covariate into the ordinal regression model. The direction of the ß coefficient for gene expression change indicates whether the change was associated with regression (negative) or progression (positive). The coefficients and the corresponding Ps generated from this model assumed a linear trend in relative gene expression. For example, for a gene with relatively increased expression in individuals whose lesions regressed, there is intermediate expression in individuals whose lesions were stable and lower expression in individuals whose lesions progressed. Conversely, genes associated with progression will also maintain this linear relationship in reverse (greater expression in individuals whose lesions progressed, intermediate in stable individuals, and lower expression in individuals whose lesions regressed).
We searched for genes associated with regression, stable, and progression using two different definitions of statistical significance. For the first definition, we used only a P cutoff of P < 0.001 (Tables 2 and 3). Our second definition was more stringent and required, in addition to a P < 0.001, a magnitude component (i.e., that the range of gene expression be at least 2-fold; Fig. 1 ).
|
|
|
| Results |
|---|
|
|
|---|
|
We analyzed these data on the log2 scale, so a 1-unit increase in the gene expression covariate was equivalent to a 2-fold increase in gene expression from baseline to 10 months. Additionally, the interpretation of coefficients in this model is slightly different than in most t test models. Although one can discuss underexpression and overexpression, for research purposes we preferred to think about overexpression in a gene trending toward regression or progression of esophageal squamous dysplasia. Figure 1 depicts the ß coefficients from the analyses in a continuum along the X axis, where zero means the three groups of individuals (those who had regression, stability in worst diagnosis, or progression) had similar changes in gene expression over the 10-month period. The more negative a coefficient, the stronger a 2-fold increase in the gene expression over the 10-month period pulls toward the regression end of the continuum. The more positive a coefficient, the stronger a 2-fold increase in the gene expression over the 10-month period pulls toward the progression end of the continuum. These coefficients are similar to point estimates, but they are not fold changes. Every value is based on what occurs to esophageal dysplasia status given a 2-fold increase in the RNA expression of a particular gene from the beginning to the end of the intervention period.
Using the Gene Ontology database, we categorized the 149 differentially expressed genes into known or probable functional categories. The most common categories included immune response (15 genes: FCER1A, CD58, HLA-DPA1, HLA-DRA, HLA-DQB1, CRP, ORMS, SERPINA7, SERPINA1, CNR2, NFATC4, NFRKB, MBP, INHBB, and CMKLR1), cell cycle (15 genes: RET, TGF-
, MET, CDC14B, CDC42, DLG3, EGR4, IL9, MAS1, NAPILI, CDC25C, SOCS1, APBB2, BLZF1, and NPPB), metabolism (15 genes: RDH11, AGPAT3, UGT2B28, GBA3, CHAT, HSDIIB1, ATP12A, DCT, ABCA1, UQCRB, ATP5J, PPARA, PNLIPRP2, TXNRD2, and ABHD6), calcium transport or calcium ion activity (10 genes: THBS1, TACSTD2, CADPS, RAMP3, CLCA2, CDH9, GRIN1, TRPC4, STC2, and ZZEF1), regulation of transcription (9 genes: ZBTB7A, ZNF480, KIAA0194, GTF2A1, PRDM10, ZNF1A1, NR4A3, ZNF42, and MYCBP2), signal transduction (7 genes: PNOC, EPHB1, BRAP, SEMA4F, GP6, NPPB, and PHEX), cytoskeleton and microtubules (5 genes: KIF5B, RDX, TUBB5, DST, and ASTN), nucleotide processing and biosynthesis (4 genes: MTAP, SLC29A1, ADCY2, and UMPS), G-coupled signaling (4 genes: F2RL2, GRPR5, BAI3, and RASL113), and apoptosis (3 genes: CASP2, BAX, and CIDEC).
We used the EASE software package to identify pathways overrepresented in the 149 differentially expressed genes. The 15 genes associated with the immune response pathway were greater than the expected number based on the total number of immune response genes on the microarray chip (P < 0.01). No other pathway was overrepresented in the 149 genes we found in our analysis. Further analysis of the immune response genes revealed that those associated with regression had functions consistent with up-regulation of immune response, including antigen presentation [HLA-DPA1, HLA-DRA, and HLA-DQB1 (14, 15)], protection of T cells from cell death [CD58 (16)], and T-cell activation [FCERA1 (17, 18); Table 4 ]. In contrast, genes associated with progression were involved in immunosuppression [CNR2 (1921) and NFATC4 (22, 23)], acute phase response and inflammation [CRP (24), ORMS (25, 26), SERPINA7, and SERPINA1 (27, 28)], or a decrease in the level of B and T cells [NFRKB (23, 29)]. Three genes identified by the EASE program [MBP (3032), INHBB, and CMKLR1 (33, 34)] had unknown functions but are speculated to play a role in immune response. All 5 of the immune genes associated with regression and 2 (CRP and MBP) of the 10 immune genes associated with progression also had gene expression ranges that exceeded 2-fold.
|
| Discussion |
|---|
|
|
|---|
Immune response was the only pathway with more differentially expressed genes than expected based on the total number of such genes examined on the microarray. Among differentially expressed immune response genes, those associated with regression involved immune stimulation, whereas those associated with progression related to immune suppression and inflammation. Differences in the underlying immune response between these two groups may be important in distinguishing mechanisms associated with regression and progression of preneoplastic lesions. Immune reactions may be of critical importance in the enhancement or impairment of tumor oncogenic capacity. Human neoplastic cells express tumor-associated antigens in ways that immunologically differentiate them from normal, and the immune system may be activated to effectively react against these tumor cells (36). Genes associated with antigen presentation, similar to those observed in regression here, produce a cascade of events that leads to activation and expansion of T cells as well as the local accumulation of immune cells. These events may result in the elimination of cells, tissue repair, and resolution of inflammation. In contrast to this focal immune response, progression of preneoplastic lesions is believed, in part, to be caused by the subversion of immune response and the appearance of a more generalized inflammatory reaction (3740). The most simplistic explanation for this phenotype is that defective antigen presentation in preneoplastic/tumor cells does not produce an adequate immune response (41).
Other groups of genes that were common in our 149 differentially expressed genes included genes related to the cell cycle (n = 15), metabolism (n = 15), and calcium transport or calcium ion activity (n = 10). Fifteen cell cyclerelated genes were differentially expressed in this analysis, including 3 genes (CDC14B, TGF-
, and RET), whose range of gene expression was at least 2-fold. Fourteen of these cell cycle genes were overexpressed in individuals whose lesions progressed. Although the exact function of many of these genes is unknown, some have been shown to play an important role in the cell cycle and/or have been implicated in the progression of cancer. TGF-
, for example, is often overexpressed in head and neck and other squamous cell carcinomas and is believed to be a potent mitogen involved in wound healing (42). Its overexpression is thought to be an early event, found in both dysplasia and normal mucosa surrounding dysplastic tissue (43, 44). Two other potent mitogens, MET and RET, which are tyrosine kinase proto-oncogenes, are also involved in a variety of cancers (45, 46).
Fifteen genes involved in endogenous or xenobiotic metabolism were found to be significantly associated with change in worst squamous diagnosis. Several of the significantly associated genes (AGPAT3, ATP5J, and ATP12A) play a role in ATP metabolism, and increased expression was associated with progression, which suggests a response to increased energy requirements.
Genes involved in calcium regulation and calcium signaling also seemed to be important in the progression of esophageal squamous dysplasia. We found significant associations with 10 calcium-related genes. A recent array study that compared ESCC tumor and normal squamous tissue found a large number of changes in calcium-binding or calcium-modulating genes (47). One notable gene associated with progression of squamous dysplasia in our study is tumor-associated calcium signal transducer 2 (TACSTD2; also known as TROP-2), which has previously been identified as a tumor antigen (48) that transduces calcium signals (49). The TROP-2 protein has been reported to be overexpressed in the sera of ESCC patients (50). The increased expression of this gene in the normal tissue of subjects without frank tumors suggests that this gene and protein may be important at an early stage of tumor formation.
The results reported here represent a unique effort in carcinogenesis research to bridge the gap in our understanding between morphology and molecular biology. A cancer prevention trial was used, in which the primary end point, regression, or progression of a surrogate cancer end point (intraepithelial neoplasia) was evaluated in relation to change in gene expression profiles from the beginning to the end of the intervention (51). Because the trial found that supplementation with selenium improved the morphologic profile among participants who started the trial with mild dysplasia, we were able to examine molecular changes (changes in RNA expression profiles) coincidental with morphologic changes. All the genes identified here are of interest because of their potential role in the natural history of esophageal squamous dysplasia. In addition, as part of future cancer prevention strategies, the genes associated with regression represent potential targets for gene or pathway enhancement and the genes associated with progression are potential targets for pathway disruption. Similarly, both regression- and progression-associated genes merit consideration as early detection markers.
There are several caveats that should be considered in evaluating this study. The tissues analyzed here were only from histologically normal esophageal tissue in persons with dysplasia elsewhere in their esophagus. We do not know if gene expression results would be the same if the premalignant lesions themselves had been analyzed. Although there were several strong design features built into this study (i.e., the end point was the paired before-versus-after supplementation morphology difference in a single individual, and the exposure was the paired before-versus-after supplementation gene expression profile in the same person), the study was only modest in size. It should thus be considered an exploratory effort for which the findings must be replicated elsewhere using different study populations and alternative study designs.
In conclusion, we examined the association between gene expression changes over a 10-month period in normal esophageal mucosal biopsy specimens and regression, stability, and progression in worst esophageal squamous lesions. We found many significant associations, including numerous genes involved in immune response.
| Acknowledgments |
|---|
| Footnotes |
|---|
Received 2/20/06. Revised 4/25/06. Accepted 5/ 4/06.
| References |
|---|
|
|
|---|
and
subunits of the mast cell receptor for immunoglobulin E are located on human chromosome band 1q23. Immunogenetics 1990;32:1836.[Medline]
production by human mononuclear cells. Int Immunopharmacol 2002;2:53743.[CrossRef][Medline]
B factorsthe distant relatives. Int J Biochem Cell Biol 2004;36:116670.[CrossRef][Medline]
1-antitrypsin: a paradigm for other SERPINs. Int J Biochem Cell Biol 1997;29:150111.[CrossRef][Medline]
qter, and trisomy 8q24.3
qter. Am J Med Genet 2002;108:3226.[Medline]
and EGFR expression in oral leukoplakia and oral submucous fibrosis by quantitative immunohistochemistry. Oncology 2001;61:28492.[Medline]
and EGFR mRNA in laryngeal carcinoma tissues. Chin Med J (Engl) 1999;112:108892.[Medline]This article has been cited by other articles:
![]() |
S. R. Mallery, J. C. Zwick, P. Pei, M. Tong, P. E. Larsen, B. S. Shumway, B. Lu, H. W. Fields, R. J. Mumper, and G. D. Stoner Topical Application of a Bioadhesive Black Raspberry Gel Modulates Gene Expression and Reduces Cyclooxygenase 2 Protein in Human Premalignant Oral Lesions Cancer Res., June 15, 2008; 68(12): 4945 - 4957. [Abstract] [Full Text] [PDF] |
||||
![]() |
V. Pagmantidis, C. Meplan, E. M van Schothorst, J. Keijer, and J. E Hesketh Supplementation of healthy volunteers with nutritionally relevant amounts of selenium increases the expression of lymphocyte protein biosynthesis genes Am. J. Clinical Nutrition, January 1, 2008; 87(1): 181 - 189. [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 |