EXPERIMENTAL STUDY

Analysis by cDNA microarrays of gene expression patterns of human adrenocortical tumors

E P Slater, S M Diehl, P Langer, B Samans1, A Ramaswamy2, A Zielke and D K Bartsch

Department of Surgery, 1Institute of Medical Biometry and Epidemiology and 2 Department of Pathology, Philipps-University Marburg, Baldingerstrasse, 35033 Marburg, Germany

(Correspondence should be addressed to E P Slater; Email: slater@med.uni-marburg.de)

Abstract

Objectives: Adrenocortical carcinoma (ACC) is a rare malignant neoplasm with extremely poor prog- nosis. The molecular mechanisms of adrenocortical tumorigenesis are still not well understood. The comparative analysis by cDNA microarrays of gene-expression patterns of benign and malignant adrenocortical tumors allows us to identify new tumor-suppressor genes and proto-oncogenes under- lying adrenocortical tumorigenesis.

Design and methods: Total RNA from fresh-frozen tissue of 10 ACC and 10 benign adrenocortical ade- nomas was isolated after histologic confirmation of neoplastic cellularity of at least 85%. The refer- ence consisted of pooled RNA of 10 normal adrenal cortex samples. Amplified RNA of tumor and reference was used to synthesize Cy3- and Cy5-fluorescently labeled cDNA in a flip-color technique. D-chips containing 11 540 DNA spots were hybridized and scanned and the images were analyzed by ImaGene 3.0 software.

Results: The comparative analysis of gene expression revealed many genes with more than fourfold expression difference between ACC and normal tissue (42 genes), cortical adenoma and normal tissue (11 genes), and ACC and cortical adenoma (21 genes) respectively. As confirmed by real- time PCR, the IGF2 gene was significantly upregulated in ACCs versus cortical adenomas and normal cortical tissue. Genes that were downregulated in adrenocortical tumors included chromogra- nin B and early growth response factor 1.

Conclusions: Comprehensive expression profiling of adrenocortical tumors by the cDNA microarray technique is a very powerful tool to elucidate the molecular steps associated with the tumorigenesis of these ill-defined neoplasms. To evaluate the role of identified genes, further detailed analyses, including correlation with clinical data, are required.

European Journal of Endocrinology 154 587-598

Introduction

Adrenal masses are a common disorder, affecting 3-7% of the population. Most turn out to be benign adrenocortical adenomas, which may be functional or nonfunctional. Much more rarely, these masses represent primary adrenal carcinoma (1). Adrenocorti- cal carcinoma (ACC) is a highly malignant tumor with an incidence of ~ 1 per 1.7 million inhabitants per year in the West. Although ACC is rare, its highly aggressive behavior and 5-year survival rate of only 10-20% urgently require the decoding of its molecular basis to develop new strategies for diagnosis and treatment (2-6). The genetic background of adrenocortical tumorigenesis is poorly characterized. In other endo- crine tissues, such as thyroid, there is conclusive evi- dence that hyperplasia and adenomas can precede

cancer. In the adrenal, there are clinical cases of either hyperplasia or adenoma associated with later development of cancer. However, only a few studies have attempted to characterize this process on a mol- ecular basis (7). Although it is unclear whether there is an adenoma-carcinoma sequence, common patterns seen in adenomas and carcinomas, and the accumu- lation of chromosomal imbalances with tumor pro- gression support the existence of an adenoma- carcinoma sequence (8). X-chromosome inactivation analysis has shown that ACCs are of monoclonal origin, whereas benign adenomas may be monoclonal or polyclonal (9-11). The evidence gathered so far shows that the transition from adrenal adenoma to car- cinoma involves a monoclonal proliferation of cells that, among other yet to be characterized defects, have undergone chromosomal duplication at the

11p15.5 locus, leading to overexpression of the insulin- like growth factor (IGF)2 gene and abrogation of expression of the CDKN1C and H19 genes (12, 13). TP53 has been shown to be involved in progression to carcinoma in a subset of patients, and it has been suggested that the frequency of adrenocorticotropic hormone (ACTH) receptor deletion might also be involved (1). Other key oncogenes and tumor suppres- sor genes remain to be identified. However, a recent study has reported that chromosomal loci 1p, 2p16, 11q13 and 17p may harbor potential tumor-suppres- sor genes, and chromosomes 4, 5 and 12 potential oncogenes associated with adrenal tumorigenesis (1). Therefore, detailed analysis of the genes involved is highly desirable. The development of the cDNA micro- array technique offers the opportunity to analyze a large number of genes, and allows comparative analysis of gene-expression profiling in benign and malignant adrenocortical tumors and the identification of tumor- suppressor genes and proto-oncogenes associated with the initiation and progression of adrenocortical tumors. Findings from these analyses might clarify adrenocortical tumorigenesis and lead to the establish- ment of new diagnostic and prognostic markers as well as the characterization of novel strategies for treatment. Therefore, we have performed a comprehensive and representative analysis of neoplastic and nonneoplastic adrenocortical tissue samples. We were particularly interested in evaluating the differential profile of 11 500 genes with established importance for develop- ment and progression of malignant diseases (1), to assess whether combinations of genes can predict malignant tumors (2) and to validate the results by real time RT-PCR of selected candidate genes (3).

Materials and methods

Normal and tumor samples

The adrenocortical tissues analyzed in this study were obtained from the collection of fresh frozen adrenal tissue of the Department of Surgery, Philipps-University Marburg, Germany, collected between 1996 and 2003. For the purpose of this study, 10 ACCs and 10 adrenocor- tical adenomas (four from Conn’s syndrome, four from Cushing’s disease and two nonfunctional adenomas), as well as 10 nonneoplastic adrenal cortical tissue samples, were evaluated. The classification was deter- mined by both conventional histologic methods (14) and Weiss score (15) where adenomas met fewer than four and carcinomas at least four of the criteria. The ethics committee of the university approved this study, and all patients participating in the study consented to sampling. Frozen tumor samples were formalin-fixed and embedded in paraffin, and sections were evaluated by hematoxylin and eosin (HE) staining regarding diag- nosis and neoplastic cellularity. Only tumor samples with neoplasticity of greater than 85% were included

in the analysis. Control samples consisted of 10 histologically confirmed normal adrenal cortices.

