Abstract

Background. Expression of long noncoding RNA (lncRNA) ELF3 antisense RNA 1 (ELF3-AS1) is observed in some cancers, while its role in hepatocellular carcinoma (HCC) is unclear. The study aimed to investigate the relationship between ELF3-AS1 and HCC based on database, bioinformatics, and statistical analysis. Methods. In this study, Kruskal–Wallis test, Wilcoxon sign-rank test, logistic regression, Kaplan–Meier method, Cox regression analysis, gene set enrichment analysis (GSEA), and immunoinfiltration analysis were used to assess the relationship between ELF3-AS1 expression and clinical characteristics of HCC patients, the relationship between ELF3-AS1 expression and prognosis of HCC patients, and the possible functions of ELF3-AS1 in HCC. Results. High expression of ELF3-AS1 in patients with HCC was related to T stage ( < 0.001), gender ( = 0.006), residual tumor ( = 0.008), histologic grade ( < 0.001), adjacent hepatic tissue inflammation ( = 0.011), AFP ( < 0.001), and vascular invasion ( = 0.028). High ELF3-AS1 expression was associated with poor overall survival (OS) ( = 0.001) and DSS ( = 0.047). ELF3-AS1 expression ( = 0.011) was independently correlated with OS in HCC patients. In the high ELF3-AS1 expression group, GPCR-radioligand binding, M phase, Class A/1 (rhodopsin-like receptors), cell cycle checkpoints, translation, mitotic metaphase and anaphase, signaling by robo receptors, keratinization, and rRNA processing were differentially enriched by GESA. ELF3-AS1 expression was associated with immune infiltrating cells. Conclusions. ELF3-AS1 expression was associated with poor prognosis in HCC. ELF3-AS1 expression was significantly associated with immune infiltration. ELF3-AS1 is a promising biomarker that can be used for the diagnosis and prognosis of HCC.

1. Introduction

Hepatocellular carcinoma (HCC) is one of the leading causes of cancer deaths in the world [1]. In the last few years, more than 700,000 people have died each year from HCC, and this number is increasing every year [2]. HCC is a complex neoplastic disease that is highly aggressive, with a pathology that includes changes in tumor cell behavior and the vascular system [3]. Nowadays, surgery is only effective for patients with early HCC [4]. Most patients are diagnosed at an advanced stage, and radiotherapy and chemotherapy are not effective [4]. The high recurrence and metastasis rates of HCC reduce patient survival [5]. Despite some advances in surgical treatment and neoadjuvant therapy, the prognosis of HCC is poor [6]. The lack of useful markers makes it difficult for clinicians to predict the clinical outcome of patients with HCC. Therefore, new markers need to be explored to identify HCC patients with poor prognosis.

The length of long noncoding RNAs (lncRNAs) is over 200 bp [7]. lncRNAs play a key role in the formation of cancer [8]. Growing evidence for the involvement of aberrant lncRNAs in the development and progression of HCC [4, 6, 9, 10]. Therefore, screening for lncRNAs that are clinically relevant to HCC is important for the early diagnosis and effective treatment of HCC.

The lncRNA ELF3 antisense RNA 1 (ELF3-AS1) is the antisense transcript for ELF3. ELF3-AS1 was used as an oncogene in bladder cancer (BLCA) [11]. ELF3-AS1 is upregulated in oral squamous cell carcinoma (OSCC) [12]. ELF3-AS1 is upregulated in osteosarcoma (OS) [13]. ELF3-AS1 is upregulated in nonsmall cell lung cancer (NSCLC) [14]. ELF3-AS1 is a prognostic marker and potential therapeutic target for gliomas [15]. However, the correlation between ELF3-AS1 and HCC has not been studied.

In this study, we compared the differences in ELF3-AS1 expression between HCC tumor tissues and normal samples, explored the correlation between ELF3-AS1 expression and clinical characteristics of HCC, and assessed the prognostic value of ELF3-AS1 in HCC. For ELF3-AS1 high and low expression groups, genomic enrichment analysis (GSEA) was performed to explore the possible functions of ELF3-AS1. The possible function of ELF3-AS1 in HCC was explored by immune infiltration analysis. The study may provide new directions for the development of diagnostic and therapeutic strategies for HCC.

2. Materials and Methods

2.1. Differential Expression of ELF3-AS1
2.1.1. Baseline Information

