Plasma Exosomal Long Noncoding RNA lnc-SLC2A12-10:1 as a Novel Diagnostic Biomarker for Gastric Cancer.

PURPOSE: Exosomes participate in cellular communications by transmitting active molecules, including long noncoding RNAs (lncRNAs) and are regarded as suitable candidates for disease diagnosis. This study aimed to identify gastric cancer (GC)-specific exosomal lncRNA and investigate the potential diagnostic value of plasma exosomal lncRNA in GC. PATIENTS AND METHODS: Exosomes from the culture media (CM) of four GC cells (GCCs) and human gastric epithelial cells were isolated. Exosomal RNA was extracted, and lncRNA microarray assay was performed to identify GC-specific exosomal lncRNAs. The expression levels of the candidate exosomal lncRNAs were validated in 120 subjects via quantitative reverse transcription PCR (qRT-PCR). The receiver operating characteristic (ROC) curve and area under curve were used to estimate the diagnostic capacity. We investigated the potential relationship between plasma exosomal lncRNA expression and the clinicopathological parameters of GC.
RESULTS: A total of 199 exosomal lncRNAs were expressed at considerable higher levels in GCCs than those in normal controls, among which the top 10 upregulated lncRNAs were selected for further validation in cell, CM, and plasma. qRT-PCR revealed that lnc-SLC2A12-10:1 was remarkably upregulated in exosomes derived from patients with GC and GCCs. The area under the ROC curve was 0.776, which was higher than the diagnostic accuracies of CEA, CA 19-9, and CA72-4. The expression level of exosomal lnc-SLC2A12-10:1 was also significantly correlated with tumor size, TNM stage, lymph node metastasis, and degree of differentiation. The postoperative expression levels of exosomal lnc-SLC2A12-10:1 were lower compared with those of preoperative levels.
CONCLUSION: Our study suggested that exosomal lnc-SLC2A12-10:1 may be a potential noninvasive biomarker for the diagnosis and prognosis monitoring of GC. Further large-scale studies are necessary to validate its performance in GC progression.

Introduction

Gastric cancer (GC) is the fifth most common cancer and the third leading cause of cancer-related death worldwide, with its highest incidence rate recorded in East Asia.1,2 In China, GC remains the second most prevalent cancer, with approximately 679,100 new cases and 498,000 deaths in 2015.3 Most patients are diagnosed at the advanced stage, and their 5-year overall survival rates are less than 25% due to the lack of sensitive and specific early diagnosis.4 Therefore, novel and reliable GC-specific biomarkers should be developed for early GC diagnosis.

Long noncoding RNAs (lncRNAs) are a group of noncoding RNAs with a length longer than 200 nucleotides and limited protein-coding potential.5 Emerging evidences indicate that lncRNAs are involved in tumour initiation and progression through the regulation of gene expression at the epigenetic, transcriptional, and post-transcriptional levels.6,7 Aberrant lncRNAs exert effects on cancer detection, and circulating lncRNAs have been reported as potential biomarkers for the diagnosis and prognosis of different cancers.8 For example, plasma lncRNA H19 has been identified in GC and may serve as a potential biomarker for GC diagnosis.9 Additionally, different plasma lncRNAs panels can be used as biomarkers for GC detection.10,11 These data indicate that plasma lncRNAs may serve as biomarkers for GC.

Exosomes, one of the extracellular vesicles with a size in the range of 50–150 nm, are secreted by many cells and involved in intercellular communication by transmitting intracellular cargos, such as nucleic acids, lipids, and proteins, including lncRNA.12 Exosomes secreted from tumour cells can be transferred to the circulation, and exosomal RNAs are correlated with the clinicopathological features of cancer.13,14 Furthermore, exosomes can protect RNAs from degradation, and the stability and longevity of RNA molecules are ideal for noninvasive tumour diagnosis.15 Hence, we supposed that GC cells (GCCs) can generate and release exosomes, and the examination of plasma exosomal lncRNAs may provide a convenient and noninvasive method for GC diagnosis and prognosis.

