Identification of Potential miRNA-mRNA Regulatory Axis of Intrahepatic cccDNA in Patients with Chronic Hepatitis B Virus Infection in the Grey Zone

authors:

avatar Wen-Ting Zhang 1 , avatar Antonio Gil Gómez 2 , avatar Macarena López Sánchez 3 , avatar Pérez-Del-Pulgar Sofia 4 , avatar Manuel Romero Gómez 2 , avatar Shuai-Shuai Gao 5 , *

Xijing Hospital of Digestive Disease, Air Force Medical University, Xi'an, China
Institute of Biomedicine of Seville, Sevilla, Spain
Clinical Bioinformatics Area, Virgen del Rocío University Hospital, Sevilla, Spain
University of Barcelona, Barcelona, Spain
Xi'an Daxing Hospital, Xi'an, China

How To Cite Zhang W, Gómez A G, Sánchez M L, Sofia P, Gómez M R, et al. Identification of Potential miRNA-mRNA Regulatory Axis of Intrahepatic cccDNA in Patients with Chronic Hepatitis B Virus Infection in the Grey Zone. Hepat Mon. 2023;23(1):e133782. https://doi.org/10.5812/hepatmon-133782.

Abstract

Background:

Intrahepatic covalently closed circular DNA (cccDNA) plays a critical role in the life cycle of the hepatitis B virus (HBV). Growing evidence suggests that microRNAs (miRNAs) may regulate cccDNA expression and contribute to the natural history of chronic hepatitis B (CHB).

Objectives:

This study aimed to investigate potential miRNA-mRNA regulatory axes of intrahepatic cccDNA in CHB-GZ patients and to identify new therapeutic targets.

Methods:

Thirteen CHB-GZ patients were included and divided into two groups based on cccDNA levels: Reference group (n = 7) with cccDNA < 1 copy/cell and control group (n = 6) with cccDNA ≥ 1 copy/cell. Transcriptome-wide miRNA and mRNA expression profiles in liver tissue were determined. Differentially expressed miRNAs (DE-miRNAs) and mRNAs (DE-mRNAs) were defined by |logFC| > log1.5 and P < 0.05. Enrichment analyses of gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were performed. Candidate miRNA-mRNA interaction pairs were acquired from miRTarBase, and a miRNA-mRNA network was constructed using Cytoscape based on the Spearman correlation coefficient (r).

Results:

We identified 19 DE-miRNAs and 340 DE-mRNAs. The most enriched GO terms were related to biological processes that regulate the virus life cycle, viral process, and viral genome replication. KEGG pathway enrichment analysis suggested that these predicted targets were associated with hepatitis B. Finally, we found a high correlation between miR-4295 and ZNF224, which suggests that they may form a potential regulatory axis of intrahepatic cccDNA in CHB-GZ patients.

Conclusions:

Our study suggests that miR-4295 and ZNF224 may be the potential regulatory axis of intrahepatic cccDNA in patients with CHB-GZ.

1. Background

Hepatitis B virus (HBV) infection remains a significant public health burden worldwide, affecting over 257 million people, and is closely associated with developing liver cirrhosis and hepatocellular carcinoma (1). HBV is a small enveloped DNA virus that establishes its genome in the nucleus of infected hepatocytes by forming the HBV DNA minichromosome, known as the covalently closed circular DNA (cccDNA) (2). The cccDNA serves as a template for the transcription of all viral RNAs, playing a critical role in the viral life cycle and contributing to the persistence of infection and relapse (3). Although effective treatment strategies based on interferon and nucleos(t)ide analogs suppress HBV replication, they do not target cccDNA, challenging the complete eradication of HBV (4). Therefore, the only way to cure chronic HBV infection is to eliminate cccDNA from infected hepatocytes.

MicroRNAs (miRNAs) are a large class of small non-coding RNAs (21 - 25 nucleotides) that regulate gene expression post-transcriptionally by targeting mRNAs (5). Growing evidence suggests that miRNAs may regulate intrahepatic cccDNA expression and play a role in the natural history of chronic hepatitis B (CHB). They can directly affect HBV replication by binding to HBV transcripts or indirectly by targeting cellular factors involved in HBV replication (6). For example, miR-548ah promotes HBV replication and expression by targeting histone deacetylase 4 (HDAC4) (7), which may be related to the deacetylation state of histone binding to cccDNA. Later, Zhang et al. showed that miR-449a directly targeted cAMP-responsive element binding protein 5 (CREB5), which, in turn, induces the expression of farnesoid X receptor α (FXRα). This factor regulates cccDNA activity in specific HBV replication systems (8). While most CHB patients belong to the five phases in the natural history of HBV defined by the EASL in 2017, some individuals do not meet any of the usual immune states and are known to be in the "grey zone (GZ)" (9).