Software used was R (version 3.6.3) (statistical analysis and visualization). R package type was basic R package. The molecule is ELF3-AS1 (ENSG00000234678). The subgroup is median. Data are RNAseq data and clinical data in level 3 HTSeq-FPKM format from the TCGA (https://portal.gdc.cancer.gov) HCC (hepatocellular carcinoma) project. Data filtering is through retaining data with clinical information. Data conversion included RNAseq data in FPKM (Fragments Per Kilobase per Million) format converted to TPM (transcripts per million reads) format and grouped according to molecular expression.

2.1.2. Unpaired Samples

Software used was R (version 3.6.3). R package type was mainly ggplot2. Molecules used were ELF3-AS1. Data are RNAseq data in level 3 HTSeq-FPKM format from the TCGA HCC project. Data filtering is none. Data transformation included FPKM format RNAseq data converted to TPM format and log2 transformed for sample-to-sample expression comparisons.

2.1.3. Paired Samples

Software used was R (version 3.6.3). R package was mainly ggplot2. Molecules included ELF3-AS1. Data are RNAseq data in level 3 HTSeq-FPKM format from the TCGA HCC project. Data filtering is through retaining paired samples. Data transformation included RNAseq data in FPKM format converted to TPM format and log2 transformed for sample-to-sample expression comparisons.

2.1.4. ROC Analysis

Software used was R (version 3.6.3). R packages used were pROC package (for analysis) || ggplot2 package. Molecule is ELF3-AS1. Clinical variables are tumor vs. normal. Data are RNAseq data and clinical data in level 3 HTSeq-FPKM format from the TCGA HCC project. Data filtering is through retaining data with clinical information. Data transformation included RNAseq data in FPKM format converted to TPM format and log2 transformed for analysis. ROC results were interpreted as false positive rate (FPR) is the horizontal coordinate and true positive rate (TPR) is the vertical coordinate.

2.2. The Relationship between ELF3-AS1 and Clinical Characteristics
2.2.1. Correlation of Gene Expression with Clinical Characteristics

Software used was R (version 3.6.3). R package mainly included ggplot2. Molecule is ELF3-AS1. Clinical variables: T stage, gender, residual tumor, histological grade, adjacent hepatic tissue inflammation, AFP, and vascular invasion. Data are RNAseq data and clinical data in level 3 HTSeq-FPKM format from the TCGA HCC project. Data filtering is through retaining data with clinical information. Data transformation included RNAseq data in FPKM format converted to TPM format and log2 transformed for analysis.

2.2.2. Logistics Analysis

Software used was R (version 3.6.3). R package type was mainly the basic package. Statistical method is the dichotomous logistic model. Dependent variable is ELF3-AS1. Type of independent variable is low high dichotomous. Data are RNAseq data in level 3 HTSeq-FPKM format from the TCGA HCC project. Data filtering is through retaining data with clinical information. Data conversion included RNAseq data in FPKM format converted to TPM format and grouped according to molecular expression.

2.3. The Relationship between ELF3-AS1 and Prognosis
2.3.1. Kaplan–Meier Method

Software used was R (version 3.6.3). R packages are survminer package and survival package. Molecule is ELF3-AS1. Subgroups are 0–50 vs. 50–100. Prognosis type is overall survival and disease specific survival. Data are RNAseq data and clinical data in level 3 HTSeq-FPKM format from the TCGA HCC project. Data filtering is through retaining data with clinical information. Data conversion included RNAseq data in FPKM format converted to TPM format and analyzed by grouping them according to molecular expression. Additional data are prognostic data from an article [16].

2.3.2. COX Regression

Software used was R (version 3.6.3). R package is survivor package. Statistical method is COX regression module. Prognosis type is overall survival. Included variables are clinical characteristics and ELF3-AS1. Data are RNAseq data and clinical data in level 3 HTSeq-FPKM format from the TCGA HCC project. Data filtering is through retaining data with clinical information. Data conversion included RNAseq data in FPKM format converted to TPM format and analyzed by grouping them according to molecular expression. Additional data are prognostic data from an article [16].

2.4. Enrichment of ELF3-AS1-Related Pathways
2.4.1. Single Gene Differential Analysis

Software used was R (version 3.6.3). R package is deseq2 [17]. Target molecule is ELF3-AS1. Low expression group is 0–50%. High expression group is 50–100%. Data are RNAseq data in level 3 HTSeq-Counts format from the TCGA HCC project. Data filtering is filtering of paracancer samples.

2.4.2. GSEA Analysis

Software used was R (version 3.6.3). R package are ggplot2 package and clusterProfiler package [18] (for GSEA analysis). Method is gene set enrichment analysis (GSEA) [19]. Species is Homo sapiens. Reference gene collection is c2.cp.v7.2.symbols.gmt (curated). Gene set database is MSigDB Collections (database hyperlink) (with detailed descriptions of individual gene sets).

2.5. The Relationship between ELF3-AS1 Expression and Immune Infiltrating Cells

Software used was R (version 3.6.3). R package is GSVA package [20]. Immunocell algorithm is ssGSEA (algorithm built into the GSVA package). Molecule is ELF3-AS1. Immune cells are 24 immune cells. Data are RNAseq data and clinical data in level 3 HTSeq-FPKM format from the TCGA HCC project. Data filtering is removal of paracancerous tissue. Expression profile data conversion included FPKM format RNAseq data converted to TPM format and log2 transformed for analysis. Other data included markers for 24 immune cells obtained from an article [21].

3. Results

3.1. Clinical Characteristics

As shown in Table 1, the characteristics of HCC in TCGA were collected. The age included 177 patients (≤60, 47.5%) and 196 patients (>60, 52.5%). The median age is 61 years, with a range of 52 to 69 years. The race included 185 white patients, 160 Asian patients, and 17 Black or African American patients. The pathologic stage included 173 stage I (49.4%), 87 stage II (24.9%), 85 stage III (24.3%), and 5 stage IV (1.4%). The tumor status included 202 tumor-free (56.9%) and 153 with tumor (43.1%). The gender included 121 females (32.4%) and 253 males (67.6%). The weight included 184 (≤70) (53.2%) and 162 (>70) (46.8%). The height included 201 (<170) (58.9%) and 140 (≥170) (41.1%). The BMI included 177 (≤25) (52.5%) and 160 (>25) (47.5%). The residual tumor included 327 R0 (94.8%), 17 R1 (4.9%), and 1 R2 (0.3%). The histologic grade included 55 G1 (14.9%), 178 G2 (48.2%), 124 G3 (33.6%), and 12 G4 (3.3%). The adjacent hepatic tissue inflammation included 118 none (49.8%), 101 mild (42.6%), and 18 severe (7.6%). The albumin included 69 (<3.5) (23%) and 231 (≥3.5) (77%). The AFP included 215 (≤400) (76.8%) and 65 (>400) (23.2%). The prothrombin time included 208 (≤4) (70%) and 89 (>4) (30%). The Child-Pugh grade included 219 A (90.9%), 21 B (8.7%), and 1 C (0.4%). The fibrosis Ishak score included 75 0 (34.9%), 31 1/2 (14.4%), 28 3/4 (13%), and 81 5/6 (37.7%). The vascular invasion included 208 no (65.4%) and 110 yes (34.6%).

3.2. ELF3-AS1 Expression Correlated with Poor Clinical Characteristics of HCC

As shown in Figure 1(a), ELF3-AS1 was highly expressed in HCC tissues (0.978 ± 0.061 vs. 2.381 ± 0.052,  < 0.001). As shown in Figure 1(b), LF3-AS1 was highly expressed in HCC tissues ( < 0.001), based on 50 HCC tissues and their matched normal liver tissues. As shown in Figure 1(c), the AUC of ELF3-AS1 was 0.904. As shown in Table 2, clinical and gene expression data were collected. As shown in Table 2, the value of T stage is 0.002, the value of pathologic stage is 0.003, the value of gender is 0.008, the value of race is 0.037, the value of residual tumor is  = 0.004, the value of histologic grade is <0.001, the value of adjacent hepatic tissue inflammation is 0.025, the value of AFP is <0.001, and the value of vascular invasion is 0.037. As shown in Figure 2 and Table 3, the value of T stage is <0.001, the value of gender is 0.006, the value of residual tumor is 0.008, the value of histologic grade is <0.001, the value of adjacent hepatic tissue inflammation is 0.011, the value of AFP is <0.001, and the value of vascular invasion is 0.028.

3.3. Relationship between ELF3-AS1 and Survival of HCC Patients

As shown in Figure 3, the expression of ELF3-AS1 was positively correlated with poor OS ( = 0.001), disease-specific survival (DSS) ( = 0.047) of HCC patients. As shown in Table 4, the results of univariate analysis showed that high ELF3-AS1 expression levels ( < 0.001) were associated with T stage ( < 0.001), pathologic stage ( < 0.001), and tumor status ( < 0.001). The results of multivariate analysis showed that ELF3-AS1 expression ( = 0.011), pathologic stage ( = 0.013), age ( = 0.013), and tumor status ( = 0.005) were independently correlated with OS in multivariate analysis. The results suggested that increased expression of ELF3-AS1 was associated with poor OS.

3.4. ELF3-AS1-Related Pathways

A dataset of 191 significant differences was enriched in ELF3-AS1 low expression phenotypes. As shown in Table 5 and Figure 4, the top 9 low -value datasets include GPCR-radioligand binding, M phase, Class A/1 (rhodopsin-like receptors), cell cycle checkpoints, translation, mitotic metaphase and anaphase, signaling by robo receptors, keratinization, and rRNA processing.

3.5. Correlation of ELF3-AS1 Expression with Immune Infiltration

As shown in Figure 5 and Table 6, ELF3-AS1 expression is negatively and significantly correlated with infiltration levels of CD8 T cells ( = 0.001), cytotoxic cells ( = 0.045), eosinophils ( = 0.002), neutrophils ( < 0.001), NK cells ( = 0.001), Tcm ( < 0.001), Th17 cells ( = 0.003), and TReg ( = 0.017) and positively correlated with that of aDC ( = 0.002), B cells ( = 0.004), macrophages ( < 0.001), NK CD56bright cells ( < 0.001), T cells ( = 0.042), TFH ( < 0.001), and Th2 cells ( < 0.001).

4. Discussion

LncRNAs play a key role in tumorigenesis and progression [22]. ZNF385D-AS2 may be a useful biomarker for prognosis in patients with HCC [9]. High LUCAT1 expression is an independent prognostic factor for HCC [6]. F11-AS1 may serve as a therapeutic target for HCC [4]. GIHCG is a biomarker that can be used to predict the prognosis of patients with HCC [10]. Therefore, the study of lncRNAs as new HCC biomarkers and therapeutic targets is very important.

ELF3-AS1 is upregulated in OSCC, OS, and NSCLC [1214]. In this study, high ELF3-AS1 expression in HCC was associated with the characteristics including T stage, gender, residual tumor, histologic grade, adjacent hepatic tissue inflammation, AFP, and vascular invasion. High ELF3-AS1 expression was associated with poor OS and DSS. ELF3-AS1 expression was independently correlated with OS in HCC patients.

ELF3-AS1 mediated BLCA tumorigenesis through enhanced ELF3-AS1/KLF8 signaling [11]. ELF3-AS1 may be regulating GLUT1 to promote the proliferation of OSCC cells [12]. ELF3-AS1 may promote OS cell proliferation by upregulating KLF12 through methylation of the miR-205 gene [13]. ELF3-AS1 promotes cancer cell invasion and migration through downregulation of miR-212 by methylation in NSCLC [14]. ELF3-AS1 is not only an important prognostic marker, but also a potential therapeutic target for glioma [15]. In this study, based on GESA, ELF3-AS1 was related to pathways including GPCR-radioligand binding, M phase, Class A/1 (rhodopsin-like receptors), cell cycle checkpoints, translation, mitotic metaphase and anaphase, signaling by robo receptors, keratinization, and rRNA processing.

Immune infiltrating cells in HCC are currently a hot topic, and a major understanding of immune infiltrating cells can support the development of new and emerging immunotherapies. The development and progression of HCC are associated with a unique immune response profile of the liver microenvironment, in which CD4+CD25+Foxp3 regulatory T cells (Tregs) play a key role through their immunosuppressive effects [23]. Immunotherapy enhances the immune response against hepatocellular carcinoma by blocking Treg activity [23]. In this study, we explored the relationship between ELF3-AS1 expression in HCC and various immune infiltrations. ELF3-AS1 expression was associated with infiltration of CD8 T cells, cytotoxic cells, eosinophils, neutrophils, NK cells, Tcm, Th17 cells, Treg, aDC, B cells, macrophages, NK CD56bright cells, T cells, TFH, and Th2 cells in HCC. These correlations may suggest potential mechanisms by which ST-AS1 inhibits the function of CD8 T cells, cytotoxic cells, eosinophils, neutrophils, NK cells, Tcm, Th17 cells, and Treg and subsequently promotes the function of aDC, B cells, macrophages, NK CD56bright cells, T cells, TFH, and Th2 cells.

Despite some limitations, this is an exploratory study to investigate the relationship between ELF3-AS1 and HCC. This study is based on database and bioinformatics analysis. We have not yet investigated in detail the direct mechanism of ELF3-AS1 involvement in HCC. Further studies are needed regarding the direct mechanism of ELF3-AS1-mediated HCC.

5. Conclusion

ELF3-AS1 was highly expressed in HCC relative to normal tissue and related to poor OS and DSS. ELF3-AS1 might participate in the development of HCC by pathways including GPCR-radioligand binding, M phase, Class A/1 (rhodopsin-like receptors), cell cycle checkpoints, translation, mitotic metaphase and anaphase, signaling by robo receptors, keratinization, and rRNA processing. ELF3-AS1 expression was associated with immune infiltrating cells. This study investigated the role of ELF3-AS1 in HCC and provided a promising biomarker for the diagnosis and prognosis of HCC.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest relevant to this study.

Authors’ Contributions

TC and YP conceived and designed the study. TC, CZ, XW, and YP analyzed the data. TC and YP wrote and revised the paper. All authors read and approved the final manuscript.

Acknowledgments

This study was supported by grants from the National Natural Science Foundation of China (grant numbers 81960431 and 81960535).