In the present study, we systematically investigated the dysregulated lncRNAs in GCC-derived exosomes according to lncRNA microarray results. Data showed that lnc-SLC2A12-10:1 was significantly upregulated in GCC exosomes and plasma of GC patients. The diagnostic potential of plasma exosomal lnc-SLC2A12-10:1 was further investigated, and its relationship with the clinicopathological features of GC patients was evaluated.

Patients and Methods

Cell Culture

Four human GCC lines (MGC803, BGC823, SGC7901, and AGS) and one normal human gastric (GES-1) cell line were purchased from the Chinese Academy of Sciences Cell Bank of Type Culture Collection. All cell lines were cultured at 37 °C with 5% CO2 in 1640/DMEM (HyClone) containing 10% exosome-depleted FBS supplemented with 100 U/mL penicillin and 100 μg/mL streptomycin (Gibco). Culture media (CM) were collected, centrifuged at 2000 g for 10 min, followed by another 12,000 g for 10 min at 4 °C, and then stored at −80 °C until exosome extraction.

Patients

Plasma samples from 60 GC patients and 60 age-matched healthy controls were collected from Henan Provincial People’s Hospital between February 2018 and March 2019. All enrolled patients were newly diagnosed GC subjects. The inclusion criteria were as follows: (1) all subjects were diagnosed by gastroscopy or histopathology and (2) did not undergo any preoperative treatment. The exclusion criteria were as follows: (1) combined with other cancer patients and (2) combined with heart, kidney, and liver dysfunction. All clinicopathological data for the GC samples, including age, sex, clinical stage, and histological grade, were obtained from the clinical and pathological records. Moreover, 20 paired GC tissues and their adjacent normal tissues were obtained from resection before other therapies were performed. This study was performed in accordance with the rules of the Declaration of Helsinki of 1975 (revised in 2013) and approved by the Ethics Committee of Henan Provincial People’s Hospital. Informed consent forms were obtained from all subjects before they participated in the study. Plasma samples were separated within 2 h after collection following a two-step centrifugation protocol (3000 g for 10 min at 4 °C, 12,000 g for 10 min at 4 °C) to thoroughly remove cell debris. The separated plasma was stored in RNase-free centrifuge tubes at −80 °C until exosome extraction.

Exosome Isolation and Identification

Plasma and CM were filtered through a 0.22 μm pore polyvinylidene fluoride filter (Millipore). Subsequently, exosomes were isolated using ExoQuick solution or ExoQuick-TC solution (System Biosciences, USA) according to the manufacturer’s protocol. The exosome pellets were resuspended in 50 μL of PBS and stored at −80 °C until further analysis. The concentration and size distribution of the exosomes were determined using nanoparticle tracking analysis (NTA) on a Zeta Potential/Particle Sizer NICOMPTM 380 ZLS analyzer (Santa Barbara, California, USA) following our previous description.16

Transmission Electron Microscopy (TEM)

For TEM, 10 μL of exosome suspension was absorbed onto carbon-coated copper grids (200 mesh) for 1 min. The samples were washed with double-distilled water and negatively stained with 2% uranyl acetate solution for 1 min. After air-drying, the samples were visualized at 87,000× in a Phillips Tecnai TEM at 80 kV.

Western Blot

Briefly, the cells or exosomes were lysed in standard RIPA buffer supplemented with protease and phosphatase inhibitor cocktails (Roche). The amount of proteins was measured with a BCA protein assay kit (Beyotime, China). The protein was solubilized with loading buffer (5×) and heated at 100°C for 10 min. Proteins were separated by SDS-PAGE and then transferred to a 0.2-μm PVDF membrane (Bio- Rad, USA). After blocking with Odyssey Blocking Buffer (Li-COR Biosciences, USA), the membrane was incubated with primary antibody (1:1000) at 4°C overnight, then incubated with IRDye 800CW or 680 secondary antibodies (1:5000, LI-COR Biosciences, USA). The Odyssey Infrared Imaging System was used to visualize targeted protein bands.