2. Objectives

The objective was to investigate the regulatory role of miRNAs on cccDNA expression in CHB-GZ patients, establish a miRNA-mRNA interaction network, and provide a theoretical basis for finding new therapeutic targets.

3. Methods

3.1. Patients and Liver Tissue Samples

Liver tissues were collected from 13 CHB-GZ patients at the Digestive Diseases Unit of the Virgen del Rocío University Hospital (HUVR). Patients with other liver diseases and viral infections were excluded from the study. The research was approved by the Ethics Committee of the Virgen Macarena-Virgen del Rocío University Hospitals (Approval ID: CEI-HUVM-VR-30/04/2019), and each patient signed an informed consent form. The biopsies were snap-frozen in liquid nitrogen and stored at -80°C until use.

3.2. RNA Extraction

Total or small RNA extractions from tissue samples were performed using the mirVana™ miRNA Isolation Kit (Ambion® Life Technologies, CA, USA) according to the manufacturer's instructions. The total RNA concentration and purity of each specimen were quantified using a NanoDrop™ 2000c spectrophotometer (Thermo Scientific). The purity was estimated by the absorbance ratio at 260 to 280 nm.

3.3. Quantification of Intrahepatic cccDNA

Intrahepatic DNA was isolated from liver tissues using the MasterPure™ DNA and RNA purification kit (Epicenter, Illumina, San Diego, CA, USA) according to the manufacturer’s instructions. The amount of purified DNA was measured using a spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA). Total intrahepatic cccDNA was determined by real-time PCR, and quantification was performed using the ViiA 7 real-time PCR system mediated by a pair of cccDNA-selective primers and probe. The primers and the probe are shown in Table 1. The quantification range is between 5 × 101 and 5 × 107 copies per reaction for intrahepatic cccDNA. The copies of the viral genome were normalized against the number of cells, estimated by the copy number of the human β-globin gene.

Table 1.

The Primers and the Probe of Intrahepatic cccDNA Quantification

cccDNASequencePosition
ForwardCCGTGTGCACTTCGCTTCA1577 – 1595
ReverseGCACAGCTTGGAGGCTTGA1884 – 1866
ProbeCATGGAGACCACCGTGAACGCCC1609 – 1631

3.4. Microarray Profiling of miRNA and mRNA Expression

Transcriptome-wide miRNA and mRNA expression profiles were determined using the GeneChip™ miRNA 4.0 Assay (Affymetrix, ThermoFisher CA, USA) and the Clariom™ S Assay, human, respectively. Sample labeling, microarray hybridization, and washing were performed as recommended by the manufacturer. Finally, the arrays were scanned using the Affymetrix GCS 3000 scanner (Affymetrix, Santa Clara, CA, USA), and signal values were evaluated using the Affymetrix® GeneChip™ Command Console Software (10).

3.5. Screening of Candidate DE-miRNAs and DE-mRNAs

Thirteen CHB-GZ patients' liver tissues were divided into a reference group (n = 7) with cccDNA < 1 copy/cell and a control group (n = 6) with cccDNA ≥ 1 copy/cell, based on the level of cccDNA intrahepatic. We used one copy per cell of cccDNA intrahepatic as the grouping standard because the mean value of cccDNA in HBV-infected patients is 0.01 to 1 copy per cell (11). Moreover, Liang et al. revealed that baseline HBV intrahepatic cccDNA > one copy/cell independently predicts liver inflammation (12). Differentially expressed miRNAs (DE-miRNAs) and mRNAs (DE-mRNAs) between the two groups were selected using the "limma” package in R software. Raw data were initially normalized using the quantile algorithm. The threshold points for DE-miRNAs and DE-mRNAs were |logFC| > log1.5 and P < 0.05.

3.6. Functional Enrichment Analysis