RNA isolation was performed as follows. Frozen con- trol and tumor tissue samples were dissected by the pathologist (A M) and homogenized in the presence of TRIzol Reagent (Invitrogen, Karlsruhe, Germany) by the manufacturer’s protocol. Total RNA was then further purified by digestion with DNase I and recovery of RNA with the RNeasy kit (Qiagen, Hilden, Germany) by the supplier’s protocol. To determine the integrity of RNA, standard RT-PCR for the amplification of 17x- hydroxylase and ß-actin was performed on 25 ng RNA with Qiagen’s OneStep RT-PCR kit by the manu- facturer’s protocol. The primer sequences were as fol- lows:

· Cyp17_forward TCTCTTGCTGCTTACCCTAG; Cyp17_reverse TCAAGGAGATGACATTGGTT (GenBank accession no. NM_000102).

· ß-actin_forward GATGATGATATCGCCGCGCTCG- TCGTC; ß-actin_reverse GTGCCTCAGGGCAGCG- GAACCGCTCA (GenBank accession no. M10277).

The reaction mixture was incubated at 50 ℃ for 30 min and 95 ℃ for 15 min, followed by 30 cycles of standard PCR (1-min denaturation at 95℃, 1-min annealing at 58 °℃ and 1-min extension at 72°℃). PCR products were visualized by ethidium bromide staining of PAGE. Only samples showing expression of Cyp17 as well as ß-actin were included in analysis. RNA (2 µg) from 10 ACC and 10 adenomas (eight functional and two non- functional adenomas) was amplified with the Messa- geAmp aRNA Kit (Ambion, Huntingdon, UK). For the reference, the 10 control samples were amplified and then pooled. Amino allyl-cDNA was synthesized with 2 µg aRNA and then labeled and purified with the CyScribe Post-Labeling Kit (Amersham Biosciences, Frei- burg, Germany). Samples were fluorescently labeled with Cy3 and Cy5 by the flip-color technique.

For gene-expression profiling, the reference and tumor samples were mixed, denatured and then hybridized to microarrays for 16 h at 56 ℃ and washed at a stringency of 0.1% SSC and 0.1% SDS. The microarray contains 11 540 DNA spots; detailed protocols and data descrip- tion of the chip are available from the website: www.im- t.uni-marburg.de. Each experiment was performed as a sandwich hybridization with two arrays. Spot intensities were extracted from a scanned image with ImaGene 3.0 Software (BioDiscovery, Los Angeles, CA, USA). For each spot, median signal and background intensities for both channels were obtained. To account for spot differences, the background-corrected ratios of the two channels were calculated and log2-transformed. To balance the fluorescence intensities for the two dyes as well as to allow comparison of expression levels across exper- iments, the raw data were standardized. We used a spatial and intensity-dependent standardization (like Yang et al. (16)) to correct for inherent bias on

each chip (the lowest scatter-plot). As each gene was measured twice in the sandwich hybridization, mean log-ratios M were calculated from replicates. If gene repli- cates differed more than the maximum of threefold and 75% of the calculated average log-ratio, or the back- ground intensity was higher than signal intensity, the spot was excluded on that array. Differentially expressed genes were selected by a fold-change difference of at least 2 and an absolute value of the t-statistic of 1.96. Prior to the cluster analysis, the expression profile of each gene was centered by subtracting the mean observed value. Average linkage hierarchic clustering was then performed for genes as well as for chips with the Euclidean distance metric as implemented in the pro- gram Genesis (17).

The microarray results were validated for three candi- date genes known to be expressed in adrenocortical tissue, including chromogranin B (CgB), early growth response gene 1 (Egr-1), and IGF2 by real-time RT- PCR analysis. For validation, 10 µg total RNA were reverse transcribed with Superscript II reverse transcrip- tase (Invitrogen) and an oligo dT15 primer, according to the manufacturer’s instructions and previously pub- lished methods (18). Real-time PCR with the LightCycler System (Roche, Mannheim, Germany) was performed in a reaction mixture of 20 ul using the QuantiTect SYBR Green PCR Kit according to the manufacturer (Qiagen). Primers designed for analysis were as follows:

· Egr1 (GenBank accession no. NM001964) for- ward CGAGCAGCCCTACGAGCACCTGAC and reverse TGCGCAGCTCAGGGGTGGGCTCTG.

· IGF2 (GenBank accession no. BC000531) for- ward CCGTGCTTCCGGACAACT and reverse GGACTGCTTCCAGGTGTCATATT.

· CgB (GenBank accession no. AL035461) forward TGCCAGTGGATAACAGGAAC and reverse TCTT- CAGGACTTGGCGGCA.

· GAPDH forward CGTCTTCACCACCATGGAGA and reverse CGGCCATCACGCCACAGTTT.

The cycle threshold values for each gene were analyzed relative to those for GAPDH.

Results

For this study, 10 adrenocortical cancers, as established by histology, and 10 adenomas, of which eight were functional, including four aldosterone-producing ade- nomas, four cortisol-producing adenomas, and two nonfunctional adenomas (incidentalomas), were selected from the tissue bank of the Department of Sur- gery, Philipps-University Marburg. The clinicopatholo- gic data of the adrenocortical tumors analyzed are summarized in Table 1. Normal adrenal cortex, adreno- cortical adenoma and ACC were analyzed histologically

Table 1 Pathologic and clinical aspects of adrenocortical cases used for profiling studies.
DesignationTissue typeAgeGenderStage*SyndromeMetastatic sites **Follow-up
A1ACC50FIINoneNoneDOD
A2ACC46FINoneNoneNED
A3ACC59FIIINoneNoneNED
A4ACC43FIVCushingLungDOD
A5ACC15MIICushingNoneDOD
A6ACC74FIIINoneNoneDURC
A7ACC52MIIINoneNoneDOD
A8ACC51FIIICushing, AGSNoneDOD
A9ACC37FIIConnNoneNED
A10ACC37FIIICushingNoneAWD
A11ACC84FIINoneNoneNED
A12ACC43MIINoneNoneDOD
A15ACC17FIIICushing, AGSNoneNED
T1Adenoma46MN.A.ConnNED
T2Adenoma45FN.A.ConnNED
T3Adenoma47FN.A.ConnNED
T4Adenoma54MN.A.ConnNED
T5Adenoma65MN.A.CushNED
T6Adenoma27FN.A.CushNED
T7Adenoma47FN.A.CushNED
T8Adenoma45FN.A.CushNED
T9Adenoma65FN.A.IncNED
T10Adenoma52FN.A.IncNED
T11Adenoma64MN.A.IncNED
T12Adenoma73FN.A.IncNED