lncRNA Microarray Analysis

Exosome pellets from 10 mL of CM were collected for lncRNA microarray assay. Total RNA was extracted using miRNeasy Micro Kit (Qiagen) and quantified with a NanoDrop 2000c spectrophotometer (Thermo Scientific). The quality of RNA was assessed by capillary electrophoresis on an Agilent 2100 Bioanalyzer (Agilent Technologies, CA). The lncRNA microarray analysis was performed by Shanghai Biotechnology Corporation using the Agilent human lncRNA V6 Microarray. The microarray data have been uploaded to the GEO database (GSE148334).

RNA Extraction and Quantitative Real-Time PCR (qRT-PCR)

The total RNA from tissues/cells was extracted using TRIzol reagent (Invitrogen), and exosomal RNA was extracted using a miRNeasy Micro Kit (QIAGEN) according to the manufacturers’ instructions. The quantity and quality of the RNA were evaluated using a NanoDrop spectrophotometer (Thermo Fisher Scientific). Purified RNA was reversely transcribed into cDNA using the PrimeScript™ RT reagent kit (Takara). Then, qRT-PCR was performed using SYBR Green assays (Takara) on an ABI 7900 system (Applied Biosystems). The reactions were incubated at 95 °C for 10 min, followed by 45 cycles of 95 °C for 5 s and 60 °C for 30 s. All experiments were conducted in triplicate, and the products were confirmed by melting curve analysis following each reaction. The level of each candidate lncRNA was normalized to that of GAPDH. The primers used for the lnc-SLC2A12-10:1 were (forward) 5′-CACTGAAACCAGCTTGCAT-3′ and (reverse) 5′-TCACTCCCTTGGCTATGAGG-3′. GAPDH: (forward) 5′-CTCTGCTCCTCCTGTTCGAC-3′ and (reverse) 5′-GCGCCCAATACGACCAAATC-3′. The relative expression levels of the target lncRNAs were calculated using the 2−ΔΔCT method.

Serum Tumour Marker Assay

The expression levels of CEA, CA19-9, and CA72-4 were determined according to the analysis requirements of Roche COBAS E602. The CEA, CA19-9, and CA72-4 reference intervals are 0–5.0 ng/mL, 0–35.0 U/mL, and 0–6.9 U/mL, respectively.

Statistical Analysis

All statistical analyses were performed with SPSS 20.0 and GraphPad Prism 7.0 Software. The relationships of the expression level of plasma exosomal lnc-SLC2A12-10:1 with clinicopathological parameters were evaluated by a Chi-square test. Other statistical significance between groups was determined by a Student’s t-test or one-way ANOVA. Receiver operating characteristic (ROC) curve and area under curve (AUC) were used to estimate the diagnostic value of each index for GC. A combined ROC was calculated on the basis of the logistic regression model. P value <0.05 was regarded as statistically significant.

Results

Characterization of Purified Exosomes

Exosomes are characterized by their typical morphology and conserved size, as well as the presence of specific protein markers. To ensure that exosomes were extracted successfully, we identified the exosomes isolated from the CM (Exo-CM) and plasma (Exo-P) through TEM, particle size analysis, and Western blot. The purified exosomes displayed a typical structure, with a diameter ranging from 80 nm to120 nm (Figure 1A). Conventional exosomal protein markers CD63 and TSG101 were detected in the exosome samples but not in the GCCs (Figure 1B). The expression of CD63 was relatively higher when exosomes were isolated using commercial kit. NTA showed the size distribution of exosomes, which was consistent with the range of exosomes reported in the literature (Figure 1C). The characterization of the exosome purified from cell CM was the same as that from plasma. Hence, these results confirm that exosomes were successfully purified from CM and plasma.

Figure 1

Identification of purified exosomes. (A) Representative transmission electron microscopy image of purified exosomes from culture medium (Exo-CM) and plasma (Exo-P) (scale bar, 100 nm). (B) Western blot analysis of CD63 and TSG101 in GC cell and exosomes. (C) NTA of the size distribution and number of exosomes.

