Detection of differentially expressed genes and association with clinicopathological features in laryngeal squamous cell carcinoma
- Authors:
- Published online on: September 18, 2012 https://doi.org/10.3892/ol.2012.920
- Pages: 1354-1360
Abstract
Introduction
Head and neck cancer is the sixth most common type of cancer in the world, accounting for more than 540,000 new cases and 271,000 mortalities each year (1). These cancers occur in the lips, oral cavity, nasal cavity, paranasal sinuses, pharynx and larynx, 90% of which are squamous cell carcinomas. They significantly affect long-term survival and the quality of life of patients. The development of novel strategies is required for prevention and early detection, to reduce cancer incidence and overcome problems associated with treatment of late-stage tumors. Improved prediction of outcome will lead to treatment decisions that prolong patients’ survival and quality of life.
Predictions concerning the outcome of head and neck cancers are currently based mostly on clinicopathological features, including tumor stage, differentiation, size and regional lymph node or distant metastasis. However, increasing numbers of studies utilize aberrant gene expression and genomic and epigenetic alterations to predict prognosis.
It has been established that tobacco smoke and alcohol consumption are the most significant risk factors for head and neck cancers. They contribute to these cancers through multiple genetic alterations, including the silencing of tumor suppressor genes and oncogene activation (2). A large body of knowledge has accumulated regarding gene alterations that are associated with the development of this deadly disease. However, greater understanding of the links between gene alteration and head and neck cancer development and progression is required.
DNA microarray profiling is an innovative technology that facilitates analysis of a great number of genes simultaneously. In this study, we performed an analysis of nearly the entire human genome in order to detect altered gene expression between primary laryngeal squamous cell carcinoma (LSCC) and adjacent normal tissues. We identified 349 genes that are differentially expressed between normal and malignant tissues, a number of which have not previously been associated with LSCC. Thus, we selected specific tumor-related genes from microarray data, which have not been reported in LSCC before. We then verified the differential expression of the seven genes in another set of LSCC tissues. In the future, these genes may be further evaluated as biomarkers or potential therapeutic targets for LSCC.
Materials and methods
Tissue samples
Between October 2007 and March 2010, tissue biopsy specimens (tumor and matched adjacent normal tissues) were collected from 36 patients (Table I) with LSCC at the Department of Otorhinolaryngology Head and Neck Surgery, Drum Tower Hospital Affiliated to Nanjing University School of Medicine (Nanjing, China). Pathological analyses confirmed the diagnosis of each patient. The tissues were immediately frozen in liquid nitrogen and stored at −80°C until use. Our institutional Human Ethics Committee approved the study. Informed consent was obtained either from the patient or the patient’s family.
RNA isolation and microarray analysis
To isolate RNA from the tissue specimens, both tumor and normal mucosae were put into liquid nitrogen and ground into powder in TRIzol reagent (Invitrogen Life Technologies, Carlsbad, CA, USA) using a rotor-stator homogenizer. Total RNA was then isolated following the manufacturer’s instructions. The integrity of the RNA was verified by visual inspection after 1% agarose gel electrophoresis; the 28S ribosomal RNA band intensity was two times that of the 18S ribosomal RNA band (3). Sample purity was ensured by an OD260/OD280 ratio >1.8, measured with a spectrophotometer.
For DNA microarray analysis, biotinylated probes were prepared using 2 μg of total RNA. Briefly, the total RNA obtained from tumor and normal tissues was mixed with 100 pmol of T7-oligo(dT)24 primer and denatured at 70°C for 10 min, then chilled on ice. The first-strand cDNA synthesis was performed using Superscript II reverse transcriptase (Life Technologies, Carlsbad, CA, USA) and the second-strand with DNA Polymerase I, E. coli DNA ligase and RNase H. The biotinylated probes were then prepared from the entire cDNA reaction using an ENZO Bioarray High Yield RNA Transcript Labeling kit (ENZO Diagnostics, Toronto, Canada).
The purified probes were incubated with 1X fragmentation buffer at 95°C for 35 min to reduce the average probe length. Hybridization was performed at 45°C for 20 h with biotinylated probes on the microarrays. The non-specific binding of these probes was removed by low stringency washes (10 times) and high stringency washes (4 times) using a GeneChip Fluidics Station 400 wash station (Agilent, San Diego, CA, USA). The positive signal was detected by incubating the microar-rays with streptavidin phycoerythrin (Molecular Probes, Camarillo, CA, USA) and scanned with a GeneArray Scanner (Hewlett-Packard, San Diego, CA, USA). The scanned data were analyzed with GeneChip Analysis Suite 3.3 (Agilent).
Semi-quantitative reverse transcription polymerase chain reaction (RT-PCR)
To confirm the differential gene expression of laryngeal cancer revealed during cDNA microarray analysis, we used a 2-step method of semi-quantitative RT-PCR starting with tissues from 32 cases of laryngeal cancer and matched normal adjacent tissues. Briefly, total RNA was first reverse transcribed into cDNA using Superscript II reverse transcriptase (Life Technologies) and then amplified in a programmable Applied Biosystems 2720 thermal cycler (Singapore). For each reaction, a 50-μl PCR mixture containing 200 μM dNTPs, 1.25 units Taq polymerase in 10X Taq polymerase buffer (Takara Bio, Inc., Shiga, Japan), and corresponding concentrations of primers (Table II) was set to an initial denaturing at 95°C for 5 min and then appropriate PCR cycles for different genes of 94°C for 1 min, annealing temperature (Table II) for 1 min, 72°C for 30 sec and a final extension at 72°C for 10 min in a programmable 2720. The PCR reactions were performed in triplicate.
The PCR-amplified gene products were visualized in a 2% (w/v) agarose gel stained with ethidium bromide. Images of resulting gels were captured with LabWorks45 (UVP, Upland, CA, USA). The genes detected by PCR were SENP1, CD109, CKS2, LAMA2, LAMA3, ITGAV, ITGB8 and β-actin (Table II). β-actin was used as the loading control and normalizing reference for each gene in these tissue samples. The primers were designed according to their GenBank sequences using the Primer 3 online tool.
Protein extraction and western blot analysis
Both LSCC and the matched adjacent normal tissues were homogenized for total cellular protein extraction using a commercial protein kit from Pierce Biotechnology (Rockford, IL, USA). The protein concentration of the homogenates was determined by a bicinchoninic acid protein assay kit (Shenergy Biocolor, Shanghai, China).
Equal amounts of the protein samples (50 μg) were separated via 10–15% sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE), followed by electrophoretic transfer onto polyvinylidene difluoride membranes (Roche Diagnostics, Indianapolis, IN, USA). These membranes were incubated with 5% non-fat milk in phosphate-buffered saline (PBS) for 2 h and then with the primary antibody at 4°C overnight. The primary antibodies CKS2 (#ab54658) and SENP1 (#ab3656) were purchased from Abcam (Cambridge, MA, USA). CD109 (#SC33115) was obtained from Santa Cruz Biotechnology, Inc. (Santa Cruz, CA, USA). The next day, the membranes were washed with PBS 3 times and then incubated with an anti-goat or anti-rabbit horseradish peroxidase-conjugated secondary antibody (Santa Cruz Biotechnology, Inc.). The immunoreactive signals were visualized using an enhanced chemiluminescence detection kit (Pierce Biotechnology) and quantified with a densitometer (Kodak Digital Science 1D Analysis Software, Rochester, NY, USA).
Statistical analysis
DNA microarray data were analyzed using the Agilent GeneChip Analysis Suite 3.3 and summarized as fold changes. The data from semi-quantitative RT-PCR and western blot analysis were summarized as percentages of controls. The differential expression levels of genes between the tumor and normal tissues were statistically analyzed with paired-sample t-tests using SPSS 16.0 software (SPSS, Inc., Chicago, IL, USA). The association of gene expression levels with clinicopathological data was statistically analyzed with an independent-samples t-test. P<0.05 was considered to indicate a statistically significant difference.
Results
Detection of differentially expressed genes between the primary LSCC and corresponding normal tissues
In this study, we first randomly selected 4 pairs of primary laryngeal cancer and corresponding normal tissues for DNA microarray analysis. We then isolated RNA from the frozen tissues and performed DNA microarray analysis in Agilent chips. We identified that 10,909 genes were differentially expressed between laryngeal cancer and the matched normal tissues in case 1; 10,223 genes in case 2; 5,730 genes in case 3; and 14,665 genes in case 4. Among these differentially expressed genes, there were 349 that were identified in all four cases, of which 112 were significantly upregulated with intensity ratios up 2.0, while 237 were downregulated with ratios down 0.5 (Table III).
Table III.Differentially expressed genes between primary laryngeal cancer and corresponding normal tissues. |
Validation of microarray data using semi-quantitative RT-PCR or western blot analysis
From the microarray data we chose 7 genes whose differential status was validated with semi-quantitative RT-PCR or western blot analysis, using source material from 32 cases of laryngeal cancer and the corresponding normal tissues. The results demonstrated that expression of these 7 genes were in accordance with the microarray data (Figs. 1 and 2). Expression levels of SENP1, CD109, CKS2, LAMA3, ITGAV and ITGB8 mRNA were all increased compared with the normal tissues, while LAMA2 mRNA was significant decreased in tumor tissues compared with normal tissues. As shown in Table IV, of the 32 laryngeal cancers, compared with normal epithelial tissues mRNA expression of SENP1 was significantly elevated in 22 cases (68.8%), LAMA3 in 23 (71.9%), CD109 in 26 (81.3%), CKS2 in 25 (78.1%), ITGAV in 22 (68.8%) and ITGB8 in 20 (62.5%), while LAMA2 was significantly less in 18 (56.3%). Western blot data showed that of these 32 laryngeal cancer tissues, compared with the corresponding normal tissues, SENP1 protein levels were markedly higher in 21 cases (65.6%), CD109 in 24 (75%) and CKS2 in 23 (71.9%; Table V).
Table V.Western blot analysis of SENP1, CD109 and CKS2 protein expression in laryngeal cancer and the corresponding adjacent normal tissues. |
Association of the expression of three genes with patient clinicopathological data
Statistical analysis revealed that protein levels of the three genes SENP1, CD109 and CKS2 were significantly different between tumor and corresponding normal tissues (P≤0.05). We examined their expression levels for associations with clinicopathological data, including age and tumor classification, stage, differentiation and lymph node metastasis (Table VI). SENP1 expression differed between stage I+II and III+IV tumors. CKS2 expression differed with tumor classification, tumor differentiation and lymph node metastasis. CD109 expression differed between glottic carcinoma and subglottic carcinoma.
Table VI.Table IV. Association of SENP1, CD109 and CKS2 expression levels with patient clinicopathological data. |
Discussion
A profile of the genes that were differentially expressed between laryngeal cancers and corresponding normal mucosae was created using cDNA microarray analysis. A total of 349 differentially expressed genes were identified in four patients, of which 112 were significantly upregulated and 237 were downregulated. We also identified certain genes that were altered in LSCC, including P63 and Mdm2. We also found other genes that were altered in human cancers, but which have not been identified before in LSCC. Thus, we selected 7 genes to study the differential mRNA and protein expression in LSCC using semi-quantitative RT-PCR and western blot analysis, respectively. The data demonstrated that, compared with the normal mucosae, 6 of the 7 genes were upregulated in laryngeal cancer and 1 was downregulated.
Associations between these genes and clinicopathological data from the patients were identified. For example, the expression of SENP1 was associated with tumor stage, while CKS2 was associated with tumor classification, differentiation and lymph node metastasis. These results imply that the detection of elevated levels of expression of SENP1 and CKS2 should be further evaluated as tumor markers for early detection or prognosis of laryngeal cancer.
The identification of genes differentially expressed between normal and malignant tissues is the first step to understanding how altered expression may contribute to tumorigenesis. These genes are likely to represent critical points of alteration in pathways that regulate the cell cycle, cell-cell adhesion and cell motility. Our current data identified a large number of genes differentially expressed in tumor tissues compared with normal tissues of the larynx. Specifically, in 4 cancer cases there were 10,909, 10,223, 5,730 and 14,665 differentially expressed genes. Only 349 genes were identified to be differentially expressed in all 4 cases. These data indicate that each patient has an individualized profile with genes that may be targeted for personalized medical treatment in future studies. However, there are also genes that are commonly altered in laryngeal cancer that may be evaluated as biomarkers for early detection and prediction of prognosis of laryngeal cancer.
Of the genes whose expression is commonly elevated in laryngeal cancer, the small ubiquitin-like modifier (SUMO) is involved in numerous cellular processes, including nuclear-cytosolic transport, transcriptional regulation, apoptosis, protein stability, response to stress and progression through the cell cycle (4). The family of sentrin/SUMO-specific proteases (SENPs) is one of a group of enzymes that process newly synthesized SUMO1s into the conjugate form and catalyze the deconjugation of SUMO-containing species to regulate the function of the SUMO protein. SENP1, a member of the SENP family, has been reported to be overexpressed in colon cancer tissues (5) and has been demonstrated to regulate androgen receptor transactivation by targeting histone deacetylase I, and induce c-Jun activity through de-SUMOylation of p300 (6). Previous studies have demonstrated that SENP1 was able to transform normal prostate epithelium into a dysplasia, and also directly modulate several oncogenic pathways in prostate cells (7,8). However, it remains unknown whether and how the expression of SENP1 plays a role in laryngeal cancer.
In the present study, we identified that SENP1 mRNA and protein were highly expressed in laryngeal cancer tissues. SENP1 expression was also statistically different between stage I+II and III+IV tumors. Together, these findings indicate that SENP1 may play a role in tumorigenesis and progression of laryngeal cancer.
CD109 is a glycosylphosphatidylinositol-linked glycoprotein which belongs to the α2 macroglobulin/C3/C4/C5 family of thioester-containing proteins. A previous study revealed that CD109 was a useful diagnostic marker for basal-like breast carcinoma (9). Another study identified that CD109 may be involved in bladder tumorigenesis and might be a potential target for cancer immunotherapy (10). Zhang et al (11) detected CD109 expression in half of lung squamous cell carcinomas cases investigated, but not in lung adenocarcinomas, large cell carcinomas or small cell carcinoma. In addition, CD109 expression was found to be upregulated in approximately 50% of esophageal squamous cell carcinoma cases. In the present study, we identified that CD109 was highly expressed in laryngeal cancer tissues, and CD109 expression differed between glottic carcinoma and subglottic carcinoma. However, the number of cases of subglottic carcinoma in this study were too few to make a definitive statement; these data require confirmation with additional studies with larger sample sizes.
CKS2, an essential component of cyclin/cyclin-dependent kinase complexes, contributes to cell cycle progression. Earlier studies identified CKS2 as a transcriptional target that was downregulated by the tumor suppressor p53 (12). Other studies demonstrated that expression of CKS1 and CKS2 were elevated in prostate cancer, while knockdown of CKS2 expression induced programmed cell death and inhibited tumorigenicity (13). Another study demonstrated that CKS2 was significantly expressed in metastasized tumors (14). Using oligomicroarray analysis and qRT-PCR, Uchikado et al (15) identified CKS2 as a gene associated with the lymph node metastasis of esophageal squamous cell carcinoma. The present study also revealed that CKS2 was highly expressed in laryngeal cancer and associated with tumor classification, tumor differentiation and lymph node metastasis.
Laminin, a basement membrane protein consisting of α, β and γ chains, plays a critical role in the maintenance of tissue structures (16). Laminin expression is a prerequisite for normal embryonic development (17). Abnormal expression of LAMA332 and its integrin receptors is a hallmark of certain types of tumor and is considered to promote the invasion of colon, breast and skin cancer cells (18). Our present study confirmed these previous findings, although its role in laryngeal cancer requires further study.
Integrins are a group of cell adhesion molecules that regulate a wide variety of dynamic cellular processes, including cell migration, phagocytosis, growth and embryonic development. The interaction of integrins with extracellular ligands is regulated from inside the cell through the short cytoplasmic α- and β-integrin tails, and transmits biochemical and mechanical signals to the cytoskeleton to change cell shape, behavior and fate (19). ITGAVB6 has a role in the inhibition of colon cancer cell apoptosis through targeting the mitochondrial pathway (20). Another study (21) demonstrated that antisense ITGAV and ITGB3 inhibited tumor vascularization and growth, but enhanced the apoptosis of tumor cells. Antisense ITGAV suppressed tumor growth more markedly than antisense ITGB3. Loss of the ITGB8 subunit resulted in abnormal blood vessel development in the yolk sac, placenta and brain (22); animals lacking the ITGB8 gene die either at mid-gestation (due to insufficient vascularization of the placenta and yolk sac) or shortly after birth with severe intra-cerebral hemorrhage (22). Our present study also demonstrated altered expression of these integrins in laryngeal cancer.
In conclusion, our study provides the first evidence that SENP1, CD109, CKS2, LAMA2, LAMA3, ITGAV and ITGB8 are differentially expressed in laryngeal cancer tissue specimens. Further study of these seven genes may aid the understanding of the multistep process of laryngeal tumorigenesis, and evaluate them as tumor biomarkers for early detection or prediction of prognosis of laryngeal cancer.
Acknowledgements
We would like to thank Dr J.D. Wu of Nanjing Medical University for his helpful advice and Mr. J.Y. Shen and Mr. L. Zhang for their technical assistance. This study was supported in part by grants from the Nanjing Technology Development Project Fund (No. 200702075), the Medical Science and Technology Development Foundation, Nanjing Department of Health (No. YKK10175), and the BenQ Medical Center Research Fund (No. SRD20100001).