*Tumors graded according to UICC.

** At time of operation.

Cush: Cushing’s disease; Inc: Incidentaloma; N.A .: not applicable; Conn: Conn’s syndrome; AGS: adrenogenital syndrome; DOD: death of disease; NED: no evidence of disease; AWD: alive with disease; DURC: death of unrelated cause.

to confirm tissue origin, lack of necrosis and, for tumors, neoplasticity of greater than 85%, before pro- ceeding with nucleic acid purification. The normal tissue displayed regularly shaped nuclei and fat con- tent. Fig. 1 shows representative examples of HE stain- ing of tissue samples used for the RNA purification. The adrenal origin of the samples was additionally ensured by the RNA expression of 17a-hydroxylase (Cyp17), in addition to ß-actin, by RT-PCR. All samples chosen for further analysis were positive for both Cyp17 and

Figure 1 HE staining of normal adrenal cortex (A), adenoma (B) and adrenocortical carcinoma (C). x 100.

A

B

C

ß-actin (data not shown). Microarray analysis demon- strated that ACCs were more dissimilar to normal adrenal than adenomas. The adenomas were more closely related to each other and to normal adrenal; this is not an entirely unexpected result given the histo- logic similarity of the tissues. In particular, 40 differen- tially expressed genes were found in adenomas in comparison to normal adrenal (Table 2), and 144 dif- ferentially expressed genes were detected in ACCs (Table 3). As shown in Fig. 2 and Table 4, more than 60 genes were found with at least a threefold change in mRNA levels. Both up- and downregulated genes were identified (Table 4). These genes, with at least threefold differences in mRNA expression, were then subjected to cluster analysis. Transcriptional profiles, which distinguish between benign and malignant adrenocortical tumors, identified several differentially expressed transcripts as demonstrated by cluster analysis (Fig. 2). The sample dendrogram revealed the similarities among the adenomas (T) and the carci- nomas (A) respectively, and clearly distinguished the two groups. Two major clusters of genes were differen- tially expressed in the carcinomas compared with the adenomas: genes that were expressed at a higher level in the carcinomas and those that were expressed at a lower level in the carcinomas. The former include the gene for IGF-2 and potential oncogenes. The latter rep- resent potential tumor-suppressor genes.

Some of the differentially expressed genes were chosen for further analysis to ensure the validity of the microarray results. As expected from the literature (1, 12, 19), the IGF-2 gene expression was 3-7-fold greater in ACCs than in adenomas. Real-time RT-PCR confirmed this result where the cycle threshold crossing point (CT) for the ACC samples was 3-5 values less, representing a maximal increase in expression of 8-32-fold in comparison to the control or adenoma samples (Fig. 3).

As one of our main interests is the identification of potential tumor-suppressor genes whose investigation may contribute to the understanding of the pathome- chanism of ACC, we decided to concentrate in our initial analysis on clearly downregulated genes. The tissue-specifically expressed gene, chromogranin B (CgB), was found to be downregulated in both adenomas (28-fold) (Table 2) and ACCs (13-fold) (Table 3). The loss of a tissue-specifically expressed gene could reflect dedifferentiation of tumor tissue. This decrease in expression was confirmed by real- time RT-PCR, which demonstrated an increase in the CT value from 4 to 5, representing at most a 16-32- fold decrease in expression, for the adenomas analyzed, and from 3 to 4, or at most a 8-16-fold decrease for the ACCs.

The Egr-1 gene was downregulated in ACC in com- parison to normal adrenal by eightfold (Fig. 3). This finding was also confirmed by real-time RT-PCR, where the CT values for the ACCs were increased by

Table 2 List of genes found to be differentially regulated in adenomas relative to normal adrenal cortical tissue in microarray analysis.
Gene nameAccession no.Fold regulationUp/down
Protein tyrosine phosphatase, receptor type, f polypeptide (PTPRF), interacting protein (liprin), alpha 2H088504.42up
Chemokine (C-C motif) ligand 2AA4251023.94up
Activity-regulated cytoskeleton-associated proteinH861173.84up
Hypothetical protein FLJ10052R548223.76up
Purkinje cell protein 4AA4528263.36up
Mucolipin 3AA1717183.29up
Neuronal pentraxin IIAA6830412.40up
Cytoplasmic FMR1 interacting protein 2H120442.38up
B-cell CLL/lymphoma 11B (zinc finger protein)H097482.27up
Chromogranin B (secretogranin 1)W3776928.46down
Myxovirus (influenza virus) resistance 1, interferon-inducible protein p78 (mouse)AA45704212.08down
ESTH114539.86down
yz80b09.s1 Soares_multiple_sclerosis_2NbHMSP Homo sapiens cDNA clone IMAGE:289337 3'N926467.60down
Chromosome 1 open reading frame 29AA4101887.54down
Insulin-like growth factor binding protein 6AA4794287.35down
Brain expressed, X-linked 1W605824.70down
Contactin 1H193154.46down
Dickkopf homolog 3 (Xenopus laevis)AA4259474.33down
Sjogren syndrome antigen B (autoantigen La)H294854.08down
Chemokine (C-C motif) ligand 2AA9385633.99down
Cathepsin HAA4872313.93down
Sodium channel, voltage-gated, type III, betaR539303.82down
Insulin-like growth factor binding protein 6AA4787243.71down
Heat shock 105 kDa/110 kDa protein 1AA4851513.71down
Homo sapiens cDNA FLJ38885 fis, clone MESAN2017417T890943.63down
Sodium channel, voltage-gated, type III, betaAA1348243.55down
Solute carrier family 40 (iron-regulated transporter), member 1T52564 | T572353.43down
Ubiquitin carboxyl-terminal esterase L1 (ubiquitin thiolesterase)AA6704383.35down
KIAA1576 proteinAA6093483.34down
Protein tyrosine phosphatase, receptor type, N polypeptide 2AA4645903.30down
Sal-like 2 (Drosophila)H232543.29down
Connective tissue growth factorAA5987943.16down
Complexin 2H099663.10down
Integral membrane protein 2CAA0342133.08down
Hypothetical protein FLJ39155R081413.07down
Calmodulin 1 (phosphorylase kinase, delta)R76554 | R762773.02down
Homo sapiens cDNA FLJ39226 fis, clone OCBBF2007232.H090782.99down
Aldo-keto reductase family 1, member C1 (dihydrodiol dehydrogenase 1; 20-alpha (3-alpha)-hydroxysteroid dehydrogenase)R931242.95down
gr04f12.s1 soares fetal liver spleen 1NFLS Homo sapiens cDNA clone IMAGE: 204335H599162.92down
Nucleolar protein 4AA4300332.87down

