|
|
||||||||
Guest Access | Sign In via User Name/Password |
|||||||||
* From the Department of Surgery (Drs. Hoang, DCunha, and Maddaus, and Ms. Casmey), Division of Cardiovascular and Thoracic Surgery, and the Thoracic Oncology Laboratory (Drs. M. Kratzke, Frizelle, and R. Kratzke), Division of Hematology, Oncology, and Transplant, Department of Medicine, University of Minnesota Medical School, Minneapolis, MN.
Correspondence to: Robert A. Kratzke, MD, Division of Hematology, Oncology, and Transplant, 111E MVAMC, 1 Veterans Dr, Minneapolis, MN 55417; e-mail: kratz003{at}tc.umn.edu
| Abstract |
|---|
|
|
|---|
Design: Complementary DNA (cDNA) microarray transcriptional profiling studies of 10 MPM cell lines and 4 MPM primary tumor specimens were performed using hierarchic clustering. To confirm microarray data, we used real-time polymerase chain reaction and immunoblotting.
Results: Cluster analysis differentiated among epithelial (E), sarcomatoid, and biphasic MPM variants. Expression profiling identified common overexpressed or underexpressed genes in MPM. Notably, matriptase messenger RNA was found to be overexpressed by 826-fold in E MPM, with protein expression subsequently confirmed by immunoblot analysis. This recently characterized trypsin-like serine protease has been implicated in tumor invasion and metastasis of E-derived cancers, but has not been described until now in MPM. We also identified other novel genes, such as insulin-like growth factor binding protein 5 and a cDNA clone similar to proteolipid MAL2.
Conclusions: Thus, further large-scale profiling of MPM may elucidate previously unrecognized molecular mechanisms by identifying novel genes that are involved in malignant transformation. Our study has now found matriptase to be one of these mesothelioma-associated genes, with potential pathogenic and therapeutic significance.
Key Words: expression profiling matriptase mesothelioma microarray real-time polymerase chain reaction
| Introduction |
|---|
|
|
|---|
Aside from asbestos exposure, other factors such as ionizing radiation or tumor DNA virus sinvian virus-40 may act synergistically in MPM pathogenesis.56 Also, several well-defined acquired genetic targets have been identified in MPM, including the 9p21 locus (p16INK4a, p14ARF) and the 22q11-q13.1 locus (NF2).2 However, the molecular mechanisms controlling the transformation of mesothelial cells remain poorly defined. This is underscored by the observation that these well-characterized etiologies incompletely account for the known incidences of MPM. About 10 to 20% of MPM occurrences have been documented in patients without previous exposure to asbestos,2 and only 60% of MPM tumors are known to contain SV-40 viral DNA.6 Accordingly, multiple active pathways are thought to be possible.
Only two studies78 to date have used the microarray technique to search for additional genes and pathways that are potentially involved in MPM biology. These investigations were limited by the extremely small number of MPM cell lines examined and by the lack of additional analysis of primary tumor tissues. To identify novel candidate genetic targets, we used hierarchic cluster analysis in our current study to compare the gene profiles of a larger sampling of MPM cell lines and primary tumors vs nonmalignant mesothelium on complementary DNA (cDNA) microarrays.
| Materials and Methods |
|---|
|
|
|---|
|
RNA Isolation
Cells growing asynchronously were lysed in reagent (Trizol; Invitrogen Life Technologies) when they reached about 80% confluence and then were processed as described previously,11 with modifications. Before RNA precipitation by isopropyl alcohol, the RNA aqueous phase was mixed with an equal volume of 70% ethanol, and total RNA was subsequently extracted on a silica gel-based membrane spin column (Qiagen; Valencia, CA) per the manufacturers instructions. On-column deoxyribonuclease digestion was performed with all samples to eliminate potential genomic DNA. RNA yield and purity were determined by spectrophotometry. Integrity was verified on 1.5% agarose-formaldehyde gels stained with ethidium bromide. Tissue specimens were processed in a similar fashion, each starting from 100 mg frozen tumor.
Microarrays
Microarray experiments were performed (MicroMax Human cDNA Microarray System II-TSA; Perkin Elmer LifeSciences; Boston, MA) according to the manufacturers instructions. Briefly, 1 µg total RNA of each MPM sample was reverse-transcribed and simultaneously labeled with fluorescein-deoxyuridinetriphosphate to produce target cDNA. CRL-9444 target cDNA was labeled with biotin-deoxyuridinetriphosphate and served as the reference "normal" sample in each microarray experiment, as validated by others.7 Target cDNA was hybridized on 4,800 gene MicroMax microarrays. The washing and detection steps were based on a sequential fluorescence detection process with horseradish peroxidase-conjugated antibodies and tyramide-linked cyanine-3 (Cy3) or cyanine-5 (Cy5) dyes.1213 Microarray slides were scanned at 532 nm (Cy3, MPM samples) and 635 nm (Cy5, reference) [ScanArray Express confocal laser scanner; Perkin Elmer]. Some microarray experiments were selected for duplicate processing to assure reproducibility.
Microarray Analysis
Each gene was represented on a microarray in duplicate, so that a corresponding raw expression ratio was defined as the mean signal intensity of the Cy3/Cy5 replicates. Expression ratios then were filtered per microarray, according to the following algorithm,14 with modifications. Briefly, the 50th percentile of all measured expression ratios was used as a positive control for each gene. Every gene measurement was divided by this synthetic positive control, assuming that it was at least 0.01. The bottom 10th percentile of expression ratios was used as a test for correct background subtraction. This was never less than the negative of the synthetic positive control. Across all 14 MPM specimens, the threshold values used to define significant relative expression changes were set at 3.0 for overexpression and 0.30 for underexpression.
In a separate but complementary approach, the raw expression ratios for each gene again were subjected to median-centered normalization, followed by K groups analysis (separate groups test), which is a combination of gene rank-ordering tests followed by a Pearson correlation analysis (Expressionist Suite, version 3.1; GeneData; San Francisco, CA). Microarray experiments were divided into three separate groups, according to histologic type. An iterative approach selected a K groups threshold value of 0.9, which best identified a subset of significant, differentially expressed, simultaneously divergent genes expressed between MPM histologic types. This gene subset then was used in agglomerative two-dimensional hierarchic cluster analysis.1516
Real-Time Polymerase Chain Reaction
To examine the reliability of our microarray data, we identified an overexpressed gene (matriptase) with a potential role in MPM pathogenesis and determined its relative expression in identical samples (n = 14) used for microarray analysis by quantitative real-time polymerase chain reaction (QRT-PCR). We designed a probe (TaqMan probe; Applied Biosystems; Foster City, CA) and primer set from the matriptase open reading frame (Primer Express software, version 2.0; Applied Biosystems). The probe and primer set for the endogenous control gene ß-glucuronidase (GUSB) was identical to that of a published sequence.17 The sequences of matriptase primers and probe are as follows: forward, 5'-GCGCTCCCTGAAGTCCTTT-3'; reverse, 5'-GTCCTGGGTCCTCTGTACTGTTTT-3'; and probe, 5'-TCACCTCAGTGGTGGCTTTCCCCA-3'. The gene specificity of sequences derived for primers and probes was confirmed by BLASTn (National Center for Biotechnology Information; Bethesda, MD) searches against the dbEST database (National Center for Biotechnology Information) and the nonredundant set of GenBank.
cDNA was generated from 1 µg total RNA (Gene-Amp system; Perkin Elmer) primed with random hexamers. In a 25-µL polymerase chain reaction (TaqMan; final volume after cDNA added), components included the following: 1x TaqMan buffer; 200 mmol/L each of four deoxynucleotidetriphosphates; 5 mmol/L magnesium chloride, 1.25 U AmpliTaq Gold, 0.5 U AmpErase uracil N-glycosylase (all from Applied Biosystems); 200 nmol/L each of forward and reverse primers; and 100 nmol/L probe (TaqMan). Then, 100 ng cDNA was added, and samples were processed in triplicate under the following cycling parameters (ABI PRISM 7900HT Sequence Detection System; Applied Biosystems): 50°C for 2 min; 95°C for 10 min; and 40 cycles at 95°C for 15 s and at 60°C for 1 min. The coefficients of variation for triplicate reactions and between assays were < 10%. The expected 79-base pair amplicon product was eluted from 2% agarose gel after electrophoresis and ethidium bromide staining, and its identity was confirmed by sequence analysis.
The relative amount of mRNA for a particular sample was represented by the threshold cycle (Ct) of amplification. The relative quantitation of matriptase transcripts was calculated by the comparative Ct method.1819 GUSB was used as an endogenous control for the normalization of sample loading, because levels of GUSB are consistent across diseased and nondiseased lung tissues.17 The GUSB Ct value was subtracted from the matriptase CT value to obtain a
Ct. Likewise, we determined the difference between the
Ct values of MPM samples for the target gene and the
Ct value of the calibrator (CRL-9444). The calibrator was chosen to be the same reference sample used in microarray hybridizations to facilitate comparison between the different techniques. The relative, normalized quantitative matriptase expression level (2
Ct) was calculated for all samples.
Immunoblotting
Tumor tissues were kept frozen while being ground to a fine powder, and cell lines were washed in phosphate-buffered saline solution before the addition of cell lysis buffer.9 Protein concentration was determined by the Bradford dye-binding protein assay (Bio-Rad Laboratories; Hercules, CA). Proteins (50 µg) were resolved by 10% sodium dodecyl sulfate-polyacrylamide gel electrophoresis, were transferred to a nitrocellulose membrane (Amersham Biosciences; Piscataway, NJ), and then were probed with a matriptase-specific antibody, using previously described protocols.20 The antimatriptase murine monoclonal antibody M3221 was obtained as a generous gift from Dr. Chen-Yong Lin (Lombardi Cancer Center, Georgetown University Medical Center; Washington, DC) and was used at a dilution of 1:2,000. Detection and visualization employed a horseradish peroxidase-labeled secondary antibody coupled to the detection system (ECL Detection System; Amersham Biosciences), followed by exposure to autoradiography film (Amersham Biosciences).
Statistical Analysis
A correlation analysis was performed to determine differences in gene expression results between microarray and QRT-PCR techniques. Matriptase expression ratios (ie, tumor specimen/reference sample ratio) calculated for samples from each technique were log-transformed and compared. The Pearson correlation coefficient (r value) was calculated for all samples. We also calculated a linear regression line. Statistical analysis was performed using a statistical software package (S-Plus 6; Insightful; Seattle, WA).
| Results |
|---|
|
|
|---|
|
|
|
Verification of Microarray Data
We validated the microarray expression data for the matriptase gene at the level of transcription. The matriptase-specific QRT-PCR primer-probe set quantitatively estimated relative amounts of matriptase mRNA transcripts. Overall, results were very similar between the two different techniques, with a correlation coefficient of 0.710 (Fig 3
). This data confirmed the reliability of our microarray strategy to identify differentially expressed genes in MPM specimens. We then extended this validation by estimating MPM matriptase protein content in seven of the MPM cell lines and in the four tissue specimens to determine whether the mRNA overexpression was reflected by a corresponding increase in proteins (Fig 4
). The M32 monoclonal antibody detected, as expected, both the single-chain and two-chain forms of matriptase (70-kd) in 8 of 11 tested samples. Thus, these results collectively revealed a positive qualitative correlation among the three different analyses of matriptase expression in our MPM specimens.
|
|
| Discussion |
|---|
|
|
|---|
To the best of our knowledge, this study represents the largest of its kind to specifically aim at finding novel mesothelioma-associated genes, which may provide for insights into mechanisms of tumorigenesis. We analyzed a total of 14 specimens of MPM cell lines and fresh-frozen tissues on a 4,800-gene cDNA microarray. Rihn and colleagues7 were the first group to employ microarrays to identify potential genes involved in pleural cell malignant transformation, but the major limitation of their study was the use of a single MPM cell line. From subsequent microarray studies, it has been acknowledged that gene-profiling precision and reliability increase with the number of samples grouped for analysis.24 Another study by Kettunen and colleagues8 had similar limitations. It investigated the expression patterns of four MPM cell lines (E, one cell line; B, three cell lines) derived from primary tumors on a 588-gene, cancer-specific array. In contrast, Gordon and colleagues2526 studied larger numbers of MPM tissue samples with oligonucleotide microarrays, but their focus was on identifying a minimal number of genes for use as potential clinical diagnostic tests and prognostic tests.
Among genes identified as significantly overexpressed across specimens on microarrays, we quantitatively validated matriptase overexpression by QRT-PCR and further confirmed our results with immunoblotting. Matriptase overexpression may be important in the pathogenesis of MPM. Matriptase has been identified as a type II transmembrane serine protease expressed in a wide variety of epithelial human cancers, including prostate cancer,27 breast cancer,28 and ovarian cancer,29 but not in mesenchymal tissues.9 It has been proposed to have multiple functions, acting as a potential activator of critical molecules associated with tumor invasion and metastasis. One study30 demonstrated that matriptase could activate hepatocyte growth factor/scatter factor, a mesenchymal cell-derived cancer cell growth and motility protein. Takeuchi and colleagues31 characterized matriptase to be an in vitro activator of prourokinase plasminogen activator, linking it to the activation of additional protease systems implicated in cancer cell invasion and metastasis. Our validated microarray expression data for matriptase are consistent with those in the literature, indicating the highest levels in E MPM, moderate levels in B MPM, and (relatively) very low levels in S MPM. B MPM shares qualities of the other two histologic types and could be expected to have intermediate matriptase expression, as found in our microarray data. Matriptase overexpression in MPM is, thus, a novel finding that we validated at both the mRNA and protein level. However, future investigations will need to better define the functional role of matriptase in the molecular pathogenesis of MPM.
We also identified additional novel genes with potential importance in molecular mechanisms relating to MPM. For example, IGFBP5 was underexpressed in MPM, but IGF-I was overexpressed across all MPM specimens. The IGF axis was previously characterized to be a complex system composed of the IGF ligands, a superfamily of IGFBPs, and IGF transmembrane receptors.3233 Dysregulation at each level in this system has been implicated in cancer growth and in the progression of multiple cell types, including prostate,34 lung,3536 and colon.37 The precise functional role of IGFBP5 in cancer remains poorly defined and is likely to be cell type-specific.38 Some studies have documented IGFBP5-induced growth stimulation of prostate cancer cells39 and breast cancer cells,40 while other studies have found that IGFBP5 inhibited cervical carcinoma41 and osteosarcoma42 cell proliferation. Thus, we speculate that IGFBP5 acts as an inhibitor of IGF-I in MPM and that underexpression of IGFBP5 may contribute to uncontrolled cellular expansion via an IGF-mediated autocrine growth loop. In addition, we identified from the discriminatory 180-gene subset a cDNA clone similar to the transport and signal transduction proteolipid MAL2 (member of the MAL protein family) that may be important, since MAL has been implicated as a molecular marker in B-cell lymphoma,43 and recently was identified in E-derived renal and thyroid carcinomas.44 Furthermore, MAL2 gene expression and protein function have been studied in hepatocellular and colorectal carcinoma cell lines.45
A potential limitation of our study is the use of a single type of reference for normal (CRL-9444), which may have introduced experimental bias. Currently, there is no consensus on what constitutes an ideal reference for any particular type of microarray experiment.46 Multiple strategies have been used successfully by others. Ross and colleagues47 used a "universal" reference RNA that was derived from an equal mix of 11 diverse human tumor cell lines. Alternatively, Holloway et al46 and Kikuchi and colleagues48 used commercially available RNA references derived from pooled, tissue-specific nonmalignant samples. A third, equally effective strategy, described, for example, by Miura and colleagues,49 used immortalized, nonmalignant cells (bronchial epithelial) that are known to transform into the cancer of interest (lung adenocarcinoma). Thus, without an absolute standard, we chose the mesothelial cell line CRL-9444 as our normal control for the following reasons: (1) any human mesothelial cell is a rational baseline since it represents the noncancerous cell of origin for malignant mesotheliomas; (2) a cell line allowed for adequate and consistent amounts of RNA to be isolated for multiple cDNA microarray experiments; (3) Rihn and colleagues7 validated its use as a normal reference in their mesothelioma cDNA microarray study; and (4) continued and widened usage of this reference mesothelial cell line in subsequent cDNA microarray studies may facilitate the comparison of data from different research groups.
In summary, transcriptional profiling of MPM cell lines and primary tumor tissues has identified novel genes like matriptase and IGFBP5 with significant differential expression. These genes may have a role in the pathogenesis of MPM, a malignancy known to be heterogeneous and composed of complex phenotypes. Furthermore, these genes may represent novel therapeutic targets. In our study, hierarchic cluster analysis yielded a set of genes that was able to discriminate among histologic MPM types. Our findings suggest that an alternative method of MPM classification may be based on distinct gene expression profiles. Clearly, additional validated, large-scale, microarray-based studies of MPM are needed to confirm these results. Functional gene studies are further required to establish any potential clinical significance.
| Acknowledgements |
|---|
| Footnotes |
|---|
Dr. Hoang was supported by grants from the Lillehei Heart Institute, the Veterans of Foreign Wars/Ladies Auxiliary Cancer Research Center Endowment Fund, and a by training grant from the National Institutes of Health (grant No. T32HL07062) that was awarded to the Division of Hematology, Oncology, and Transplant, University of Minnesota. Dr. R. Kratzke was supported, in part, by grants from the National Institutes of Health (grant No. R21CA83689) and from the Veterans Affairs Research Service.
Received for publication June 3, 2003. Accepted for publication October 6, 2003.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
B. A. Whitson, B. A. Jacobson, S. Frizelle, M. R. Patel, D. Yee, M. A. Maddaus, and R. A. Kratzke Effects of insulin-like growth factor-1 receptor inhibition in mesothelioma. Ann. Thorac. Surg., September 1, 2006; 82(3): 996 - 1001. [Abstract] [Full Text] [PDF] |
||||
![]() |
K.-i. Kiyomiya, M.-S. Lee, I-C. Tseng, H. Zuo, R. J. Barndt, M. D. Johnson, R. B. Dickson, and C.-Y. Lin Matriptase activation and shedding with HAI-1 is induced by steroid sex hormones in human prostate cancer cells, but not in breast cancer cells Am J Physiol Cell Physiol, July 1, 2006; 291(1): C40 - C49. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. List, R. Szabo, A. Molinolo, B. S. Nielsen, and T. H. Bugge Delineation of Matriptase Protein Expression by Enzymatic Gene Trapping Suggests Diverging Roles in Barrier Function, Hair Formation, and Squamous Cell Carcinogenesis Am. J. Pathol., May 1, 2006; 168(5): 1513 - 1525. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. List, R. Szabo, A. Molinolo, V. Sriuranpong, V. Redeye, T. Murdock, B. Burke, B. S. Nielsen, J. S. Gutkind, and T. H. Bugge Deregulated matriptase causes ras-independent multistage carcinogenesis and promotes ras-mediated malignant transformation Genes & Dev., August 15, 2005; 19(16): 1934 - 1950. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. D. Hoang, X. Zhang, P. D. Scott, T. J. Guillaume, M. A. Maddaus, D. Yee, and R. A. Kratzke Selective Activation of Insulin Receptor Substrate-1 and -2 in Pleural Mesothelioma Cells: Association with Distinct Malignant Phenotypes Cancer Res., October 15, 2004; 64(20): 7479 - 7485. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |