Study of the G Protein Nucleolar 2 Value in Liver Hepatocellular Carcinoma Treatment and Prognosis
Total Page:16
File Type:pdf, Size:1020Kb
Hindawi BioMed Research International Volume 2021, Article ID 4873678, 12 pages https://doi.org/10.1155/2021/4873678 Research Article Study of the G Protein Nucleolar 2 Value in Liver Hepatocellular Carcinoma Treatment and Prognosis Yiwei Dong,1 Qianqian Cai,2 Lisheng Fu,1 Haojie Liu,1 Mingzhe Ma,3 and Xingzhong Wu 1 1Department of Biochemistry and Molecular Biology, School of Basic Medical Sciences, Fudan University, Key Lab of Glycoconjugate Research, Ministry of Public Health, Shanghai 200032, China 2Shanghai Key Laboratory of Molecular Imaging, Shanghai University of Medicine and Health Sciences, Shanghai 201318, China 3Department of Gastric Surgery, Fudan University Shanghai Cancer Center, Shanghai 200025, China Correspondence should be addressed to Xingzhong Wu; [email protected] Received 10 May 2021; Accepted 29 June 2021; Published 20 July 2021 Academic Editor: Tao Huang Copyright © 2021 Yiwei Dong et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. LIHC (liver hepatocellular carcinoma) mostly occurs in patients with chronic liver disease. It is primarily induced by a vicious cycle of liver injury, inflammation, and regeneration that usually last for decades. The G protein nucleolar 2 (GNL2), as a protein- encoding gene, is also known as NGP1, Nog2, Nug2, Ngp-1, and HUMAUANTIG. Few reports are shown towards the specific biological function of GNL2. Meanwhile, it is still unclear whether it is related to the pathogenesis of carcinoma up to date. Here, our study attempts to validate the role and function of GNL2 in LIHC via multiple databases and functional assays. After analysis of gene expression profile from The Cancer Genome Atlas (TCGA) database, GNL2 was largely heightened in LIHC, and its overexpression displayed a close relationship with different stages and poor prognosis of carcinoma. After enrichment analysis, the data revealed that the genes coexpressed with GNL2 probably participated in ribosome biosynthesis which was essential for unrestricted growth of carcinoma. Cell functional assays presented that GNL2 knockdown by siRNA in LIHC cells MHCC97-H and SMCC-7721 greatly reduced cell proliferation, migration, and invasion ability. All in all, these findings capitulated that GNL2 could be a promising treatment target and prognosis biomarker for LIHC. 1. Background been the most commonly used biomarker in LIHC man- agement [4, 13, 14]. However, LIHC is classified as a Liver cancer is composed of primary carcinoma and meta- complicated disease accompanied by several risk factors static carcinoma [1, 2]. As far as the most common caused pathogenic mechanisms. Thus, to characterize inducers of carcinoma-associated mortality in the world, LIHC only by a single biomarker is not feasible. To explore primary liver carcinoma ranks second, thus becoming a the pathogenesis of LIHC and uncover the candidate bio- great public health challenge [3–5]. Primary liver carci- markers become urgently needed. noma consists of liver hepatocellular carcinoma (LIHC), As a protein-encoding gene, the G protein nucleolar 2 intrahepatic cholangiocarcinoma (iCCA), and other rarely (GNL2) is also named by nucleostemin (NS), which is essen- seen neoplasms [6, 7]. Mostly, LIHC occurs in patients tial for stem cell growth and development [15, 16]. GNL2 is with chronic liver disease, and it is generally caused by a highly expressed in tissues, such as testis and bone marrow vicious cycle of liver damage, inflammation, and regenera- [17, 18]. Presently, there are no researches on the particular tion that usually span several decades [8–11]. Like other biological function of GNL2, and little is known about its malignant neoplasms, the widely generated biological alter- relationship with carcinoma pathogenesis. ation in LIHC pathogenesis comprises activated oncogene Here, our study is aimed at investigating whether there is and inactivated neoplasm inhibitor genes [12]. Since its an association of GNL2 expression with the prognosis in discovery over 60 years ago, alpha-fetoprotein (AFP) has LIHC patients. In addition, our study also wants to validate 2 BioMed Research International 10 8 6 value P 10 4 –Log 2 0 –4 –2 0 2 4 Log2 fold change Figure 1: Volcano plot illustrated the DEGs existing in liver carcinoma samples and normal ones based on TCGA dataset. Genes with upregulation were presented as the red dot, while those with downregulation were shown as the blue dot. its inhibitory effect on LIHC cell proliferation, migration, and annotation enrichment, we utilized the online tool Enrichr invasion by in vitro experiments. (https://maayanlab.cloud/Enrichr/). Enrichment annota- tions were presented as Kyoto Encyclopedia of Genes and 2. Material and Methods Genomes (KEGG) pathway and Gene Ontology (GO) bio- logical process. Data were plotted and visualized using an fi ff 2.1. Identi cation of the Di erentially Expressed Genes online platform http://www.bioinformatics.com.cn. (DEGs) in TCGA Database. The Cancer Genome Atlas (TCGA) database, as the National Cancer Institute (NCI) 2.5. Cell Culture and Transfection. The human LIHC cell and the National Human Genome Institute (NHGRI) sup- lines MHCC97-H and SMCC-7721 were obtained from the ported the project, was employed to obtain mRNA expres- China Center for Type Culture Collection in 2020 and sion profile of liver carcinoma and normal samples. It could authenticated by short tandem repeat analysis. All cells were offer all-inclusive maps of the main alterations of genomic routinely cultured in DMEM medium (Thermo Fisher Scien- in multiple sorts of carcinomas. ∣log fold change ðFCÞ ∣ >1 2 tific, Inc., USA) supplemented with 10% FBS (Thermo Fisher and P value <0.05 were considered to be parameters for iden- Scientific, Inc., USA) in an incubator with 5% CO at 37°C. tifying the DEGs. The protein-protein network (PPI) of the 2 SiRNA targeting GNL2 (si-GNL2) and a negative control selected DEGs was constructed by the Search Tool for the siRNA (si-NC) were acquired from GeneCopoeia, Inc. Retrieval of Interacting Genes (STRING) database. (USA). When MHCC97-H and SMCC-7721 cells confluence 2.2. Gene Expression Profiling Interactive Analysis (GEPIA) reached up to 80% on the following day, they were trans- Database. GEPIA is an interacting online server for the anal- fected with indicated siRNAs utilizing Lipofectamine 2000 ffi ysis of RNA sequencing data of tumorous and healthy sam- (Invitrogen, USA) as the manual described. Reduced e - ples from the Genotype-Tissue Expression (GTEx) and ciency of GNL2 was determined by qRT-PCR at 48 h post- TCGA datasets. The URL of GEPIA is http://gepia.cancer- transfection. The siRNAs used were as follows: si-GNL2-1, ′ ′ pku.cn. In this study, we used GEPIA to determine GNL2 5 -CACGTGTGATTAAGCAGTCATCATT-3 ; si-GNL2-2, expression in different stages of liver carcinoma and then 5′-CCATACAAAGTTGTCATGAAGCAAA-3′; si-NC, 5′ perform patient survival analysis and identify similar genes. -UUCUCCGAACGUGUCACGUTT-3′. 2.3. Kaplan-Meier Plotter Database. The Kaplan-Meier plot- 2.6. Quantitative Reverse-Transcription PCR (qRT-PCR) ter (https://kmplot.com/analysis/), functioning as a web tool, Assay. Overall RNA was harvested from indicated cells utiliz- is used to validate survival biomarkers in the light of meta- ing TRIzol reagent (Invitrogen, USA). qRT-PCR was carried analysis. The information of gene expression profile and out on a Bio-Rad Single Color Real-Time PCR system (Bio- various survival data comprising overall survival (OS), Rad Laboratories, Inc., USA) with SYBR-Green Real-Time relapse-free survival (RFS), and disease-specific survival PCR Master Mix (Toyobo Life Science, Japan). β-Actin was (DSS) come from the well-known public databases, including used as an internal control. GEO (only Affymetrix microarrays), EGA, and TCGA. Kaplan-Meier plotter was taken to generate the survival curves of GNL2, and log rank P values, as well as hazard ratio 2.7. Cell Proliferation Assay. 2,000 of si-NC- and si-GNL2- (HR) with 95% confidence intervals, were calculated. transfected cells were plated in per well of 96-well plates, followed by maintaining 0, 24, 48, and 72 h at 37°C. 10 μlof 2.4. Biological Process and Pathway Enrichment Analysis. The CCK-8 reagent was added into each well at the indicated top 1000 genes coexpressed with GNL2 (sorted by Pearson time, followed by another 2 h incubation. The wavelength correlation coefficient) in LICH-tumor and LICH-normal at 450 nm was read in a Varioskan Flash spectrophotometer dataset were identified by GEPIA. To assess the gene- (Thermo Fisher Scientific, Inc., USA). BioMed Research International 3 GPM6A ICM5 TBXA2R CLEC4M GRPR DKKL1 CHRM3 MARCO SGCD ADRA1A STAB2 WDFY2 CDC37L1 TSPO PVALB DBH SP2 ADRA2B CRHBP TUBDNASE1L3 NAAA ECM1 UCN2 TRIB1 NTF3 MAT1A CIDEB GNL2 MAP2K1 INS-IGF2 SPN ALB ZFP36 CXCL12 MMP14 HIP1 TBX21 PLG IL7R ETS2 EGR1 LGALS1ADAMTS13 POSTN LCAT CCL23 DUSP3 APOF CSRNP1 FCN2 TSLP DSG1 RCAN1 FCN3 CFP LIFR C1RL MED1 C1R CSTA MASP1 C8A IL5RA GHR COLEC10 C6 GYS2 CPED1 CYP4A11 ALDH8A1 TTC36 FMO1 CYP2C8 ALDH2 ALDH6A1 ADH4 ACADSB UGT3A1 CNDP1 LRAT GLYATL1 ACADS ERLIN1 ETFDH NAT2 BCO2 MRO DCDC2 MAP7D2 GPR20 MRPL43 VIPR1 KIF13B GDPD3 RPL27A BMPER PTH1R KLC3 RPS25 GDF2 PNPT1 CLEC4G SENP7 PFIA1 CLEC1B PIAS1 CASKIN1 CPEB3 Figure 2: PPI network of the DEGs. Nodes displayed proteins, and edges exhibited the interactions between proteins. 2.8. Transwell Assay. For cell migration assay, Matrigel needed with 10% FBS, followed by maintaining 24 h at 37°C. Then, not to be precoated in the Transwell chamber (Corning, NY). migrated cells underside were rinsed by PBS, fixed by metha- The Transwell upper chamber was added 50,000 si-NC- or si- nol, and stained by DAPI. Positive cells were photographed GNL2-transfected cells suspended in serum-free DMEM and counted under a microscope. Five random fields were medium, and the lower chamber was added DMEM medium selected and counted the average.