We used the “ClusterProfiler” package in R to perform functional enrichment analysis on the gene ontology (GO) (13) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) (14) pathway to determine the roles played by DE-mRNAs. The GO analysis included three categories: Biological process (BP), cellular component (CC), and molecular function (MF). Statistically significant differences for the enrichment analysis were defined as P < 0.05.

3.7. Construction of miRNA-mRNA Network

We used miRTarBase (15) to investigate the interaction of candidate DE-miRNAs with DE-mRNAs. The miRNA-mRNA network was constructed and analyzed using Cytoscape software (16) (R package). We calculated the Spearman correlation coefficient (r) between miRNA-mRNA target genes. A strong correlation was considered when r < -0.7 and P < 0.01.

3.8. Statistical Analysis

Quantitative data were expressed as mean ± standard deviation. We used SPSS (version 25.0) and R software (version 4.0.5) for statistical analyses. Spearman correlation coefficients were calculated to evaluate the correlations. We considered P < 0.05 to be statistically significant.

4. Results

4.1. Heterogeneity Test and Screening of Candidate DE-miRNAs and DE-mRNAs

The robust multichip average (RMA) normalization method was used to standardize the data of DE-miRNAs and DE-mRNAs. Then, the “limma" package of the Bioconductor analysis tool was used to make preprocessing and normalized box plots. In comparing CHB-GZ patients with intrahepatic cccDNA < 1 copy/cell and ≥ 1 copy/cell, 19 DE-miRNAs were selected, including 12 down-regulated miRNAs and seven up-regulated. Among them, miR-4295 and miR-4271 were down-regulated, while miR-1281, miR-8071, and miR-575 were up-regulated. A total of 340 DE-mRNAs were obtained, among which 160 were down-regulated and 180 were up-regulated. Among them, ZNF224 and SMYD5 were up-regulated, whereas KPNA2, GPN2, TMEM184B, and POLD3 were down-regulated. The identification of DE-miRNAs and DE-mRNAs is shown in Figure 1.

Identification of DE-miRNAs and DE-mRNAs. A, Box plot of preprocessed miRNA data. B, Box plot of normalized miRNA data. C, Volcano plot of DE-miRNAs. D, Heatmap of DE-miRNAs. E, Box plot of preprocessed mRNA data. F, Box plot of normalized mRNA data. G, Volcano plot of DE-mRNAs. H, Heatmap of DE-mRNAs.
Identification of DE-miRNAs and DE-mRNAs. A, Box plot of preprocessed miRNA data. B, Box plot of normalized miRNA data. C, Volcano plot of DE-miRNAs. D, Heatmap of DE-miRNAs. E, Box plot of preprocessed mRNA data. F, Box plot of normalized mRNA data. G, Volcano plot of DE-mRNAs. H, Heatmap of DE-mRNAs.

4.2. Functional and Pathway Enrichment Analysis

The GO terms were significantly enriched in the biological process of the virus, such as the regulation of the viral life cycle, viral process, viral genome replication, and regulation of the kinase I-kappa B/NF-kappa B. The most enriched MF category in the GO terms was "immune receptor activity," and in the CC category, it was the "outer side of the plasma membrane." KEGG pathway enrichment analysis revealed that hepatitis B and influenza A infections were significantly enriched pathways (Figure 2).

GO enrichment analysis (BP, MF, CC) and KEGG pathway enrichment analysis.
GO enrichment analysis (BP, MF, CC) and KEGG pathway enrichment analysis.

4.3. Construction of miRNA-mRNA Regulatory Network

We used the miRTarBase database to evaluate the correlation between the DE-miRNAs and the corresponding target DE-mRNAs. It is well known that target genes are inversely regulated by miRNAs. The higher the correlation coefficient, the greater the mutual regulation effect (17). Based on this theory, we calculated the Pearson correlation coefficient between 19 DE-miRNAs and 340 DE-mRNAs and chose those predicted miRNAs that correlated negatively with target genes (r < 0) and the interaction pairs that exist in miRTarBase for further research (Figure 3A). Subsequently, we selected six pairs of miRNA-mRNA interactions (miR-4295 - ZNF224, miR-4271 - SMYD5, miR-1281 - KPNA2, miR-1281 - GPN2, miR-8071 - TMEM184B, and miR-575 - POLD3) with P < 0.05. The correlation of the above six significant interaction pairs is shown in Figure 3B. Finally, a significantly strong negative correlation was found between the miR-4295-ZNF224 interaction pair with a cut-off value of r < -0.7 and P < 0.01.

