Integrated analysis of high-throughput sequencing-based lncRNA-mediated ceRNA network in Hepatic Alveolar Echinococcosis
========================================================================================================================

* Zhen Liu
* Chang-zhen Shang
* Jin-peng Wang
* Zhi-gang Gai
* Fu-cai Ma
* Pan Xia
* Yan Wang
* Xiao Yang
* Hai-hong Zhu

## Abstract

**Background** Numerous studies have indicated that long non-coding RNAs (lncRNAs) can modulate the expression of target gene mRNAs by adsorbing microRNAs (miRNAs). The lncRNA–miRNA–mRNA ceRNA network has been theorized to play an indispensable role in many types of tumors, and has been garnering increasing attention. However, the role of the lncRNA-associated ceRNA regulatory network in Hepatic Alveolar Echinococcosis (HAE) remains unclear and requires further exploration.

**Methods** In this study, high-throughput sequencing was performed on lesion tissues and adjacent tissues from three patients with Hepatic Alveolar Echinococcosis (HAE) to identify differentially expressed RNAs. We utilized Cytoscape (version 3.10.1) to construct the lncRNA-miRNA-mRNA ceRNA network based on the interactions from the miRcode, miRTarBase, miRDB, and TargetScan databases, and identified hub lncRNAs from within the ceRNA network. Through the use of the “clusterProfiler” package in R, we performed Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway annotations for the DEGs (Differentially Expressed Genes) within the ceRNA network. Concurrently, we utilized these DEGs to construct a protein-protein interaction network (PPI). Finally, an analysis was conducted on the PCBP1-AS1-miR-20b-5p/CAPRIN2 axis within the ceRNA network.

**Results** In HAE, a total of 979 differentially expressed lncRNAs (DELncRNAs) and 870 differentially expressed mRNAs (DEmRNAs) were identified. An HAE-specific ceRNA network comprising 11 lncRNAs, 21 miRNAs, and 56 mRNAs was established, and analysis of this network led to the construction of a sub-network associated with hub lncRNAs. GO and KEGG pathway analyses indicated that the HAE-specific ceRNA network is related to molecular functions and pathways associated with cancer. Subsequent experiments using qPCR and dual-luciferase assays validated the interactions between PCBP1-AS1 and miR-20b-5p, as well as between miR-20b-5p and CAPRIN2. Analysis of the target gene in relation to clinical characteristics of HAE patients suggested that the PCBP1-AS1-miR-20b-5p/CAPRIN2 axis may influence the development of HAE.

**Conclusion** In this study, we described the gene regulation within the lncRNA-miRNA-mRNA ceRNA network during the development of Hepatic Alveolar Echinococcosis (HAE), which contributes to a deeper exploration of the molecular mechanisms underlying HAE. Additionally, we discovered that PCBP1-AS1 may regulate the expression of CAPRIN2 by adsorbing miR-20b-5p, affecting the onset and progression of HAE. PCBP1-AS1 could potentially serve as a useful target for the diagnosis and treatment of HAE.

Keywords
*   High-throughput sequencing
*   Hepatic Alveolar Echinococcosis
*   long non-coding RNA
*   competitive endogenous RNA network
*   lnc-PCBP1-AS1-miR20b-5p/CAPRIN2 axis

## Introduction

Hepatic alveolar echinococcosis (HAE) is a zoonotic disease caused by the larval of the multilocular echinococcus, accounting for 5-10% of the total incidence of echinococcosis[1]. The growth of HAE is slow and the disease often remains asymptomatic for an extended period, approximately 5 to 15 years. Its growth pattern resembles that of a malignant condition, capable of infiltrating or invading surrounding tissues and organs, and it may even metastasize to distant organs[2]. Despite the advancements made in the diagnosis and treatment of HAE, the therapeutic outcomes and prognosis are still less than satisfactory[3]. Research conducted through BLAST sequence analysis of the Echinococcus genome has indicated that approximately one-third of the genes may be unique characteristics of the species and form the basis of its biological traits, which are specific to Echinococcosis. These genes could potentially aid in refining diagnostic methods and identifying novel therapeutic targets[4]. Therefore, this study preliminarily explores the biological mechanisms and pathogenesis of HAE at the genetic level, in an effort to identify new therapeutic targets to improve treatment outcomes and prognosis.

In recent years, advances in sequencing technologies have enabled more in-depth genomic research. It is now known that less than 2% of the human genome encodes proteins. Among the various types of non-coding RNAs (ncRNAs), long non-coding RNAs (lncRNAs) have attracted significant attention. These are defined as long chains of non-coding RNA that exceed 200 base pairs in length and do not have the function of encoding proteins[5][6]. However, through their diverse mechanisms, lncRNAs play functional roles in various cell types, and many lncRNAs possess both tumor-suppressive and oncogenic capabilities[7][8]. Additionally, microRNAs (miRNAs) are endogenous RNA molecules approximately 22 base pairs in length that can bind to target RNAs, promoting the degradation of mRNA or inhibiting protein translation[9][10], they play a crucial role in regulating gene expression and are involved in numerous biological processes through their network of regulatory interactions[11][12][13].

In 2011, Salmena and colleagues first proposed the competitive endogenous RNA (ceRNA) hypothesis. This hypothesis suggests that ceRNAs interact with target miRNAs through miRNA response elements (MREs), thereby regulating the transcriptome on a large scale, Simultaneously, this concept of ceRNAs and their networks offers a novel perspective on the intricate interactions of RNA molecules within the cell, enhancing our comprehension of gene regulation and the potential for targeted therapies in the treatment of diseases[14]. In recent years, a growing body of research has indicated that lncRNAs can act as ceRNAs, regulating the expression of target mRNAs by competitively sharing miRNAs. They play a significant role in various types of tumors, Such as in hepatocellular carcinoma[15], pancreatic cancer[16], gastric cancer[17], and colorectal cancer[18].However, to date, there have been no reports on the study of lncRNAs and ceRNA networks in HAE, and the expression profile of lncRNAs and the associated ceRNA regulatory networks in HAE require further clarification. Therefore, it is particularly important to further understand the role of lncRNAs in the development of HAE and to investigate the pathogenic mechanisms involving ceRNA networks in HAE.

In this study, we performed high-throughput sequencing on lesion tissues and adjacent tissues from three HAE patients to identify differentially expressed lncRNAs and mRNAs. Through comprehensive analysis, we constructed a ceRNA network based on the interactions between lncRNA-miRNA and miRNA-mRNA, which will help enhance our understanding of the pathogenesis of HAE, identify new therapeutic targets, improve treatment outcomes, and refine prognosis. Furthermore, our analysis of the Lnc-PCBP1-AS1-miR-20b-5p/CAPRIN2 axis revealed that overexpressed Lnc-PCBP1-AS1 may play a key role in the development and progression of HAE by adsorbing miR-20b-5p and thereby targeting the expression of CAPRIN2.

## Methods

### Patients and Samples

This study collected 13 cases of HAE patients who met the surgical resection criteria at Qinghai Provincial People’s Hospital from September 2021 to October 2023. The collected tissues included lesion center tissue, tissue <1cm around the lesion, and tissue >1cm around the lesion. The enrolled patients had no previous history of liver surgery and were not complicated with hepatic cystic echinococcosis, other parasitic diseases, liver benign or malignant tumors, liver abscess, liver metastatic cancer, or other space-occupying liver diseases. Both preoperative imaging and postoperative pathology diagnosed hepatic alveolar echinococcosis. After the specimens were removed, the required tissue samples were quickly cut and placed in RNA-free freezing tubes, then promptly stored in liquid nitrogen at −80°C for subsequent RNA extraction. This study was approved by the Ethics Committee of Qinghai Provincial People’s Hospital (approval number 2021-32). All participants voluntarily joined the project and signed an informed consent form. Table 1 presents the clinical data of the 13 HAE patients.