lncRNA Microarray Assay for Identification of GC-Specific Exosomal lncRNA

Exosomal RNA from the CM of four GCCs and GES-1 were purified, and lncRNA microarray assay was conducted for the identification of GC-specific exosomal lncRNAs. After the combined analysis, 199 exosomal lncRNAs were found to be upregulated (fold change of ≥2 and P ≤ 0.05) in CM from GCCs compared with those in GES-1 (Figure 2A). The heatmap showed the different expression levels of the 199 upregulated lncRNAs, and we preferentially selected the top 10 upregulated lncRNAs for further verification (Figure 2B).

Figure 2

Screening of upregulated GC-specific exosomal lncRNAs according to lncRNA microarray assay. (A) A total of 199 lncRNAs were upregulated in CM-derived exosome from four GCCs compared with GES-1. (B) Heatmap result of upregulated lncRNAs based on lncRNA microarray analysis.

Verification of Upregulated Exosomal lncRNAs in GC Specimens and GCCs

To confirm the robustness of the lncRNA microarray assay, we selected the top 10 upregulated lncRNAs for further validation in cell and exosomes from CM and plasma. qRT-PCR showed that the top 10 upregulated lncRNAs increased in GCCs and GCC CM-derived exosome, which was consistent with the microarray results (data not shown). Among these upregulated lncRNAs, one particular lncRNA (lnc-SLC2A12-10:1) was significantly upregulated in the plasma exosomes and tissues of the patients with GC compared with controls (Figure 3). These results suggested that exosomal lnc-SLC2A12-10:1 may be a potential noninvasive biomarker for GC diagnosis.

Figure 3

Relative expression levels of lnc-SLC2A12-10:1 in different samples. (A) Expression levels of lnc-SLC2A12-10:1 in GC cell and control GES cell. (B) lnc-SLC2A12-10:1 expression levels in exosomes derived from GC cell and control GES cell. (C) lnc-SLC2A12-10:1 expression levels in plasma exosome of GC patients (n=60) and healthy controls (n=60). (D) lnc-SLC2A12-10:1 expression levels in 20 paired GC tissues and corresponding adjacent normal tissues.

Expression Pattern and Diagnostic Accuracy of Plasma Exosomal lnc-SLC2A12-10:1 in Patients with GC

To further investigate whether lnc-SLC2A12-10:1 can be used as a biomarker for GC diagnosis, we analyzed the lnc-SLC2A12-10:1 expression levels in 60 patients with GC and 60 healthy controls by qRT-PCR. The exosomal lnc-SLC2A12-10:1 expression level was significantly upregulated in GC patients compared with healthy controls. The ROC curve and AUC were also generated to evaluate the diagnostic accuracy of lnc-SLC2A12-10:1 and compared with the conventional tumour marker (CEA, CA19-9, CA72-4). The AUC value of exosomal lnc-SLC2A12-10:1 was 0.776 (95% confidence interval [CI]: 0.694–0.858, P < 0.001), which was significantly higher than those of CEA, CA 19–9, and CA 72–4, of 0.677 (95% CI: 0.582–0.772, P < 0.001), 0.660 (95% CI: 0.563–0.758, P < 0.01), and 0.633 (95% CI: 0.533–0.732, P < 0.05), respectively (Figure 4A). At the optimal cut-off value, the sensitivity and specificity of exosomal lnc-SLC2A12-10:1 were 78.3% and 75.0%, respectively. The AUC value increased to 0.851 (95% CI 0.780–0.922, P < 0.001) when combined with CEA, CA 19–9, and CA 72–4. The combinative diagnostic capability was better than that of each marker alone (Figure 4B). These results showed that exosomal lnc-SLC2A12-10:1 may be an appropriate diagnostic marker for GC.

Figure 4