A, Diagram of miRNA-mRNA interaction; B, Correlation of six significant miRNA-mRNA interaction pairs: miR-4295/ZNF224, miR-4271/SMYD5, miR-1281/KPNA2, miR-1281/GPN2, miR-8071/TMEM184B, and miR-575/POLD3.
A, Diagram of miRNA-mRNA interaction; B, Correlation of six significant miRNA-mRNA interaction pairs: miR-4295/ZNF224, miR-4271/SMYD5, miR-1281/KPNA2, miR-1281/GPN2, miR-8071/TMEM184B, and miR-575/POLD3.

5. Discussion

The miRNA-mRNA regulatory network is associated with the initiation and progression of various human diseases (17-19). However, the miRNA-mRNA regulatory network that regulates intrahepatic cccDNA of CHB-GZ patients is still unknown. In this work, we aimed to identify several miRNAs, their target genes, and miRNA-mRNA regulatory network that may be involved in regulating the expression of cccDNA in patients with CHB-GZ through a series of bioinformatic analyses. These analyses provide us with a theoretical basis for locating new therapeutic targets.

A total of 19 DE-miRNAs, including 12 down-regulated miRNAs and seven up-regulated miRNAs, and 340 DE-mRNAs, including 160 down-regulated mRNAs and 180 up-regulated mRNAs, were screened. A series of bioinformatic analyses established a possible miRNA-mRNA regulatory network. A total of six pairs of miRNA-mRNA (miR-4295 - ZNF224, miR-4271 - SMYD5, miR-1281 - KPNA2, miR-1281 - GPN2, miR-8071 - TMEM184B, and miR-575 - POLD3) were obtained, among which miR-4295 - ZNF224 showed a relatively high correlation.

ZNF224 is a Kruppel-associated box zinc finger protein (KRAB-ZFP) family member, the largest family of transcriptional regulators in higher vertebrates. It is characterized by an N-terminal KRAB domain and a C-terminal array of DNA-binding zinc fingers (20). Despite their numerical abundance, the functions of KRAB-ZFP have not been clearly defined for a long time. However, cumulative data indicate that some processes involve imprinting, cell differentiation, metabolic control, and gender dimorphism (20). The ZNF224 gene is differentially expressed in human tissues and cells. According to the distinctive cellular environments, it exhibits different essential functions in transcription or post-transcriptional gene regulation (21). ZNF224 was initially identified as a transcriptional repressor involved in the specific silencing of genes through chromatin-modifying activities, the recruitment of the corepressor KRAB-associated protein 1 (KAP1), and the protein arginine methyltransferase 5 (PRMT5) in the promoter of its target genes (21). In carcinogenesis, ZNF224 was assigned a dual role of tumor suppressor and oncogene, thus causing apoptosis or stimulating cell proliferation (22). Its function as an oncogene was discovered in bladder cancer, breast cancer, and chronic lymphocytic leukemia (22). However, ZNF224 in chronic myeloid leukemia is a tumor suppressor that can promote apoptotic gene expression and suppress anti-apoptotic gene expression to improve drug resistance (23).

MiR-4295 is a newly discovered miRNA located on chromosome 10q25.2. It is a unique miRNA with many biological functions (24). MiR-4295 is widely expressed in vivo, participating in the biological processes of multiple cancers. At the same time, it promotes cancer by affecting seven cellular activities, including cell proliferation, cell migration, cell cycle, cell invasion, epithelial-mesenchymal transition (EMT), and apoptosis (24). For example, miR-4295 promotes the growth of bladder cancer cells by targeting the B-cell translocation gene 1 (BTG1) (25). It contributes to cell invasion and proliferation of pancreatic ductal adenocarcinoma through glypican-5 (GPC5) (26). It can also promote cell migration, proliferation, and invasion of osteosarcoma by targeting interferon regulatory factor 1 (IRF1) (27). Although miR-4295 is dysregulated in various cancers, it has also been found to inhibit cancer (24). The expression of miR-4295 is related to prognosis and can be affected by numerous factors connected to the therapeutic effects of various drugs (24). The target genes of miR-4295 include TP63, USP28, RUNX3, GPC5, CDKN1A, BTG1, NPTX1, LRIG1, and PTPN14 (25, 28, 29). In addition, downstream genes of miR-4295 also include SOS2, EREG, VEGF, VEGF-A, FLT1, FLT2, and FOXF1 (30, 31). Three signaling pathways related to miR-4295, including N-myc/miR-4295/RUNX3, miR-4295/GPC5/Wnt/β-catenin, and miR-4295/EGFR/PI3K/Akt, have been identified (24). The Wnt/β-catenin pathways are well-known in hepatitis B infection and hepatocellular carcinoma.