View this table:
[Table 1](http://medrxiv.org/content/early/2024/07/03/2024.07.02.24309861/T1)

Table 1 Clinicopathological characteristics of 13 HAE patients.

### RNA Extraction and Sequencing

Total RNA extraction was performed using TRIzol reagent (Thermo Scientific, USA). The RNA samples were assessed using 1% agarose gel electrophoresis, NanoPhotometer® spectrophotometer (Implen, CA, USA), Qubit® 3.0 Fluorometer (Implen, CA, USA), and Agilent 2100 RNA Nano 6000 Assay Kit (Implen, CA, USA). Three micrograms of total RNA from both the central lesion tissue and surrounding tissue of HAE patients were used as starting material for library construction. After library construction, initial quantification was done using Qubit 3.0, and the quality of the libraries was evaluated using the Agilent 2100, Bio-RAD CFX 96 real-time PCR instrument, and Bio-RAD KIT iQ SYBR GRN to ensure library quality. Clustering was followed by paired-end sequencing using the Illumina 4000 high-throughput sequencing platform (HiSeq/MiSeq). Uniquely mapped reads were used for gene read count and FPKM (Fragments Per Kilobase of transcript per Million mapped reads) calculation.

### Differentially Expressed Analysis

The R software package DESeq2 was utilized to calculate the differential expression of transcripts using a negative binomial distribution test. The cutoff criteria for identifying differentially expressed genes (DEGs) were set at a P-value <0.05 and a Fold Change(FC)>2/<0.5.

### Construction of a ceRNA Regulatory Network

The construction of the ceRNA network involved using the miRcode database ([http://www.mircode.org/](http://www.mircode.org/)) to match differentially expressed lncRNAs (FC >2 or <0.5, *P* <0.05) with miRNAs. Based on the miRWalk database ([http://mirwalk.umm.uni-heidelberg.de/](http://mirwalk.umm.uni-heidelberg.de/)), three bioinformatics algorithms—miRTarBase, miRDB, and TargetScan—were employed to predict the target mRNAs of miRNAs (FC >2 or <0.5, *P* <0.05). The predicted genes were then intersected with the genes obtained from sequencing. Following data integration, we constructed a comprehensive lncRNA-miRNA-mRNA regulatory network, which was visualized using Cytoscape (version v3.10.1).

Additionally, sub-networks were reconstructed using the BinGO plugin in Cytoscape[19].

### Functional Enrichment and Protein–Protein Interaction Analysis

To gain a deeper understanding of the mechanisms by which differentially expressed mRNAs (DEmRNAs) contribute to the development of HAE, we utilized the “clusterProfiler” package in R software for functional annotation analysis. We also employed the ggplot2 package in R to create a bubble chart. *P* <0.05 was considered a threshold for significance. The protein-protein interaction (PPI) network was established using the STRING online search tool with a combined score of ≥0.4. The PPI network was visualized using Cytoscape (version v3.10.1).

### Real-Time Fluorescent Quantitative PCR (qRT-PCR)

We determined the relative expression levels of selected lncRNAs, miRNAs, and mRNAs using qRT-PCR. This study used GAPDH and U6 as internal references for qRT-PCR, and the 2^-ΔΔCt method was employed to calculate the relative expression levels of the selected genes. The primer sequences for lncRNAs, mRNAs, and miRNAs, as well as the internal references, are listed in Table 2 and Table 3.

View this table:
[Table 2](http://medrxiv.org/content/early/2024/07/03/2024.07.02.24309861/T2)

Table 2 LncRNA and mRNA primer sequences.

View this table:
[Table 3](http://medrxiv.org/content/early/2024/07/03/2024.07.02.24309861/T3)

Table 3 miRNA primer sequences.

### Cell culture and transfection

The 293T human embryonic kidney cell line was sourced from Promega, and the cells were cultured under conditions of DMEM medium supplemented with 10% FBS and 1% PS, at 37°C in a humidified atmosphere with 5% CO2. The overexpression plasmids for hsa-miR-20b-5p and the control plasmids were purchased from GeneCopoeia, while the wild-type and mutant plasmids for PCBP1-AS1 and CAPRIN2 were obtained from Universal Bio. Transfection of 293T human embryonic kidney cells was carried out using Lipofectamine 3000 reagent (Invitrogen, Shanghai, China), with 0.5μg of plasmid, 0.05μM of miRNA, and 1μL of P3000™ reagent. Subsequently, the Dual-Glo® Luciferase assay system was used to measure the luminescence of both firefly and Renilla luciferases. Luciferase activity was represented by the ratio of firefly to Renilla luciferase luminescence intensity. The experiment was conducted independently in triplicate.

### Statistical analysis

Data were statistically analyzed using SPSS 28.0 software (IBM, Armonk, NY, USA) or GraphPad Prism 9.0 (GraphPad Software, San Diego, CA, USA). Non-parametric data by group were analyzed using the Friedman rank-sum test, and experimental data are presented in the form of *M* (*P*25, *P*75). Student’s t-test or one-way analysis of variance (ANOVA) was used to analyze differences between two experimental groups or among more than two groups, respectively. Data were expressed as the mean ± standard deviation (SD). *P* <0.05 was considered statistically significant.

## Results

### Screening for the Coding Potential of LncRNAs

We sequenced lncRNAs from the central lesion tissue and surrounding tissue of three patients with human HAE. Since lncRNAs are RNAs that do not encode proteins, we predicted their coding potential to differentiate RNAs with encoding functions from those without encoding potential in the sequencing data. To distinguish lncRNAs, we utilized a comprehensive set of the most widely applied coding potential analysis methods, including: CPC analysis[20], CNCI analysis [21], pfam protein domain analysis [22], and PLEK analysis [23]. Specifically, CPC analysis identified 9,046 lncRNAs, CNCI analysis identified 35,119 lncRNAs, pfam protein domain analysis identified 37,943 lncRNAs, and PLEK analysis identified 2,602 lncRNAs. By integrating the predictions from these four methods, we obtained a total of 26,365 lncRNAs. (Figure 1)

![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/07/03/2024.07.02.24309861/F1.medium.gif)

[Figure 1:](http://medrxiv.org/content/early/2024/07/03/2024.07.02.24309861/F1)

Figure 1: 
Venn diagram showing the results of four coding potential prediction methods.

### LncRNA expression profiles in HAE

To date, tens of thousands of lncRNA genes have been identified in humans and single-celled eukaryotic organisms. Based on the relative positional relationship between lncRNAs and protein-coding genes, lncRNAs are classified into intergenic lncRNA (lincRNA), intronic lncRNA, antisense lncRNA, sense lncRNA, and bidirectional lncRNA[24]. According to the lncRNA sequencing results, there are 9,022 intronic lncRNAs (denoted as ‘i’), 9,057 sense overlapping lncRNAs (denoted as ‘o’), 5,733 intergenic lncRNAs (denoted as ‘u’), and 2,553 antisense lncRNAs (denoted as ‘x’), with the majority containing two or more exons (Figure 2).

![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/07/03/2024.07.02.24309861/F2.medium.gif)

[Figure 2:](http://medrxiv.org/content/early/2024/07/03/2024.07.02.24309861/F2)

Figure 2: 
Characteristics of LncRNAs in HAE sequencing. (A)Distribution bar chart of lncRNA types. (B) Bar chart of the number of exons contained in lncRNA.

### The selection of DELncRNAs in HAE

To investigate the expression of LncRNAs in patients with HAE, high-throughput sequencing technology was employed, with a P-value <0.05 and a fold change (FC) >2 or FC <0.5 as the criteria for screening differentially expressed LncRNAs (DElncRNAs). A total of 979 DElncRNAs were identified (591 upregulated and 388 downregulated). These differentially expressed LncRNAs are displayed as a heatmap (Figure 3A) and a volcano plot (Figure 3B), they exhibit good discrimination and clustering between the central lesion tissue group and the surrounding lesion tissue group. The top 10 up-regulated and down-regulated DElncRNAs with their names, FoldChange, Log2FoldChange and P-Value are listed in Table 4. All DElncRNAs with their names, Fold Change, Log2FoldChange and P-Value are listed in Supplementary Table 1.

![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/07/03/2024.07.02.24309861/F3.medium.gif)

[Figure 3:](http://medrxiv.org/content/early/2024/07/03/2024.07.02.24309861/F3)

Figure 3: 
Identification of DELncRNAs in HAE sequencing from lesion tissue and surrounding tissue.

(A) Heatmap of DElncRNAs. (B) Volcano plot of DElncRNAs. Red represents upregulated genes, blue represents downregulated genes, and gray represents genes without significant differences.

View this table:
[Table 4](http://medrxiv.org/content/early/2024/07/03/2024.07.02.24309861/T4)

Table 4 Top 10 up-regulated and down-regulated lncRNAs.

### DEmRNAs in HAE

Most lncRNAs can act as miRNA sponges, indirectly regulating downstream target genes. To explore the potential “ceRNA” mechanism of lncRNAs in HAE and to reveal the differentially expressed mRNAs in the same tissues, our research group also conducted mRNA sequencing. A total of 870 mRNAs (743 upregulated and 127 downregulated) were identified as differentially expressed RNAs (DERNAs) in HAE (*P* <0.05, FC >2 / <0.5). These differentially expressed mRNAs are displayed as a heatmap (Figure 4A) and a volcano plot (Figure 4B). The top 10 up-regulated and down-regulated DEmRNAs with their names, FoldChange, Log2FoldChange and P-Value are listed in Table 5. All DEmRNAs with their names, Fold Change, Log2FoldChange and P-Value are listed in Supplementary Table 2.

![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/07/03/2024.07.02.24309861/F4.medium.gif)

[Figure 4:](http://medrxiv.org/content/early/2024/07/03/2024.07.02.24309861/F4)

Figure 4: 
Identification of DEmRNAs in HAE sequencing from lesion tissue and surrounding tissue. (A) Heatmap of DEmRNAs. (B) Volcano plot of DEmRNAs. Red represents upregulated genes, blue represents downregulated genes, and gray represents genes without significant differences.

View this table:
[Table 5](http://medrxiv.org/content/early/2024/07/03/2024.07.02.24309861/T5)

Table 5 Top 10 up-regulated and down-regulated mRNAs.

### Construction of the ceRNA Network in HAE

As previously mentioned, lncRNAs exert their functions primarily by adsorbing miRNAs. Therefore, based on high-throughput sequencing analysis, we constructed a comprehensive lncRNA-miRNA-mRNA network to study the roles of differentially expressed lncRNAs and mRNAs in HAE. Through the interactions of lncRNA-miRNA and miRNA-mRNA, we obtained a total of 11 lncRNAs, 21 miRNAs, and 56 mRNAs. The constructed lncRNA-miRNA-mRNA ceRNA regulatory network includes 88 nodes and 183 edges (Figure 5A). It is well known that lncRNAs and mRNAs have co-expression patterns in the ceRNA network. Thus, we selected a hub lncRNA (degree ≥16) and its connected miRNAs and mRNAs to reconstruct a sub-network. As shown in (Figure 5B), the PCBP1-AS1-miRNA-mRNA sub-network consists of 1 lncRNA node, 16 miRNA nodes, 51 mRNA nodes, and 91 edges. The DElncRNA and its matching miRNA and 21 miRNAs and their matching DEmRNAs can be found in (Table 6 and Table7). The degree of the ceRNA network were listed in Supplementary Table 3.

![Figure 5:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/07/03/2024.07.02.24309861/F5.medium.gif)

[Figure 5:](http://medrxiv.org/content/early/2024/07/03/2024.07.02.24309861/F5)

Figure 5: 
Construction of lncRNA-miRNA-mRNA network in HAE. (A) lncRNA-miRNA-mRNA ceRNA network in HAE. (B) The lncRNA PCBP1-AS1 sub-network. Red nodes indicate up-regulated RNAs while blue nodes indicate down-regulated RNAs. lncRNAs are denoted by rectangles, miRNAs by light green diamonds, and mRNAs by ellipses.

View this table:
[Table 6](http://medrxiv.org/content/early/2024/07/03/2024.07.02.24309861/T6)

Table 6 DElncRNAs in the ceRNA network were targeted by miRNAs.

View this table:
[Table 7](http://medrxiv.org/content/early/2024/07/03/2024.07.02.24309861/T7)

Table 7 The 21 miRNAs with their target DEmRNAs in the ceRNA network.

### Functional Enrichment Analysis and PPI Network Construction

GO and KEGG enrichment analyses indicate that the DEmRNAs in the ceRNA network are significantly associated with the positive regulation of canonical Wnt signaling pathway, histone lysine methylation, and the positive regulation of the Wnt signaling pathway, among other biological processes (BP) (Figure 6A). The most enriched Cellular Component (CC) was the cellular microtubule (Figure 6B). The enrichment of Molecular Functions (MF) was primarily related to histone methyltransferase activity (H3-K9 specific), histone-lysine N-methyltransferase activity and microtubule binding (Figure 6C). DEmRNAs related to the ceRNA network were significantly enriched in 5 KEGG pathways, namely adherens junction, platinum drug resistance, insulin resistance, HIF-1 signaling pathway, and toxoplasmosis (Figure 6D). The top 10 results of the functional enrichment analysis are listed in the Supplementary Table 4. The enrichment analysis suggests that the HAE-specific ceRNA network may participate in the pathogenesis and development of HAE by regulating these biological processes and pathways.

![Figure 6:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/07/03/2024.07.02.24309861/F6.medium.gif)

[Figure 6:](http://medrxiv.org/content/early/2024/07/03/2024.07.02.24309861/F6)

Figure 6: 
GO and KEGG enrichment analysis of DEmRNAs in the ceRNA network. (Top10). (A) Bubble Plot of BP. (B) Bubble Plot of CC. (C) Bubble Plot of MF. (D) Bubble Plot of KEGG.

Additionally, further analysis was conducted on the interacting mRNAs within the ceRNA network, leading to the construction of a Protein-Protein Interaction (PPI) network. The PPI network comprises 19 nodes and 17 edges (Figure 7). These mRNAs are likely to play a significant role in the disease occurrence and progression of HAE and serve as subjects for subsequent research. Notably, genes with high network degree and comprehensive scores within the network, such as YAP1, ZEB1, PAFAH1B1, TJP1 and YWHA, were found to be primarily enriched in biological processes like “ Regulation of developmental process” and “Positive regulation of developmental process”. This suggests that these genes are crucial for the development and progression of HAE.

![Figure 7:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/07/03/2024.07.02.24309861/F7.medium.gif)

[Figure 7:](http://medrxiv.org/content/early/2024/07/03/2024.07.02.24309861/F7)

Figure 7: 
The PPI network of DEmRNAs within the ceRNA network. Circles represent DEmRNAs, lines represent interactions between encoded proteins, gradient from dark red to light yellow based on degree.

### Validation of screened lncRNA-PCBP1-AS1, miR-20b-5p and CAPRIN2 by qRT-PCR

We utilized qRT-PCR to analyze the expression levels of PCBP1-AS1, miR-20b-5p, and CAPRIN2 in 13 HAE lesion tissues, as well as in tissues less than 1cm and more than 1cm surrounding the lesion. The qRT-PCR results indicated that PCBP1-AS1 and CAPRIN2 were highly expressed (*P* <0.05) in the lesion tissues (Figure 8A and Figure 8B), which is consistent with the results from high-throughput sequencing; whereas miR-20b-5p exhibited low expression(*P* <0.05) in the HAE lesion tissues (Figure 8C).

![Figure 8:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/07/03/2024.07.02.24309861/F8.medium.gif)

[Figure 8:](http://medrxiv.org/content/early/2024/07/03/2024.07.02.24309861/F8)

Figure 8: 
Expression of Candidate gene in HAE tissues assessed by qRT-PCR.(A) PCBP1-AS1 expression in HAE;(B) CAPRIN2 expression in HAE;(C) miR-20b-5p expression in HAE. The data are shown as the Median and interquartile range (IQR). n=13. **P* <0.05, ***P* <0.01, \***|*P* <0.001.

### miR-20b-5p is directly regulated by PCBP1-AS1 and interacts with CAPRIN2

In the ceRNA network, there is a targeted binding relationship among PCBP1-AS1, miR-20b-5p, and CAPRIN2. The sequencing data show a consistent expression trend for PCBP1-AS1 and CAPRIN2, and qRT-PCR indicates that both PCBP1-AS1 and CAPRIN2 are overexpressed in the lesion tissues, while miR-20b-5p is underexpressed in the same tissues. Subsequently, we used the RNAhybrid ([https://bibiserv.cebitec.uni-bielefeld.de/rnahybrid](https://bibiserv.cebitec.uni-bielefeld.de/rnahybrid)) and Targetscan databases([https://www.targetscan.org/](https://www.targetscan.org/)) to predict potential binding sites between PCBP1-AS1 and miR-20b-5p, as well as between miR-20b-5p and CAPRIN2 (Figure 9A and Figure 9B), and calculated the minimum free energy of these binding sites using RNAhybrid (Figure 9C). In order to further confirm the predicted results, we transfected them into 293T cells. The results showed that miR-20b-5p could decrease the fluorescence values in cells transfected with PCBP1-AS1-WT and CAPRIN2-WT. In contrast, there was no significant change in fluorescence values in cells transfected with PCBP1-AS1-MUT and CAPRIN2-MUT (Figure 9D and Figure 9E). These results suggest that PCBP1-AS1 directly targets miR-20b-5p, and miR-20b-5p interacts with CAPRIN2, with miR-20b-5p expression being negatively correlated with that of PCBP1-AS1 and CAPRIN2.

![Figure 9:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/07/03/2024.07.02.24309861/F9.medium.gif)

[Figure 9:](http://medrxiv.org/content/early/2024/07/03/2024.07.02.24309861/F9)

Figure 9: 
(A) The binding sites between miR-20b-5p and PCBP1-AS1;(B) The binding sites between miR-20b-5p and CAPRIN2;(C) The interaction between THAP9-AS1 and miR-484 and the minimum free energy (mfe) was predicted with RNAhybrid software;(D) Results of the dual-luciferase assay for miR-20b-5p and PCBP1-AS1;(E) Results of the dual-luciferase assay for miR-20b-5p and CAPRIN2. The data are shown as the Mean ± standard deviation. (n=3). \***|*P* <0.001, \**\*|\**P* <0.0001.

### Analysis of CAPRIN2 expression in HAE lesions and Clinicopathological characteristics

Through qRT-PCR and dual-luciferase reporter assays, the targeted regulatory relationship among PCBP1-AS1, miR-20b-5p, and CAPRIN2 has been demonstrated. Specifically, PCBP1-AS1 can affect the expression of the CAPRIN2 gene by adsorbing miR-20b-5p. To further investigate the association between this regulatory relationship and the clinical characteristics of HAE patients, we compared the expression levels of CAPRIN2 with various clinical parameters of the patients. The analysis revealed that overexpression of CAPRIN2 was not significantly correlated with the patient’s sex, age, multiple intrahepatic lesions, or distant metastasis of HAE (*P* > 0.05), but it was significantly correlated with the size of HAE lesions and intrahepatic vascular invasion (*P* <0.05) (Figure 10), indicating that patients with overexpressed CAPRIN2 may have a faster disease progression rate than those with low CAPRIN2 expression. This could imply that the expression level of CAPRIN2 could serve as a potential biomarker for predicting the prognosis of HAE patients.

![Figure 10](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/07/03/2024.07.02.24309861/F10.medium.gif)

[Figure 10](http://medrxiv.org/content/early/2024/07/03/2024.07.02.24309861/F10)

Figure 10 
The expression of CAPRIN2 correlated significantly with Clinicopathological characteristics.

These findings suggest that the regulatory axis of PCBP1-AS1-miR-20b-5p/CAPRIN2 may play a key role in controlling the progression and vascular invasion of HAE, which is of great significance for understanding the molecular mechanisms of HAE and developing new therapeutic strategies. Future research can further explore the specific mechanisms of action of this molecular axis and assess its feasibility as a potential therapeutic target.

## Discussion

HAE is a zoonotic disease with severe clinical manifestations and poor prognosis. Currently, radical surgery combined with albendazole medication remains the most ideal treatment option for HAE patients[25][26]. In recent years, numerous studies have gradually revealed the ability of LncRNAs to regulate the expression of target genes through various pathways, playing a significant role in promoting the occurrence, development, and metastasis of tumors. Moreover, the ceRNA hypothesis has proposed an innovative regulatory mechanism, where lncRNAs act as ceRNAs and form extensive regulatory networks by competitively binding to endogenous miRNAs. The imbalance of such ceRNA regulatory networks can significantly affect the biological functions of an organism, leading to the onset of diseases. Additionally, multiple studies have confirmed that ceRNA network play a crucial role in the development of various tumors. However, in HAE, the ceRNA network mediated by lncRNAs has not yet been comprehensively studied, and many ceRNA molecules still await discovery and understanding.

In this study, we utilized high-throughput sequencing technology to precisely identify DElncRNAs and DEmRNAs in HAE lesion tissues and their adjacent tissues. We attempted to screen for ceRNA molecules with specific functions using the hierarchical action model of lncRNA-miRNA-mRNA. By integrating the interactions between miRNAs, DEmRNAs, and DElncRNAs, we constructed an HAE-characteristic ceRNA regulatory network, which includes 89 key nodes and 183 interacting edges. Furthermore, we identified hub LncRNAs and built a sub-network centered around PCBP1-AS1.

Furthermore, we conducted GO and KEGG enrichment analyses to delve deeper into the biological pathways and functions involved in the ceRNA network for DEmRNAs. The GO enrichment analysis revealed that DEmRNAs in the ceRNA network are primarily involved in the positive regulation of the Wnt signaling pathway and histone lysine methylation, key biological processes that may be significantly associated with the occurrence of HAE. At the molecular level, the DEmRNAs are mainly associated with histone methyltransferase activity (specifically H3-K9) and histone lysine N-methyltransferase activity. Notably, DNA methyltransferase 3A (DNMT3A) is related to histone methyltransferase activity, and DNA methyltransferases are enzymes responsible for methylating DNA in mammals, which can lead to gene silencing[27]. Moreover, a multitude of histone lysine methyltransferases catalyze various lysine methylation events, decorating the core histones, and their aberrant expression is frequently observed in cancers, developmental disorders, and other pathologies[28]. In simple terms, HAE may be related to abnormalities in the intracellular epigenetic regulatory mechanisms, which affect gene expression and cellular function.

The KEGG pathway enrichment analysis results further revealed the roles of DEmRNAs involved in the ceRNA network across multiple biological pathways, such as adherens junction, platinum drug resistance, insulin resistance, HIF-1 signaling pathway, and toxoplasmosis. The enrichment of these pathways suggests that the pathogenesis of HAE may involve several complex biological processes. It is particularly noteworthy that adherens junctions play a critical role in maintaining the integrity of migrating cell groups, facilitating cell coordination, and cell rearrangement[29]. They are also implicated in the collective invasion of cancer cells[30], and adherens junctions mediate adhesive forces between adjacent cells[31], initiating key intracellular signals that regulate major cellular functions, including proliferation, differentiation, and migration[32]. Interestingly, the pathogenesis of HAE may be related to the disruption of Kupffer cells[33]and metabolic pathways[34], indicating a significant association between adherens junctions and HAE.

To systematically analyze the relationships and functions of DEmRNAs involved in the ceRNA network in HAE, we constructed a PPI network. Through this analysis, we identified several genes with high comprehensive scores, such as YAP1, ZEB1, PAFAH1B1, TJP1 and YWHAZ, which are primarily enriched in biological processes related to the regulation of developmental processes. Simultaneously, we observed that YAP1 and YWHAZ have a higher number of nodes in the PPI network, which typically indicates that DEmRNAs with more nodes play a more central role within the network. Research has indicated that YAP1 is one of the most critical effectors of the Hippo pathway and interacts with other pro-oncogenic pathways, promoting cancer development in various ways, including the promotion of malignant phenotypes, expansion of cancer stem cells, and enhancement of cancer cell drug resistance[35][36]. Additionally, YWHAZ has been closely associated with the progression of and as a therapeutic target for various types of cancer, including hepatocellular carcinoma[37], gastric cancer[38], and pancreatic cancer[39]. All in all, these genes are highly relevant to the occurrence and progression of related cancers and merit further research to clarify their roles in HAE.

Compelling evidence has shown that lncRNAs can compete for binding sites on miRNAs, thereby regulating the expression of target mRNAs. Within the ceRNA network, we discovered that Lnc-PCBP1-AS1 plays a pivotal role, which may influence variable targets through 16 lncRNA-associated miRNAs. Interestingly, a number of miRNAs, including miR-363[40], miR-125a-5p[41], miR-193a[42], miR-4465[43], miR-4319[44], miR-125b-5p[45], miR-22-3p[46], miR-107[47], miR-17-5p[48], miR-20b-5p[49], and miR-24-3p[50], have all been associated with the initiation and progression of hepatocellular carcinoma. Furthermore, miR-20b-5p can be sponged by the lncRNA WWOX-AS1, which upregulates WWOX and thereby inhibits the progression of hepatocellular carcinoma[51]. The lncRNA MALAT1 can also function as a ceRNA, regulating TXNIP by sequestering miR-20b-5p[52]. Consequently, the predicted miRNAs can be utilized to identify lncRNAs that play a significant role in gene regulation in HAE, thereby enhancing our comprehension of the ceRNA network.

Within the ceRNA network, PCBP1, miR-20b-5p, and CAPRIN2 exhibit a high degree of intramodular connectivity, suggesting that they may play a pivotal role in regulating gene expression and biological functions. So, they have been selected as candidate genes for further research. It has been reported that PCBP1-AS1 is closely associated with the proliferation and metastasis of hepatocellular carcinoma[53]. Additionally, the dysregulation of miR-20b-5p has been confirmed in hepatocellular carcinoma, where it negatively regulates CPEB3, affecting the progression of the disease[54]. Concurrently, our research has discovered that PCBP1-AS1 shares miR-20b-5p response elements with CAPRIN2 and regulates the expression of CAPRIN2 by spongeing miR-20b-5p. Also,Experiments have validated that CAPRIN2 is indeed a genuine target of miR-20b-5p.Besides, the target gene CAPRIN2 plays a crucial role in the Wnt signaling pathway[55], abnormal Wnt signaling is the basis for a variety of human diseases[56], and dysregulation of the Wnt signaling pathway has been proven to be associated with a range of cancers[57]. Moreover, numerous preclinical experiments have indicated that inhibiting Wnt signaling can prevent the growth and survival of cancer cells[58]. The aforementioned GO enrichment results indicate that DEmRNAs involved in the ceRNA network in HAE are primarily engaged in the regulation of the Wnt signaling pathway. This suggests that the PCBP1-AS1-miR-20b-5p/CAPRIN2 axis may play a key role in HAE through the Wnt signaling pathway. Concurrently, we have analyzed the target gene CAPRIN2 in relation to clinical characteristics of patients and found that overexpression of CAPRIN2 is significantly associated with the size of lesions and the extent of vascular invasion in the liver in HAE. This suggests that PCBP1-AS1 might control CAPRIN2 expression in HAE by binding miR-20b-5p, affecting the HAE start and spread. However, the specific mechanisms by which the PCBP1-AS1-miR-20b-5p/CAPRIN2 axis modulates the development of HAE require further experimental validation. Considering that PCBP1-AS1 is a hub lncRNA in the ceRNA network, it could potentially serve as a target for HAE diagnosis and treatment. Overall, our research enhances the understanding of the regulatory mechanisms of the PCBP1-AS1-miR-20b-5p/CAPRIN2 axis in HAE, which may provide valuable insights to advance HAE research.

In conclusion, we have constructed a ceRNA network using lncRNA-miRNA and miRNA-mRNA interactions in HAE, laying the groundwork for further investigation into the regulatory mechanisms of HAE. However, this study still has some limitations, such as a lack of extensive experiments and a small sample size. Subsequently, we will expand our collection of clinical data based on these findings and conduct a more detailed analysis of the pathogenesis of HAE to clarify the roles of these differentially expressed genes within the ceRNA regulatory network. Equally, we will also focus on how the lnc-PCBP1-AS1-miR-20b-5p/CAPRIN2 axis regulates the progression of HAE, its functional and mechanistic aspects, which could lay the foundation for the diagnosis and treatment of HAE.

## Author contributions

All authors contributed to the article and approved the submitted version.

## Funding

This research was supported in part by the Qinghai Provincial Science and Technology Department under the project “Research on differential expression of LncRNA and circRNA and related signaling pathways in hepatic alveolar echinococcosis.” (Project No. 2022-ZJ-747) and “KunLun talents High-end Innovation and Entrepreneurship Talent Program” of Qinghai Province (Youth Talent character [2021] No. 13).

## Data availability

The experimental data has been deposited in the Gene Expression Omnibus (GEO) database under the accession number GSE261626.

## 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.

## Supporting information

Supplemental Data 1 [[supplements/309861_file02.docx]](pending:yes)

Supplemental Data 2 [[supplements/309861_file03.docx]](pending:yes)

Supplemental Data 3 [[supplements/309861_file04.docx]](pending:yes)

Supplemental Data 4 [[supplements/309861_file05.docx]](pending:yes)

*   Received July 2, 2024.
*   Revision received July 2, 2024.
*   Accepted July 3, 2024.


*   © 2024, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/)

## Reference

1.  [1].McManus DP, Zhang W, Li J, Bartley PB. Echinococcosis. Lancet. 2003 Oct 18;362(9392):1295-304.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(03)14573-4&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=14575976&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F03%2F2024.07.02.24309861.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000186036900024&link_type=ISI) 

2.  [2].Xu X, Qian X, Gao C, Pang Y, Zhou H, Zhu L, Wang Z, Pang M, Wu D, Yu W, Kong F, Shi D, Guo Y, Su X, Hu W, Yan J, Feng X, Fan H. Advances in the pharmacological treatment of hepatic alveolar echinococcosis: From laboratory to clinic. Front Microbiol. 2022 Aug 8; 13:953846.
    
    

3.  [3].Nunnari G, Pinzone MR, Gruttadauria S, Celesia BM, Madeddu G, Malaguarnera G, Pavone P, Cappellani A, Cacopardo B. Hepatic echinococcosis: clinical and therapeutic aspects. World J Gastroenterol. 2012 Apr 7;18(13):1448–58.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3748/wjg.v18.i13.1448&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22509076&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F03%2F2024.07.02.24309861.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000302501200004&link_type=ISI) 

4.  [4].Zheng H, Zhang W, Zhang L, Zhang Z, Li J, Lu G, Zhu Y, Wang Y, Huang Y, Liu J, Kang H, Chen J, Wang L, Chen A, Yu S, Gao Z, Jin L, Gu W, Wang Z, Zhao L, Shi B, Wen H, Lin R, Jones MK, Brejova B, Vinar T, Zhao G, McManus DP, Chen Z, Zhou Y, Wang S. The genome of the hydatid tapeworm Echinococcus granulosus. Nat Genet. 2013 Oct;45(10):1168–75.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.2757&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24013640&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F03%2F2024.07.02.24309861.atom) 

5.  [5].Kopp F, Mendell JT. Functional Classification and Experimental Dissection of Long Noncoding RNAs. Cell. 2018 Jan 25;172(3):393–407.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2018.01.011&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29373828&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F03%2F2024.07.02.24309861.atom) 

6.  [6].Djebali S, Davis CA, Merkel A, Dobin A, Lassmann T, Mortazavi A, Tanzer A, Lagarde J, Lin W, Schlesinger F, Xue C, Marinov GK, Khatun J, Williams BA, Zaleski C, Rozowsky J, Röder M, Kokocinski F, Abdelhamid RF, Alioto T, Antoshechkin I, Baer MT, Bar NS, Batut P, Bell K, Bell I, Chakrabortty S, Chen X, Chr ast J, Curado J, Derrien T, Drenkow J, Dumais E, Dumais J, Duttagupta R, Falconnet E, Fastuca M, Fejes-Toth K, Ferreira P, Foissac S, Fullwood MJ, Gao H, Gonzalez D, Gordon A, Gunawardena H, Howald C, Jha S, Johnson R, Kapranov P, King B, Kingswood C, Luo OJ, Park E, Persaud K, Preall JB, Ribeca P, Risk B, Robyr D, Sammeth M, Schaffer L, See LH, Shahab A, Skancke J, Suzuki AM, Takahashi H, Tilgner H, Trout D, Walters N, Wang H, Wrobel J, Yu Y, Ruan X, Hayashizaki Y, Harrow J, Gerstein M, Hubbard T, Reymond A, Antonarakis SE, Hannon G, Giddings MC, Ruan Y, Wold B, Carninci P, Guigó R, Gingeras TR. Landscape of transcription in human cells. Nature. 2012 Sep 6;489(7414):101-8.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature11233&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22955620&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F03%2F2024.07.02.24309861.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000308347000043&link_type=ISI) 

7.  [7].Slack FJ, Chinnaiyan AM. The Role of Non-coding RNAs in Oncology. Cell. 2019 Nov 14;179(5):1033–1055.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/J.CELL.2019.10.017&link_type=DOI) 

8.  [8].Herman AB, Tsitsipatis D, Gorospe M. Integrated lncRNA function upon genomic and epigenomic regulation. Mol Cell. 2022 Jun 16;82(12):2252–2266.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.molcel.2022.05.027&link_type=DOI) 

9.  [9].Grimson A, Farh KK, Johnston WK, Garrett-Engele P, Lim LP, Bartel DP. MicroRNA targeting specificity in mammals: determinants beyond seed pairing. Mol Cell. 2007 Jul 6;27(1):91–105.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.molcel.2007.06.017&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17612493&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F03%2F2024.07.02.24309861.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000248088500009&link_type=ISI) 

10. [10].Baek D, Villén J, Shin C, Camargo FD, Gygi SP, Bartel DP. The impact of microRNAs on protein output. Nature. 2008 Sep 4;455(7209):64-71.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature07242&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18668037&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F03%2F2024.07.02.24309861.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000258890200037&link_type=ISI) 

11. [11].Zhang H, Li Y, Lai M. The microRNA network and tumor metastasis. Oncogene. 2010 Feb 18;29(7):937–48.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/onc.2009.406&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19935707&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F03%2F2024.07.02.24309861.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000274604400001&link_type=ISI) 

12. [12].Thomson DW, Dinger ME. Endogenous microRNA sponges: evidence and controversy. Nat Rev Genet. 2016 May;17(5):272–83.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrg.2016.20&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27040487&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F03%2F2024.07.02.24309861.atom) 

13. [13].Tay Y, Rinn J, Pandolfi PP. The multilayered complexity of ceRNA crosstalk and competition. Nature. 2014 Jan 16;505(7483):344-52.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature12986&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24429633&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F03%2F2024.07.02.24309861.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000329621800032&link_type=ISI) 

14. [14].Salmena L, Poliseno L, Tay Y, Kats L, Pandolfi PP. A ceRNA hypothesis: the Rosetta Stone of a hidden RNA language? Cell. 2011 Aug 5;146(3):353–8.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2011.07.014&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21802130&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F03%2F2024.07.02.24309861.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000293570500010&link_type=ISI) 

15. [15].Rong Z, Wang Z, Wang X, Qin C, Geng W. Molecular interplay between linc01134 and YY1 dictates hepatocellular carcinoma progression. J Exp Clin Cancer Res. 2020 Apr 9;39(1):61.
    
    

16. [16].Li N, Yang G, Luo L, Ling L, Wang X, Shi L, Lan J, Jia X, Zhang Q, Long Z, Liu J, Hu W, He Z, Liu H, Liu W, Zheng G. lncRNA THAP9-AS1 Promotes Pancreatic Ductal Adenocarcinoma Growth and Leads to a Poor Clinical Outcome via Sponging miR-484 and Interacting with YAP. Clin Cancer Res. 2020 Apr 1;26(7):1736–1748.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTA6ImNsaW5jYW5yZXMiO3M6NToicmVzaWQiO3M6OToiMjYvNy8xNzM2IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDcvMDMvMjAyNC4wNy4wMi4yNDMwOTg2MS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

17. [17].Yang XZ, Cheng TT, He QJ, Lei ZY, Chi J, Tang Z, Liao QX, Zhang H, Zeng LS, Cui SZ. LINC01133 as ceRNA inhibits gastric cancer progression by sponging miR-106a-3p to regulate APC expression and the Wnt/β-catenin pathway. Mol Cancer. 2018 Aug 22;17(1):126.
    
    

18. [18].He F, Song Z, Chen H, Chen Z, Yang P, Li W, Yang Z, Zhang T, Wang F, Wei J, Wei F, Wang Q, Cao J. Long noncoding RNA PVT1-214 promotes proliferation and invasion of colorectal cancer by stabilizing Lin28 and interacting with miR-128. Oncogene. 2019 Jan;38(2):164–179.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41388-018-0432-8&link_type=DOI) 

19. [19].Maere S, Heymans K, Kuiper M. BiNGO: a Cytoscape plugin to assess overrepresentation of gene ontology categories in biological networks. Bioinformatics.2005 Aug 15;21(16):3448–9.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/bti551&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15972284&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F03%2F2024.07.02.24309861.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000231360600024&link_type=ISI) 

20. [20].Kong L, Zhang Y, Ye ZQ, Liu XQ, Zhao SQ, Wei L, Gao G. CPC: assess the protein-coding potential of transcripts using sequence features and support vector machine. Nucleic Acids Res. 2007 Jul;35(Web Server issue): W345-9.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkm391&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17631615&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F03%2F2024.07.02.24309861.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000255311500065&link_type=ISI) 

21. [21].Sun L, Luo H, Bu D, Zhao G, Yu K, Zhang C, Liu Y, Chen R, Zhao Y. Utilizing sequence intrinsic composition to classify protein-coding and long non-coding transcripts. Nucleic Acids Res. 2013 Sep;41(17): e166.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkt646&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23892401&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F03%2F2024.07.02.24309861.atom) 

22. [22].Finn RD, Bateman A, Clements J, Coggill P, Eberhardt RY, Eddy SR, Heger A, Hetherington K, Holm L, Mistry J, Sonnhammer EL, Tate J, Punta M. Pfam: the protein families database. Nucleic Acids Res. 2014 Jan;42(Database issue): D222-30.
    
    

23. [23].Li A, Zhang J, Zhou Z. PLEK: a tool for predicting long non-coding RNAs and messenger RNAs based on an improved k-mer scheme. BMC Bioinformatics. 2014 Sep 19;15(1):311.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1471-2105-15-311&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25239089&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F03%2F2024.07.02.24309861.atom) 

24. [24].Ponting CP, Oliver PL, Reik W. Evolution and functions of long noncoding RNAs. Cell. 2009 Feb 20;136(4):629–41.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2009.02.006&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19239885&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F03%2F2024.07.02.24309861.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000263688200015&link_type=ISI) 

25. [25].Maimaitinijiati Y, Meng Y, Chen X. Is long-term follow-up without surgical treatment a valid option for hepatic alveolar echinococcosis? World J Gastroenterol. 2022 Jun 28;28(24):2775–2777.
    
    

26. [26].Nunnari G, Pinzone MR, Gruttadauria S, Celesia BM, Madeddu G, Malaguarnera G, Pavone P, Cappellani A, Cacopardo B. Hepatic echinococcosis: clinical and therapeutic aspects. World J Gastroenterol. 2012 Apr 7;18(13):1448–58.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3748/wjg.v18.i13.1448&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22509076&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F03%2F2024.07.02.24309861.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000302501200004&link_type=ISI) 