ROC curve analysis for predicting lnc-SLC2A12-10:1 as a GC diagnosis biomarker. (A) ROC curve analysis of plasma CEA, CA19-9, CA72-4, and exosomal lnc-SLC2A12-10:1 in GC patients versus healthy controls. (B) ROC curve analysis of the combination of CEA, CA19-9, or CA72-4 and exosomal lnc-SLC2A12-10:1 in GC patients versus healthy controls, respectively.

Correlation of Plasma Exosomal lnc-SLC2A12-10:1 with Clinicopathological Characteristics

The association between exosomal lnc-SLC2A12-10:1 expression levels and the clinicopathological features of GC patients was also analysed. A total of 60 GC patients were classified into relatively high and low expression groups according to the cut-off value of the ROC curve. Statistical analysis showed that lnc-SLC2A12-10:1 expression levels were positively associated with tumour size (P = 0.035), TNM stage (P = 0.018), lymph node metastasis (P = 0.032), and differentiation (P = 0.017) of GC. Conversely, the correlation of the lnc-SLC2A12-10:1 expression with other clinical features, such as age, gender, CEA, CA19-9, and CA72-4, was not significant (Table 1).

Table 1

Correlation Between lnc-SLC2A12-10:1 Expression Levels and Clinicopathologic Features of GC Patients

Characteristics	n	High Expression	Low Expression	p-value
Gender				0.884
Male	38	20	18
Female	22	12	10
Age, years				0.782
<60	20	11	9
≥60	40	20	20
Tumor size, cm				0.035*
<5	24	8	16
≥5	36	23	13
Differentiation				0.017*
Poor +Moderate	38	27	11
Well	22	8	14
TNM stage				0.018*
T1+T2	25	9	16
T3+T4	35	24	11
Lymph node metastasis				0.032*
No	28	11	17
Yes	32	22	10
CEA, ng/mL				0.981
<5	39	21	18
≥5	21	12	9
CA19-9, U/mL				0.794
<35	36	18	18
≥35	24	13	11
CA72-4, U/mL				0.653
<6.9	41	19	22
≥6.9	19	10	9

Notes: Statistical analyses were carried out using Pearson χ2 test. *p < 0.05.

Abbreviations: CEA, carcinoembryonic antigen; CA19-9, carbohydrate antigen19-9; CA72-4, carbohydrate antigen 72–4.

Dynamic Monitoring of Plasma Exosomal lnc-SLC2A12-10:1 in GC Patients

To determine whether the expression levels of lnc-SLC2A12-10:1 are correlated with tumour load and revert to a normal state after surgery, we investigated the differences among the lnc-SLC2A12-10:1 levels of 12 paired preoperative/postoperative cases. The lnc-SLC2A12-10:1 expression levels significantly decreased 10 days after surgery (P < 0.05, Figure 5). The results indicated that exosomal lnc-SLC2A12-10:1 may be used for dynamic monitoring during the treatment of GC patients.

Figure 5

lnc-SLC2A12-10:1 expression levels in 12 paired GC patients before and after surgery.

Discussion

Early diagnosis and prompt treatment can substantially improve the overall survival rate of GC patients.17 At present, the main methods of GC diagnosis include X-ray, gastroscopy, CT, and laboratory examinations. Endoscopy followed by pathological examination is the gold standard for GC, but it is invasive and uncomfortable, leading to a limited use of widespread screening in China.18 Biomarkers in plasma and serum provide a convenient and noninvasive method for tumour diagnosis. However, the current serum tumour biomarkers used for GC diagnosis, such as CEA, CA 19–9, and CA 72–4, exhibit low sensitivity and specificity, particularly at the early stage of cancer.19 As a promising alternative to liquid biopsy, circulating tumour-derived exosomes have gained interest in the fields of noninvasive cancer diagnosis and treatment response monitoring.20

Tumour-derived ncRNAs, including lncRNAs, can be secreted into body fluids, such as blood, urine, milk, and saliva via exosomes.13,21 Exosomal lncRNAs are more stable and concentrated than general lncRNAs in peripheral circulation due to the protection of exosomal membranaceous structures.22 Tumour-derived exosomal RNAs can also accurately reflect the changes in cancer cells during tumour progression.23 Thus, circulating exosomal lncRNAs have emerged as reliable and efficient hallmarks for cancer diagnosis and prognosis monitoring.