Although no research has revealed the relationship between miR-4295 and ZNF224 to the HBV infection process, the enrichment function of the GO and KEGG pathway is primarily related to the biological process of the hepatitis B virus. More enriched GO terms were related to regulating the viral life cycle, the viral process, and the viral genome replication. Enrichment analysis of the KEGG pathway suggested that these predicted targets were related to hepatitis B and influenza A. According to an analysis, annual influenza vaccination reduces the risk of hospitalization and death in CHB patients, suggesting possible impacts of influenza-targeted therapy on HBV infection (32).

5.1. Limitations

While we conducted a thorough analysis and identified several potential miRNA-mRNA axes that may regulate intrahepatic cccDNA in CHB-GZ patients, we acknowledge some limitations to this study. Firstly, the sample size was relatively small, consisting of only 13 liver tissue samples. Additionally, the miRNA-mRNA regulatory axis we identified has yet to be validated by in vivo or in vitro experiments. Therefore, larger clinical samples and corresponding experiments are required to confirm our results in the future. Long-term follow-up is also necessary to observe disease progression.

5.2. Conclusions

In summary, we have identified a potential miRNA-mRNA regulatory axis related to intrahepatic cccDNA in patients with CHB-GZ. Our study suggests that decreased expression of miR-4295 in CHB-GZ patients with a low level of cccDNA is associated with increased expression of ZNF224, which is involved in the biological process of HBV replication and the regulation of its life cycle. These findings could provide new potential targets for the treatment of CHB infection. However, larger clinical samples and further experiments are needed to verify our findings.