27. [27].Fuks F, Hurd PJ, Deplus R, Kouzarides T. The DNA methyltransferases associate with HP1 and the SUV39H1 histone methyltransferase. Nucleic Acids Res. 2003 May 1;31(9):2305–12.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkg332&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12711675&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F03%2F2024.07.02.24309861.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000182627700014&link_type=ISI) 

28. [28].Husmann D, Gozani O. Histone lysine methyltransferases in biology and disease. Nat Struct Mol Biol. 2019 Oct;26(10):880–889.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41594-019-0298-7&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31582846&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F03%2F2024.07.02.24309861.atom) 

29. [29].Peglion F, Llense F, Etienne-Manneville S. Adherens junction treadmilling during collective migration. Nat Cell Biol. 2014 Jul;16(7):639–51.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ncb2985&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24929360&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F03%2F2024.07.02.24309861.atom) 

30. [30].Friedl P, Gilmour D. Collective cell migration in morphogenesis, regeneration and cancer. Nat Rev Mol Cell Biol. 2009 Jul;10(7):445–57.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrm2720&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19546857&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F03%2F2024.07.02.24309861.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000267266000009&link_type=ISI) 

31. [31].Smutny M, Yap AS. Neighborly relations: cadherins and mechanotransduction. J Cell Biol. 2010 Jun 28;189(7):1075–7.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamNiIjtzOjU6InJlc2lkIjtzOjEwOiIxODkvNy8xMDc1IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDcvMDMvMjAyNC4wNy4wMi4yNDMwOTg2MS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