4, representing a maximal decrease in expression of up to 16-fold, Thus, for each of these candidate genes, the differential gene expression was confirmed by real-time RT-PCR, and the findings were comparable to those of the microarray analysis (Fig. 3).

Discussion

DNA microarray technology allows comprehensive examination of the transcriptional profile of tumors. It is rapidly being applied to various problems in pathology and oncology, such as tumor classification, and it is a useful gene discovery tool to complement other similar technologies. In this study, we used DNA microarrays

to generate transcriptional profiles of benign and malig- nant adrenocortical tumors with various hormonal secretion profiles (Table 1). We demonstrated for the first time that these profiles can distinguish normal and benign tissues from malignant tumors and identify differentially expressed genes that may help to explain the pathogenesis of the disease and have diagnostic and therapeutic implications. Surprisingly, only 11 genes were found to be differentially expressed more than fourfold in comparing adenomas to normal tissue (Table 2). A comparison of the ACCs to normal tissue resulted in a list of 42 genes that were differentially expressed more than fourfold (Table 3), and the com- parison ACCs to adenomas resulted in 21 genes that were more than fourfold differentially expressed.

Table 3 List of genes found to be differentially regulated in ACC relative to normal adrenal cortical tissue in the microarray analysis.
Gene nameAccession no.Fold regulationUp/down
Dopachrome tautomerase (dopachrome delta-isomerase. tyrosine-related protein 2)AA4785536.14up
Human neuropeptide Y receptor Y1 (NPYY1)R438176.11up
Sodium channel, voltage-gated, type III, betaR539305.78up
Fibronectin 1R626124.43up
Chemokine (C-C motif) ligand 2AA9385634.30up
Sjogren syndrome antigen B (autoantigen La)H294854.29up
Fibronectin 1R626124.28up
Insulin-like growth factor 2 (somatomedin A)N545964.16up
GA binding protein transcription factor. alpha subunit 60 kDaH962414.12up
Hypothetical protein MGC5306H472574.07up
ESTN597723.93up
Hydroxy-delta-5-steroid dehydrogenase. 3 betaand steroid delta-isomerase 1R688033.84up
Dickkopf homolog 3 (Xenopus laevis)AA4259473.77up
Aldehyde dehydrogenase 1 family. member A3AA4552353.69up
Clone IMAGE:120162 mRNA sequenceT952743.67up
CDC28 protein kinase regulatory subunit 2.AA3978133.51up
Chemokine (C-C motif) ligand 2AA4251023.27up
ESTH163893.24up
Insulin-like growth factor binding protein 3AA5986013.22up
Immediate early response 3AA4577053.21up
Dentatorubral-pallidoluysian atrophy (atrophin-1)H086433.20up
Ectodermal-neural cortex (with BTB-like domain)H721223.18up
Clone DKFZp434E235H225633.18up
yq51h12.s1 Soares fetal liver spleen 1NFLS cDNA clone IMAGE:199367 3' ..R956913.11up
Dopa decarboxylase (aromatic L-amino acid decarboxylase)AA7026403.11up
Glutaminyl-peptide cyclotransferase (glutaminyl cyclase)AA2821343.03up
Fibulin 2AA4528403.02up
c-src tyrosine kinaseAA0797752.98up
Flotillin 1AA4881752.92up
Fibroblast growth factor 14AA4000472.89up
RWD domain containing 1AA4874992.89up
Hypothetical protein FLJ37078R442922.84up
Sodium channel. nonvoltage-gated 1 alphaAA4589822.84up
Troponin I, skeletal, fastAA1813342.82up
Membrane-associated tyrosine- and threonine-specific cdc2-inhibitory kinaseAA4780662.73up
Insulin-like growth factor 2 (somatomedin A)N545962.63up
ParathymosinR115262.53up
Cadherin 13. H-cadherin (heart)R417872.50up
Adrenergic. alpha-2A -. receptorT486922.49up
Glucocorticoid receptor DNA binding factor 1N722762.42up
CDC28 protein kinase regulatory subunit 2AA0100652.41up
Signal recognition particle receptor, B subunitR956932.37up
Cytochrome c oxidase subunit VIa polypeptide 1AA4822432.35up
Ubiquitin-like 1 (sentrin)AA4886262.33up
Collagen, type V, alpha 1R756352.30up
Fibroblast growth factor receptor 1 (fms-related tyrosine kinase 2. Pfeiffer syndrome)R548462.27up
Splicing factor. arginine/serine-rich 2. interacting proteinR911712.25up
ESTT975922.23up
Lymphocyte-specific protein 1T906322.18up
Likely ortholog of rat vacuole membrane protein 1AA4853732.14up
Dishevelled. dsh homolog 2 (Drosophila)R383252.13up
Kruppel-like zinc finger protein GLIS2R438262.12up
Dehydrogenase/reductase (SDR family) member 7H871442.12up
Insulin-like growth factor 2 (somatomedin A)N746232.11up
Hypothetical protein BC013949W697412.11up
Dipeptidylpeptidase 4 (CD26. adenosine deaminase complexing protein 2)W702342.11up
NucleoredoxinT642162.10up
Zinc finger protein 161AA2326472.10up
Hypothetical protein from clone 643T534042.09up
BTG family. member 3N524962.07up
Enoyl Coenzyme A hydratase domain containing 1AA1735732.06up
Glyceraldehyde-3-phosphate dehydrogenaseh529502.03up
ESTH484672.03up
Prefoldin 5AA4464532.03up
Table 3 Continued
Gene nameAccession no.Fold regulationUp/down
Integrin, alpha 7T609262.02up
ThymopoietinT639802.02up
Myxovirus (influenza virus) resistance 1, interferon-inducible protein p78 (mouse)AA45704213.60down
Chromogranin B (secretogranin 1)W3776912.72down
Actin binding LIM protein 1AA40660111.43down
ESTH114539.23down
Ferredoxin 1AA1873498.96down
FBJ murine osteosarcoma viral oncogene homolog BT619488.32down
yz80b09.s1 Soares_multiple_sclerosis_2NbHMSP Homo sapiens cDNA cloneN926468.17down
IMAGE:289337 3'.
Early growth response 1AA4865337.54down
ESTN689937.54down
Fatty acid binding protein 5 (psoriasis-associated)T600757.49down
Aldehyde dehydrogenase 1 family. member A1AA6641017.35down
zb50h07.s1 Soares_fetal_lung_NbHL19W Homo sapiens cDNA cloneN936867.08down
IMAGE:307069 3'
Chromosome 1 open reading frame 29AA4101886.96down
Contactin 1H193156.79down
Microsomal glutathione S-transferase 1AA4959366.58down
Glucose phosphate isomeraseAA4011116.54down
Regucalcin (senescence marker protein-30)H051406.28down
Interferon-induced protein with tetratricopeptide repeats 1AA4896405.04down
Chemokine (C-C motif) ligand 15R966265.03down
Solute carrier family 16 (monocarboxylic acid transporters) member 9W164244.87down
DKFZP586B1621 proteinH205434.85down
High mobility group AT-hook 1AI0424044.80down
Serine (or cysteine) proteinase inhibitor. clade G (C1 inhibitor). member 1.AA4814384.64down
(angioedema. hereditary)
Mitogen-activated protein kinase kinase kinase 5AA1508284.52down
Protein tyrosine phosphatase, receptor type, f polypeptide (PTPRF), interacting protein (liprin), alpha 2H088504.48down
Vascular cell adhesion molecule 1H070724.47down
Ubiquitin carboxyl-terminal esterase L1 (ubiquitin thiolesterase)AA6704384.31down
v-jun sarcoma virus 17 oncogene homolog (avian)W961554.30down
Activity-regulated cytoskeleton-associated proteinH861174.26down
Heat shock 105 kDa/110 kDa protein 1AA4851514.21down
Glutathione S-transferase theta 1H998134.21down
Insulin-like growth factor binding protein 6AA4794284.08down
Atp-binding cassette. sub-family B (MDR/TAP). member 1AA4598243.95down
Insulin-like growth factor binding protein 6AA4787243.89down
Zinc finger protein 145 (Kruppel-like, expressed in promyelocytic leukemia)AA1016323.88down
Fibulin 1AA1348713.85down
Chromosome 9 open reading frame 97AA2353883.84down
KIAA1268 proteinT649563.82down
Solute carrier family 26 (sulfate transporter), member 2W152633.80down
Major histocompatibility complex, class II, DP beta 1AA4865323.80down
CD163 antigenAA4016933.79down
Glutathione peroxidase 3 (plasma)AA6641803.79down
MAX proteinT894963.76down
Insulin-like growth factor 1 receptorH133003.75down
Arachidonate 5-lipoxygenase-activating proteinT496523.74down
Interferon-induced transmembrane protein 1 (9-27)AA0583233.66down
Chromosome 10 open reading frame 10N552693.65down
Scavenger receptor class B, member 1AA4438993.61down
Connective tissue growth factorAA5987943.59down
ESTT911003.57down
Discoidin domain receptor family, member 2AA2437493.56down
Homo sapiens cDNA FLJ40165 fis. clone TESTI2015962.H228543.50down
Solute carrier family 40 (iron-regulated transporter), member 1T52564 | T572353.44down
Hypothetical protein FLJ39155R081413.41down
Homo sapiens cDNA FLJ38885 fis. clone MESAN2017417.T890943.40down
Plasminogen-likeT675493.38down
Periostin, osteoblast specific factorAA5986533.37down
ATP-binding cassette, subfamily C (CFTR/MRP), member 3AA4298953.37down
Chemokine-like factor superfamily 3AA4865613.36down
ALEX3 proteinN544563.33down
Table 3 Continued
Gene nameAccession no.Fold regulationUp/down
Interferon, alpha-inducible protein (clone IFI-6-16)AA4484783.29down
Complement component 7AA5984783.27down
KIAA1576 proteinAA6093483.23down
Sodium channel, voltage-gated, type III, betaAA1348243.11down
Retinoic acid receptor responder (tazarotene induced) 2AA4820673.07down
Complexin 2H099663.05down
Glutaredoxin (thioltransferase)AA2911633.03down
ATP-binding cassette. subfamily B (MDR/TAP), member 1AA1359583.03down
AE binding protein 1AA4904622.99down
Hypothetical protein FLJ20637AA4874622.98down
Fibronectin type III domain containing 5N549012.92down
Glutathione S-transferase A4AA1523472.91down
G protein-coupled receptor 116T639712.90down
Calmodulin 1 (phosphorylase kinase, delta)R76554 | R762772.88down
Major histocompatibility complex; class Il; DR beta 5AA4857392.86down
Scm-like with four mbt domains 1AA4005122.26down
Notch homolog 3 (Drosophila)T635112.22down
Peroxisomal lon proteaseT671382.04down

The gene for chromogranin B (CgB) was found to be downregulated in both adenomas (28-fold) (Table 2) and ACCs (13-fold) (Table 3). This finding was confirmed by real-time PCR (Fig. 3). Chromogranins are represen- tative proteins contained in endocrine cells of various organs, including some ductal cells of the breast. Hom- ology between the BRCA1 protein (1214-1223) and the chromogranins has been detected, suggesting that chromogranin may play the role of a tumor suppressor, like BRCA1 (20). Our results suggest that loss of CgB may be an early event in adrenocortical tumorigenesis. In patients with lymph-node-negative primary invasive ductal breast carcinoma, CgB-negative tumors demon- strated a significantly poorer prognosis than in patients with CgB-positive tumors. In univariate analysis, a sig- nificantly increased risk of disease progression and death was present in patients with CgA-poor and CgB- poor tumors respectively (21). It was concluded that the CgB immunostaining pattern of the primary tumor can distinguish patients with increased risk of death in patients with sporadic medullary thyroid carcinoma. The problem of whether CgB represents a prognostic factor in ACCs requires a large-scale study.

Another interesting potential tumor suppressor is the Egr-1 gene. As validated by real-time PCR, it was down- regulated in ACC in comparison to normal adrenal cortex tissue and adrenal adenoma by eightfold and threefold respectively (Fig. 3). Egr-1 is a transcription factor that has previously been suggested to be a master switch regulating inflammatory parameters (22). Egr-1 is also known as nerve growth factor induced-A (NGFI-A), Krox-24, ZIF268, ETR103 and TIS8. It is a phosphorylated zinc-finger-containing transcription factor often associated with cell growth stimulation (23). The gene for Egr-1 is located on the q31.1 ‘cytokine cluster’ region of chromosome 5 in man, and is rapidly and transiently induced by a