References

  • 1.

    Lanini S, Ustianowski A, Pisapia R, Zumla A, Ippolito G. Viral hepatitis: Etiology, epidemiology, transmission, diagnostics, treatment, and prevention. Infect Dis Clin North Am. 2019;33(4):1045-62. [PubMed ID: 31668190]. https://doi.org/10.1016/j.idc.2019.08.004.

  • 2.

    Lythgoe KA, Lumley SF, Pellis L, McKeating JA, Matthews PC. Estimating hepatitis B virus cccDNA persistence in chronic infection. Virus Evol. 2021;7(1):veaa063. [PubMed ID: 33732502]. [PubMed Central ID: PMC7947180]. https://doi.org/10.1093/ve/veaa063.

  • 3.

    Cheng ST, Hu JL, Ren JH, Yu HB, Zhong S, Wai Wong VK, et al. Dicoumarol, an NQO1 inhibitor, blocks cccDNA transcription by promoting degradation of HBx. J Hepatol. 2021;74(3):522-34. [PubMed ID: 32987030]. https://doi.org/10.1016/j.jhep.2020.09.019.

  • 4.

    Xia Y, Guo H. Hepatitis B virus cccDNA: Formation, regulation and therapeutic potential. Antiviral Res. 2020;180:104824. [PubMed ID: 32450266]. [PubMed Central ID: PMC7387223]. https://doi.org/10.1016/j.antiviral.2020.104824.

  • 5.

    O'Brien J, Hayder H, Zayed Y, Peng C. Overview of MicroRNA biogenesis, mechanisms of actions, and circulation. Front Endocrinol (Lausanne). 2018;9:402. [PubMed ID: 30123182]. [PubMed Central ID: PMC6085463]. https://doi.org/10.3389/fendo.2018.00402.

  • 6.

    Dandri M. Epigenetic modulation in chronic hepatitis B virus infection. Semin Immunopathol. 2020;42(2):173-85. [PubMed ID: 32185454]. [PubMed Central ID: PMC7174266]. https://doi.org/10.1007/s00281-020-00780-6.

  • 7.

    Xing T, Zhu J, Xian J, Li A, Wang X, Wang W, et al. miRNA-548ah promotes the replication and expression of hepatitis B virus by targeting histone deacetylase 4. Life Sci. 2019;219:199-208. [PubMed ID: 30615846]. https://doi.org/10.1016/j.lfs.2018.12.057.

  • 8.

    Zhang X, Liu H, Xie Z, Deng W, Wu C, Qin B, et al. Epigenetically regulated miR-449a enhances hepatitis B virus replication by targeting cAMP-responsive element binding protein 5 and modulating hepatocytes phenotype. Sci Rep. 2016;6:25389. [PubMed ID: 27138288]. [PubMed Central ID: PMC4853741]. https://doi.org/10.1038/srep25389.

  • 9.

    Yao K, Liu J, Wang J, Yan X, Xia J, Yang Y, et al. Distribution and clinical characteristics of patients with chronic hepatitis B virus infection in the grey zone. J Viral Hepat. 2021;28(7):1025-33. [PubMed ID: 33797145]. https://doi.org/10.1111/jvh.13511.

  • 10.

    Reza AMMT, Cho S, Choi Y, Hong K, Kim J. Microarray profiling of miRNA and mRNA expression in Rag2 knockout and wild-type mouse spleens. Scientific Data. 2018;5(1). https://doi.org/10.1038/sdata.2017.199.

  • 11.

    Lucifora J, Protzer U. Attacking hepatitis B virus cccDNA--The holy grail to hepatitis B cure. J Hepatol. 2016;64(1 Suppl):S41-8. [PubMed ID: 27084036]. https://doi.org/10.1016/j.jhep.2016.02.009.

  • 12.

    Liang LB, Zhu X, Yan LB, Du LY, Liu C, Liao J, et al. Quantitative intrahepatic HBV cccDNA correlates with histological liver inflammation in chronic hepatitis B virus infection. Int J Infect Dis. 2016;52:77-82. [PubMed ID: 27686728]. https://doi.org/10.1016/j.ijid.2016.09.022.

  • 13.

    Gene Ontology C. Gene ontology consortium: Going forward. Nucleic Acids Res. 2015;43(Database issue):D1049-56. [PubMed ID: 25428369]. [PubMed Central ID: PMC4383973]. https://doi.org/10.1093/nar/gku1179.

  • 14.

    Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27-30. [PubMed ID: 10592173]. [PubMed Central ID: PMC102409]. https://doi.org/10.1093/nar/28.1.27.

  • 15.

    Huang HY, Lin YC, Cui S, Huang Y, Tang Y, Xu J, et al. miRTarBase update 2022: an informative resource for experimentally validated miRNA-target interactions. Nucleic Acids Res. 2022;50(D1):D222-30. [PubMed ID: 34850920]. [PubMed Central ID: PMC8728135]. https://doi.org/10.1093/nar/gkab1079.

  • 16.

    Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498-504. [PubMed ID: 14597658]. [PubMed Central ID: PMC403769]. https://doi.org/10.1101/gr.1239303.

  • 17.

    Lou W, Liu J, Ding B, Chen D, Xu L, Ding J, et al. Identification of potential miRNA-mRNA regulatory network contributing to pathogenesis of HBV-related HCC. J Transl Med. 2019;17(1):7. [PubMed ID: 30602391]. [PubMed Central ID: PMC6317219]. https://doi.org/10.1186/s12967-018-1761-7.

  • 18.

    Huang DP, Zeng YH, Yuan WQ, Huang XF, Chen SQ, Wang MY, et al. Bioinformatics analyses of potential miRNA-mRNA regulatory axis in HBV-related hepatocellular carcinoma. Int J Med Sci. 2021;18(2):335-46. [PubMed ID: 33390802]. [PubMed Central ID: PMC7757140]. https://doi.org/10.7150/ijms.50126.

  • 19.

    Ma S, Xie Z, Zhang L, Yang Y, Jiang H, Ouyang X, et al. Identification of a Potential miRNA–mRNA Regulatory Network Associated With the Prognosis of HBV-ACLF. Front Mol Biosci. 2021;8. https://doi.org/10.3389/fmolb.2021.657631.

  • 20.

    Ecco G, Imbeault M, Trono D. KRAB zinc finger proteins. Development. 2017;144(15):2719-29. [PubMed ID: 28765213]. [PubMed Central ID: PMC7117961]. https://doi.org/10.1242/dev.132605.

  • 21.

    Cesaro E, Sodaro G, Montano G, Grosso M, Lupo A, Costanzo P. The complex role of the ZNF224 transcription factor in cancer. Adv Protein Chem Struct Biol. 2017;107:191-222. [PubMed ID: 28215224]. https://doi.org/10.1016/bs.apcsb.2016.11.003.

  • 22.

    Sobocinska J, Molenda S, Machnik M, Oleksiewicz U. KRAB-ZFP transcriptional regulators acting as oncogenes and tumor suppressors: An overview. Int J Mol Sci. 2021;22(4). [PubMed ID: 33672287]. [PubMed Central ID: PMC7926519]. https://doi.org/10.3390/ijms22042212.

  • 23.

    Sodaro G, Cesaro E, Montano G, Blasio G, Fiorentino F, Romano S, et al. Role of ZNF224 in c-Myc repression and imatinib responsiveness in chronic myeloid leukemia. Oncotarget. 2018;9(3):3417-31. [PubMed ID: 29423056]. [PubMed Central ID: PMC5790473]. https://doi.org/10.18632/oncotarget.23283.

  • 24.

    Li Y, Zhang Y, Zou Y, Duan S. The paradoxical roles of miR-4295 in human cancer: Implications in pathogenesis and personalized medicine. Genes Dis. 2022;9(3):638-47. [PubMed ID: 35782974]. [PubMed Central ID: PMC9243315]. https://doi.org/10.1016/j.gendis.2020.09.007.

  • 25.

    Nan YH, Wang J, Wang Y, Sun PH, Han YP, Fan L, et al. MiR-4295 promotes cell growth in bladder cancer by targeting BTG1. Am J Transl Res. 2016;8(11):4892-901. [PubMed ID: 27904689]. [PubMed Central ID: PMC5126331].

  • 26.

    Yuan Q, Zhang Y, Li J, Cao G, Yang W. High expression of microRNA-4295 contributes to cell proliferation and invasion of pancreatic ductal adenocarcinoma by the down-regulation of Glypican-5. Biochem Biophys Res Commun. 2018;497(1):73-9. [PubMed ID: 29407175]. https://doi.org/10.1016/j.bbrc.2018.02.023.

  • 27.

    Cheng JP, Huang B, Duan JH, Yi KJ, Zhuang ZL. miR-4295 promotes cell proliferation, migration and invasion of osteosarcoma through targeting interferon regulatory factor 1. Oncol Lett. 2020;20(5):260. [PubMed ID: 32989394]. [PubMed Central ID: PMC7517570]. https://doi.org/10.3892/ol.2020.12123.

  • 28.

    Jin H, Xu J, Guo X, Huang H, Li J, Peng M, et al. XIAP RING domain mediates miR-4295 expression and subsequently inhibiting p63alpha protein translation and promoting transformation of bladder epithelial cells. Oncotarget. 2016;7(35):56540-57. [PubMed ID: 27447744]. [PubMed Central ID: PMC5302933]. https://doi.org/10.18632/oncotarget.10645.

  • 29.

    Yan R, Li K, Yuan DW, Wang HN, Zhang Y, Dang CX, et al. Downregulation of microRNA-4295 enhances cisplatin-induced gastric cancer cell apoptosis through the EGFR/PI3K/Akt signaling pathway by targeting LRIG1. Int J Oncol. 2018;53(6):2566-78. [PubMed ID: 30320337]. [PubMed Central ID: PMC6203147]. https://doi.org/10.3892/ijo.2018.4595.

  • 30.

    Zhang L, Wang H, Wang C. Persistence of smoking induced non-small cell lung carcinogenesis by decreasing ERBB pathway-related microRNA expression. Thorac Cancer. 2019;10(4):890-7. [PubMed ID: 30868748]. [PubMed Central ID: PMC6449224]. https://doi.org/10.1111/1759-7714.13020.

  • 31.

    Liang G, Duan C, He J, Ma W, Dai X. PTPN14, a target gene of miR-4295, restricts the growth and invasion of osteosarcoma cells through inactivation of YAP1 signalling. Clin Exp Pharmacol Physiol. 2020;47(7):1301-10. [PubMed ID: 32141101]. https://doi.org/10.1111/1440-1681.13296.

  • 32.

    Su FH, Huang YL, Sung FC, Su CT, Hsu WH, Chang SN, et al. Annual influenza vaccination reduces total hospitalization in patients with chronic hepatitis B virus infection: A population-based analysis. Vaccine. 2016;34(1):120-7. [PubMed ID: 26614589]. https://doi.org/10.1016/j.vaccine.2015.10.129.