32. [32].Harris TJ, Tepass U. Adherens junctions: from molecules to morphogenesis. Nat Rev Mol Cell Biol. 2010 Jul;11(7):502–14.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrm2927&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20571587&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F03%2F2024.07.02.24309861.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000280076000015&link_type=ISI) 

33. [33].Liu Y, Tian F, Shan J, Gao J, Li B, Lv J, Zhou X, Cai X, Wen H, Ma X. Kupffer Cells: Important Participant of Hepatic Alveolar Echinococcosis. Front Cell Infect Microbiol. 2020 Jan 29; 10:8.
    
    

34. [34].Lin C, Chen Z, Zhang L, Wei Z, Cheng KK, Liu Y, Shen G, Fan H, Dong J. Deciphering the metabolic perturbation in hepatic alveolar echinococcosis: a 1H NMR-based metabolomics study. Parasit Vectors. 2019 Jun 13;12(1):300.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13071-019-3554-0&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31196218&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F03%2F2024.07.02.24309861.atom) 

35. [35].Johnson R, Halder G. The two faces of Hippo: targeting the Hippo pathway for regenerative medicine and cancer treatment. Nat Rev Drug Discov. 2014 Jan;13(1):63–79.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrd4161&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24336504&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F03%2F2024.07.02.24309861.atom) 