In recent years, numerous studies have demonstrated that circulating lncRNAs could function as biomarkers in several cancers, including GC. For instance, serum exosomal lncRNA CRNDE-h can effectively distinguish CRC patients from benign colorectal diseases and NC subjects, with a significantly high AUC value of 0.892, and is significantly correlated with aggressive tumour behavior and poor prognosis.15 A recent study identified three upregulated serum exosomal lncRNAs (PCAT‐1, UBC1, and SNHG16) and designed a three-lncRNA panel, which offers an easy and fast noninvasive approach for the diagnosis and recurrence prediction of bladder cancer.24 In addition, serum exosomal lncRNAs ENSG00000258332.1 and LINC00635 were used for hepatocellular carcinoma diagnosis and prognosis.25 For GC research, serum exosomal lncRNA HOTTIP is typically increased in GC, and exosomal HOTTIP overexpression is an independent prognostic factor in GC patients.26 Early GC-specific exosomal lncRNAs (lncUEGC1) is identified through RNA sequencing and have been proved a convincing biomarker for early-stage GC.27 Consistent with previous studies, our study first identified GC-specific exosomal lncRNA through lncRNA microarray assay. Then, we primitively validated that exosomal lnc-SLC2A12-10:1 was significantly upregulated in the plasma of GC patients and had a certain diagnostic value. Cell CM was selected instead of direct plasma samples because acquiring a considerable amount of total RNA that can be used for microarray is difficult due to the limited exosomes in the plasma, and the exosomal lncRNAs derived from GCCs may be exclusive to GC.

Among the top 10 upregulated lncRNAs, lnc-SLC2A12-10:1 was the uniquely upregulated lncRNA in the plasma-derived exosome of GC patients. Thus, we focused on this special exosomal lncRNA in our subsequent work. We further evaluated the diagnostic performance of exosomal lnc-SLC2A12-10:1 for GC detection. Plasma exosomal lnc-SLC2A12-10:1 exhibited an AUC value of 0.776 in distinguishing GC from healthy volunteers in our cohort of subjects, which was better than those of traditional serum tumour markers (CEA, CA 19–9, and CA 72–4). Furthermore, we found that the combination of CEA, CA 19–9, CA 72–4, and lnc-SLC2A12-10:1 would more accurately discriminate between cases and controls. The AUC of the combination can be increased to 0.851, and the diagnostic efficacy was similar to those of previous studies.27,28 In the present study, the expression level of exosomal lnc-SLC2A12-10:1 was also significantly correlated with tumour size, TNM stage, lymph node metastasis, and differentiation. The exosomal lnc-SLC2A12-10:1 expression level significantly decreased 10 days after surgery compared with before surgery. Despite the relatively small sample size, the results of this study suggested that exosomal lnc-SLC2A12-10:1 may be a useful prognostic indicator.

lnc-SLC2A12-10:1 is located on chromosome 6, with a length of 701 nt. To the best of our knowledge, lnc-SLC2A12-10:1 has not been investigated before. Other studies have demonstrated that lncRNAs are involved in tumour initiation and progression through gene expression regulation.6 Understanding the function of plasma lncRNAs in tumorigenesis would facilitate their clinical application. Our preliminary tests showed that lnc-SLC2A12-10:1 may be a carcinogenic factor that can promote GCC proliferation and invasion in vitro (data not shown).

Taken together, our study suggested that exosomal lnc-SLC2A12-10:1 may be a potential noninvasive biomarker for GC diagnosis and prognosis monitoring. Of course, we recognized that this is a preliminary report on the clinical values of exosomal lnc-SLC2A12-10:1 in GC. Further studies with large-scale clinical samples, dynamic monitoring, and long-term follow up are necessary to support our conclusions. Functional analysis is also required to elucidate the possible regulating mechanism of exosomal lncRNAs in GC progression.