Skip to main content

Iron-sulphur cluster biogenesis factor LYRM4 is a novel prognostic biomarker associated with immune infiltrates in hepatocellular carcinoma

Abstract

Background

LYRM4 is necessary to maintain the stability and activity of the human cysteine desulfurase complex NFS1-LYRM4-ACP. The existing experimental results indicate that cancer cells rely on the high expression of NFS1. However, the role of LYRM4 in liver hepatocellular carcinoma (LIHC) remains unclear.

Methods

In this study, we combined bioinformatics analysis and clinical specimens to evaluate the mRNA, protein expression, and gene regulatory network of LYRM4 in LIHC. Furthermore, we detected the activity of several classical iron-sulphur proteins in LIHC cell lines through UV-vis spectrophotometry.

Results

The mRNA and protein levels of LYRM4 were upregulated in LIHC. Subsequent analysis revealed that the LYRM4 mRNA expression was related to various clinical stratifications, prognosis, and survival of LIHC patients. In addition, the mRNA expression of LYRM4 was significantly associated with ALT, tumour thrombus, and encapsulation of HBV-related LIHC patients. IHC results confirmed that LYRM4 was highly expressed in LIHC tissues and showed that the expression of LYRM4 protein in LIHC was significantly correlated with age and serum low-density lipoprotein (LDL) and triglyceride (TG) content. In particular, the mRNA expression of key iron- sulphur proteins POLD1 and PRIM2 was significantly overexpressed and correlated with poor prognosis in LIHC patients. Compared with hepatocytes, the activities of mitochondrial complex I and aconitate hydratase (ACO2) in LIHC cell lines were significantly increased. These results indicated that the iron-sulphur cluster (ISC) biosynthesis was significantly elevated in LIHC, leading to ISC-dependent metabolic reprogramming. Changes in the activity of ISC-dependent proteins may also occur in paracancerous tissues. Further analysis of the biological interaction and gene regulation networks of LYRM4 suggested that these genes were mainly involved in the citric acid cycle and oxidative phosphorylation. Finally, LYRM4 expression in LIHC was significantly positively correlated with the infiltrating levels of six immune cell types, and both factors were strongly associated with prognosis.

Conclusion

LYRM4 could be a novel prognostic biomarker and molecular target for LIHC therapy. In particular, the potential regulatory networks of LYRM4 overexpression in LIHC provide a scientific basis for future research on the role of the ISC assembly mechanism and LYRM4-mediated sulphur transfer routes in carcinogenesis.

Background

LIHC is the fifth most common malignancy and the second most frequent cause of cancer-related mortality worldwide [1, 2]. Recent studies from several countries have reported a significant age-specific increase in LIHC development among people over 75 years old [2, 3]. Based on annual projections, the World Health Organization predicts that more than 1 million patients with LIHC will die by 2030 [4]. LIHC accounts for approximately 85%–90% of all primary liver malignancies. The rising incidence of LIHC is correlated with the increasing prevalence of hepatitis B virus (HBV), hepatitis C virus (HCV) infections and alcohol abuse, particularly in East Asia [5, 6]. Although therapies with sorafenib and regorafenib provide a modest survival benefit, overall anti-tumor responses are limited [7], with a 5-year survival rate of only 18% for LIHC [8]. Therefore, it is necessary to screen novel biomarkers for early diagnosis and identify new therapeutic targets for LIHC by screening the changes in the gene function network related to the formation and progression of LIHC.

Interestingly, traditional methods focus more on the LIHC itself to stratify and screen prognostic gene sets. In contrast, Gong et al. [9] developed a new method to identify LIHC subtypes and prognostic gene sets based on the activity changes of immunologic and hallmark gene sets in tumour and non-tumour tissues, since the activity changes of these gene sets in adjacent non-tumour tissues may substantially impact prognosis by affecting the proliferation of the remaining liver cells and colonisation of circulating tumour cells after treatment measures such as hepatectomy, so as to improve patient outcomes with personalised medicine and accurate assessment of prognosis.

Iron-sulphur (Fe-S) clusters are ubiquitous in all life forms. They have become integral parts of many physiological functions, including energy conversion, protein translation, DNA replication, and repair [10]. The human mitochondrial cysteine desulfurase complex consists of a catalytic subunit (NFS1), an LYR protein (LYRM4, which is also an alternative name to ISD11), and an acyl carrier protein (ACP), in which the C-terminus of LYRM4 is essential for forming a stable and active complex with NFS1 [11,12,13]. The cysteine desulfurase complex further interacts with mitochondrial scaffold protein ISCU, frataxin (FXN) (is a kinetic activator of persulfide transfer) and ferredoxin FDX2 (allows persulfide reduction into sulphide), forming the “core ISC complex” for de novo Fe-S cluster synthesis [14]. It also acts as a sulphur donor for other processes, such as tRNA synthesis [15,16,17] and molybdenum cofactor biosynthesis [17, 18]. In addition, the LYR superfamily forms lock-and-key interactions with acyl-ACP, which is associated with fatty acid biosynthesis to coordinate the expression, Fe/S cofactor maturation, and activity of the respiratory complexes [19,20,21].

A recent report showed that NFS1 is overexpressed in most well-differentiated lung adenocarcinomas because incipient lung tumours survive in a high oxygen environment. NFS1 activity is particularly important for maintaining iron-sulphur co-factors present in at least 48 human cell-essential enzymes. Thus, overexpression of NFS1 in adenocarcinomas leads to resistance to high oxygen tension and protects cells from undergoing ferroptosis in response to oxidative damage [22]. Consistent with this, our analysis of RNA-Seq data from 159 HBV-related LIHC patients showed that NFS1 mRNA was overexpressed in tumour tissues compared to paired adjacent liver tissues (p ˂ 0.05) [5]. Shi et al. [23] used RNA-interference techniques to demonstrate that human LYRM4 is essential for both ISC biogenesis and maintenance of normal cellular iron homeostasis in HeLa cells. However, there have been no reports on the differential expression of NFS1 and its auxiliary protein LYRM4 in tumour tissues.

Liver cells contain approximately two thousand mitochondria, which are the energy and metabolic factories of cells [24]. Therefore, it is necessary to evaluate the expression of LYRM4 in pan-cancer and to investigate whether abnormal expression of LYRM4 is closely related to the tumorigenesis and development of LIHC. In the present study, we investigated the expression profile, prognostic significance, and gene regulation network of LYRM4 in LIHC based on different open databases by utilizing multi-dimensional analysis methods. Moreover, we used IHC to detect the protein expression level of LYRM4 in LIHC and paired adjacent normal tissues and analysed the correlation between their characteristics. In addition, co-expression analysis was performed between LYRM4 and the genes involved in mitochondrial ISC biosynthesis in tumour and non-tumour tissues. The changes in ISC-dependent key iron-sulphur proteins and their relationship with the prognosis of LIHC patients were also analysed. Finally, we explored the correlation between LYRM4 expression and immune infiltration. These findings will aid in revealing the potential pathogenesis mechanism of LYRM4 in LIHC and show that LYRM4 may be a useful biomarker and target in the treatment of LIHC.

Methods

Reagents and antibodies

Horseradish peroxidase (HRP)-conjugated anti-rabbit and anti-mouse immunoglobulin G, and monoclonal antibody against β-Actin were obtained from Beyotime (Shanghai, China). BCA protein assay kit, SuperSignal West Femto maximum sensitivity substrate, and pierce protease and phosphatase inhibitor mini tablets, EDTA-free were obtained from Thermo Scientific (Waltham, MA). Rabbit polyclonal to LYRM4 was purchased from Abcam (Cambridge, UK). Streptavidin-biotin-peroxidase complex (SABC) and 3,3′-Diaminobenzidine tetrahydrochloride (DAB) were purchased from ZSGB-BIO (Beijing, China).

Tissue microarrays (TMA) and IHC staining

We collected paraffin tissue specimens from 92 patients diagnosed with LIHC from the First Affiliated Hospital of Wenzhou Medical University from May 2018 to October 2019, and detailed clinical data and laboratory examination results of all patients were tracked and recorded. The inclusion criteria were as follows: (1) the pathologic diagnosis was LIHC; (2) no chemotherapy or radiation therapy was administered before surgery; (3) detailed hospitalisation information; and (4) patients who volunteered to participate in this study and had lived in Wenzhou for more than 5 years. The exclusion criteria were as follows: (1) basic metabolic and infectious diseases; (2) history of other tumours; (3) radiotherapy and/or chemotherapy; and (4) inability to obtain complete clinical data.

All collected formalin-fixed, paraffin-embedded LIHC tissues were prepared to construct TMA, as described in a previous study [25]. IHC staining was performed following previously described protocols, with minor modifications [26]. The TMA slides were incubated at 4 °C overnight with primary antibodies against LYRM4 (1:100) or PBS as a negative control, followed by incubation with biotin-labelled goat anti-rabbit IgG at 37 °C for 30 min, and then incubated with SABC at room temperature for 30 min. Signals were visualised with DAB, counterstained with hematoxylin, dehydrated in ethanol, cleared in xylene, and then fixed.