36. [36].Pan D. The hippo signaling pathway in development and cancer. Dev Cell. 2010 Oct 19;19(4):491–505.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.devcel.2010.09.011&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20951342&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F03%2F2024.07.02.24309861.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000283522700004&link_type=ISI) 

37. [37].Zhao JF, Zhao Q, Hu H, Liao JZ, Lin JS, Xia C, Chang Y, Liu J, Guo AY, He XX. The ASH1-miR-375-YWHAZ Signaling Axis Regulates Tumor Properties in Hepatocellular Carcinoma. Mol Ther Nucleic Acids. 2018 Jun 1; 11:538–553.
    
    

38. [38].Liu J, Song S, Lin S, Zhang M, Du Y, Zhang D, Xu W, Wang H. Circ-SERPINE2 promotes the development of gastric carcinoma by sponging miR-375 and modulating YWHAZ. Cell Prolif. 2019 Jul;52(4): e12648.
    
    

39. [39].Cao F, Jiang Y, Chang L, Du H, Chang D, Pan C, Huang X, Yu D, Zhang M, Fan Y, Bian X, Li K. High-throughput functional screen identifies YWHAZ as a key regulator of pancreatic cancer metastasis. Cell Death Dis. 2023 Jul 14;14(7):431.
    
    

40. [40].Zhou P, Huang G, Zhao Y, Zhong D, Xu Z, Zeng Y, Zhang Y, Li S, He F. MicroRNA-363-mediated downregulation of S1PR1 suppresses the proliferation of hepatocellular carcinoma cells. Cell Signal. 2014 Jun;26(6):1347–54.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F03%2F2024.07.02.24309861.atom) 

