- 1Department of Laboratory Medicine, Affiliated Hospital of Nantong University, Nantong, China
- 2Research Center of Clinical Medicine, Affiliated Hospital of Nantong University, Nantong, China
- 3Medical School of Nantong University, Nantong University, Nantong, China
- 4Department of Medical Ultrasonics, Affiliated Hospital of Nantong University, Nantong, China
Background: It has been reported that long non-coding RNAs (lncRNAs) can be regarded as a biomarker and had particular clinical significance for early screening and gastric cancer (GC) diagnosis. Therefore, this study aimed to investigate whether serum HCP5 could be a new diagnostic biomarker.
Methods: Filtered out the HCP5 from the GEO database. The specificity of HCP5 was verified by real-time fluorescence quantitative PCR (qRT-PCR), and then the stability of HCP5 was verified by room temperature storage and repeated freeze-thaw experiments. Meanwhile, the accuracy of HCP5 was verified by agarose gel electrophoresis (AGE) and Sanger sequencing. Simultaneously, the expression level of serum HCP5 was detected by qRT-PCR in 98 patients with primary gastric cancer, 21 gastritis patients, 82 healthy donors, and multiple cancer types. Then, the methodology analysis was carried on. Moreover, receiver operating characteristic (ROC) was used to evaluate its diagnostic efficiency.
Results: qRT-PCR method had good repeatability and stability in detecting HCP5. The expression level of HCP5 in the serum of gastric cancer patients was remarkably higher than that of healthy controls, and it could distinguish gastritis patients from healthy donors. Besides, the expression of HCP5 was increased dramatically in MKN-45 and MGC-803. The FISH assay showed that HCP5 was mainly distributed in the cytoplasm of MKN-45 and BGC-823 cells. When HCP5 was combined with existing tumor markers, the diagnostic efficiency of HCP5 was the best, and the combined diagnosis of carcinoembryonic antigen (CEA), carbohydrate antigen199 (CA199), and HCP5 can significantly improve the diagnostic sensitivity. Besides, compared with the expression levels of thyroid cancer (THCA), colorectal cancer (CRC), and breast cancer (BRCA), serum HCP5 in gastric cancer was the most specific. Moreover, the high expression of serum HCP5 was related to differentiation, lymph node metastasis, and nerve invasion. The term of serum HCP5 after the operation was significantly lower than that of patients with primary gastric cancer.
Conclusion: Serum HCP5 can be used as a potential biomarker of non-invasive fluid biopsy, which had a unique value in the early diagnosis, development, and prognosis of gastric cancer.
Introduction
Gastric cancer (GC) is one of the most common malignant tumors globally, with high mortality and morbidity (1). As shown by data, 50% of gastrointestinal cancer cases have occurred in China, of which gastric cancer accounts for the most significant proportion, and the 5-year overall survival rate is meager, which was reported <35% in 2013-2015 (2). Among several factors, helicobacter pylori infection was the most common one, including other high-risk factors like genetic susceptibility, diet, drinking, smoking, etc. (3). At this stage, gastric cancer treatment mainly relies on surgical techniques, traditional radiotherapy, chemotherapy, and neoadjuvant therapies. However, since GC is rarely diagnosed early, it has usually advanced to most patients’ late stages. Furthermore, traditional surgery, radiotherapy, and chemotherapy are often not as effective as early-stage disease, malignant invasion, and metastasis later appeared (4). That means it is urgent to find a screening index with high specificity and sensitivity. Studies have shown that molecular analysis of non-invasive body fluids may help diagnose GC (5). In recent years, more and more literature has demonstrated that long non-coding RNA (lncRNA) played an essential role in the pathogenesis and process of GC (6–8). Hence, lncRNA has become the focus of attention in recent years.
LncRNA is a non-coding RNA with a length of more than 200nt, lacking the potential of coding peptides (9). Increasing evidence has shown that lncRNA can regulate genes by directly binding to genes, participating in translation inhibition, splicing modification, and messenger RNA (mRNA) degradation, and using microRNA (miRNA) as a competitive endogenous RNA (ceRNA) to prevent mRNA degradation and stabilize mRNA (10, 11). It plays a vital role in the regulation of splicing, transcription control, and post-transcriptional processing (12). Moreover, lncRNA may act as tumor suppressor genes and oncogenes, affecting cell proliferation (13), apoptosis (14), differentiation (15), metastasis (16), DNA damage (17), angiogenesis (18), and immune response (19), and so on. Simultaneously, the existing literature has confirmed that lncRNAs also exist in serum, plasma, urine, and exosomes of various cancer types, indicating that they can also be used as classic markers for liquid biopsy. For example, FAM83H-AS1 and lncRNA-ATB were significantly high in serum breast cancer patients, showing the potential ability to monitor breast cancer progression and staging (20). Besides, lncRNA GAS5 resulted in a specific reduction of plasma expression in coronary artery disease (CAD) patients, which could be used as a specific biomarker for diagnosing CAD (21). What’s more, the up-regulation of LncRNA XIST transported by serum extracellular vesicles was related to the progression of colorectal cancer (CRC) (22). It was also reported HYMA1, OTX2-AS1 linc00477, and loc100506688 in urine exosomes, which were found in bladder cancer, could potentially serve as biomarkers and therapeutic targets (23). Consequently, lncRNAs could be used as molecular markers for tumor diagnosis and new targets for tumor treatment shortly (24).
Histocompatibility leukocyte antigen complex P5 (HCP5) is a momentous lncRNA located between MHC class I polypeptide related sequence A (MICA) and MHC class I polypeptide related sequence B (MICB) genes in the MHC class I chain-related gene (MHCI) region, involving many autoimmune diseases and malignant tumors (25). Moreover, the abnormal expression of HCP5 is closely related to cell proliferation, migration, invasion, apoptosis, lymphatic metastasis, and drug resistance in various cancers (26–29). Hence, HCP5 is expected to become a biomarker and therapeutic target in multiple cancers (25). Functional tests verified that HCP5 could function as an oncogene in osteosarcoma and was activated by the transcription factor 1 (SP1), convincing that SP1 induced the up-regulation of HCP5, which impacted the development of osteosarcoma (30). Moreover, Wu et al. (31) observed that overexpression of HCP5 could interact with miR-3619-5p, which upregulated peroxisome proliferative activated receptor gamma (PPARG) coactivator 1 alpha (PPARGC1A) through miR-3619-5p/the AMP-activated protein kinase (AMPK)/transcription complex peroxisome proliferator-activated receptor (PPAR) coactivator-1α (PGC1α)/enhancer-binding protein beta (CEBPB) axis, regulating 5-Fluorouracil (5-Fu) and oxaliplatin resistance in gastric cancer. This study surveyed 98 GC patients, 82 healthy donors, 21 gastritis patients, and 56 patients from diverse cancers. A diagnostic model consisting of serum HCP5 and existing indicators carcinoembryonic antigen (CEA) and carbohydrate antigen199 (CA199) was established to improve diagnostic efficacy. Besides, the role of HCP5 in the dynamic monitoring of tumors in GC patients has also been discovered.
Methods
Serum Samples Collection
From January 2016 to January 2021, Serum samples of 98 GC patients, 21 gastritis patients, 19 thyroid cancer (THCA) patients, 20 colorectal cancer (CRC) patients, 17 breast cancer (BRCA) patients, and 82 healthy donors were collected at the clinical laboratory of Affiliated Hospital of Nantong University (Nantong, China). All serum samples were stored in RNase-free test tubes at -80°C for later use. All samples mentioned above were collected under the ethics of the World Medical Association, and informed consent was obtained for experiments on human subjects. Besides, this study was approved by the Human Research Ethics Committee of Nantong University Hospital. (Ethics review report number: 2018-L055).
RNA Extraction
According to the serum extraction kit (BioTeke, Wuxi, China), total RNA was extracted from 300µl serum. Then, 900µl lysis buffer was added, mixed by pipetting, vortexed for 1min, and incubated at 15-30°C for 5min. Adding 180μl chloroform in them and let them stand. After then, the serum samples were centrifuged at 12000rpm for 10min in a 4°C low-temperature high-speed centrifuge, followed by colorless water phases and organic phases. The colorless aqueous phase was then absorbed and transferred into a 1.5ml RNase-free tube, added to 700μl 70% ethanol, mixed upside down, and moved to the adsorption column. Centrifuged at 12000rpm for 1min and discarded the liquid. They were then adding 500μl deproteinized liquid, centrifuged at 12000rpm for 1min. After cleaning with 700μl rinsing solution, leaving for 2min and the liquid was transferred to a new 1.5ml RNase-free tube. Then, added 30μl RNase-free water in a warm bath, standing at room temperature for 2min, centrifuged at 12000rpm, and twice in 1min. After centrifugation, put it in the refrigerator at -80°C for later use.
Complementary DNA (cDNA) Synthesis and Real-Time Fluorescence Quantitative Polymerase Chain Reaction (qRT-PCR)
The reverse transcription kit (Thermo Fisher Scientific) was used, adding 4µl 5×Reaction Buffer, 2µl 10mM deoxyribonucleotide triphosphate (dNTPs), 1µl Oligo (dT) Primer, 1µl RNase inhibitor (20U/µl), and 1µl Reverse Transcriptase (200U/µl) to form the reaction system, which was reverse-transcribed into cDNA with the total 11µl RNA solution. Then, we incubated the 20µl total system at 25°C for 5 minutes, 42°C for 60 minutes, 70°C for 5 minutes and stored the synthesized cDNA at -80°C when not in use. After that, qRT-PCR assays were performed on the Roche Lightcycler 480 (Roche, Switzerland), using SYBR Green I mix (Roche) as the fluorescent dyes for a total of 20µl system. The primers used in this article were all synthesized by RiboBio Corporation (Suzhou, China). The primer sequences are as follows: HCP5, 5’-TGAGAGCAGGACAGGAAAA-3’ (forward) and 5’-CCAACCAGACCCTAAGTGA-3’ (reverse); 18S ribosomal RNA (18S): 5’-CGCTCGCTCC TCTCCTACTT-3’ (forward) and 5’-CGGGTTGGTTTTGATCTGATAA-3’ (reverse).
The PCR cycling program included that activating the enzyme at 95°C for 10min, denaturation at 95°C for 15s, annealing at 60°C for 30s, then collecting fluorescence information at 80°C for a total of 45 cycles. The 2-ΔΔCt method was used to calculate the relative expression of the target gene. ΔΔCt means the difference between the experimental group and the control group.
Cell Culture
Human gastric epithelial cell (GES-1) and gastric cancer cell line (MKN-45, MGC-803, BGC-823, HGC-27) were purchased from the Cell Bank of the Chinese Academy of Sciences (Shanghai, China). Cells were cultured in RPMI-1640 (Corning, USA) containing 10% fetal bovine serum (FBS) (Gibco, USA) and 1% penicillin-streptomycin in a humid condition of 5% CO2 at 37°C.
Cell Secretion Assay
The 1, 3, 5, and 7 days cell supernatants of MKN-45, MGC-803, and GES-1 were extracted. Furthermore, the expression level of lncRNA was later detected.
Linear Verification Experiment
To verify the linearity of lncRNA, a healthy donor’s serum cDNA was diluted by 10, 102, and 103 times, respectively.
Stability Verification Experiment
The mixed serum samples were placed at room temperature for 0, 6, 12, 18, and 24 h. Besides, repeated freezing and thawing for 0, 1, 3, 5, and 10 times to verify the stability of lncRNA.
Agarose Gel Electrophoresis (AGE)
50×TAE 1:49 diluted to 1×, configured 2% gels, and waited for use. 1µl the loading buffer and 5µl products were added, then electrophoresed at 110-120V for about 30min. The accuracy of the product was verified by sequencing analysis.
Cell Passage
Firstly, the original culture medium was discarded and washed 1-2 times with phosphate buffer saline (PBS). Secondly, cells in the T25 culture flasks were digested for 1min by adding about 1ml 0.25% trypsin-EDTA (Gibco, USA). Next, adding 2-3ml of the complete medium to culture flasks to stop digestion, rinsed, and suspended until all the cells were flushed down. They were then transferred to 15ml centrifuge tubes and centrifuged at 1000 rpm for 5min. Adding about 3ml complete medium to infiltrate the bottom of T25 culture flasks. While centrifugation was finished, we discarded the supernatant, adding about 4ml complete medium to suspend and evenly 2ml suspension to each bottle, mixed well, and then incubated in the condition with 5% CO2 at 37°C.
Cell Cryopreservation
The original culture medium was discarded, rinsed 1-2 times with PBS, added approximately 1ml 0.25% trypsin-EDTA to digest cells in the T25 culture bottle for 1min. 2-3ml of complete medium was added to terminate digestion, transferred to a 15ml centrifuge tube and centrifuged for 5min at 1000rpm. Then, the suspension was transferred to the freezer tube and quickly placed in a refrigerator at -80°C for later use.
Fluorescence In Situ Hybridization (FISH)
According to instructions of the manufacturer, the fluorescence in situ hybridization kit (RiboBio, Guangzhou, China) was used to repair and penetrate cells. 200μl prehybridization solution was incubated at 37°C for 30min, followed by 20μmol FISH probes/hybridization buffer for overnight incubation at 37°C. Rinsed each well with washing solution I, II, III and PBS in turn. 4’,6-Diamidino-2’-phenylindole (DAPI) was then dyed in the dark for 10min, and photos were taken under a fluorescence microscope. The HCP5 Fish probe was designed and synthesized by RiboBio (RiboBio, Guangzhou, China).
Data Analysis
All statistics were analyzed using SPSS 20.0 software (IBM SPSS Statistics, Chicago, Illinois, USA) and GraphPad Prism V.8.00 software (GraphPad Software, La Jolla, USA). Meanwhile, all data were expressed as the mean ± standard deviation (SD) of values obtained in three independent experiments. According to the actual situation, independent sample t-test, one-way analysis of variance (ANOVA), and χ2 test were used for statistical analysis. The receiver operating characteristic curve (ROC) and the area under the ROC curve (AUC) were used to evaluate the diagnostic value of HCP5. Besides, bivariate logistic regression was used to analyze the diagnostic value of HCP5, CEA, and CA199. The ROC curve and AUC were obtained through non-parametric analysis. Youden index (Youden index = sensitivity + specificity - 1) was used as the expression value of serum HCP5. The overall survival (OS) was estimated by using the Kaplan-Meier method. All values of p<0.05 were considered to reach a statistically significant difference.
Results
The Innovative Detection Method of Serum HCP5
LncRNA, as a classic non-coding RNA, played a vital role in liquid biopsies as a biomarker. The literature confirmed that 18S had good linearity and stability compared with the reference genes glyceraldehyde 3-phosphate dehydrogenase, U6 small nuclear RNA (U6), β-actin (ACTB), and tubulin (TUB) in gastric cancer (32). Therefore, 18S was selected as the endogenous reference gene in this assay. To verify whether HCP5 was suitable for clinical laboratory analysis, we performed a methodological assessment of HCP5 detected by the qRT-PCR method. To verify its linearity, we built a continuous 10-fold dilution of cDNA and diluted the original serum cDNA concentrations 10, 102 and 103 times. Different concentrations of serum HCP5 were detected by qRT-PCR, and data analysis was performed (Figure 1A). At the same time, 18S was analyzed, as shown in Figure 1B. It could be seen that R2 of the standard curve of HCP5 was 0.960, and the regression equation was y=-0.901x+23.57, indicating that qRT-PCR was an effective method for detecting different concentrations of serum HCP5. In addition, we selected the mixture of serums for HCP5 precision determination. The results showed that the inter-and intra-assay coefficients of variation (CV) for HCP5 were 2.99% and 2.11%, respectively (Table 1). Secondly, to verify the stability of HCP5 in liquid biopsies, we placed the mixed standard human serum samples at room temperatures of 0, 6, 12, 18, and 24h. Similarly, we performed frozen-thawed cycles for another set of mixed samples 0, 1, 3, 5, and 10 times. All of these samples were stored in RNase-free tubes and kept at -80°C for later use. As shown in Figures 1C, D, when environmental conditions changed, the level of serum HCP5 remained significantly unchanged, indicating that HCP5 had good stability. Besides, the smooth unimodal melting curve also indicated the high specificity of amplified product serum HCP5 (Figure 1E). By analyzing the PCR products through AGE, the result showed that the band was 122bp and the accuracy was verified (Figure 1F). Sanger sequencing showed that the PCR product sequence was consistent with HCP5, with a sequence size of 122bp, further verifying the accuracy of the PCR method (Figure 1G). In this section, from all of these results, we found that the detection of HCP5 by qRT-PCR was susceptible and specific.
Figure 1 Methodology evaluation of HCP5 in GC serum samples. (A, B) Standard curves of serum HCP5 and 18S in a ten-fold serial dilution to show the linearity. (C, D) Stability of HCP5 under room temperature incubation time or multiple freeze-thaw cycles. Data were presented as raw Ct value (n=3). (E) The specificity of PCR products by the melting curve. (F) The validation of PCR products by the agarose gel electrophoresis. (G) The product sequence was verified by Sanger sequencing.
The Expression Level and Diagnostic Efficacy of Serum HCP5 in Gastric Cancer
To verify whether HCP5 can be used as an emerging GC diagnostic marker. We detected the expression of serum HCP5 in 98 GC patients and 82 healthy donors by qRT-PCR. The results showed that the expression level of serum HCP5 in GC patients was significantly higher than that in healthy donors (Figure 2A). In addition, HCP5 can clearly distinguish patients with gastritis and gastric cancer (Figure 2B). Besides, CEA and CA199 are usually widely used for screening and auxiliary diagnosis of GC. As shown in Figures 2C, D, the expression levels of serum CEA and CA199 in GC patients were also higher than those in healthy controls. The ROC curves showed that the area under the curve of serum HCP5 detected by qRT-PCR was 0.818 (95% CI:0.757-0.880 P<0.001) (Figure 2E). Under this cut-off value, the sensitivity was 80%, and the specificity was 70% (Table 2). Moreover, the clinical indicators CEA and CA199 area under the curve were 0.725 (95% CI:0.650-0.799, P<0.001) and 0.687 (95% CI:0.610-0.764, P<0.001), respectively (Figures 2F, G). Besides, we also analyzed the expression level of HCP5 and the comprehensive diagnostic analysis between CEA and CA199. As mentioned above, compared with the separate detection, the three combined diagnoses provided the high AUC 0.870 (95% CI:0.819-0.921, P<0.001) in the distinction between gastric cancer and healthy donors (Figure 2H). Then, through logistic regression analysis, the diagnostic efficiency of several serum biomarkers was discussed. Besides, compared with single detection, combined detection can better distinguish GC patients from healthy donors, especially when compared with the combined diagnosis model of HCP5, CEA, and CA199, the combination of which can increase sensitivity and specificity, which could reach 81% and 79%, respectively (Table 2). The above data showed that the single diagnosis of HCP5 had a significant AUC, and the combined diagnosis of HCP5, CEA, and CA199 can improve the diagnosis efficiency of GC.
Figure 2 Concentrations of serum HCP5, CEA, and CA199 in GC cases and diagnostic efficacy. (A) Detection of serum HCP5 in GC patients (n=98), healthy donors (n=82). (B) Detection of serum HCP5 in GC (n=98) and gastritis (n=21) patients. (C) The concentrations of serum CEA. (D) The concentrations of serum CA199. *P < 0.05, **P < 0.01, ****P < 0.001, NS means no statistically difference. (E) The ROC curve of serum HCP5 in GC cases. (F) The ROC curve of serum CEA in GC cases. (G) The ROC curve of serum CA199 in GC cases. (H) The diagnostic efficacy of the combined diagnosis.
Table 2 Combination of serum HCP5, CEA and CA199 levels significantly improves the diagnostic sensitivity between GC patients and healthy donors.
Origin of HCP5 in GC Samples
To understand the overall expression level of HCP5 in GC, we collected 20 pairs of tissue samples. The result showed that HCP5 was significantly highly expressed in GC tissues (Figure 3A). Because of the high expression of HCP5 in serum GC, we assumed whether HCP5 was secreted by GC cells. Then, we detected the expression level of HCP5 in two GC cell lines (MKN-45, MGC-803) and normal control cells GES-1 for 1, 3, 5, and 7 days, respectively (Figure 3B). The results showed that compared with normal cell lines, HCP5 increased in MKN-45 and MGC-803 over time, especially in MKN-45. Hence, we speculated that the high expression of HCP5 in peripheral blood might be derived from the secretion of some tumor cells. Furthermore, we detected the subcellular localization of lncRNA HCP5, and the results showed that HCP5 was mainly distributed in the cytoplasm of MKN-45 and BGC-823 cells (Figure 3C).
Figure 3 The origin of serum HCP5 in GC cases. (A) The expression level of HCP5 in 20 pairs of GC tissues. (B) HCP5 was secreted into the culture medium by MKN-45 and MGC-803 cells in a time-dependent manner. *P < 0.05, **P < 0.01, ****P < 0.001. (C) The FISH assay of HCP5 in MKN-45 and BGC-823 cells. Scale bars: 20μm.
Diagnostic Efficacy of Serum HCP5 in Various Cancer Types and Distinguishing Between Gastritis and Healthy Donors
To verify the specificity of HCP5 expression levels, serum expression levels of 19 thyroid cancer patients, 20 colorectal cancer patients, and 17 breast cancer patients were examined. The results showed that HCP5 expression was not statistically significant not only in thyroid cancer but also in colorectal cancer, and in breast cancer had a little significance (Figures 4A–C). At the same time, we detected expression levels of serum HCP5, CEA, and CA199 in 21 patients with gastritis patients and 21 healthy donors, respectively (Figures 4D–F). Collectively, the results indicated that serum HCP5 could distinguish gastritis patients from healthy donors. Thus, it can be seen HCP5 can be regarded as a specific biomarker for GC diagnosis.
Figure 4 The specificity of serum HCP5 in GC and differentiating from gastritis and healthy donors. (A) Detection of serum HCP5 in 19 thyroid cancer patients. (B) Detection of serum HCP5 in 20 colorectal cancer patients. (C) Detection of serum HCP5 in 17 breast cancer patients. (D) The expression levels of HCP5 in gastritis patients (n=21) and healthy donors (n=21). (E) The concentrations of CEA in gastritis patients (n=21) and healthy donors (n=21). (F) The concentrations of CA199 in gastritis patients (n=21) and healthy donors (n=21). *P < 0.05, **P < 0.01, NS means no statistically difference.
Correlation Between Serum HCP5 Expression and GC Clinicopathological Parameters
The clinicopathological parameters of 98 patients with GC were summarized in Table 3. Based on the cut-off value of serum HCP5, 98 GC patients were divided into the high-value and low-value groups. Chi-square test showed that although we did not find any correlation with serum HCP5 expression in pathological parameters such as gender, age, tumor size, TNM stage, distant metastasis, and CEA. Nevertheless, the high expression of serum HCP5 was correlated with differentiation (P<0.05), lymph node metastasis (P<0.05), and nerve invasion (P<0.05), indicating that serum HCP5 expression was correlated with some GC clinicopathological parameters.
Tumor Dynamic Monitoring of Serum HCP5 in GC Patients
Due to previous reports and our current research, circulating lncRNAs may be secreted by tumor cells, and they would return to an average level after surgery. To verify whether the expression of serum HCP5 was related to tumor dynamic monitoring, we compared serum HCP5 in unpaired samples, including 98 patients with primary gastric cancer, 46 patients with surgical treatment, and 57 patients with tumor recurrence. It was found that the serum HCP5 expression of patients after treatment was significantly lower than that of GC patients and relapsed patients before treatment (Figure 5A). Meanwhile, we investigated the differences in serum HCP5 levels in 15 pairs of GC patients before and after surgery. The results showed the expression level of HCP5 after gastrectomy was significantly lower than that of patients with primary GC, which showed that HCP5 could be used for dynamic monitoring (Figure 5B). Besides, the survival curve revealed that the overall survival rate of patients with high HCP5 expression was significantly lower than that of patients with low HCP5 expression (P=0.0079) (Figure 5C). To sum up, we believed that serum HCP5 could be regarded as a new marker for early diagnosis of GC and tumor dynamic monitoring.
Figure 5 Serum HCP5 in monitoring tumor dynamics in GC patients. (A) Detection of serum HCP5 expression in GC pre-operation patients (n=98), post-operation patients (n=46) and recurrence patients (n=57). (B) Altered expression of serum HCP5 in 15 paired samples preoperatively and postoperatively. (C) The survival curve of patients with GC. *P < 0.05, **P < 0.01, ****P < 0.001, NS means no statistically difference.
Discussion
Gastric cancer is the common cancer with high mortality and morbidity. Recently, the treatment effect and prediction were still poor and low. Surgical resection was still the primary treatment at present. lncRNA can be divided into sense, antisense, intronic, bidirectional transcripts, intergenic and enhancer RNAs (33). Current studies have shown that lncRNA carried out strict regulation and abnormal expression during the development of different cancers, which can be divided into (1) chromatin regulation; (2) regulation of histone modification; (3) regulation of DNA methylation; (4) interaction with transcription factors; (5) regulation of genomic tissue through enhancer cycle; (6) post-transcriptional regulation (34). For example, kinds of literature have shown that NEAT1 affects chromatin remodeling, increased histone acetylation level, promoted aldehyde dehydrogenase 1 family (ALDH1) and c-Myc expression, and improved drug resistance and tumor dryness 5-FU (35). Zhang et al. (36) found that lnc-LALC bonded to enhancer of zeste homolog 2 (EZH2), recruited DNA methyltransferase (DNMT) to Leucine zipper putative tumor suppressor 1 (LZTS1) promoter, and changed the expression of LZTS1. Another literature reported that in ischemic heart disease, the combination of lncCIRBIL and BCL2-related transcription factor1 (Bclaf1) had a protective effect on I/R damage and was a potential target for treating ischemic heart disease (37). What’s more, it also has been convinced that the chromatin cyclization between enhancers E1 and E2 and lncRNA promoter played a co-regulatory role in prostate cancer (38). Furthermore, Wang et al. (39) found that linc00336 can act as an endogenous sponge of miRNA to inhibit lung cancer hypertrophy. Simultaneously, more and more pieces of literature have confirmed that lncRNA can be used as a potential biomarker for tumor screening and prognosis monitoring (20–24). Inspired by this idea, this study was devoted to finding suitable tumor biomarkers.
HCP5, also named P5-1, is located on chromosome 6p21.33. Previous experiments have shown that HCP5 can participate in regulating various tumors, such as HCP5, which was mediated by transforming growth factor β (TGF β) and regulated by recombinant SMAD family member 3 (SMAD3) transcription, could promote the growth and metastasis of lung adenocarcinoma (LUAD) tumors (40). Moreover, HCP5 can encourage the development of cervical cancer by inhibiting microRNA-15a regulating metastasis-associated in colon cancer-1 (MACC1) (30). The above studies illustrated that HCP5 was a regulator in the process of tumor development, and it may become a predictor of tumor diagnosis and treatment. At present, real-time fluorescence quantitative PCR (41), microarray (42), and next-generation sequencing (NGS) (43) are often used for quantitative analysis of ncRNA. qRT-PCR was considered a reliable method for detecting non-coding RNA gene expression because of its high sensitivity, speed, and efficiency, as the small sample requirements and ease of interpretation of the results were significant advantages (44). Consequently, in this study, we chose the qRT-PCR method to evaluate serum HCP5 expression, and 18S was selected as the internal reference. The results showed that the expression of HCP5 in GC patients by this method had good linearity, stability, specificity, and repeatability. Besides, large serum samples of gastric cancer detected by qRT-PCR showed that HCP5 could not only significantly distinguish between GC patients and healthy donors, but also could distinguish GC patients from gastritis patients well. The diagnostic efficiency of HCP5 alone was substantially higher than that of CEA and CA199, and the combined diagnosis of HCP5, CEA, and CA199 can improve the diagnostic efficiency of gastric cancer. Then, we detected the expression level of serum HCP5 in thyroid, colorectal, and breast cancer, found that its expression in the serum GC was the most specific. At the same time, we also detected the expression level of serum HCP5 in gastritis and normal controls. Compared with the typical clinical diagnostic indicators CEA and CA199, the results showed that serum HCP5 could better distinguish gastritis patients from healthy donors. From the above experiments, it can be seen that HCP5 showed a remarkable advantage as a specific biomarker for the diagnosis of GC. Besides, clinicopathological parameters also showed that the high expression of serum HCP5 in GC was notably correlated with differentiation (P<0.05), lymph node metastasis (P<0.05), and nerve invasion (P<0.05), indicating that serum HCP5 expression was correlated with some GC clinicopathological parameters. Moreover, the statistics of 20 pairs of GC tissues showed that HCP5 had a conspicuous trend of high expression. What’s more, previous literature reported a pronounced correlation between the expression of circulating biomarkers in blood and cells (45). To verify whether HCP5 was secreted by tumor cells, we also examined the expression of HCP5 in GES-1, MGC-803, and MKN-45. We found that the expression of HCP5 in MGC-803 and MKN-45 was dramatically increased, especially in MKN-45. It was speculated that HCP5 might come from the secretion of tumor cells and was mainly related to tumor metastasis. Meanwhile, the FISH assay showed that HCP5 was mainly present in the cytoplasm of MKN-45 and BGC-823. Besides, by comparing the expression of serum HCP5 in patients with primary gastric cancer, patients undergoing surgery, and patients with postoperative recurrence, we found that serum HCP5 rebounded in patients with tumor recurrence and it had the ability to monitor tumor dynamics. Meanwhile, we also detected the expression level of HCP5 in 15 pairs of GC patients before and after an operation. The results showed that the expression level of HCP5 after gastrectomy was strikingly lower than that in patients with primary GC, which indicated that HCP5 could be a specific biomarker and played an important role in dynamic monitoring. The survival curve showed that the OS rate of patients with high serum HCP5 was significantly lower than that of patients with low expression of serum HCP5. It showed a positive correlation between a high level of HCP5 and poor survival and prognosis of GC.
In short, our study showed that the increase of serum HCP5 could significantly distinguish between patients with primary gastric cancer and healthy controls, and the combined diagnosis of HCP5, CEA, and CA199 had high diagnostic efficiency. At present, gastritis is mainly divided into gastroscopy and laboratory examination, which included gastric juice analysis, pepsinogen test, serum gastrin test, and Helicobacter pylori test. There may be no significant difference in laboratory indexes among different types and locations of Lauren tumors, which had some limitations (46). For Helicobacter pylori testing, studies have shown that bacteria could be eliminated in atrophic gastritis development, so it also had certain limitations (47). Hence, it was indispensable to find other specific serum biomarkers combined with other preventive measures. Compared with the traditional clinical indexes, serum HCP5 can better distinguish gastritis from healthy donors. Thus, it can be concluded that serum HCP5 may be a potential biomarker in dynamic monitoring of gastric cancer and tumor. Nevertheless, there were some limitations in this study. The samples of this study were only limited to one hospital in China. If conditions permit, serum samples from patients in multiple hospitals could be collected for further study. In addition, due to the limited sample size of patients with gastritis, it was necessary to increase further gastritis samples and long-term observation of patients with gastritis to obtain new insights into the potential mechanism of the diagnostic and prognostic value of HCP5.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Ethics Statement
The studies involving human participants were reviewed and approved by the ethics committee of the Affiliated Hospital of Nantong University (ethical review report number: 2018-L055). The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.
Author Contributions
SQ designed the experiment, conducted the data analysis, made charts, and wrote the article. LY performed the experiment and conducted the data collection. SK gave some suggestions and carried out the revision of the article. YX gave advice for some charts and participated in the final revision. BL and SJ provided resources and guidance for the paper. All authors contributed to the article and approved the submitted version.
Funding
This project was supported by grants from the National Natural Science Foundation of China [grant number: 81871720, 82072363], and Jiangsu Province “science, education, and health” key discipline [grant number: ZDXKB2016011].
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2021.684531/full#supplementary-material
Supplementary Figure 1 | (A) The ROC curve of the poorly, moderately poorly groups and well, medium-well groups. (B) The ROC curve of the lymph node metastasis positive and negative groups. (C) The ROC curve of the nerve invasion positive and negative groups.
Supplementary Figure 2 | (A) Serum HCP5 expression of differentiation. (B) Serum HCP5 expression of lymph node metastasis. (C) Serum HCP5 expression of nerve invasion.
Supplementary Figure 3 | The time-dependent studies on HGC-27, BGC-823, and GES-1 cells.
Supplementary Figure 4 | (A) The ROC curve of the TIII-TIV group. (B) The ROC curve of the TI-TII group.
Glossary
GC: gastric cancer
lncRNA: long non-coding RNA
mRNA: messenger RNA
miRNA: microRNA
ceRNA: competitive endogenous RNA
CAD: coronary artery disease
CRC: colorectal cancer
HCP5: Histocompatibility leukocyte antigen complex P
MICA: MHC class I polypeptide related sequence A
MICB: MHC class I polypeptide related sequence B
MHCI: MHC class I chain-related gene
SP1: the transcription factor 1
PPARG: peroxisome proliferative activated receptor gamma
PPARGC1A: PPARG coactivator 1 alpha
AMPK: the AMP-activated protein kinase
PPAR: peroxisome proliferator-activated receptor
PGC1α: PPAR coactivator-1α
CEBPB: enhancer-binding protein β
5-Fu: 5-Fluorouracil
CEA: carcinoembryonic antigen
CA199: carbohydrate antigen199
THCA: thyroid cancer
BRCA: breast cancer
cDNA: Complementary DNA
qRT-PCR: real-time fluorescence quantitative polymerase chain reaction
dNTPs: deoxyribonucleotide triphosphate
18S: 18S ribosomal RNA
FBS: fetal bovine serum
AGE: Agarose Gel Electrophoresis
PBS: phosphate buffer saline
FISH: fluorescence in situ hybridization
DAPI: 4’’,6-Diamidino-2’-phenylindole
SD: standard deviation
ANOVA: one-way analysis of variance
ROC: receiver operating characteristic curve
AUC: area under the ROC curve
OS: overall survival
U6: U6 small nuclear RNA
ACTB: β-actin
TUB: tubulin
CV: coefficients of variation
STAD: stomach adenocarcinoma
ALDH1: aldehyde dehydrogenase 1 family
EZH2: enhancer of zeste homolog 2
DNMT: DNA methyltransferase
LZTS1: Leucine zipper putative tumor suppressor 1
Bclaf1: BCL2-related transcription factor1
TGF β: transforming growth factor β
SMAD3: recombinant SMAD family member 3
LUAD: Lung adenocarcinoma
MACC1: metastasis-associated in colon cancer-1
NGS: next-generation sequencing
References
1. Sung H, Ferlay J, Siegel R, Laversanne M, Soerjomataram I, Jemal A, et al. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA: Cancer J Clin (2021) 70:313. doi: 10.3322/caac.21660
2. Zeng H, Chen W, Zheng R, Zhang S, Ji JS, Zou X, et al. Changing Cancer Survival in China During 2003-15: A Pooled Analysis of 17 Population-Based Cancer Registries. Lancet Glob Health (2018) 6(5):e555–67. doi: 10.1016/S2214-109X(18)30127-X
3. Cai Z, Liu Q. Understanding the Global Cancer Statistics 2018: Implications for Cancer Control. Sci China Life Sci (2019) 87:162–8. doi: 10.1007/s11427-019-9816-1
4. Portuondo J, Tran Cao H, da Costa W, Sada Y, Harris A, Massarweh N. Treatment Approach, Hospital Practice Patterns, and Receipt of Multimodality Therapy as Measures of Quality for Locally Advanced Gastric Cancer. J Surg Oncol (2021) 123:1724–35. doi: 10.1002/jso.26460
5. Yamamoto H, Watanabe Y, Sato Y, Maehata T, Itoh F. Non-Invasive Early Molecular Detection of Gastric Cancers. Cancers (2020) 12(10):2880. doi: 10.3390/cancers12102880
6. Zhang F, Wang H, Yu J, Yao X, Yang S, Li W, et al. Lncrna CRNDE Attenuates Chemoresistance in Gastric Cancer Via SRSF6-Regulated Alternative Splicing of PICALM. Mol Cancer (2021) 20(1):6. doi: 10.1186/s12943-020-01299-y
7. Shuai Y, Ma Z, Liu W, Yu T, Yan C, Jiang H, et al. TEAD4 Modulated LncRNA Mnx1-AS1 Contributes to Gastric Cancer Progression Partly Through Suppressing BTG2 and Activating BCL2. Mol Cancer (2020) 19(1):6. doi: 10.1186/s12943-019-1104-1
8. Wang C, Zhu C, Xu J, Wang M, Zhao W, Liu Q, et al. The Lncrna UCA1 Promotes Proliferation, Migration, Immune Escape and Inhibits Apoptosis in Gastric Cancer by Sponging Anti-Tumor MiRNAs. Mol Cancer (2019) 18(1):115. doi: 10.1186/s12943-019-1032-0
9. Hung T, Chang H. Long Noncoding RNA in Genome Regulation: Prospects and Mechanisms. RNA Biol (2010) 7(5):582–5. doi: 10.4161/rna.7.5.13216
10. Sun W, Lv J, Duan L, Lin R, Li Y, Li S, et al. Long Noncoding RNA H19 Promotes Vascular Remodeling by Sponging let-7a to Upregulate the Expression of Cyclin D1. Biochem Biophys Res Commun (2019) 508(4):1038–42. doi: 10.1016/j.bbrc.2018.11.185
11. Svitich O, Sobolev V, Gankovskaya L, Zhigalkina P, Zverev V. The Role of Regulatory RNAs (miRNAs) in Asthma. Allergol Immunopathol (2018) 46(2):201–5. doi: 10.1016/j.aller.2017.09.015
12. Mercer T, Dinger M, Mattick J. Long Non-Coding RNAs: Insights Into Functions. Nat Rev Genet (2009) 10(3):155–9. doi: 10.1038/nrg2521
13. Kretz M, Webster D, Flockhart R, Lee C, Zehnder A, Lopez-Pajares V, et al. Suppression of Progenitor Differentiation Requires the Long Noncoding RNA Ancr. Genes Dev (2012) 26(4):338–43. doi: 10.1101/gad.182121.111
14. Zhang J, Sui S, Wu H, Zhang J, Zhang X, Xu S, et al. The Transcriptional Landscape of lncRNAs Reveals the Oncogenic Function of LINC00511 in ER-Negative Breast Cancer. Cell Death Dis (2019) 10(8):599. doi: 10.1038/s41419-019-1835-3
15. Wang P, Xue Y, Han Y, Lin L, Wu C, Xu S, et al. The STAT3-Binding Long Noncoding RNA lnc-DC Controls Human Dendritic Cell Differentiation. Science (New York NY) (2014) 344(6181):310–3. doi: 10.1126/science.1251456
16. Zhang H, Liu L, Wang Y, Zhao G, Xie R, Liu C, et al. KLF8 Involves in TGF-Beta-induced EMT and Promotes Invasion and Migration in Gastric Cancer Cells. J Cancer Res Clin Oncol (2013) 139(6):1033–42. doi: 10.1007/s00432-012-1363-3
17. Shahabi S, Kumaran V, Castillo J, Cong Z, Nandagopal G, Mullen D, et al. Linc00261 Is An Epigenetically Regulated Tumor Suppressor Essential for Activation of the DNA Damage Response. Cancer Res (2019) 79(12):3050–62. doi: 10.1158/0008-5472.Can-18-2034
18. Yu B, Wang S. Angio-Lncrs: LncRNAs That Regulate Angiogenesis and Vascular Disease. Theranostics (2018) 8(13):3654–75. doi: 10.7150/thno.26024
19. Wu M, Fu P, Qu L, Liu J, Lin A. Long Noncoding RNAs, New Critical Regulators in Cancer Immunity. Front Oncol (2020) 10:550987. doi: 10.3389/fonc.2020.550987
20. El-Ashmawy N, Hussien F, El-Feky O, Hamouda S, Al-Ashmawy G. Serum LncRNA-ATB and FAM83H-AS1 as Diagnostic/Prognostic Non-Invasive Biomarkers for Breast Cancer. Life Sci (2020) 259:118193. doi: 10.1016/j.lfs.2020.118193
21. Yin Q, Wu A, Liu M. Plasma Long Non-Coding RNA (LncRNA) GAS5 Is a New Biomarker for Coronary Artery Disease. Med Sci Monitor Int Med J Exp Clin Res (2017) 23:6042–8. doi: 10.12659/msm.907118
22. Yu J, Dong W, Liang J. Extracellular Vesicle-Transported Long Non-Coding RNA (LncRNA) X Inactive-Specific Transcript (XIST) in Serum Is a Potential Novel Biomarker for Colorectal Cancer Diagnosis. Med Sci Monitor Int Med J Exp Clin Res (2020) 26:e924448. doi: 10.12659/msm.924448
23. Berrondo C, Flax J, Kucherov V, Siebert A, Osinski T, Rosenberg A, et al. Expression of the Long Non-Coding RNA Hotair Correlates With Disease Progression in Bladder Cancer and Is Contained in Bladder Cancer Patient Urinary Exosomes. PloS One (2016) 11(1):e0147236. doi: 10.1371/journal.pone.0147236
24. Kopp F, Mendell J. Functional Classification and Experimental Dissection of Long Noncoding RNAs. Cell (2018) 172(3):393–407. doi: 10.1016/j.cell.2018.01.011
25. Zou Y, Chen B. Long Non-Coding RNA HCP5 in Cancer. Clinica Chimica acta; Int J Clin Chem (2021) 512:33–9. doi: 10.1016/j.cca.2020.11.015
26. Liang L, Xu J, Wang M, Xu G, Zhang N, Wang G, et al. Lncrna HCP5 Promotes Follicular Thyroid Carcinoma Progression Via miRNAs Sponge. Cell Death Dis (2018) 9(3):372. doi: 10.1038/s41419-018-0382-7
27. Wang L, Luan T, Zhou S, Lin J, Yang Y, Liu W, et al. Lncrna HCP5 Promotes Triple Negative Breast Cancer Progression as a ceRNA to Regulate BIRC3 by Sponging Mir-219a-5p. Cancer Med (2019) 8(9):4389–403. doi: 10.1002/cam4.2335
28. Wu J, Chen H, Ye M, Wang B, Zhang Y, Sheng J, et al. Long Noncoding RNA HCP5 Contributes to Cisplatin Resistance in Human Triple-Negative Breast Cancer Via Regulation of PTEN Expression. Biomed Pharmacother = Biomed Pharmacother (2019) 115:108869. doi: 10.1016/j.biopha.2019.108869
29. Liu Y, Wang J, Dong L, Xia L, Zhu H, Li Z, et al. Long Noncoding Rna HCP5 Regulates Pancreatic Cancer Gemcitabine (Gem) Resistance by Sponging Hsa-Mir-214-3p To Target Hdgf. OncoTargets Ther (2019) 12:8207–16. doi: 10.2147/ott.S222703
30. Zhao W, Li L. SP1-Induced Upregulation of Long Non-Coding RNA HCP5 Promotes the Development of Osteosarcoma. Pathol Res Pract (2019) 215(3):439–45. doi: 10.1016/j.prp.2018.12.006
31. Wu H, Liu B, Chen Z, Li G, Zhang Z. MSC-Induced Lncrna HCP5 Drove Fatty Acid Oxidation Through Mir-3619-5p/AMPK/PGC1α/CEBPB Axis to Promote Stemness and Chemo-Resistance of Gastric Cancer. Cell Death Dis (2020) 11(4):233. doi: 10.1038/s41419-020-2426-z
32. Grace M, McLeland C, Gagliardi S, Smith J, Jackson W, Blakely W. Development and Assessment of a Quantitative Reverse transcription-PCR Assay for Simultaneous Measurement of Four Amplicons. Clin Chem (2003) 49: (9):1467–75. doi: 10.1373/49.9.1467
33. Beermann J, Piccoli M, Viereck J, Thum T. Non-Coding RNAs in Development and Disease: Background, Mechanisms, and Therapeutic Approaches. Physiol Rev (2016) 96(4):1297–325. doi: 10.1152/physrev.00041.2015
34. Schmitz S, Grote P, Herrmann B. Mechanisms of Long Noncoding RNA Function in Development and Disease. Cell Mol Life Sci CMLS (2016) 73(13):2491–509. doi: 10.1007/s00018-016-2174-5
35. Zhu Y, Hu H, Yuan Z, Zhang Q, Xiong H, Hu Z, et al. Lncrna NEAT1 Remodels Chromatin to Promote the 5-Fu Resistance by Maintaining Colorectal Cancer Stemness. Cell Death Dis (2020) 11(11):962. doi: 10.1038/s41419-020-03164-8
36. Zhang C, Wang L, Jin C, Zhou J, Peng C, Wang Y, et al. Long Non-Coding RNA lnc-LALC Facilitates Colorectal Cancer Liver Metastasis Via Epigenetically Silencing LZTS1. Cell Death Dis (2021) 12(2):224. doi: 10.1038/s41419-021-03461-w
37. Zhang Y, Zhang X, Cai B, Li Y, Jiang Y, Fu X, et al. The Long Noncoding RNA lncCIRBIL Disrupts the Nuclear Translocation of Bclaf1 Alleviating Cardiac Ischemia-Reperfusion Injury. Nat Commun (2021) 12(1):522. doi: 10.1038/s41467-020-20844-3
38. Chen Z, Song X, Li Q, Xie L, Guo T, Su T, et al. Tmprss2androgen Receptor-Activated Enhancers Simultaneously Regulate Oncogene and lncRNA in Prostate Cancer. Cells (2019) 8(8):864. doi: 10.3390/cells8080864
39. Wang M, Mao C, Ouyang L, Liu Y, Lai W, Liu N, et al. Long Noncoding RNA LINC00336 Inhibits Ferroptosis in Lung Cancer by Functioning as a Competing Endogenous RNA. Cell Death Differ (2019) 26(11):2329–43. doi: 10.1038/s41418-019-0304-y
40. Jiang L, Wang R, Fang L, Ge X, Chen L, Zhou M, et al. HCP5 is a SMAD3-Responsive Long Non-Coding RNA That Promotes Lung Adenocarcinoma Metastasis Via miR-203/SNAI Axis. Theranostics (2019) 9(9):2460–74. doi: 10.7150/thno.31097
41. Benes V, Castoldi M. Expression Profiling of microRNA Using Real-Time Quantitative PCR, How to Use It and What Is Available. Methods (San Diego Calif) (2010) 50(4):244–9. doi: 10.1016/j.ymeth.2010.01.026
42. Qin Y, Hou Y, Liu S, Zhu P, Wan X, Zhao M, et al. A Novel Long Non-Coding RNA Lnc030 Maintains Breast Cancer Stem Cell Stemness by Stabilizing Sqle mRNA and Increasing Cholesterol Synthesis. Advanced Sci (Weinheim Baden Wurttemberg Germany) (2021) 8(2):2002232. doi: 10.1002/advs.202002232
43. Sikora M, Marycz K, Smieszek A. Small and Long Non-Coding RNAs as Functional Regulators of Bone Homeostasis, Acting Alone or Cooperatively. Mol Ther Nucleic Acids (2020) 21:792–803. doi: 10.1016/j.omtn.2020.07.017
44. Anfossi S, Babayan A, Pantel K, Calin G. Clinical Utility of Circulating non-Coding RNAs - An Update. Nat Rev Clin Oncol (2018) 15(9):541–63. doi: 10.1038/s41571-018-0035-x
45. Pritchard C, Kroh E, Wood B, Arroyo J, Dougherty K, Miyaji M, et al. Blood Cell Origin of Circulating microRNAs: A Cautionary Note for Cancer Biomarker Studies. Cancer Prev Res (Philadelphia Pa) (2012) 5(3):492–7. doi: 10.1158/1940-6207.Capr-11-0370
46. Bornschein J, Selgrad M, Wex T, Kuester D, Malfertheiner P. Serological Assessment of Gastric Mucosal Atrophy in Gastric Cancer. BMC Gastroenterol (2012) 12:10. doi: 10.1186/1471-230x-12-10
Keywords: gastric cancer, GC, HCP5, biomarker, diagnosis, gastritis
Citation: Qin S, Yang L, Kong S, Xu Y, Liang B and Ju S (2021) LncRNA HCP5 : A Potential Biomarker for Diagnosing Gastric Cancer. Front. Oncol. 11:684531. doi: 10.3389/fonc.2021.684531
Received: 23 March 2021; Accepted: 24 May 2021;
Published: 18 June 2021.
Edited by:
Jaw-Yuan Wang, Kaohsiung Medical University Hospital, TaiwanReviewed by:
Matthew Brendan O’Rourke, University of Technology Sydney, AustraliaPeng-Chan Lin, National Cheng Kung University, Taiwan
De-Chuan Chan, Tri-Service General Hospital, Taiwan
Copyright © 2021 Qin, Yang, Kong, Xu, Liang and Ju. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Bo Liang, 281692191@qq.com; Shaoqing Ju, jsq814@hotmail.com
†These authors have contributed equally to this work