large number of stressful stimuli as well as growth fac- tors and cytokines. In fact, Egr-1 has been proposed as a master switch of gene expression underlying coordi- nated responses to various types of injury (23).

Apoptosis-regulating genes, such as TP53 and ANp73a, are known to increase the expression of Egr-1 (24, 25). Egr-1 has previously been implicated in the development and maintenance of prostate cancer (26, 27). Thus, further detailed, large-scale studies are warranted to clarify the role of Egr-1 in ACC.

One of the most interesting overexpressed genes is IGF2 on chromosome 11p15. Like other groups, we found a 3-7-fold overexpression of IGF2 in ACCs com- pared with adenomas and normal adrenal cortex (1, 13, 19). IGF-I and IGF-II are polypeptides involved in metabolism, growth and cell differentiation. They are synthesized in various tissues and have endocrine and auto/paracrine mechanisms of action depending on the tissue origin (2, 3). Both peptides are normally produced in adrenocotical cells (4, 19, 28). Through its action on steroidogenesis enzymes, IGF-I maintains the differentiated functions of the cell (19, 29). The precise role of IGF-II in mature adrenocortical gland is less clear. The IGF2 gene induction has previously been described as one of the most significant differences between ACCs and adrenal adenomas (1). Genetic altera- tions involving the 11p15 locus are very common in malignant tumors, but are found only in rare adrenal adenomas (13). The fact that we found this as well is further verification of our results. However, in adreno- cortical tumors associated with Beckwith-Wiedemann syndrome, allelic losses at the 11p15 locus have been found in both adenomas and carcinomas as well as in familial carcinomas. Current data suggest that abnorm- alities in structure and/or expression of the IGF-II gene are a late event in the multistep tumorigenesis of sporadic adrenocortical neoplasms.

Figure 2 Hierarchic cluster analysis of differentially expressed genes in adrenal cortical adenomas and carcinomas. Both genes and individual tumors are clustered in this diagram. Green indicates downregulation of a gene in a given tumor; red indicates upregulation relative to the mean. Gray spots indicate missing data. Table 4 provides the fold regulation for each gene. The numbers at the top ident- ify the individual tumors. Brackets above the tumors and on the right side of gene list indicate the clustering.

3.0

1:1

-3.0

2

AA487231

AA17 17 18

R54848

RD3124

R026 12

R62612

H96241

T95274

R68803

R95001

T97592

AA485373

AA598053

R43910

R41787

N54596

AA488626

N54690

T53404

R11520

H47257

H12044

AA452820

AA282134

NO3680

H14208

AA484106

aa459824

T98298

R56774

AA+13899

R16838

R96626

AA482067

H00813

TO 1948

T60075

AA486533

AA437705

T49652

AA480001

AA478553

H15533

N59772

AI042404

N30000

H23137

R43817

T89490

AA401111

AA405030

AA40 1693

H16389

AA291103

AA484691

AA064180

AA455235

AA134871

AA406001

AA604101

AA152347

N59799

AA429895

H20543

AA187349

TO1100

₩15203

R41972

H28681

AA150828

H05140

H22854

H19203

AA484526

Several other genes were over- or underexpressed in ACCs compared with adrenocortical adenomas and normal adrenal cortex tissues, but we have not yet validated the expression pattern by real-time PCR. The presence of a gene on the list in Table 4 does not indicate that the gene product is either necessary or sufficient for causing ACC, but only that it is expressed as part of the complex pattern of gene expression that occurs during the initiation of the disease.

During the course of our studies, three other groups performed similar analyses of adrenocortical tissue (19, 30, 31). In the search for reliable markers for the

clinical management of adrenal tumors, de Fraipont et al. (31) designed an adrenal-specific microarray. They identified two clusters of genes, the IGF2 cluster and the steroidogenesis cluster, which, in combination, provide a good predictor of malignancy. As in our studies, they also confirmed the finding from the earlier analysis by Giordano et al. (19) demonstrating upregu- lation of IGF2 expression (30). Adrenal hyperplasia was analyzed in a similar fashion by Bourdeau et al. (30).

In summary, our findings indicate that microarray analysis can distinguish between ACC and adenomas by the differential expression of a set of the genes