41. [41].Kim JK, Noh JH, Jung KH, Eun JW, Bae HJ, Kim MG, Chang YG, Shen Q, Park WS, Lee JY, Borlak J, Nam SW. Sirtuin7 oncogenic potential in human hepatocellular carcinoma and its regulation by the tumor suppressors MiR-125a-5p and MiR-125b. Hepatology. 2013 Mar;57(3):1055–67.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/hep.26101&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23079745&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F03%2F2024.07.02.24309861.atom) 

42. [42].Salvi A, Conde I, Abeni E, Arici B, Grossi I, Specchia C, Portolani N, Barlati S, De Petro G. Effects of miR-193a and sorafenib on hepatocellular carcinoma cells. Mol Cancer. 2013 Dec 13; 12:162.
    
    

43. [43].Bu WJ, Fang Z, Li WL, Wang X, Dong MJ, Tao QY, Zhang L, Xu YQ. LINC00240 sponges miR-4465 to promote proliferation, migration, and invasion of hepatocellular carcinoma cells via HGF/c-MET signaling pathway. Eur Rev Med Pharmacol Sci. 2020 Oct;24(20):10452–10461.
    
    

44. [44].Han S, Shi Y, Sun L, Liu Z, Song T,  Liu Q. MiR-4319 induced an inhibition of epithelial-mesenchymal transition and prevented cancer stemness of HCC through targeting FOXQ1. Int J Biol Sci. 2019 Nov 15;15(13):2936-2947.
    
    

