miRNA expression profile in primary gastric cancers and paired lymph node metastases indicates that miR-10a plays a role in metastasis from primary gastric cancer to lymph nodes
- Authors:
- Published online on: December 6, 2011 https://doi.org/10.3892/etm.2011.411
- Pages: 351-356
Abstract
Introduction
Gastric cancer is one of the most common malignant tumors of the digestive system, and its early diagnosis rate is still low in China, where gastric cancer is typically diagnosed at a more advanced stage; metastasis is typically present at the time of diagnosis. However, conventional strategies based on radical surgery for the treatment of gastric cancer are not yet satisfactory. Therefore, investigation of the mechanisms of metastasis and recurrence of gastric cancer, and exploration of the methods for the effective prediction, prevention and treatment of tumor metastasis are attracting increased attention in both basic and clinical cancer research.
microRNAs (miRNAs), discovered in 1993, are a class of non-coding RNAs 19–25 nt in length that regulate gene expression at the post-transcriptional level (1,2). These miRNAs regulate complicated cellular processes such as differentiation, proliferation and apoptosis (3) and are widely involved in tumorigenesis, tumor progression and metastasis (4,5). Recent research suggests that miRNAs are involved in tumor invasion and metastasis, in the context of breast cancer, colorectal cancer, esophageal squamous cell carcinoma, and pancreatic endocrine tumors (6,7). Few studies have reported the role of miRNAs in gastric cancer owing to its relatively low incidence in Europe and America. Therefore, profiling should be performed in different types of gastric cancer to identify differentially expressed miRNAs. In the present study, we investigated differences in miRNA profiling between primary gastric cancer and paired lymph node metastasis, and then detected the differentially expressed miRNAs in 4 cell lines, 33 gastric cancer specimens and corresponding adjacent gastric mucosa tissue, in an effort to identify the miRNAs involved in gastric cancer lymph node metastasis. Furthermore, our data provide evidence regarding the effects of miRNAs on metastasis involved in the molecular pathogenesis of gastric cancer.
Materials and methods
Specimens
Specimens used in the miRNA array
Cases of full-thickness gastric wall with invasion and lymph nodes with tumor metastasis were assessed intraoperatively; specimens from the primary tumor were retained. The specimens were embedded in optimal cutting temperature compound (OCT) and stored at −80°C until use. Experimental samples were selected to pathologically confirm that invasive adenocarcinoma of the stomach lesions was of full thickness and to verify the existence of lymph node metastasis. In total, 5 pairs of specimens (each of which included primary tumor and paired lymph node metastasis tissues) were selected, and their clinicopathological characteristics are summarized in Table I.
Specimens used for validation
Surgical specimens were obtained from 33 randomly selected patients with primary tumors among those who underwent D2 radical resection of gastric cancer at ZhongShan Hospital, Fudan University (Shanghai, China) from May, 2006 to December, 2007. All eligible cases had not received preoperative tumor-related treatment, were pathologically confirmed as gastric adenocarcinoma, and displayed full-thickness gastric wall invasion. Clinicopathological statistical data are shown in Table IV. Another 10 cases were randomly selected from the 33 cases mentioned above; normal gastric mucosa tissues were taken from the tumor marginal zone (>5 cm from the tumor). These samples were mixed in proper proportion after RNA extraction and used as the control group (indicated below with an N). Random sampling was performed with SPSS software. Ethical approval was obtained from the ZhongShan Hospital Research Ethics Committee.
Gastric cancer cell lines
Human gastric cancer cell lines, AGS and SGC7901, were purchased from the Cell Bank of the Chinese Academy of Sciences. The human gastric mucosal cell line GES-1 and moderately differentiated tubular adenocarcinoma cell line MKN-28 were a gift from the Institute of Digestive Disease, Shanghai Jiaotong University. All cells mentioned above were cultured in RPMI-1640 medium supplemented with 10% fetal bovine serum in 5% CO2 at 37°C.
Methods
Laser capture microdissection (LCM)
The collected tissues embedded in OCT were cut in 8-μm sections at −20°C, after being stained with hematoxylin and eosin (H&E). The stained slides were placed into a Veritas Microdissection Instrument (Veritas, Mountain View, CA). Areas that were enriched in tumor cells (10 mm2/slides) were selected and captured using ultraviolet laser cutting according to the manufacturer’s recommended protocol. The tissue samples were immediately placed in a microcentrifuge tube containing lysis/binding buffer (Ambion, Inc., USA), which was vortex-mixed and stored upside down at −20°C. RNA concentration and integrity were determined by UV absorption and agarose gel electrophoresis, respectively.
miRNA array
miRCURY LNA™ microRNA arrays (Agilent Technologies, Santa Clara, CA), which contain more than 1,700 capture probes, were used to compare the expression profiles of 5 paired primary gastric cancer tissues and metastatic lymph nodes collected from the same case. A difference in expression was defined as a change of >1.5-fold in the corresponding level of miRNA detection between paired specimens. The Student’s t-test was used to analyze whether there were significant differences between paired samples.
qPCR assessment of miRNA (8)
miRNA genes were selected from the Sanger Center miRNA Registry (http://www.mirbase.org). The sequence-specific primers for miRNAs and endogenous control U6 snRNA were synthesized by Sangon (Shanghai, China) (Table II). The reverse-transcription reaction was carried out with miRNA-specific stem-loop RT primer and reverse transcriptase (Fermantas, USA) according to the manufacturer’s instructions. Real-time quantification polymerase chain reaction (qPCR) analysis was carried out after cDNA was synthesized. Reaction conditions were as follows: 95°C for 10 min; followed by 40 cycles of 95°C for 15 sec, 60°C for 15 sec, and 72°C for 45 sec for the amplification. The threshold cycle (Ct) is defined as the fractional cycle number at which the fluorescence passes the fixed threshold. The gene expression Ct values of miRNAs from each sample were calculated by normalization to the internal control U6 snRNA. All experiments were repeated in triplicate.
Bioinformatic analysis of miRNA
miRGen (http://www.diana.pcbi.upenn.edu/miRGen/v3/miRGen.html) was employed for the prediction of miR-10a target genes. Functional analysis of the predicted targets was performed using MAS3.0 Bioinformatics Resources (http://bioinfo.capitalbio.com/mas3/).
Statistical analysis
miRNA expression in tumor samples was detected by real-time quantitative PCR and then subjected to N normalization. The relationship between changes in miRNA expression and clinicopathologic characteristics was analyzed using one-way ANOVA. Two-tailed P<0.05 was judged to be significant. All data and statistics processing were performed with Excel 2007 (Microsoft, USA) and SPSS 16.0 (SPSS, Chicago, IL).
Results
Results of the miRNA array
Considering the expression level of miRNAs in primary gastric cancer as baseline, there were 76 upregulated and 106 downregulated miRNAs associated with lymph node metastasis, among which four statistically significant miRNAs were included: three downregulated (miR-24-1*, miR-510 and miR-1284) and one upregulated miRNA (miR-10a) (Table III).
Differentially expressed miRNAs in primary gastric tumors
We examined the two most obviously changed miRNAs (miR-10a and miR-510) in 33 cases of primary gastric cancer and corresponding normal tissue (N). Each detected miRNA exhibited a good amplification curve; the melting curve had a single peak. The results suggested that the products of the real-time PCR were specific, and the experimental results were reliable.
Expression levels of miR-10a and miR-510 were measured in 33 cases of primary gastric cancer and its corresponding normal tissue. The average ΔCt (9) values of the tumor tissue were 12.48±1.91 and 12.78±1.76, respectively. The average ΔCt values of the corresponding adjacent normal tissue were 14.52±0.60 and 13.11±0.80, respectively. The expression of miR-10a was higher in gastric cancer than that in the adjacent normal mucosa (P<0.05) (Fig. 1).
Clinicopathologic characteristics
We investigated the relationship between expression of miR-10a and miR-510 and each of six clinicopathological characteristics (patient age, gender, tumor differentiation, major tumor histological classification, lymphatic invasion and lymph node metastasis). The results showed that only miR-10a expression was significantly related to node metastasis (P=0.047); this measure was not related to the state of lymphatic invasion (P=0.169) (Fig. 2 and Table IV).
Expression of miR-10a in the cell lines
The expression of miR-10a varied in all four cell lines. The expression of miR-10a in the MKN28, SGC7901 and AGS cell lines was 15-fold, 11.7-fold and 10.1-fold higher than that in the GES-1 cell line (Fig. 3). As MKN28 and SGC7901 were established from lymph node metastasis, the expression of miR-10a was higher in these two metastasis original cell lines than that in the AGS cell line which was established from a primary tumor.
Bioinformatic analysis of target genes for miR-10a
Potential target genes of miR-10a were predicted using miRGen; the corresponding results from the PicTar and TargetScanS databases were chosen as candidate targets; 118 potential target genes were obtained. KEGG analysis showed that the predicted target genes involved many signal pathways involved in tumorigenesis and tumor metastasis, such as the Wnt signal pathway and that initiated by cell adhesion molecules (CAMs) (Table V).
Table V.The most significant 5 pathways identified in the functional analysis of the miR-10a target genes. |
Discussion
Few research efforts to date have examined the lymph node metastasis of gastric cancer using high-throughput methods (10,11). These publications have focused mainly on the genomic and proteomic differences between primary gastric cancer patients with metastasis and those without metastasis. Samples taken from different parts of the same person may help to avoid false-positive results on the miRNA array (caused by background differences among genes). In order to improve the specificity of the chip-based results, we chose to investigate primary gastric cancer and paired lymph node metastasis sites in the same patient. As the ratio of lymph node metastasis increases with the increased depth of invasion (12,13), the specimens selected in the study all displayed full-thickness invasion.
The results obtained identified four differentially expressed miRNAs. miR-10a was upregulated while miR-24-1*, miR-510 and miR-1284 were downregulated in lymph node metastatic sites as compared to the primary gastric cancer. We hypothesized that if some relationship exists between the level of miRNA expression in the primary tumor sites and metastasis, it is quite possible that miRNAs contribute to the process of gastric cancer lymph node metastasis and may have the ability to serve as indicators in lymph node metastasis detection. In order to prove the hypothesis, we investigated the expression level of miR-10a and miR-510, which were the most obviously changed in the miRNA array, in the gastric cancer primary sites of 33 patients. The results showed that the expression level of miR-10a was affected by the existence of lymphatic metastasis, while there was no obvious relationship with the state of lymphatic invasion. Furthermore, miR-10a was upregulated in gastric cancer at both the tissue and cell levels; more obvious differences were noted at the cell level. Noticeably, the expression of miR-10a in the two metastasis original cell lines (MKN28 and SGC7910) was higher than that of the AGS cell line which was established from a primary tumor. This find also indicated that miR-10a was closely associated with lymph node metastasis. To assess the possibility of using miR-10a as a biomarker in lymph node metastasis, further research with more samples is needed.
miR-10a belongs to the miR-10 family; the other member is miR-10b. These two members differ from each other in only one base. Recently, several reports found that miR-10 expression was aberrantly increased in several types of cancers (14–17), upregulated in glioma, liver cancer, colon cancer, melanoma and breast cancer; downregulated in acute or chronic lymphocytic leukemia and head and neck squamous-cell carcinoma. However, studies concerning the function of miR-10 are still lacking. Fortunately, miR-10 was proven to play a pivotal role in tumor invasion and metastasis in these studies. Ma et al (18) found that miR-10b was significantly correlated with breast tumor cell metastasis. They also demonstrated the ability of grafted tumors in mice to invade locally and transfer to distant sites and this was associated with upregulated expression of miR-10b in breast cancer cells lacking metastatic ability. In subsequent studies, the group effectively inhibited the distant metastasis of grafted mouse tumors using miR-10b antagomirs (chemically modified anti-miRNA oligonucleotides) (19). These results stimulated related research. Weiss et al (20) found that miR-10a was highly expressed in pancreatic cancer cell lines with high metastatic ability. Experimental manipulation of miR-10a expression altered the ability of pancreatic cancer cells to metastasize and invade. This pattern was also observed in zebrafish transplanted with pancreatic cancer. In gastric cancer, research on miR-10 is still limited to describing changes in miR-10 expression, with no elucidation of the detailed mechanisms. Li et al (21) found a seven-miRNA signature in gastric cancer (including miR-10b), which had a close relationship with disease-free survival and overall survival rate. Another study excluded a role for miR-10 family members (22). Although there is no report on whether miR-10 contributes to the mechanism of gastric cancer metastasis, the relationship of miR-10a to gastric cancer lymph node metastasis demonstrated in our research provides promising evidence for subsequent research.
miRNAs function post-transcriptionally during biological processes; therefore, definitive identification of miRNA target genes is essential. For this reason, we used bioinformatic methods to predict miR-10a target genes, and then analyzed the possible mechanisms of its role in gastric tumorigenesis and metastasis. We chose the most commonly used target gene-predicting databases, PicTar (23) and TargetScanS (24), which are available on the miRGen website (25). To narrow the scope of this hunt for target genes and improve the specificity of the prediction, we chose the results common to both databases and subjected these to functional analysis. The analysis of KEGG (26) metabolic pathways showed that these predicted targets took part in many pathways involved in tumorigenesis, invasion and metastasis. For example, the Wnt signaling pathway plays an important role in organism development and gastrointestinal tumorigenesis (27), and cell adhesion molecules (CAMs) localized at the cell surface take part in intercellular and extracellular interactions.
These findings strongly suggest that miR-10a is closely correlated with gastric tumorigenesis and metastasis. A possible mechanism may be that it functions through interaction with target genes involved in tumorigenesis- and metastasis-related pathways.
Acknowledgements
This study is supported by the Key Project of the Chinese Ministry of Education (no. 108015). We thank other members of our department for their valuable comments and help.