Table 4 List of Genes found to be differentially regulated in ACC relative to adenomas in the microarray analysisa.
Gene nameAccession no.Fold regulationUp/down
Cathepsin HAA4872313.81up
Mucolipin 3AA1717184.34up
Fibroblast growth factor receptor 1 (fms-related tyrosine kinase 2,R548464.41up
Pfeiffer syndrome)
Aldo-keto reductase family 1, member C1 (dihydrodiol dehydrogenase 1; 20-alpha (3-alpha)-hydroxysteroid dehydrogenase)R931243.45up
Fibronectin 1R626123.73up
Fibronectin 1R626123.51up
GAbinding protein transcription factor,alpha subunit 60 kDaH962413.12up
Homo sapiens clone IMAGE:120162 mRNA sequenceT952743.24up
Hydroxy-delta-5-steroid dehydrogenase, 3 beta- and steroid deltaisomerase 1R688035.67down
yq51h12.s1 soares fetal liver spleen 1NFLS Homo sapiens cDNA cloneR956914.02up
IMAGE: 199367 3', mRNA sequence
ESTT975923.26up
Likely ortholog of rat vacuole membrane protein 1AA4853732.89up
Periostin, osteoblast specific factorAA5986533.38up
yg22a04.s1 Soares infant brain 1NIB Homo sapiens cDNA cloneR439103.09up
IMAGE:32 962 3', mRNA sequence
Cadherin 13; H cadherin (heart)R417873.57up
Insulin-like growth factor 2 (somatomedin A)N545966.47up
Ubiquitin-like 1 (sentrin)AA4886262.96up
Insulin-like growth factor 2 (somatomedin A)N545962.91up
Hypothetical protein from clone 643T534042.98up
ParathymosinR115262.89up
Hypothetical protein MGC5306H472572.93down
Cytoplasmic FMR1 interacting protein 2H120446.03down
Purkinje cell protein 4AA4528263.10down
Glutaminyl-peptide cyclotransferase (glutaminyl cyclase)AA2821343.08down
zb50h07.s1 Soares_fetal_lung_NbHL19W Homo sapiens cDNA cloneN936869.65down
IMAGE:307069 3', mRNA sequence
ParalemminH142083.55down
Hyaluronoglucosaminidase 1AA4641963.41down
ATP-binding cassette, sub-family B (MDR/TAP), member 1AA4598244.57down
Putative membrane proteinT982983.22down
Bone morphogenetic protein 1R567743.18down
Scavenger receptor class B, member 1AA4438996.23down
Cytochrome p450, family 17, subfamily A, polypeptide 1R168382.87down
Chemokine (C-C motif) ligand 15R966264.73down
Retinoic acid receptor responder(tazarotene induced)2AA4820673.39down
Glutathione S-transferase theta 1H998132.95down
FBJ murine osteosarcoma viral oncogene homolog BT619486.82down
Fatty acid binding protein 5 (psoriasis-associated)T600753.38down
Early growth response 1AA4865333.26down
Immediate early response 3AA4577053.04down
Arachidonate 5-lipoxygenase-activating proteinT496523.54down
Chemokine-like factor super family 3AA4865613.19down
Dopachrome tautomerase(dopachrome delta-isomerase, tyrosinerelated protein 2)AA4785532.88down
Fibroblast growth factor 11H155333.62down
ESTN597726.40down
High mobility group AT-hook 1A10424043.73down
Glutathione S-transferase A3N300963.48down
Likely ortholog of mouse cancer related gene- liver 2H231373.13down
Human neuropeptide Y receptor Y1 (NPYY1)R438173.60down
MAX proteinT894963.10down
Glucose phosphate isomeraseAA4011118.53down
Microsomal glutathione S-transferase 1AA4959366.87down
CD 163 antigenAA4016934.24down
Transcribed sequence with moderate similarity to protein pir:A36563H163893.28down
(H.sapiens) A36563 mannose receptor precursor - human
Glutaredoxin (thioltransferase)AA2911633.09down
AdlicanAA4646913.50down
Glutathione peroxidase 3 (plasma)AA6641803.08down
Aldehyde dehydrogenase 1 family, member A3AA4552353.64down
Fibulin 1AA1348712.88down
Actin binding LIM protein 1AA4066015.91down
Aldehyde dehydrogenase 1 family, member A1AA6641015.25down
Table 4 Continued
Gene nameAccession no.Fold regulationUp/down
Glutathione S-transferase A4AA1523473.49down
Hypothetical protein MGC35366N597993.07down
ATP-binding cassette, subfamily C (CFTR/MRP), member 3AA4298955.69down
DKFZP586B1621 proteinH205433.14down
Ferredoxin 1AA1873499.33down
ESTT911002.94down
Solute carrier family 26 (sulfate transporter), member 2W152635.84down
DKFZp586P1124R419723.10down
Delta-notch-like EGF repeat-containing transmembraneH286812.95down
Mitogen-activated protein kinase kinase kinase 5AA1508282.93down
Regucalcin (senescence marker protein-30)H051406.18down
Homo sapiens cDNA FLJ40165fis, clone TESTI2015962H228543.87down
Peroxiredoxin 3H192033.28down
Interleukin 1 receptor, type 1AA4645262.84down

ªListed in order of appearance in hierarchic cluster, Fig. 2.

Figure 3 Real-time RT-PCR. Shown is a comparison of the results of the microarray analysis and quantitative PCR experiments docu- menting the expression of the indicated genes. The plot shows the change in the level of expression found in the microarray analysis and the change in the cycle-threshold value from control found in real-time RT-PCR for each gene relative to GAPDH. ACC: adrenocortical carcinoma.

IGF2

CgB

egr1

10

10

log2 expression (relative to reference)

5

5

cycles minus reference

0

0

-5

-5

-10

10

T6

T8

T10 T12

A4

A7

A9

A15

T6

T8

T10

T12

A4

A7

A9

A15

T6

T8

T10

T12

A4

A7

A9

A15

Adenomas

ACCs

Adenomas

ACCs

Adenomas

ACCs

analyzed. Further analysis of the IGF2 gene validated the experimental design. Underexpression of the CgB gene was found in both types of neoplasms, and the Egr-1 gene was downregulated in the ACCs. Further detailed analyses are warranted to elucidate the role of these genes in adrenal tumorigenesis.

Acknowledgements

We thank Michael Krause and Martin Eilers for assist- ance with the microarrays, Serdar Sel for help with the real-time PCR, and Brunhilde Chaloupka for excel- lent technical assistance.

References

1 Sidhu S, Gicquel C, Bambach CP, Campbell P, Magarey C, Robinson BG & Delbridge LW. Clinical and molecular aspects of adrenocortical tumourigenesis. Australian and New Zealand Journal of Surgery 2003 73 727-738.

2 Wooten MD & King DK. Adrenal cortical carcinoma. Cancer 1993 72 3145-3155.

3 Luton JP, Cerdas S, Billaud L, Thomas G, Guilhaume B, Bertagna X, Laudat MH, Louvel A, Chapus Y, Blondeau P, Bonnin A & Bricaire H. Clinical features of adrenocortical carci- noma, prognostic factors, and the effect of mitotane therapy. New England Journal of Medicine 1990 322 1195-1201.

4 Soreide JA, Braband K & Thoresen SO. Adrenal cortical carcinoma in Norway, 1970-1984. World Journal of Surgery 1992 16 663-668.

5 Copeland PM. The incidentally discovered adrenal mass. Annals of Internal Medicine 1983 98 940-945.

6 Lipsett MB, Hertz R & Ross GT. Clinical and pathophysiologic aspects of adrenocortical carcinoma. American Journal of Medicine 1963 35 374-383.

7 Stratakis CA. Genetics of adrenocortical tumors: gatekeepers, landscapers and conductors in symphony. Trends in Endocrinology and Metabolism 2003 14 404-410.