45. [45].Lin XH, Li DP, Liu ZY, Zhang S, Tang WQ, Chen RX, Weng SQ, Tseng YJ, Xue RY, Dong L. Six immune-related promising biomarkers may promote hepatocellular carcinoma prognosis: a bioinformatics analysis and experimental validation. Cancer Cell Int. 2023 Mar 23;23(1):52.
    
    

46. [46].Cui S, Chen Y, Guo Y, Wang X, Chen D. Hsa-miR-22-3p inhibits liver cancer cell EMT and cell migration/ invasion by indirectly regulating SPRY2. PLoS One. 2023 Feb 7;18(2): e0281536.
    
    

47. [47].Chen HA, Li CC, Lin YJ, Wang TF, Chen MC, Su YH, Yeh YL, Padma VV, Liao PH, Huang CY. Hsa-miR-107 regulates chemosensitivity and inhibits tumor growth in hepatocellular carcinoma cells. Aging (Albany NY). 2021 Apr 26;13(8):12046–12057.
    
    

48. [48].Yang F, Yin Y, Wang F, Wang Y, Zhang L, Tang Y, Sun S. miR-17-5p Promotes migration of human hepatocellular carcinoma cells through the p38 mitogen-activated protein kinase-heat shock protein 27 pathway. Hepatology. 2010 May;51(5):1614–23.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/hep.23566&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20209605&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F03%2F2024.07.02.24309861.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000277261400019&link_type=ISI) 