Evaluation of IHC staining

IHC images were reviewed and scored independently by two pathologists who had no prior knowledge of the clinicopathological features of the TMA specimens. The relative protein expression of LYRM4 was analysed by calculating the integrated optical density per stained area (IOD/area, mean OD value, MOD) using Image-Pro Plus 6.0 (Media Cybernetics, Rockville, MD, USA) according to a previously described protocol [27]. Specific scoring criteria were as follows: H-score < 7 was categorised as the low expression group, and H-score ≥ 7 was categorised as the high expression group.

Cell culture

The human LIHC cell lines HCCLM3, MHCC97H, HepG2 and hepatocytes L-02 cells were obtained from the Cell Bank of the Chinese Academy of Sciences (Shanghai, China). All cells were cultured in high-glucose DMEM (GIBCO, Waltham, MA, USA) containing 10% FBS (fetal bovine serum) (GIBCO, USA) and antibiotics (100 U/ml penicillin and 100 μg/ml streptomycin) (GIBCO, USA), incubated at 37 °C in a humidified incubator with 5% CO2. Cell lines were routinely tested and confirmed to be mycoplasma free during this study.

Western blot analysis

Proteins from whole-cell was extracted with 1% Triton X-100 lysis buffer supplemented with protease and phosphatase inhibitors. Protein concentrations of the extracts were determined by the BCA assay kit. 40 μg total protein from each sample was separated by ExpressPlus PAGE Gel (10 × 8, 4–20%, 15 wells) (Genscript, Nanjing, China) and transferred onto 0.22 µm PVDF membrane (Bio-Rad, Hercules, CA, USA) using a semidry transfer system (Bio-Rad, USA). Blots were blocked at room temperature for 15 min in QuickBlock Blocking Buffer for Western Blot (Beyotime, Shanghai, China) on a shaker, and then incubated with primary antibodies specific to LYRM4 (1:2000) and β-actin (1:5000) overnight at 4 ºC, respectively. The membrane was washed in TBST for 3 × 10 min and then incubated with HRP-conjugated anti-rabbit (1:5000) and anti-mouse (1:20000) IgG secondary antibody at room temperature for 1.5 h while gently agitating. Signals were detected with SuperSignal West Femto maximum sensitivity substrate according to the manufacturer`s protocol.

Isolation and purification of mitochondria

A mitochondria extraction kit (Solarbio, Beijing, China) was used to isolate and purify the mitochondria as described previously [28, 29].

Mitochondrial complex I activity assay

Complex I activity was measured as previously described [30, 31].

Mitochondrial aconitate hydratase activity assay

The prepared mitochondrial protein (5 µL; 2 µg/µL) was added to the assay buffer containing Tris (pH 8.0; 100 mM), cis-aconitic acid (1 mM), MgCl2 (5 mM), and isocitrate dehydrogenase (4 U/ml). The mixture was incubated at 37 °C for 2 min. After reading the baseline at 340 nm, NADP+ was added, and the increase in absorbance was recorded at 340 nm for 80 s. The enzymatic activity of each mitochondrial enzyme was corrected using citrate synthase (CS) activity.

Oncomine database analysis

The Oncomine 4.5 database (https://www.oncomine.org/resource/login.html) is the largest oncogene chip database and integrated data mining platform, including 715 gene expression datasets and clinical data from 86,733 cancer tissues and normal tissues [32]. The retrieval condition was set as: Analysis Type/Cancer vs. Normal Analysis, Cancer Type/Liver Cancer, Dataset Filters/Data Type/mRNA or DNA, Sample Filters/Sample Type/Clinical Specime. Then, the mRNA expression and DNA copy number of LYRM4 in LIHC were analyzed by the Oncomine database. The results are displayed with the fold change over 2. Student's t-test was used to analyze differences in the expression of LYRM4 between normal controls and cancer specimens. A cutoff of p < 0.05 was considered statistically significant.

GEPIA2 analysis

The GEPIA2 database (http://gepia2.cancer-pku.cn/#index), developed by Peking University in China, is based on TCGA and GTEx database which includes RNA sequencing and expression data from 33 malignant tumour, 8587 normal tissue, and 9736 tumour samples [33]. The expression differences of LYRM4 between tumours and normal tissues in distinct types of cancer were determined within the GEPIA2 database. We also used GEPIA2 to analyze the correlation between LYRM4 expression and pathologic staging and survival of LIHC patients, including Overall Survival and Disease Free Survival, in LIHC. The p, p(HR) or Pr(> F) values from a log-rank test were included in the plot.

UALCAN database analysis

UALCAN (http://ualcan.path.uab.edu/) is a comprehensive and user-friendly web resource to perform in-depth analyses of TCGA gene expression data [34]. In this study, we used UALCAN to analyze the expression differences of LYRM4 and its promoter methylation in LIHC, as well as in various tumor subgroups based on age, individual cancer stages, tumor grade and other clinical stratified criteria. Student’s t-test was used to generate p-values; after Bonferroni correction for multiple measures, p was still < 0.05, which was statistically significant.

STRING analysis

The String database (https://string-db.org/) is a web-based platform for analyzing the interactions between genes or proteins, including direct physical interactions between proteins and indirect functional correlations between proteins [35]. It integrates experimental data, abstracts from PubMed and data from other databases. In addition, it provides bioinformatics-based predictions. We used STRING to search the network map of proteins associated with LYRM4. The database key setting parameter was set to medium confidence 0.4 and the maximal number of interactors was 50.

GeneMANIA analysis

GeneMANIA (http://genemania.org/) is a commonly used interactive website for constructing PPI networks and predicting the function of genes. The website features several bioinformatic methods, including gene co-expression, co-localization, pathway, physical interactions and predicted [36]. We used GeneMANIA to construct the gene networks and predict the function of LYRM4 and enriched genes in LIHC identified by GSEA (miRNA 495 and transcription factor GGAANCGGAANY_UNKNOWN).

LinkedOmics analysis

The LinkedOmics database (http://www.linkedomics.org/login.php) is a comprehensive and user-friendly web interface with integrated multi-omics data and clinical data from 11,158 patients (32 cancer types) from TCGA project [37]. In this study, LinkedOmics was utilized to study genes differentially expressed in correlation with LYRM4 in LIHC. Pearson’s correlation coefficient was applied in statistical analysis of the result produced by LinkedOmics. GSEA was utilized to perform various enrichment analyses, including GO, KEGG pathways, kinase-target, miRNA-target and transcription factor-target. The result of LinkFinder was ranked based on criteria of FDR < 0.05 and 500 simulations.

TargetScan analysis

TargetScan (http://www.targetscan.org/vert_72/) is a flexible, user-friendly web interface that predicts biological targets of miRNAs by searching for the presence of conserved 8mer, 7mer, and 6mer sites that match the seed region of each miRNA [38]. We used TargetScan to analyse miRNAs differentially expressed in connection with LYRM4 and predict the biological target genes of these miRNAs.

FunRich analysis

FunRich (http://www.funrich.org/) is an independent programming tool for functional enrichment and interaction network analysis of genes and proteins [39]. In this study, we used Funrich tool for the gene sets that were enriched in the target network of LYRM4 and miRNA enrichment analysis, including Biological process, Cellular component, Molecular function and Biological pathway.

TIMER2.0 analysis

TIMER2.0 (http://timer.cistrome.org/) is a comprehensive resource for the systematic investigation of immune infiltrates across diverse cancer types [40]. We first analysed the correlations between the expression of LYRM4 and the abundance of the six immune cell types in LIHC using Spearman tests (tumour purity adjusted). The survival module was used to draw Kaplan-Meier plots for immune infiltrates and LYRM4 to determine survival differences. Finally, the expression of LYRM4 in various tumours was analysed using the TIMER database, and the results were analysed statistically using the Wilcoxon rank sum test. Statistical significance was set at p < 0.05.

Statistical analysis

All statistical analyses were performed using SPSS 21.0 (IBM, Armonk, NY, USA). Data conforming to a normal distribution are presented as mean ± SD. Comparisons between two groups were performed using an independent Student’s t-test. Measurement data that did not conform to the normal distribution were represented by M (P25, P75), and the comparison between the two groups was performed using the rank sum test in the IHC results. The χ2 test was performed to evaluate the relationship between clinicopathological features and LYRM4 expression in the IHC results. The Wilcoxon signed-rank test was used to analyse the significance of protein and mRNA expression in adjacent normal and cancer tissues of LIHC and HBV-Related LIHC. Statistical significance was set at p < 0.05.

Results

LYRM4 expression and its correlation with clinicopathological parameters of LIHC

To determine differences in LYRM4 mRNA expression in tumour and paired normal tissues, the LYRM4 gene expression profile across diverse cancer types was determined using the Oncomine, UALCAN, and GEPIA2 databases. LYRM4 expression was higher in lymphoid neoplasms, such as diffuse large B-cell lymphoma (DLBC), leukaemia, sarcoma, thymoma (THYM), LIHC, colorectal cancer, gastric cancer, lung cancer, melanoma, and lymphoma, compared to normal samples. The mRNA level of LYRM4 was significantly decreased in kidney chromophobe (KICH), brain, and CNS cancer tissues compared with matched normal tissues (Additional file 1: Figure S1a–b). The difference in LYRM4 expression in LIHC was further evaluated using the Oncomine and GEPIA2 databases. The results showed that the mRNA expression and DNA copy number variation (CNV) of LYRM4 were significantly higher in LIHC tissues than in normal tissues (p ˂ 0.05) (Fig. 1a-c, e, f and Additional file 1: Figure S1c). Consistently, LYRM4 mRNA levels were also significantly upregulated in HBV-related LIHC patients older than 54 years (Fig. 1g) [5].

Fig. 1
figure1

Expression levels of LYRM4 mRNA and DNA copy number in LIHC tissue samples. ad Box plots showing LYRM4 mRNA expression levels in Roessler Liver 2, Chen Liver, Wurmbach Liver, and Mas Liver datasets, respectively (Oncomine database). ef Box plots showing LYRM4 copy number in The Cancer Genome Atlas (TCGA) Liver and Guichard Liver datasets, respectively. The fold change, associated p values, and overexpression gene rank were based on Oncomine analysis. g The mRNA expression of LYRM4 in HBV-related LIHC and paired non-tumour liver tissues was investigated using RNA-seq (n = 79). These figure results were obtained from the studies of Gao et al. [5], and all the LIHC patients are over 54 years old. The top and bottom black dots in the figure represent the maximum and minimum values of samples, respectively

Further subgroup analysis of multiple clinicopathological features of LYRM4 in LIHC reliably indicated high levels of LYRM4 transcription. The transcription level of LYRM4 was significantly higher in LIHC patients than that in healthy people in subgroup analyses based on age (normal vs. age [21–40 y], normal vs. age [41–60 y], normal vs. age [61–80 y], and normal vs. age [81-100 y], p ˂ 0.001). Moreover, in the sex, weight, race, individual cancer stage, tumour grade, nodal metastasis, and TP53 mutation status subgroup analyses, the expression of LYRM4 mRNA was also significantly higher in LIHC patients (Fig. 2).

Fig. 2
figure2

Relative expression of LYRM4 mRNA in subgroups of patients with LIHC, classified according to sample types (a), sex (b), age (c), weight (d), race (e), individual cancer stages (f), tumour grade (g), nodal metastasis status (h), and TP53 mutation status (i) (UALCAN database). Data are shown as mean ± SE. *, p < 0.05; **, p < 0.01; ***, p < 0.001. The asterisk indicates a significant difference between the two sets of data

In addition, we analysed the association between LYRM4 mRNA expression and various clinicopathological factors of HBV-related LIHC using paired tumours and adjacent liver tissues from 157 patients (Additional file 2: Table S1) [5]. As shown in Table 1, Pearson chi-square test analyses revealed that the mRNA expression of LYRM4 was significantly associated with tumour thrombus (p ˂ 0.05). Multivariate analysis showed that the mRNA expression of LYRM4 was related to ALT (glutamic-pyruvic transaminase), tumour thrombus, and encapsulation of HBV-related LIHC patients (p ˂ 0.05).

Table 1 Association between LYRM4 mRNA expression and various clinicopathological factors of HBV-related LIHC patients

Promoter methylation levels of LYRM4 in LIHC

To evaluate the promoter methylation levels of LYRM4 in LIHC, we used the UALCAN database to analyse the correlation between LYRM4 expression and promoter DNA methylation. No statistical difference in the promoter methylation level of LYRM4 was observed in LIHC tissues compared to normal tissues. However, further subgroup analysis of the LYRM4 promoter methylation profile revealed great significance according to race, tumour grade, and nodal metastasis status in LIHC patients (Fig. 3).

Fig. 3
figure3

LYRM4 promoter methylation profile of LIHC patient subgroups based on sample types (a), race (b), tumour grade (c), and nodal metastasis status (d) (UALCAN database). Data are shown as mean ± SE. *, p < 0.05; **, p < 0.01; ***, p < 0.001. The asterisk indicates a significant difference between the two sets of data

Prognostic potential of LYRM4 in LIHC

To further assess whether high LYRM4 mRNA expression levels are associated with clinical and pathologic parameters of LIHC patients, we used the GEPIA2 database to analyse the prognostic value of LYRM4 in LIHC. According to the Kaplan-Meier survival analysis, LIHC patients with high LYRM4 expression had lower overall survival rates than those with low LYRM4 expression (Fig. 4a, p < 0.01). However, LYRM4 mRNA expression seemed to have little impact on the disease/progression-free survival of LIHC patients (Fig. 4b).

Fig. 4
figure4

LYRM4 is a poor prognostic factor for LIHC progression and patient survival (GEPIA2 database). a Kaplan-Meier overall survival (OS) curve of LIHC patients (n = 364) based on LYRM4 mRNA expression level. b Kaplan-Meier disease free survival analysis was performed on LIHC patients according to the expression level of LYRM4

Elevated protein expression of LYRM4 contributes to reprogramming of iron-sulphur proteins-related metabolic pathways in LIHC cells

To further explore the clinical significance of LYRM4 in LIHC, LYRM4 protein expression was analysed using the Human Protein Atlas database. As shown in Fig. 5a, LYRM4 was not detected in the normal liver and showed weak to medium staining in LIHC. To increase the credibility of the results, IHC staining was performed in TMA containing 92 archived paraffin-embedded LIHC specimens (Additional file 3: Figure S2–4). Consistently, protein analysis on TMA slides confirmed that LYRM4 was dramatically overexpressed in LIHC tissues compared to adjacent normal tissues (Fig. 5b-c). The association between LYRM4 protein expression and clinicopathological features of LIHC was analysed using the chi-square test or rank sum test (Additional file 4: Table S2). As shown in Table 2, these results suggested that the expression of LYRM4 was significantly associated with age (p = 0.0410), serum LDL (p = 0.0340), and TG content (p = 0.0130). No significant relationship was found between LYRM4 protein expression and variables such as alcohol consumption. However, drinking habits dramatically increased the expression level (mean) of LYRM4 in both LIHC and paracancerous tissues (Additional file 5: Figure S14). Additionally, we found increased LYRM4 in both HBV-related (MHCC97H) and negative (HepG2) LIHC cell lines (Fig. 5d and Additional file 3: Figure S5) compared to that in hepatocytes L-02 cells.

Fig. 5
figure5

Protein expression of LYRM4 and Fe-S cluster biosynthesis were significantly increased in LIHC and affected LIHC patient prognosis. a The protein expression levels of LYRM4 were elevated in LIHC (The Human Protein Atlas database). b Representative immunohistochemistry images of LYRM4 in tissue arrays containing human LIHC specimens and matched adjacent normal tissues. Regions in squares are magnified 4× in bottom panels. c Summary statistics of H-score based on only intact and paired specimens (n = 92). d Western blot analysis of LYRM4 in normal hepatocytes and three LIHC cell lines. eg Box plots showing NDUFS8 (e), POLD1 (f), and PRIM2 (g) mRNA expression levels in LIHC (GEPIA2). hi Kaplan-Meier diagram of the relationship between POLD1 (h) and PRIM2 (i) gene expression and survival in LIHC patients (GEPIA2). j Protein expression of mitochondrial aconitase ACO2 in LIHC increased (The Human Protein Atlas database). kl Activities of mitochondrial complex I and ACO2 were detected in the human LIHC cell lines and hepatocytes, and the histogram represented the relative activity of control (hepatocytes). Data are presented as means ± SD (n = 3). *, p < 0.05; **, p < 0.01

Table 2 Association between LYRM4 protein expression and various clinicopathological factors in LIHC patients

It is well known that ISC assembly proteins are encoded by the gene cluster iscSUA-hscBA-fdx in Escherichia coli and are highly conserved among aerobic organisms from bacteria to humans [41]. To determine whether similar transcriptional mechanisms operate in LIHC cells, a co-expression analysis was performed between LYRM4 and the 19 proteins involved in mitochondrial Fe/S protein assembly and export at transcript levels in LIHC using the GEPIA2 database. As shown in Table 3, in the liver, the expression of the LYRM4 gene was significantly positively correlated with the expression of 16 other genes, excluding NFS1, NUBPL, and MFRN2, among which the co-expression correlation coefficient (R) of four genes (NFU1, ISCU, ISCA1, and ACP1) were all greater than 0.7 (p = 0). However, in LIHC normal tissues and LIHC tumour tissues, the co-expression correlation between LYRM4 and these genes was sequentially reduced, or even negatively correlated.

Table 3 The co-expression analysis was performed between LYRM4 and the proteins involved in mitochondrial Fe/S protein assembly and export at transcript levels in LIHC using GEPIA2 web

Cysteine desulfurase activity is crucial for maintaining the iron-sulphur co-factors present in many cell-essential proteins [13, 22]. Furthermore, we found that the mRNA expression of key iron-sulphur proteins NADH dehydrogenase [ubiquinone] iron-sulphur protein 8 (NDUFS8), DNA polymerase delta catalytic subunit (POLD1), and DNA primase large subunit (PRIM2) were significantly elevated (Fig. 5e–g). Moreover, the high mRNA expression of POLD1 and PRIM2 was correlated with poor prognosis in LIHC patients (Fig. 5h–i). In addition, we explored the expression patterns of the classical [4Fe-4S] protein ACO2 in the mitochondria of LIHC patients using the Human Protein Atlas. As shown in Fig. 5j, low protein expression of ACO2 was observed in normal hepatocytes, while medium to high protein expression of ACO2 was observed in LIHC tissues. Compared with hepatocytes, the activities of mitochondrial complex I and ACO2 in LIHC cell lines were also significantly elevated (Fig. 5k–l and Additional file 3: Figure S6–13). These results indicated that Fe-S cluster biosynthesis was significantly enhanced in LIHC patients and affected the prognosis of LIHC patients.

Biological interaction network of LYRM4

In order to better understand the molecular function of LYRM4 and the signaling pathways in which LYRM4 is involved, 50 interactors in correlation with LYRM4 were determined using STRING database (Fig. 6a). In addition, protein-protein interaction network analysis via GeneMANIA revealed 20 genes were enriched in LYRM4 network. The enriched gene set plays a role in mitochondrial electron transport, NADH to ubiquinone, mitochondrial ATP synthesis coupled electron transport, mitochondrial respiratory chain complex I and NADH dehydrogenase activity (Fig. 6b). Furthermore, GO enrichment and KEGG pathway analyses for biological interaction network of LYRM4 were performed via FunRich tool. Abundantly enriched GO terms indicated that these genes mainly encode proteins located in the mitochondrion and are involving in energy pathways, metabolism, and oxidoreductase activity (Fig. 6c–e). Similarly, KEGG pathway analysis showed enrichment in respiratory electron transport, ATP synthesis by chemiosmotic coupling, and citric acid cycle signaling pathways (Fig. 6f).

Fig. 6
figure6

Biological interaction network of LYRM4 and enrichment analysis of gene sets enriched in the target network of LYRM4 (STRING and GeneMANIA databases). a Protein network associated with LYRM4. b Protein-protein interaction (PPI) network and functional analysis of gene sets enriched in the target network of LYRM4. Colours of the network edges indicate the bioinformatics methods applied: Physical Interactions, Co-expression, Predicted, Co-localization, Pathway, Genetic Interactions, Shared protein domains. Colours of the network nodes indicate the biological functions of enriched gene sets. c Biological processes. d Cellular components. e Molecular functions. f KEGG pathway analysis

Enrichment analysis of LYRM4 functional networks in LIHC

GO and KEGG pathway analyses of co-expression genes correlated with LYRM4 in LIHC

To analyze the co-expression genes associated with LYRM4 in LIHC, we used LinkedOmics to analyze mRNA sequencing data from 371 LIHC patients in the TCGA. A Pearson test was conducted to analyse connections among LYRM4 and genes differentially expressed in LIHC (Fig. 7a). The top 50 significant genes positively or negatively correlated with LYRM4 were shown in the heatmaps (Fig. 7b-c). These results indicate that the expression of LYRM4 is closely correlated with a wide range of genes at transcriptome level in LIHC. Three individual genes (RPS10, TOMM6 and WDR46) showed the strongest positive correlation with the expression of LYRM4 (Additional file 6: Figure S15a–c). This result reflects changes in the component of the 40S ribosomal subunit and mitochondrial protein import, as well as the rRNA modification process in the nucleus and cytosol. Gene Set Enrichment Analysis (GSEA) results demonstrated that differentially expressed genes related to LYRM4 were mainly located in the ribosomes, cytosolic parts and mitochondrial protein complexes, where they primarily participate in protein localization to endoplasmic reticulum, translational initiation and ribonucleoprotein complex biogenesis. They act as structural constituents of ribosome and rRNA binding (Fig. 8a–c and Additional files 7, 8, 9: Table S3–5). KEGG pathway analysis showed enrichment in the ribosome, spliceosome, oxidative phosphorylation and proteasome pathways (Fig. 8d–e and Additional file 10: Table S6).

Fig. 7
figure7

Genes differentially expressed in correlation with LYRM4 in LIHC (LinkedOmics database). a Pearson test was used to analyse correlations between LYRM4 and genes differentially expressed in LIHC. bc Heatmaps showing genes (top 50) that are positively and negatively correlated with LYRM4 in LIHC. Red indicates positively-correlated genes and blue indicates negatively-correlated genes

Fig. 8
figure8

Significantly enriched GO annotations and KEGG pathways of LYRM4 in LIHC. a Cellular components. b Biological processes. c Molecular functions. d KEGG pathway analysis. The orange column represents the Leading Edge Number, and the black represents the false discovery rate (FDR). The FDR from GSEA in the figure is 0. e KEGG pathway annotations of the oxidative phosphorylation pathway. The nodes marked in red are associated with the Leading Edge Genes. FDR < 0.05 was considered statistically significant

LYRM4 networks of kinase, miRNA or transcription factor targets in LIHC

To further screen the targets of LYRM4 in LIHC, we listed the 5 most significant LYRM4-related genes in the kinase, miRNA and transcription factor target networks in LIHC (Additional files 11, 12, 13: Table S7–9). The results indicated that the most leading edge numbers in kinase target, miRNA target and transcription factor target were Kinase EGFR, AGCACTT, MIR-93, MIR-302A, MIR-302B, MIR-302C, MIR-302D, MIR-372, MIR-373, MIR-520E, MIR-520A, MIR-526B, MIR-520B, MIR-520C, MIR-520D and SCGGAAGY-V$ELK1-02, respectively (Table 4). Since FDR value of the Kinase-target networks of LYRM4 in LIHC was greater than 0.05, we did not carry out subsequent protein-protein interaction (PPI) network construction. The PPI networks constructed by GeneMANIA revealed that gene sets related to miRNA 495 were mainly involved in the regulation of mRNA binding, protein K48-linked deubiquitination, polysome, and maintenance of protein location (Additional file 14: Figure S16). The gene set rich in transcription factor PAX4 is mainly responsible for organ development, mesenchymal cell differentiation and sequence-specific DNA binding (Additional file 15: Figure S17).

Table 4 Kinase, miRNA and transcription factor-target networks of LYRM4 in LIHC (LinkedOmics)

miRNAs related to LYRM4

According to cumulative weighted context++ score, the top 5 among 1570 miRNA families associated with LYRM4 were miR-24-3p, miR-133a-3p.1, miR-125-5p, miR-193a-5p and miR-302c-3p.2/520-3p. The conserved sites in human LYRM4 mRNA 3ˊUTR targeted by miRNA families were broadly conserved among vertebrates showed in Fig. 9a. To predict the functional enrichment information of the identified 1570 miRNAs, GO and KEGG pathway analyses were performed using Funrich database. Biological processes were significantly enriched in regulation of nucleobase, nucleoside, nucleotide and nucleic acid metabolism, signal transduction and cell communication. The main enriched Cell components were the nucleus, cytoplasm, and golgi aparatus. In the Molecular function category, transcription factor activity, GTPase activity, and transcription regulator activity were primarily enriched. In terms of KEGG pathway analysis, TRAIL signaling pathway, proteoglycan syndecan-mediated signaling events, glypican pathway and sphingosine 1-phosphate (S1P) pathway were enriched (Fig. 9b–e).

Fig. 9
figure9

Prediction and enrichment analysis of the identified miRNAs related to LYRM4 (Targetscan and Funrich). a Conserved sites in LYRM4 for microRNA families broadly conserved among vertebrates. be GO functional and KEGG pathway enrichment analysis of 1,570 identified microRNAs

Immune infiltrates in correlation with LYRM4 in LIHC

LIHC is a common malignant tumour with a poor prognosis. The immune microenvironment is poorly characterised [7]. The expression of LYRM4 increased significantly with increasing malignancy of LIHC (Fig. 2f–g). Therefore, we explored the relationship between LYRM4 expression and immune infiltration in LIHC. The results showed that the correlation between LYRM4 expression in LIHC and the infiltration levels of different types of immune cells (B cells, CD4+ T cells, CD8+ T cells, neutrophils, macrophages, and dendritic cells) was statistically significant (p < 0.05, Fig. 10a). According to the cumulative survival analysis (Fig. 10b), patients with high LYRM4 expression and high abundance of the T cell CD4+ Th2-XCELL subtype had significantly lower survival rates than those with high LYRM4 expression and low abundance of the T cell CD4+ Th2-XCELL subtype (HR = 1.66, p = 0.0376). The high LYRM4 expression and high abundance of macrophage-TIMER subtype patients exhibited significantly lower survival rates in LIHC patients with high LYRM4 expression and low abundance of macrophage-TIMER subtype (HR = 1.57, p = 0.0427) (Fig. 10c). In addition, the combination of other subtypes and abundance of immune macrophage infiltration with different expression levels of LYRM4 notably affected the survival of LIHC patients (p < 0.05, Additional file 16: Figure S18d). Finally, the expression of LYRM4 in pan-cancer was analysed using the TIMER2.0 database. As shown in Fig. 10d, LYRM4 mRNA was significantly overexpressed in colon adenocarcinoma (COAD), head and neck cancer (HNSC), kidney renal papillary cell carcinoma (KIRP), LIHC, lung squamous cell carcinoma (LUSC), and prostate adenocarcinoma (PRAD; p < 0.001).

Fig. 10
figure10

Associations between mRNA expression of LYRM4 and immune infiltration in LIHC (TIMER&TIMER2.0). a Correlation between LYRM4 expression and the abundance of immune infiltrates (B cells, CD4+ T cells, CD8+ T cells, Neutrophils, Macrophages, and Dendritic cells). Correlation r or purity adjusted r value and p value are presented in the figure. The blue line is the fitted curve. bc The effect of immune cells infiltration levels in correlation with LYRM4 expression on the prognosis of LIHC. d LYRM4 mRNA expression in different types of cancers. *, p < 0.05; **, p < 0.01; ***, p < 0.001

Discussion

Differential expression of the proteins required for ISC biosynthesis has been reported to be implications for cellular metabolism and adaptation to cellular stress [22, 42]. Alvarez et al. [43] reported that hypoxia rescues Frataxin loss by restoring ISC biogenesis. In particular, Alvarez et al. [22] found that cancer cells depend on high levels of the ISC biosynthetic enzyme NFS1. Therefore, it is necessary to investigate whether the NFS1-interacting protein LYRM4 is differentially expressed in tumours and its role in carcinogenesis. In this study, we performed bioinformatics analysis using different open databases to acquire detailed information about the potential functions of LYRM4 in LIHC and its regulatory network. We provided the first evidence that the protein expression of LYRM4 is significantly upregulated in LIHC tissues and cell lines (Fig. 5). Thus, this study laid the foundation for future experimental studies on LIHC.

AFP is the most widely used serological indicator for the diagnosis of LIHC worldwide. However, approximately 30%–40% of LIHC patients are AFP-negative, and about 20%–50% of patients with chronic hepatitis or cirrhosis are AFP-positive [44]. A combination of the triple biomarkers AFP, AFP-L3, and DCP was adopted for the early diagnosis of LIHC patients according to the Milan criteria and Child-Pugh class A. The results showed that the 5-year cumulative overall survival (OS) and disease-free survival (DFS) rates (with statistical significance) in postoperative triple-negative and triple-positive patients were 85.3%, 44.2%, 61.7%, and 35.7%, respectively [45]. In addition, the OS of LIHC patients was not significantly different between the high and low AFP expression groups (GEPIA2, data not shown). Therefore, it is necessary to screen for new biomarkers for early LIHC diagnosis and poor patient prognosis.

In this study, we found that LYRM4 mRNA levels and CNVs were significantly higher in LIHC tissues than in normal tissues (Figs. 1, 2, and Additional file 1: Figure S1). However, overexpression of LYRM4 CNVs was found in only 10% of LIHC samples, based on CNVs (Fig. 1e–f). In particular, the expression level of LYRM4 was positively correlated with the degree of tumour progression and was undifferentiated in LIHC (Fig. 2f–g). We speculate that the overexpression of LYRM4 in LIHC may be closely related to the significantly reduced promoter DNA methylation levels (Fig. 3b–d).

Interestingly, in the TP53 mutation status subgroup analysis, the transcription level of LYRM4 was significantly higher in LIHC patients than in healthy people and LIHC patients with no p53 mutation (Fig. 2i). Therefore, these results suggest that TP53 mutant may regulate LYRM4 expression at the transcriptional level. Funauchi et al. [46] reported that ISCU expression was markedly decreased in the majority of LIHC tissues, and its reduced expression was significantly correlated with p53 mutation. However, in the GEPIA2 database, the mRNA expression level of ISCU in LIHC tissues was obviously higher than that in adjacent normal tissues, although the difference was not statistically significant (Additional file 17: Figure S19a). Additionally, analysis of RNA-Seq data from 159 HBV-related LIHC patients showed that ISCU mRNA was overexpressed in tumour tissues compared to paired adjacent liver tissues (Additional file 17: Figure S19b), and the Spearman’s correlation coefficient of the ISCU mRNA-protein pair was 0.472 (p = 8.07E-10) [5]. These results suggest that molecular genetic background, population, whether HBV infection and whether alcoholic fatty liver disease progression of LIHC (Additional file 5: Figure S14) may be the key factors that lead to the abnormal expression of LYRM4 and ISCU in LIHC.

To further investigate the relationship between LYRM4 protein expression and progression of LIHC. We used IHC analysis to indicate that patients with advanced age and higher serum levels of LDL and TG had higher expression of LYRM4 protein (Table 2). We speculated that high expression of LYRM4 may lead to cellular metabolic reprogramming and lipid metabolism disorders by modulating Fe-S cluster assembly, regulating the tricarboxylic acid cycle and β-oxidation [47, 48], which would increase the burden on the liver and thus adversely affect the prognosis of LIHC patients [49, 50]. Unfortunately, we did not find any other significant differences between the clinicopathological and LYRM4 expression. Consistently, Gao et al. [5] reported that Spearman’s correlation coefficient for the LYRM4 mRNA-protein pair was 0.12832 (p = 0.11238) in HBV-related LIHC. These results suggest that the high expression of LYRM4 protein in LIHC may be regulated by post-transcriptional or post-translational modification regulation. It is worthy of further multicentre and large sample size verification.

In conclusion, our data suggest that LYRM4 could serve as a potential diagnostic biomarker for LIHC. Additionally, it may be a useful prognostic marker for LIHC patients, as high expression of LYRM4 could lead to poor prognosis of LIHC (Fig. 4).

Since ISCs are essential protein cofactors with diverse physiological functions, variations in their abundance may result in changes in various downstream signalling pathways [13]. Clinical mutations of NFS1 (R72Q) [51] and LYRM4 (R68L) [52] could lead to deficits in ISC, thus reducing electron transport chain complex (ETC) activity and mitochondrial respiration. In this study, a biological interaction network of LYRM4-associated genes was constructed, and LYRM4 was involved in the regulation of mitochondrial function and oxidoreductase activity (Fig. 6). These findings are consistent with those of previous reports [10, 52]. It should be noted that the co-expression association between LYRM4 and 19 genes involved in mitochondrial Fe/S protein assembly and export became gradually disordered with the onset and progression of LIHC (Table 3). Most importantly, we found that the activities of mitochondrial complex I and ACO2 in LIHC were significantly increased. Moreover, the high expression of POLD1 and PRIM2 was related to poor clinical outcomes in LIHC (Fig. 5). Therefore, these results indicate that elevated expression of LYRM4 was significantly associated with ISC-specific metabolic reprogramming in LIHC. The regulation of ISC-specific gene set expression was gradually disturbed in paracancerous tissues, which provides a new perspective for preventing the recurrence of LIHC. Meanwhile, we can also select 1-3 iron-sulphur proteins (PRIM2, NFU1, or POLD1), which are highly expressed (Human Protein Atlas) and exhibit a large co-expression correlation coefficient with LYRM4 (Table 3), as a combined biomarkers for the early diagnosis of LIHC with LYRM4, thus effectively improving the specificity and reliability of the diagnostic results.

As shown in Table 4, we found that LYRM4 in LIHC is linked to a network of miRNAs, including miR-495, miR-18A, miR-18B, and miR-153. Yang et al. [53] found that the upregulation of miR-664, miR-485-3p, and miR-495 leads to a reduced MAT1A expression in LIHC and enhanced tumorigenesis. Xia et al. [54] reported that miR-153 plays a critical role in suppressing epithelial-to-mesenchymal transition (EMT) and LIHC progression by directly downregulating Snail expression. Therefore, we hypothesised that LYRM4 and these miRNAs may have synergistic effects in LIHC. In addition, LYRM4 may downregulate the expression of MAT1A through miR-495, thereby significantly reducing the methylation levels of several proto-oncogene promoter regions. Further studies are needed to confirm this hypothesis.

The EMT signalling pathway plays a critical role in embryological development, cancer progression, and metastasis [55]. A previous study showed that the transcription factor Elk-1 is a downstream target of Erk1/2 MAP kinases, regulating the migration of many types of cancer via its downstream targets including c-fos, EGR1, and PKCα [56]. Zhang et al. [57] demonstrated that upregulation of transcription factor PAX4 directly downregulated miR-144/451 expression and promoted migration and invasion in epithelial cancer cells. It is should also be noted that the mRNA expression of LYRM4 was significantly associated with tumour thrombus and encapsulation of HBV-related LIHC patients (Table 1). Therefore, we believe that LYRM4 may promote metastasis by regulating the expression of the transcription factors Elk-1 and PAX4 in LIHC (Table 4), which remains to be verified experimentally. Taken together, it is worth exploring whether LYRM4 overexpression mediates malignant transformation and metastasis of LIHC through regulatory metabolic reprogramming and mitochondrial retrograde signalling pathways.

According to reports, the HIF-1α/miR-210 signalling axis regulates mitochondrial free radicals in response to hypoxia in cancer cells by targeting downregulation of ISCU expression, which correlates with poor prognosis in breast cancer and HNSCC [58] and good prognosis in renal cancer [59]. Based on the TargetScan database, miR-193a-5p is the most likely target in miRNA families to regulate LYRM4 expression (Fig. 9a). Li et al. [60] found that miR-193a-5p was under-expressed in LIHC and inhibited the development of LIHC by targeting SPOCK1. Further functional enrichment of 1570 miRNAs identified as correlated with LYRM4 revealed terms closely related to the proliferation and apoptosis of LIHC (Fig. 9b–e). Therefore, it will be an interesting project to study the interaction between miR-193a-5p and LYRM4 in LIHC and its role in the development of LIHC.

Tumour-infiltrating immune cells are highly correlated with prognosis and the identification of immunotherapy targets in LIHC [61]. Our results suggest that immune-infiltrated CD4+ T cells and macrophages substantially affect the prognosis of LIHC and may have a prognostic value related to LYRM4 in LIHC (Fig. 10a–c and Additional file 16: Figure S18). Further studies are needed to verify these findings.

Conclusions

In summary, we analysed multi-omics data and clinical data from different open databases. IHC was used to explore the correlation between LYRM4 expression and clinical parameters of LIHC patients, demonstrating that LYRM4 plays an important role in hepatocarcinogenesis. In particular, alterations in ISC-specific metabolism in LIHC have been identified. Thus, LYRM4 may be a novel drug target and prognostic biomarker for LIHC. The data acquisition and analysis adopted in this study is time- and cost-efficient, providing a scientific basis for subsequent functional studies. However, functional analysis and biochemical data on the enzymatic activities of ISC biogenesis, tRNA thiolation, and molybdenum cofactor biosynthesis are needed in future studies to confirm the relationship between the overexpression of LYRM4 and the occurrence and development of LIHC in single ethnic groups with or without HBV infection or alcohol abuse.

Availability of data and materials

The datasets generated for this study are available on request to the corresponding author.

Abbreviations

ISC:

Iron-sulphur cluster

LIHC:

Liver hepatocellular carcinoma

IHC:

Immunohistochemical

LDL:

Low-density lipoprotein

TG:

Triglyceride

HBV:

Hepatitis B virus

HCV:

Hepatitis C virus

ACP:

Acyl carrier protein

FXN:

Frataxin

HRP:

Horseradish peroxidase

SABC:

Streptavidin-biotin-peroxidase complex

DAB:

3,3′-Diaminobenzidine tetrahydrochloride

TMA:

Tissue microarrays

DLBC:

Lymphoid neoplasm diffuse large B-cell lymphoma

THYM:

Thymoma

KICH:

Kidney chromophobe

CNV:

Copy number variation

PPI:

Protein-protein interaction

S1P:

Sphingosine 1-phosphate

COAD:

Colon adenocarcinoma

HNSC:

Head and neck cancer

KIRP:

Kidney renal papillary cell carcinoma

LUSC:

Lung squamous cell carcinoma

PRAD:

Prostate adenocarcinoma

DFS:

Disease-free survival

OS:

Overall survival

ETC:

Electron transport chain complex

EMT:

Epithelial-to-mesenchymal transition

FDR:

False discovery rate

AFP:

Alpha fetoprotein

AFP-L3:

Lens culinaris agglutinin-reactive fraction of AFP

DCP:

Des-γ-carboxy prothrombin

NDUFS8:

NADH dehydrogenase [ubiquinone] iron-sulfur protein 8

POLD1:

DNA polymerase delta catalytic subunit

PRIM2:

DNA primase large subunit

MAT1A:

Methionine adenosyltransferase 1A

Elk-1:

ETS domain-containing protein

PAX4:

Paired box gene 4

References

  1. 1.

    Moris D, Chakedis J, Sun SH, Spolverato G, Tsilimigras DI, Ntanasis-Stathopoulos I, Spartalis E, Pawlik TM. Management, outcomes, and prognostic factors of ruptured hepatocellular carcinoma: a systematic review. J Surg Oncol. 2018;117(3):341–53.

    PubMed  Article  Google Scholar 

  2. 2.

    Cho E, Cho HA, Jun CH, Kim HJ, Cho SB, Choi SK. A review of hepatocellular carcinoma in elderly patients focused on management and outcomes. Vivo. 2019;33(5):1411–20.

    CAS  Article  Google Scholar 

  3. 3.

    El-Serag HB, Rudolph KL. Hepatocellular carcinoma: epidemiology and molecular carcinogenesis. Gastroenterology. 2007;132(7):2557–76.

    CAS  PubMed  Article  Google Scholar 

  4. 4.

    Villanueva A. Hepatocellular carcinoma. N Engl J Med. 2019;380(15):1450–62.

    CAS  PubMed  Article  Google Scholar 

  5. 5.

    Gao Q, Zhu H, Dong L, Shi W, Chen R, Song Z, Huang C, Li J, Dong X, Zhou Y, et al. Integrated proteogenomic characterization of HBV-related hepatocellular carcinoma. Cell. 2019;179(2):561–77.

    CAS  PubMed  Article  Google Scholar 

  6. 6.

    Jiang J, Srivastava S, Seim G, Pavlova NN, King B, Zou L, Zhang C, Zhong M, Feng H, Kapur R, et al. Promoter demethylation of the asparagine synthetase gene is required for ATF4-dependent adaptation to asparagine depletion. J Biol Chem. 2019;294(49):18674–84.

    PubMed  PubMed Central  Article  Google Scholar 

  7. 7.

    Zhang Q, He Y, Luo N, Patel SJ, Han Y, Gao R, Modak M, Carotta S, Haslinger C, Kind D, et al. Landscape and dynamics of single immune cells in hepatocellular carcinoma. Cell. 2019;179(4):829–45.

    CAS  PubMed  Article  Google Scholar 

  8. 8.

    Jemal A, Ward EM, Johnson CJ, Cronin KA, Ma J, Ryerson B, Mariotto A, Lake AJ, Wilson R, Sherman RL, et al. Annual report to the nation on the status of cancer, 1975–2014, featuring survival. J Natl Cancer Inst. 2017. https://doi.org/10.1093/jnci/djx030.

    Article  PubMed  PubMed Central  Google Scholar 

  9. 9.

    Gong J, Li R, Chen Y, Zhuo Z, Chen S, Cao J, Zhang Q, Chong Y, Hu B. HCC subtypes based on the activity changes of immunologic and hallmark gene sets in tumor and nontumor tissues. Brief Bioinform. 2021. https://doi.org/10.1093/bib/bbaa427.

    Article  PubMed  Google Scholar 

  10. 10.

    Saha PP, Srivastava S, Kumar SKP, Sinha D, D’Silva P. Mapping key residues of ISD11 critical for NFS1-ISD11 subcomplex stability: implications in the development of mitochondrial disorder, COXPD19. J Biol Chem. 2015;290(43):25876–90.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  11. 11.

    Friemel M, Marelja Z, Li K, Leimkühler S. The N-terminus of iron-sulfur cluster assembly factor ISD11 is crucial for subcellular targeting and interaction with l-cysteine desulfurase NFS1. Biochemistry. 2017;56(12):1797–808.

    CAS  PubMed  Article  Google Scholar 

  12. 12.

    Cai K, Frederick RO, Tonelli M, Markley JL. Mitochondrial cysteine desulfurase and isd11 coexpressed in Escherichia coli yield complex containing acyl carrier protein. ACS Chem Biol. 2017;12(4):918–21.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  13. 13.

    Boniecki MT, Freibert SA, Mühlenhoff U, Lill R, Cygler M. Structure and functional dynamics of the mitochondrial Fe/S cluster synthesis complex. Nat Commun. 2017;8(1):1287.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  14. 14.

    Gervason S, Larkem D, Mansour AB, Botzanowski T, Müller CS, Pecqueur L, Le Pavec G, Delaunay-Moisan A, Brun O, Agramunt J, et al. Physiologically relevant reconstitution of iron-sulfur cluster biosynthesis uncovers persulfide-processing functions of ferredoxin-2 and frataxin. Nat Commun. 2019;10(1):3566.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  15. 15.

    Nakai Y, Umeda N, Suzuki T, Nakai M, Hayashi H, Watanabe K, Kagamiyama H. Yeast Nfs1p is involved in thio-modification of both mitochondrial and cytoplasmic tRNAs. J Biol Chem. 2004;279(13):12363–8.

    CAS  PubMed  Article  Google Scholar 

  16. 16.

    Black KA, Dos Santos PC. Shared-intermediates in the biosynthesis of thio-cofactors: Mechanism and functions of cysteine desulfurases and sulfur acceptors. Biochimica et biophysica acta. 2015;1853(6):1470–80.

    CAS  PubMed  Article  Google Scholar 

  17. 17.

    Leimkühler S, Bühning M, Beilschmidt L. Shared sulfur mobilization routes for tRNA thiolation and molybdenum cofactor biosynthesis in prokaryotes and eukaryotes. Biomolecules. 2017;7(1):5.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  18. 18.

    Marelja Z, Stöcklein W, Nimtz M, Leimkühler S. A novel role for human Nfs1 in the cytoplasm: Nfs1 acts as a sulfur donor for MOCS3, a protein involved in molybdenum cofactor biosynthesis. J Biol Chem. 2008;283(37):25178–85.

    CAS  PubMed  Article  Google Scholar 

  19. 19.

    Angerer H, Schönborn S, Gorka J, Bahr U, Karas M, Wittig I, Heidler J, Hoffmann J, Morgner N, Zickermann V. Acyl modification and binding of mitochondrial ACP to multiprotein complexes. Biochim Biophys Acta Mol Cell Res. 2017;10(10):9.

    Google Scholar 

  20. 20.

    Cory SA, Van Vranken JG, Brignole EJ, Patra S, Winge DR, Drennan CL, Rutter J, Barondeau DP. Structure of human Fe-S assembly subcomplex reveals unexpected cysteine desulfurase architecture and acyl-ACP-ISD11 interactions. Proc Natl Acad Sci U S A. 2017;114(27):E5325–34.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  21. 21.

    Van Vranken JG, Nowinski SM, Clowers KJ, Jeong MY, Ouyang Y, Berg JA, Gygi JP, Gygi SP, Winge DR, Rutter J. ACP acylation is an Acetyl-CoA-Dependent modification required for electron transport chain assembly. Mol Cell. 2018;71(4):567–80.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  22. 22.

    Alvarez SW, Sviderskiy VO, Terzi EM, Papagiannakopoulos T, Moreira AL, Adams S, Sabatini DM, Birsoy K, Possemato R. NFS1 undergoes positive selection in lung tumours and protects cells from ferroptosis. Nature. 2017;551(7682):639–43.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  23. 23.

    Shi Y, Ghosh MC, Tong WH, Rouault TA. Human ISD11 is essential for both iron-sulfur cluster assembly and maintenance of normal cellular iron homeostasis. Hum Mol Genet. 2009;18(16):3014–25.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  24. 24.

    Yan J, Jiang J, He L, Chen L. Mitochondrial superoxide/hydrogen peroxide: an emerging therapeutic target for metabolic diseases. Free Radic Biol Med. 2020;152:33–42.

    CAS  PubMed  Article  Google Scholar 

  25. 25.

    Nocito A, Bubendorf L, Tinner EM, Süess K, Wagner U, Forster T, Kononen J, Fijan A, Bruderer J, Schmid U, et al. Microarrays of bladder cancer tissue are highly representative of proliferation index and histological grade. J Pathol. 2001;194(3):349–57.

    CAS  PubMed  Article  Google Scholar 

  26. 26.

    Zhou C, Sun H, Zheng C, Gao J, Fu Q, Hu N, Shao X, Zhou Y, Xiong J, Nie K, et al. Oncogenic HSP60 regulates mitochondrial oxidative phosphorylation to support Erk1/2 activation during pancreatic cancer cell growth. Cell Death Dis. 2018;9(2):161.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  27. 27.

    Liang J, Cao R, Zhang Y, Xia Y, Zheng Y, Li X, Wang L, Yang W, Lu Z. PKM2 dephosphorylation by Cdc25A promotes the Warburg effect and tumorigenesis. Nat Commun. 2016;7:12431.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  28. 28.

    Dai W, Wang X, Teng H, Li C, Wang B, Wang J. Celastrol inhibits microglial pyroptosis and attenuates inflammatory reaction in acute spinal cord injury rats. Int Immunopharmacol. 2019;66:215–23.

    CAS  PubMed  Article  Google Scholar 

  29. 29.

    Xue X, Yang JY, He Y, Wang LR, Liu P, Yu LS, Bi GH, Zhu MM, Liu YY, Xiang RW, et al. Aggregated single-walled carbon nanotubes attenuate the behavioural and neurochemical effects of methamphetamine in mice. Nat Nanotechnol. 2016;11(7):613–20.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  30. 30.

    Bao Y, Ding S, Cheng J, Liu Y, Wang B, Xu H, Shen Y, Lyu J. Carnosine inhibits the proliferation of human cervical gland carcinoma cells through inhibiting both mitochondrial bioenergetics and glycolysis pathways and retarding cell cycle progression. Integrative Cancer Therapies. 2018;17(1):80–91.

    CAS  PubMed  Article  Google Scholar 

  31. 31.

    Liu Y, Huang Y, Zhang J, Pei C, Hu J, Lyu J, Shen Y. TIMMDC1 knockdown inhibits growth and metastasis of gastric cancer cells through metabolic inhibition and AKT/GSK3β/β-catenin signaling pathway. Int J Biol Sci. 2018;14(10):1256–67.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  32. 32.

    Rhodes DR, Kalyana-Sundaram S, Mahavisno V, Varambally R, Yu J, Briggs BB, Barrette TR, Anstet MJ, Kincead-Beal C, Kulkarni P, et al. Oncomine 3.0: genes, pathways, and networks in a collection of 18,000 cancer gene expression profiles. Neoplasia. 2007;9(2):166–80.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  33. 33.

    Tang Z, Li C, Kang B, Gao G, Zhang Z. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 2017;45(W1):W98–102.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  34. 34.

    Chandrashekar DS, Bashel B, Balasubramanya SAH, Creighton CJ, Ponce-Rodriguez I, Chakravarthi B, Varambally S. UALCAN: a portal for facilitating tumor subgroup gene expression and survival analyses. Neoplasia. 2017;19(8):649–58.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  35. 35.

    von Mering C, Huynen M, Jaeggi D, Schmidt S, Bork P, Snel B. STRING: a database of predicted functional associations between proteins. Nucleic Acids Res. 2003;31(1):258–61.

    Article  CAS  Google Scholar 

  36. 36.

    Warde-Farley D, Donaldson SL, Comes O, Zuberi K, Badrawi R, Chao P, Franz M, Grouios C, Kazi F, Lopes CT, et al. The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function. Nucleic Acids Res. 2010;38:W214.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  37. 37.

    Vasaikar SV, Straub P, Wang J, Zhang B. LinkedOmics: analyzing multi-omics data within and across 32 cancer types. Nucleic Acids Res. 2018;46(D1):D956–63.

    CAS  PubMed  Article  Google Scholar 

  38. 38.

    Agarwal V, Bell GW, Nam JW, Bartel DP. Predicting effective microRNA target sites in mammalian mRNAs. Elife. 2015;4:e05005.

    PubMed Central  Article  PubMed  Google Scholar 

  39. 39.

    Pathan M, Keerthikumar S, Ang CS, Gangoda L, Quek CY, Williamson NA, Mouradov D, Sieber OM, Simpson RJ, Salim A, et al. FunRich: an open access standalone functional enrichment and interaction network analysis tool. Proteomics. 2015;15(15):2597–601.

    CAS  PubMed  Article  Google Scholar 

  40. 40.

    Li T, Fan J, Wang B, Traugh N, Chen Q, Liu JS, Li B, Liu XS. TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Res. 2017;77(21):e108–10.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  41. 41.

    Tan G, Cheng Z, Pang Y, Landry AP, Li J, Lu J, Ding H. Copper binding in IscA inhibits iron-sulphur cluster assembly in Escherichia coli. Mol Microbiol. 2014;93(4):629–44.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  42. 42.

    Chan SY, Zhang YY, Hemann C, Mahoney CE, Zweier JL, Loscalzo J. MicroRNA-210 controls mitochondrial metabolism during hypoxia by repressing the iron-sulfur cluster assembly proteins ISCU1/2. Cell Metab. 2009;10(4):273–84.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  43. 43.

    Ast T, Meisel JD, Patra S, Wang H, Grange RMH, Kim SH, Calvo SE, Orefice LL, Nagashima F, Ichinose F, et al. Hypoxia rescues frataxin loss by restoring iron sulfur cluster biogenesis. Cell. 2019;177(6):1507–21.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  44. 44.

    Kim KH, Lee SY, Hwang H, Lee JY, Ji ES, An HJ, Kim JY, Yoo JS. Direct monitoring of fucosylated glycopeptides of alpha-fetoprotein in human serum for early hepatocellular carcinoma by liquid chromatography-tandem mass spectrometry with immunoprecipitation. Proteomics Clin Appl. 2018;12(6):27.

    Article  CAS  Google Scholar 

  45. 45.

    Ryu T, Takami Y, Wada Y, Tateishi M, Matsushima H, Mikagi K, Saitsu H. Double- and triple-positive tumor markers predict early recurrence and poor survival in patients with hepatocellular carcinoma within the milan criteria and Child-Pugh class A. J Gastrointest Surg. 2017;21(6):957–66.

    PubMed  Article  Google Scholar 

  46. 46.

    Funauchi Y, Tanikawa C, Yi Lo PH, Mori J, Daigo Y, Takano A, Miyagi Y, Okawa A, Nakamura Y, Matsuda K: Regulation of iron homeostasis by the p53-ISCU pathway. Sci Rep 2015, 5(16497).

  47. 47.

    Lill R, Freibert SA. Mechanisms of mitochondrial iron-sulfur protein biogenesis. Annual Rev Biochem. 2020;89:471–99.

    CAS  Article  Google Scholar 

  48. 48.

    Tamir S, Paddock ML, Darash-Yahana-Baram M, Holt SH, Sohn YS, Agranat L, Michaeli D, Stofleth JT, Lipper CH, Morcos F, et al. Structure-function analysis of NEET proteins uncovers their role as key regulators of iron and ROS homeostasis in health and disease. Biochimica et Biophysica Acta. 2015;1853(6):1294–315.

    CAS  PubMed  Article  Google Scholar 

  49. 49.

    Li Z, Guan M, Lin Y, Cui X, Zhang Y, Zhao Z, Zhu J. Aberrant lipid metabolism in hepatocellular carcinoma revealed by liver lipidomics. Int J Mol Sci. 2017;18(12):2550.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  50. 50.

    Sun Y, Zhang N, Ding YL, Yu LJ, Cai J, Ma D, Yang W, Lu WK, Niu JL. Effect of lipid metabolism disorder on liver function in patients with malignant tumors after chemotherapy: a case-control study. Lipids Health Dis. 2019;18(1):108.

    PubMed  PubMed Central  Article  Google Scholar 

  51. 51.

    Farhan SM, Wang J, Robinson JF, Lahiry P, Siu VM, Prasad C, Kronick JB, Ramsay DA, Rupar CA, Hegele RA. Exome sequencing identifies NFS1 deficiency in a novel Fe-S cluster disease, infantile mitochondrial complex II/III deficiency. Mol Genet Genomic Med. 2014;2(1):73–80.

    CAS  PubMed  Article  Google Scholar 

  52. 52.

    Lim SC, Friemel M, Marum JE, Tucker EJ, Bruno DL, Riley LG, Christodoulou J, Kirk EP, Boneh A, DeGennaro CM, et al. Mutations in LYRM4, encoding iron-sulfur cluster biogenesis factor ISD11, cause deficiency of multiple respiratory chain complexes. Hum Mol Genet. 2013;22(22):4460–73.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  53. 53.

    Yang H, Cho ME, Li TW, Peng H, Ko KS, Mato JM, Lu SC. MicroRNAs regulate methionine adenosyltransferase 1A expression in hepatocellular carcinoma. J Clin Invest. 2013;123(1):285–98.

    CAS  PubMed  Article  Google Scholar 

  54. 54.

    Xia W, Ma X, Li X, Dong H, Yi J, Zeng W, Yang Z. miR-153 inhibits epithelial-to-mesenchymal transition in hepatocellular carcinoma by targeting snail. Oncol Rep. 2015;34(2):655–62.

    CAS  PubMed  Article  Google Scholar 

  55. 55.

    Iwatsuki M, Mimori K, Yokobori T, Ishi H, Beppu T, Nakamori S, Baba H, Mori M. Epithelial-mesenchymal transition in cancer development and its clinical significance. Cancer Sci. 2010;101(2):293–9.

    CAS  PubMed  Article  Google Scholar 

  56. 56.

    Ranjan A, Iyer SV, Ward C, Link T, Diaz FJ, Dhar A, Tawfik OW, Weinman SA, Azuma Y, Izumi T, et al. MTBP inhibits the Erk1/2-Elk-1 signaling in hepatocellular carcinoma. Oncotarget. 2018;9(30):21429–43.

    PubMed  PubMed Central  Article  Google Scholar 

  57. 57.

    Zhang J, Qin X, Sun Q, Guo H, Wu X, Xie F, Xu Q, Yan M, Liu J, Han Z, et al. Transcriptional control of PAX4-regulated miR-144/451 modulates metastasis by suppressing ADAMs expression. Oncogene. 2015;34(25):3283–95.

    CAS  PubMed  Article  Google Scholar 

  58. 58.

    Favaro E, Ramachandran A, McCormick R, Gee H, Blancher C, Crosby M, Devlin C, Blick C, Buffa F, Li JL, et al. MicroRNA-210 regulates mitochondrial free radical response to hypoxia and krebs cycle in cancer cells by targeting iron sulfur cluster protein ISCU. PLoS One. 2010;5(4):0010345.

    Article  CAS  Google Scholar 

  59. 59.

    McCormick RI, Blick C, Ragoussis J, Schoedel J, Mole DR, Young AC, Selby PJ, Banks RE, Harris AL. miR-210 is a target of hypoxia-inducible factors 1 and 2 in renal cancer, regulates ISCU and correlates with good prognosis. Br J Cancer. 2013;108(5):1133–42.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  60. 60.

    Li P, Xiao Z, Luo J, Zhang Y, Lin L. MiR-139-5p, miR-940 and miR-193a-5p inhibit the growth of hepatocellular carcinoma by targeting SPOCK1. J Cell Mol Med. 2019;23(4):2475–88.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  61. 61.

    Rohr-Udilova N, Klinglmüller F, Schulte-Hermann R, Stift J, Herac M, Salzmann M, Finotello F, Timelthaler G, Oberhuber G, Pinter M, et al. Deviations of the immune cell landscape between healthy liver and hepatocellular carcinoma. Sci Rep. 2018;8(1):6220.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

Download references

Acknowledgements

We would like to thank Dr. Lirui Liu (The University of Shenzhen) for English language editing.

Funding

This study was supported by the National Natural Science Foundation of China (No. 31071531, 81671124, 31870775 and 31772642), Hunan Provincial Natural Science Foundation (2019JJ40134), Scientific Research Project of Hunan Education Department (18A098), Hunan Provincial Graduate Research and Innovation Project (CX2018B404), Key Discipline of Zhejiang Province in Medical Technology (first class, category A), Double first-class construction project of Hunan Agricultural University (SYL201802003), and China Postdoctoral Science Foundation (2019T120705).

Author information

Affiliations

Authors

Contributions

ZW, GT, and JL: designed and supervised the study. YP, YL, JZ, JHL, and QZ: participated in data analysis and figure preparation. YP, YC and XY: collected the patient samples and performed the study. YP, XY, and HZ: performed the research, analyzed and interpreted the data. YP, and TX: wrote and revised the manuscript. JF, JL, and ZW: reviewed the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Jianxin Lyu or Zheng Wang.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no conflicts of interest.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Figure S1.

LYRM4 expression levels in different human cancers and LYRM4 overexpression level in LIHC is negatively correlated with the degree of malignancy of tumor cells.

Additional file 2: Table S1.

Clinicopathologic Data, related to Table 1.

Additional file 3: Figure S2–13.

Original data.

Additional file 4: Table S2.

Clinicopathologic information of immunohistochemical samples, related to Table 2.

Additional file 5: Figure S14.

Scatterplot of alcohol consumption correlation to LYRM4 protein expression in LIHC patients identified by immunohistochemical staining.

Additional file 6: Figure S15.

Gene expression correlation analysis for LYRM4, RPS10, TOMM6, and WDR46 (LinkedOmics database).

Additional file 7: Table S3.

Significantly enriched GO annotations (biological processes) of LYRM4 in LIHC (LinkedOmics).

Additional file 8: Table S4.

Significantly enriched GO annotations (cellular components) of LYRM4 in LIHC (LinkedOmics).

Additional file 9: Table S5.

Significantly enriched GO annotations (molecular functions) of LYRM4 in LIHC (LinkedOmics).

Additional file 10: Table S6.

Significantly enriched KEGG pathway annotations of LYRM4 in LIHC (LinkedOmics).

Additional file 11: Table S7.

Significantly enriched kinase-target networks of LYRM4 in LIHC (LinkedOmics).

Additional file 12: Table S8.

Significantly enriched miRNA-target networks of LYRM4 in LIHC (LinkedOmics).

Additional file 13: Table S9.

Significantly enriched transcription factor-target networks of LYRM4 in LIHC (LinkedOmics).

Additional file 14: Figure S16.

Protein-protein interaction (PPI) network of miRNA 495-target networks (GeneMANIA database).

Additional file 15: Figure S17.

PPI network of transcription factor PAX4-target networks (GeneMANIA database).

Additional file 16: Figure S18.

The effect of LYRM4 expression in correlation with the abundance of immune infiltrates on the prognosis of LIHC (TIMER2.0 database).

Additional file 17: Figure S19.

IscU mRNA expression levels in LIHC.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Pang, Y., Tan, G., Yang, X. et al. Iron-sulphur cluster biogenesis factor LYRM4 is a novel prognostic biomarker associated with immune infiltrates in hepatocellular carcinoma. Cancer Cell Int 21, 463 (2021). https://doi.org/10.1186/s12935-021-02131-3

Download citation

Keywords

  • LYRM4
  • Hepatocellular carcinoma
  • Prognosis
  • Diagnostic biomarker
  • Functional network analysis