8 Russell AJ, Sibbald J, Haak H, Keith WN & McNicol AM. Increas- ing genome instability in adrenocortical carcinoma progression with involvement of chromosomes 3, 9 and X at the adenoma stage. British Journal of Cancer 1999 81 684-689.

9 Beuschlein F, Reincke M, Karl M, Travis WD, Jaursch-Hancke C, Abdelhamid S, Chrousos GP & Allolio B. Clonal composition of

human adrenocortical neoplasms. Cancer Research 1994 54 4927-4932.

10 Gicquel C, Leblond-Francillard M, Bertagna X, Louvel A, Chapula Y, Luton J-P, Girard F & Le Bouc Y. Clonal analysis of human adrenocortical carcinomas and secreting adenomas. Clini- cal Endocrinology 1994 40 465-477.

11 Kjellman M, Larsson C & Bäckdahl M. Genetic background of adrenocortical tumor development. World Journal of Surgery 2001 25 948-956.

12 Gicquel C, Raffin-Sanson M-L, Gaston V, Bertagna X, Plouin P-F, Schlumberger M, Louvel A, Luton J-P & Le Bouc Y. Structural and functional abnormalities at 11p15 are associated with the malignant phenotype in sporadic adrenocortical tumors: study on a series of 82 tumors. Journal of Clinical Endocrinology and Metabolism 1997 82 2559-2565.

13 Gicquel C, Bertagna X, Schneid H, Grancillard-Leblond M, Luton J-P, Girard F & Le Bouc Y. Rearrangements at the 11p15 locus and overexpression of insulin-like growth factor-II gene in sporadic adre- nocortical tumors. Journal of Clinical Endocrinology and Metabolism 1994 78 1444-1453.

14 Lack E. Tumors of the adrenal gland and extra-adrenal paraganglia. In Atlas of Tumor Pathology, 3rd Series, Fascicle 19, pp 125-151. Ed. J Rosai. Washington, DC: Armed Forces Institute of Pathology, 1997.

15 Weiss LM. Comparative histologic study of 43 metastasizing and nonmetastasizing adrenocortical tumors. American Journal of Sur- gery and Pathology 1984 8 163-169.

16 Yang YH, Dudoit S, Luu P, Lin DM, Peng V, Ngai J & Speed TP. Nor- malization for cDNA microarray data: a robust composite method addressing single and multiple slide systematic variation. Nucleic Acids Research 2002 30 e15.

17 Sturn A, Quackenbush J & Trajanoski Z. Genesis: cluster analysis of microarray data. Bioinformatics 2002 18 207-208.

18 Higuchi R, Fockler C, Dollinger G & Watson R. Kinetic PCR ana- lysis: real-time monitoring of DNA amplification reactions. Biotechnology 1993 11 1026-1030.

19 Giordano TJ, Thomas DG, Kuick R, Lizyness M, Misek DE, Smith AL, Sanders D, Aljundi RT, Gauger PG, Thompson NW, Taylor JMG & Hanash SM. Distinct transcriptional profiles of adrenocortical tumors uncovered by DNA microarray analysis. American Journal of Pathology 2003 162 521-531.

20 Yoshida R, Ohuchi N & Kimura N. Clinicopathological study of chromogranin A, B and BRCA1 expression in node-negative breast carcinoma. Oncology Reports 2002 9 1363-1367.

21 Scopsi L, Sampietro G, Boracchi P & Collini P. Argyrophilia and chromogranin A and B immunostaining in patients with sporadic medullary thyroid carcinoma. A critical appraisal of their prog- nostic utility. Journal of Pathology 1998 184 414-419.

22 Yan SF, Fujita T, Lu J, Okada K, Shan ZY, Mackman N, Pinsky DJ & Stern DM. Egr-1, a master switch coordinating upregulation of divergent gene families underlying ischemic stress. Nature Medi- cine 2000 6 1355-1361.

23 Gashler A & Sukhatme VP. Early growth response protein 1 (Egr-1): prototype of a zinc-finger family of transcription factors. Progress in Nucleic Acid Research Molecular Biology 1995 50 191-224.

24 Kartasheva NN, Lenz-Bauer C, Hartmann O, Schäfer H, Eilers M & Dobbelstein M. ANp73 can modulate the expression of various genes in a p53-independent fashion. Oncogene 2003 22 8246-8254.

25 Madden SL, Galella EA, Zhu J, Bertelsen AH & Beaudry GA. SAGE transcript profiles for p53-dependent growth regulation. Oncogene 1997 15 1079-1085.

26 Abdulkadir SA, Qu Z, Garabedian E, Song SK, Peters TJ, Svaren J, Carbone JM, Naughton CK, Catalona WJ, Ackerman JJ, Gordon JI, Humphrey PA & Milbrandt J. Impaired prostate tumorigenesis in Egr1-deficient mice. Nature Medicine 2001 7 101-107.

27 Adamson ED & Mercola D. Egr1 transcription factor: multiple roles in prostate tumor cell growth and survival. Tumour Biology 2002 23 93-102.

28 Ji B, Bi Y, Simeone D, Mortensen RM & Logsdon CD. Human pan- creatic acinar cells lack functional responses to cholecystokinin and gastrin. Gastroenterology 2001 121 1380-1390.

29 Ji B, Chen X-Q, Misek DE, Kuick R, Hanash S, Ernst S, Najarian R & Logsdon CD. Pancreatic gene expression during the initiation of acute pancreatitis: identification of EGR-1 as a key regulator. Physiological Genomics 2003 14 59-72.

30 Bourdeau I, Antonini SR, Lacroix A, Kirschner LS, Matyakhina L, Lorang D, Libutti SK & Stratakis CA. Gene array analysis of macronodular adrenal hyperplasia confirms clinical heterogen- eity and identifies several candidate genes as molecular mediators. Oncogene 2004 23 1575-1585.

31 de Fraipont F, El Atifi M, Cherradi N, Le Moigne G, Defaye G, Houlgatte R, Bertherat J, Bertagna X, Plouin P-F, Baudin E, Berger F, Gicquel C, Chabre O & Feige J-J. Gene expression profiling of human adrenocortical tumors using complementary deoxyri- bonucleic acid microarrays identifies several candidate genes as markers of malignancy. Journal of Clinical Endocrinology and Metab- olism 2005 90 1819-1829.