49. [49].Li Z, Wu L, Tan W, Zhang K, Lin Q, Zhu J, Tu C, Lv X, Jiang C. MiR-20b-5p promotes hepatocellular carcinoma cell proliferation, migration and invasion by down-regulating CPEB3. Ann Hepatol. 2021 Jul-Aug; 23:100345.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.aohep.2021.100345&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33812045&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F03%2F2024.07.02.24309861.atom) 

50. [50].Fan JC, Zeng F, Le YG, Xin L. LncRNA CASC2 inhibited the viability and induced the apoptosis of hepatocellular carcinoma cells through regulating miR-24-3p. J Cell Biochem. 2018 Aug;119(8):6391–6397.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/jcb.26479&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29091305&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F03%2F2024.07.02.24309861.atom) 

51. [51].Xu D, Liu X, Wu J, Wang Y, Zhou K, Chen W, Chen J, Chen C, Chen L. LncRNA WWOX-AS1 sponges miR-20b-5p in hepatocellular carcinoma and represses its progression by upregulating WWOX. Cancer Biol Ther. 2020 Oct 2;21(10):927–936.
    
    

52. [52].Li JZ, Ye LH, Wang DH, Zhang HC, Li TY, Liu ZQ, Dai EH, Li MR. The identify role and molecular mechanism of the MALAT1/hsa-mir-20b-5p/TXNIP axis in liver inflammation caused by CHB in patients with chronic HBV infection complicated with NAFLD. Virus Res. 2021 Jun; 298:198405.
    
    

53. [53].Luo T, Gao Y, Zhangyuan G, Xu X, Xue C, Jin L, Zhang W, Zhu C, Sun B, Qin X. lncRNA PCBP1-AS1 Aggravates the Progression of Hepatocellular Carcinoma via Regulating PCBP1/PRL-3/AKT Pathway. Cancer Manag Res. 2020 Jul 6; 12:5395–5408.
    
    

54. [54].Li Z, Wu L, Tan W, Zhang K, Lin Q, Zhu J, Tu C, Lv X, Jiang C. MiR-20b-5p promotes hepatocellular carcinoma cell proliferation, migration and invasion by down-regulating CPEB3. Ann Hepatol. 2021 Jul-Aug; 23:100345.
    
    

55. [55].Ding Y, Xi Y, Chen T, Wang JY, Tao DL, Wu ZL, Li YP, Li C, Zeng R, Li L. Caprin-2 enhances canonical Wnt signaling through regulating LRP5/6 phosphorylation. J Cell Biol. 2008 Sep 8;182(5):865–72.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamNiIjtzOjU6InJlc2lkIjtzOjk6IjE4Mi81Lzg2NSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzA3LzAzLzIwMjQuMDcuMDIuMjQzMDk4NjEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

56. [56].Clevers H, Nusse R. Wnt/β-catenin signaling and disease. Cell. 2012 Jun 8;149(6):1192–205.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2012.05.012&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22682243&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F03%2F2024.07.02.24309861.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000305119600006&link_type=ISI) 

57. [57].Zhan T, Rindtorff N, Boutros M. Wnt signaling in cancer. Oncogene. 2017 Mar;36(11):1461–1473.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/onc.2016.304&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27617575&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F03%2F2024.07.02.24309861.atom) 

58. [58].Polakis P. Wnt signaling in cancer. Cold Spring Harb Perspect Biol. 2012 May 1;4(5): a008052.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTE6ImNzaHBlcnNwZWN0IjtzOjU6InJlc2lkIjtzOjExOiI0LzUvYTAwODA1MiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzA3LzAzLzIwMjQuMDcuMDIuMjQzMDk4NjEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9)