1. Basic Medical School, Southwest Medical University, Luzhou 646000, Sichuan Province, China.
2. Key Laboratory of Epigenetics and Oncology, the Research Center for Preclinical Medicine, Southwest Medical University, Luzhou 646000, Sichuan Province, China.
3. Department of Pathology, the Affiliated Huaian No. 1 People's Hospital of Nanjing Medical University, Huai'an 223300, Jiangsu Province, China.
4. School of Medicine, Baylor College of Medicine, Houston 77030, Texas, USA.
5. Human Genome Sequencing Center, Baylor College of Medicine, Houston 77030, Texas, USA.
6. Department of Biochemistry, School of Life Sciences, Central South University, Changsha 410013, Hunan Province, China.
7. Department of Medical Technology, Faculty of Associated Medical Sciences, Chiang Mai University, Chiang Mai 50200, Thailand.
#Equal contributions to this work.
Furin is a proprotein convertase that activates different kinds of regulatory proteins, including SARS-CoV-2 spike protein which contains an additional furin-specific cleavage site. It is essential in predicting cancer patients' susceptibility to SARS-CoV-2 and the disease outcomes due to varying furin expressions in tumor tissues. In this study, we analyzed furin's expression, methylation, mutation rate, functional enrichment, survival rate and COVID-19 outcomes in normal and cancer tissues using online databases, and our IHC. As a result, furin presented with biased expression profiles in normal tissues, showing 12.25-fold higher than ACE2 in the lungs. The furin expression in tumors were significantly increased in ESCA and TGCT, and decreased in DLBC and THYM, indicating furin may play critical mechanistic functions in COVID-19 viral entry into cells in these cancer patients. Line with furin over/downexpression, furin promoter hypo-/hyper-methylation may be the regulatory cause of disease and lead to pathogenesis of ESCA and THYM. Furthermore, presence of FURIN-201 isoform with functional domains (P_proprotein, Peptidase_S8 and S8_pro-domain) is highest in all cancer types in comparison to other isoforms, demonstrating its use in tumorigenesis and SARS-Cov-2 entry into tumor tissues. Furin mutation frequency was highest in UCES, and its mutation might elevate ACE2 expression in LUAD and UCEC, reduce ACE2 expression in COAD, elevate HSPA5 expression in PAAD, and elevate TMPRSS2 expression in BRCA. These results showed that furin mutations mostly increased expression of ACE2, HSPA5, and TMPRSS2 in certain cancers, indicating furin mutations might facilitate COVID-19 cell entry in cancer patients. In addition, high expression of furin was significantly inversely correlated with long overall survival (OS) in LGG and correlated with long OS in COAD and KIRC, indicating that it could be used as a favorable prognostic marker for cancer patients' survival. GO and KEGG demonstrated that furin was mostly enriched in genes for metabolic and biosynthetic processes, retinal dehydrogenase activity, tRNA methyltransferase activity, and genes involving COVID-19, further supporting its role in COVID-19 and cancer metabolism. Moreover, Cordycepin (CD) inhibited furin expression in a dosage dependent manner. Altogether, furin's high expression might not only implies increased susceptibility to SARS-CoV-2 and higher severity of COVID-19 symptoms in cancer patients, but also it highlights the need for cancer treatment and therapy during the COVID-19 pandemic. CD might have a potential to develop an anti-SARS-CoV-2 drug through inhibiting furin expression.
Keywords: COVID-19, furin, malignant cancers, SARS-CoV-2, susceptibility, cordycepin (CD)
Furin, a paired basic amino acid cleaving enzyme, also known as FUR (FES upstream region), PACE (paired basic amino acid cleaving enzyme), PCSK3 (proprotein convertase subtilisin/kexin type 3), and SPC1 (subtilisin-like proprotein convertase 1) (OMIM: 136950), is cytogenetically located at 15q26.1 and was reported as a proprotein convertase that activates different kinds of regulatory proteins in the constitutive exocytic and endocytic pathway . Proprotein convertases are proteases that are synthesized as precursor proteins and require limited proteolysis to convert into the mature and bioactive proteins. Hendy et al. in 1995 first reported that furin is responsible for the physiologic processing of proparathyroid hormone to parathyroid hormone (PTH) ; Dubois et al. found that a biologically active transforming growth factor β1 (TGFB1) protein was cleaved by furin from pro-TGFB1 . Knockout of furin in transgenic mice is shown to be embryonic lethal thus, by conditional knocking out furin in T cells, Pesu et al. in 2008 showed that it allowed for normal T-cell development but impaired the regulatory function in effector T cells, which consequently generated less TGFB1 . Thus, furin may be vital in maintaining peripheral tolerance due to its nonredundant function in TGFB1 production. Furthermore, furin and other proprotein convertases have been widely reported to play proteolytic regulation of cell entry of viruses, including coronaviruses [5, 6].
Coronavirus disease 2019 (COVID-19) has quickly spread since December 2019 and the cases are rising worldwide [7-9]. At the end of August of 2021, the number of diagnosed cases worldwide is nearly 214 million and confirmed deaths are over 4.4 million (https://coronavirus.jhu.edu/). The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), responsible for COVID-19, contains an additional furin-specific cleavage site, a PRRA insertion at the S1/S2 boundary in the spike protein [10-12]. This specific cleavage site has been shown to promote viral infectivity and syncytia formation, but is not included in SARS-CoV or other coronaviruses [11, 13-15], supporting the hypothesis that this enzyme is indispensably involved in SARS-CoV-2 infection and COVID-19 pathogenesis . This is further substantiated by data that showed the loss of furin cleavage site weakens SARS-CoV-2 pathogenesis , whereas O-glycosylation modulates furin cleavage of the spike protein of SARS-CoV-2 , demonstrating a critical role for furin in SARS-CoV-2 pathogenesis and prevention. Inhibition of furin level in normal cells might help fight the viral infection, namely that of SARS-CoV-2. As a note, another protease, transmembrane serine protease 2 (TMPRSS2), also facilitates SARS-CoV-2 entry by cleaving and activating viral glycoproteins for viral uptake [18, 19].
As for COVID-19 viral infections, the critical event for the entry of genetic material into the host cell lies in the activation of the spike protein by host proteases to allow binding to host receptors. With the previous points in mind, the usage of furin inhibitors for COVID-19 therapy is urgently needed [20, 21], and the 3C-like protease (3CLpro) inhibitors for COVID-19 treatment may be the best option .
Since proprotein convertases require proteolysis to convert immature proteins to mature and bioactive proteins. Regulating furin activity disclosed its importance in processing cancer-related substrates and showed that high furin activity promotes tumorigenesis [6, 23]. For example, loss of furin in T cells suppresses mammary tumorigenesis of triple-negative breast cancer , likely by impairing its substrates proIGF1R and proIR processing . Furin also plays a critical role in KRAS and BRAF-related ERK/MAPK pathway activation and tumorigenesis of colorectal cancer . Furin's expression and its role in SARS-CoV-2 infected cancer patients is still unclear [6, 27]. Thus, it is critical to predict the cancer patients' susceptibility to SARS-CoV-2 infections and the disease outcomes by evaluating furin expression in different cancer tissues. In this study, we conducted expression profile analyses for furin in relation to COVID-19 in different types of normal tissues and cancer tissues to determine its potential role as a targeted therapeutic marker [20, 28, 29].
The furin gene and protein expression of the normal tissues were conducted in the Human Protein Atlas (HPA) (https://www.proteinatlas.org/ENSG00000140564-FURIN/tissue) [30-32]. The mRNA levels of furin in the HPA for single cell type (https://www.proteinatlas.org/ENSG00000140564-FURIN/celltype) and cell lines (v20.www.proteinatlas.org/ENSG00000140564-FURIN/cell) were also conducted respectively. The furin gene and protein expression levels in tumor tissues were evaluated in the HPA (v20.www.proteinatlas.org/ENSG00000140564-FURIN/pathology). The expression levels of the human angiotensin converting enzyme 2 (ACE2) gene in the lung normal tissues and lung cancer tissue were assessed in HPA, (v20.www.proteinatlas.org/ENSG00000140564-ACE2/tissue) and (v20.www.proteinatlas.org/ENSG00000140564-ACE2/pathology) respectively. The furin expression in multiple tumor tissues and corresponding normal control tissues through The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) databases were conducted using Gene Expression Profiling Interactive Analysis (GEPIA 2) (http://gepia2.cancer-pku.cn/#analysis) . Multiple gene comparison for furin and TMPRSS2 in tumor tissues and the corresponded TCGA normal and GTEx data were also evaluated using GEPIA 2.
Immunohistochemistry (IHC) was performed in Chinese lung and breast cancer tissues using the furin antibody [31, 32, 34]. Briefly, the de-paraffinized and re-hydrated tissue sections were treated in the 10 μM sodium citrate buffer at 95 °C for 12 minutes for antigen retrieval. The treated tissue slides were further incubated in the solution with 3% H2O2 for quenching the endogenous peroxidase. After blocking with 5% bovine serum albumin (BSA), the tissue slides were incubated with the primary antibody by overnight at 4 °C (rabbit monoclonal antibodies, 1:200; cat #: ab183495, Cambridge, CB2 0AX, UK) and then the biotin-conjugated secondary antibody (cat #: SP-9000, ZSGB-Bio, CN) for 60 minutes at room temperature. The bound secondary antibody was visualized by sequentially incubating with the Streptavidin-conjugated horseradish peroxidase (HRP) and its substrate diaminobenzidine (ZLI-9017, ZSGB-Bio, CN).
Cordycepin (CD) was purchased from Chengdu Must Bio-Technology Co.Ltd (Chengdu, Sichuan, P. R. China). Western blotting analysis for furin expression was performed in lung cancer cell line H1975 and breast cancer cell line BT549 with CD treatments (0, 10 µm, 20 µm, 40 µm) for 24 hours. The furin antibody for western blotting was used as same as IHC (rabbit monoclonal antibodies, cat #: ab183495, Cambridge, CB2 0AX, UK). The semi-quantitative RT-PCR for furin expression was performed in breast cancer cell line BT549 with CD treatments (0, 10 µm, 20 µm, 40 µm) for 24 hours. RT-PCR primers for furin were as follows, 5'-tgtggtgtaggtgtggccta-3' (RT-furin-L), 5'-gctgatggacagcgtgtaga-3' (RT-furin-R). GAPDH was served as internal control. All experiments were repeated at three times.
Homologs of furin in humans (NP_002560.1 in protein and NM_002569.4 in gene in GenBank, Ensembl ID: ENSG00000140564.10) and others were conducted from the National Center for Biotechnology Information (NCBI) program (https://www.ncbi.nlm.nih.gov/homologene?Db=homologene&Cmd=Retrieve&list_uids=1930) and were previously described [35, 36].
By analyzing GEPIA2 database (http://gepia2.cancer-pku.cn/#isoform) , we conducted the large datasets of TCGA and GTEx to determine furin isoform usage/distribution and domain structures in different tumor tissues.
The DNA methylation status in the furin promoter from the multiple tumor patients and the association between the furin expression and its promoter methylation in the normal and cancerous tissues, were performed through the DNA methylation interactive visualization database (DNMIVD) (http://18.104.22.168/dnmivd/query_gene/?cancer=ESCA&gene=furin)  and UALCAN (http://ualcan.path.uab.edu/index.html) .
Gene mutation module for furin, were conducted by TIMER2.0 (http://timer.comp-genomics.org/), and compared to the gene expression of furin, ACE2, HSPA5 and TMPRSS2, which are related to furin mutation status in TCGA database, Genes for furin, ACE2, HSPA5 and TMPRSS2 were determined to be important for SARS-CoV-2 infection.
Using fragments per kilobase of exon model per million reads mapped (FPKM) of the furin expression levels in patients, two expressed groups were classified and the correlations between furin expression and patient overall survival (OS) with the group cutoff at the median were evaluated in multiple cancer cohorts using GEPIA 2 (http://gepia2.cancer-pku.cn/#survival) to create Kaplan-Meier curves [33, 35, 39, 40]. Comparison of the survival contribution of furin, TMPRSS2, and ACE2 genes in multiple cancer types were estimated using the Mantel-Cox test, which was also conducted by GEPIA.
The data from Gene Ontology (GO), pathway of Kyoto Encyclopedia of Genes and Genomes (KEGG), and diseases/drugs of the genes co-expressed with furin were analyzed through the Enrichr database (https://maayanlab.cloud/Enrichr/enrich) . The p-value <0.01 was used as the cut off criterion.
The furin mRNA expression profiles were performed from the consensus datasets of HPA, GTEx and FANTOM5. The values for furin expression in the salivary gland was found to be highest at 337.9 NX, followed by the placenta (158.9), liver (118.4), pancreas (22.5), bone marrow (19.9); the lungs had the eighth highest expression (9.5). The B-cells was found to be lowest at 1.5 (Figure 1A). The protein levels of furin, based on the scores of high, medium, low, not detected, showed seven tissues with high expression including cerebral cortex, cerebellum, hippocampus, salivary gland, pancreas, kidney, placenta; sixteen tissues with medium expression; eleven tissues with low expression; ten tissues with no detected expression; and one tissue (soft tissue) showed low or no detected expression in different samples (Figure 1B). The expressed furin mRNA levels were also presented alongside high expression of protein. Expression profile for the furin mRNA in normal single cells showed the hepatocytes had the highest (115), and there was no expression in early spermatids, late spermatids, distal tubular cells, which are germ cells or epithelial cells (Figure 1C). Expression profile for the furin mRNA in cell lines showed the EFO-21 had the highest (46.8) and the daudi had the lowest (0.5) (Figure 1D). Thus, these findings revealed biased furin expression profiles for both mRNA and protein in different tissues, cells, and cell lines.
The furin expression in normal tissues in humans. A. The furin mRNA expression profiles in normal tissues. B. The furin protein expression profiles in normal tissues. C. The furin mRNA expression profiles in normal single cells. D. The furin mRNA expression profiles in cell lines. NX, consensus normalized expression.
ACE2 was reported to be the functional receptor for SARS-CoV-2, playing a crucial role in viral entry into human cells mainly through the lungs [18, 42]. Thus, comparison between furin and ACE2 mRNA expression were conducted in normal lungs to reveal possible correlations. We found that ACE2 expression value is 0.8 NX and furin expression value to be 12.25-fold higher (9.5/0.8=12.25) (Figure 2C), showing a correlation between ACE2 and furin expression levels and furin's potentially important role in COVID-19 pathogenesis in normal lung tissue.
Expression profile for the furin in cancer tissues from HPA database showed low cancer specificity in mRNA from 17 types of cancers (Figure 2A) and protein from 20 types of cancers (Figure 2B). Specifically, in terms of mRNA levels, liver cancer (LIHC) has the highest with 82.8 FPKM and lung cancer has the fourth with 44.6 FPKM, whereas glioma is lowest with 14.4. In terms of protein levels, thyroid cancer is the highest. Moderate cytoplasmic expression was positivity observed in lung, liver, prostate, and urothelial cancers. Only few nuclear staining was observed in several cases and the remaining cancer tissues were weakly or negative stained (Figure 2B, and data not shown).
Comparison between ACE2 and furin mRNA levels were conducted in 994 samples from the TCGA dataset. We found that ACE2 levels is 0.9 FPKM and furin levels is 44.6 FPKM, and thus levels for furin are 49.56-fold higher than that of ACE2 (Figure 2D). Further comparison between furin levels in normal and cancerous lung tissues found that furin levels are 4-fold (49.56/12.25=4.0) increased in cancerous lung tissues than that of the matched normal tissues (Figure 2E). This data implies that furin might play a critical role for SARS-CoV-2 entry and COVID-19 pathogenesis in lungs, particularly in lung cancer tissues.
The furin expressions in cancer tissues in humans. A. Expression profiles for the furin mRNA in cancer tissues. B. Expression profiles for the furin protein in cancer tissues. C. Relative mRNA levels for furin and ACE2 expression in lung normal tissues. D. Relative mRNA levels for furin and ACE2 expression in lung cancer tissues. E. Relative mRNA levels of furin in lung normal and tumor tissues in FPKM, Fragments Per Kilobase of exon per Million reads. The p values < 0.01.
The cytoplasmic and membranous sub-localization for receptors, particularly for membranous sub-localization, may facilitate the entry by virus. Thus to determine the spatial distribution of furin at a subcellular level in the human tissues, we did IHC in lung and breast tumor tissues, and the representative results are shown in Figure 4. IHC staining revealed the presence of cytoplasmic and membranous furin signals with moderate intensity in both breast cancer tissues (Figure 3A) and lung cancer tissues (Figure 3B). This cytoplasmic and membranous subcellular localization for furin in human tissues implies furin's possible role in viral infection.
Previous findings showed that furin levels have a 4-fold increase in lung cancer tissues than that of the matched normal tissues (Figure 2E). Thus, we further compared furin mRNA profile in tumor samples of 33 types of cancers, including the lungs, and their corresponding normal tissues. The results showed that all types of cancer tissues expressed furin, and the highest levels were found in LUAD (Figure 4A). Although most cancer tissues showed an increase of furin (Figure 4A), the furin expressions were significantly upregulated only in 2 types of cancers, esophageal carcinoma (ESCA) and testicular germ cell tumors (TGCT) (Figure 4A in red, 4B of left two panels, p<0.01). On the other hand, the furin expression levels were significantly decreased in other 2 types of cancers, including lymphoid neoplasm diffuse large B-cell lymphoma (DLBC) and thymoma (THYM) (Figure 4A in green, 4B of right two panels, p<0.01). Altogether, these findings imply that the furin might play a role in SARS-Cov-2 uptake into cells in some tumor tissues of certain cancer patients.
The furin cellular localization in cancer tissues of the lungs and breast. A, B. Representative staining for breast cancer patients. A 51-year-old Chinese female of infiltrating ductal carcinoma with tumor size of 3×1.5×1 cm (Her 2+, ER-, PR-). C, D. Representative staining for lung cancer patients. A 75-year-old Chinese male of medium-highly differentiated squamous cell carcinoma with tumor size of 2.5×1×0.8 cm (EGFR 2+). BCa, breast cancer, LCa, lung cancer. B, D. Enlarged images from A&C respectively.
To further understand whether methylation modification affects furin expression, the TCGA and GEO databases by DNMIVD was applied to evaluate the furin promoter methylation in ESCA, TGCT, DLBC, THYM tissues and their matched normal tissues. We analyzed these 4 cancers because furin expressions were significantly changed due to oncogenesis. Although the furin promoter methylation in TGCT and DLBC and their matched tissues are currently not available, the promoter methylation statuses in both ESCA and THYM were decreased in comparison to the matched normal tissues respectively (Figure 5A, p≤0.01; 5C, p=0.391), demonstrating furin expression and promoter methylation are significantly inversely correlated in ESCA, but not in THYM. Furthermore, using Spearman and Pearson correlations analysis, we disclosed an inverse correlation between the furin expression and its promoter methylation for ESCA and THYM tissues (Figure 5B & 5D). Thus, furin promoter methylations may be a regulatory mechanism for furin expression in the pathogenesis of ESCA and THYM patients.
ACE2 isoforms expressed in the airway epithelium have been reported to differentially contribute to host susceptibility to SARS-CoV-2 ; isoforms of other SARS-CoV-2 receptors or entry proteins might also play similar roles. Thus, we used GEPIA2 database to understand furin isoform prevalence and its structures in 31 types of tumor tissues, and the results are shown in Figure 6. From Figure 6A&B, we noticed that seven isoforms in total are expressed in tumor tissues but with differentially expressed levels (Figure 6A), and isoform ENST00000610579.4 (FURIN-201) usage is the highest in all cancer types, followed by ENST00000268171.7 (FURIN-001) as the second highest; other isoform levels were found to be very low or none (Figure 6B). Further isoform structure prediction showed that FURIN-001, FURIN-201 and FURIN-202 have the domains P_proprotein, Peptidase_S8 and S8_pro-domain, which encoded the same 794 amino acids (Figure 6C). FURIN-005 lacks all of P_proprotein and part of Peptidase_S8, whereas FURIN-003 lacks all of P_proprotein and Peptidase_S8, and part of S8_pro-domain (Figure 6C), demonstrating the functional roles for these isoforms for FURIN-001, FURIN-201 and FURIN-202 in tumorigenesis and SARS-CoV-2 entry into different tumor tissues, particularly FURIN-201.
Both furin and TMPRSS2 are proteinases [2, 3, 44] which are essential for SARS-CoV-2 proteolytic activation . Comparisons between furin and TMPRSS2 in different tumor tissues were conducted in TCGA normal datasets, and results are shown in Figure 7A. From Figure 7A, we disclosed that the furin expression is higher than that of TMPRSS2 in most of tumor tissues, except in PRAD, demonstrating that furin might facilitate tumorigenesis and COVID-19 viral entry in most types of cancers.
Homologs and conservation of the furin proteins in thirteen different species also revealed furin is highly conserved, demonstrating it might play a role in COVID-19 viral entry in other species as well (Figure 7B).
The furin expressions in tumor tissues and matched normal tissues. A. Expression profiles for furin in 33 types of tumor tissues (red files) and the matched normal tissues (green files) (TCGA normal and GTEx data). Tissue-wise expression is used in profiles. Green in DLBC and THYM indicates dowregulated expression whereas red in ESCA and TGCA indicates upregulated expression in tumor tissues compared with the matched normal tissues. B. Expression profiles for furin in four tumor tissues and their matched normal tissues (TCGA normal and GTEx data) (*: p<0.01). The tumor had the color red, whereas the normal tissue had the color grad. Tissue-wise expression is used in box plots. Right panel shows the full name of cancer types.
The furin promoter methylation in ESCA and THYM tumor tissues, and their corresponding normal tissues. A. The promoter methylation for regulating furin expression in ESCA. B. Pearson analysis for correlation between furin mRNA expression and its methylation in ESCA samples. C. The methylation of the promoter to regulate furin expression in THYM. D. Pearson analysis for correlation between furin mRNA expression and its methylation in THYM samples. Pearson correlation between methylation of furin gene promoter and FPKM was calculated.
Furin isoform usage and its structures in different types of cancers. A, B. Isoform usage for furin. The profiles for the expression distribution (violin plot) are presented in panel A, and isoform usage (bar plot) are presented in panel B. The x axis shows isoforms, and the Y axis shows cancer types. C. Isoform structures for furin are also shown. Multiple isoforms and visualized domains are presented in an interactive plot. Please note two isoforms are missing: ENST00000558794.1, ENST00000560018.1.
Tumor tissues usually contain gene mutations which causes malignant and recurrent after therapy. Furin mutations in cancers might affected or correlated other receptor expression, thus we conducted by TIMER2.0 analysis to look at the co-relationship between mutant furin and other SARS-CoV-2 receptor proteins. We firstly found that, in 23 types of cancers, UCES has the highest mutant frequency (29/531), whereas LGG had the lowest (1/525) (Figure 8A). Further analysis of the furin mutations' co-relation on ACE2, furin, HSPA5 and TMPRSS2 expression were conducted. We found that furin might elevate ACE2 expression in both LUAD (Figure 9B, C) and UCEC (Figure 9B, D); might reduce ACE2 expression in COAD (Figure 8B); might elevate HSPA5 expression in PAAD (Figure 9B, E); and might elevate TMPRSS2 expression in BRCA (Figure 8B, F). These results indicated that furin mutations might mostly increase expression in ACE2, followed by HSPA5 and TMPRSS2 in some cancers (Figure 8C-F), indicating furin mutations' possible regulatory roles for COVID-19 viral entry in cancers by increase of some receptor's expression. Of course, the patient number with furin mutants are less in PAAD which may hardly draw the solid conclusion. Further study for more samples including mutant patients should be conducted.
Given we found that furin expression is higher in some tumor tissues, clinical relationship between furin expression and OS outcomes was conducted. The results are shown in Figure 9, which shows high expression of furin is significant inversely correlated with long OS in LGG (Figure 9A, p<0.05), and inversely correlated, though not significantly, with long OS in LUSC (Figure 8B, p=0.055) and PCPG (Figure 9C, p=0.079). High expression of furin is significant correlated with long OS in COAD (Figure 9D, p>0.049) and KIRC (Figure 9E, p=0.0004). Thus, furin expression may be a favorable prognostic marker in above five types of cancer patient's survival. Given that furin and TMPRSS2 are both proteinases and ACE2 was reported as the main viral receptor for COVID-19, comparison between OS and TMPRSS2, ACE2 and furin was conducted and found that only KIRC and LGG show the same patterns of OS (Figure 8F), demonstrating their ubiquity in different cancers.
Expression comparisons between furin and TMPRSS2 in tumors and their corresponding normal TCGA and GTEx data (A). In this panel, an interactive heatmap is used where “T” represents tumor tissues and “N” normal tissues. Homologs and conservations of the furin proteins are presented in thirteen of the different species (B).
Genes co-expressed with furin were also determined. Co-expression analysis using GEPIA 2 revealed the top 100 genes most similar to furin in TCGA tumors, which has a similar expression pattern to other cancer types (Supplementary Table 1). The GO analyses from the Enrichr database (375 genes) are shown in Figure 10A-F. The results showed that we identified the following biological processes, including the alpha-amino acid metabolic process (GO: 1901605), aromatic amino acid family catabolic process (GO: 0009074), glucan biosynthetic process (GO: 0009250), glycogen biosynthetic process (GO: 0005978), kynurenine metabolic process (GO: 0070189), lipoprotein catabolic process (GO: 0042159), mitochondrion organization (GO: 0007005), molybdopterin cofactor biosynthetic process (GO: 0032324), molybdopterin cofactor metabolic process (GO: 0043545), NAD metabolic process (GO: 0019674), NADH metabolic process (GO: 0006734) NADP metabolic process (GO: 0006739), and organic cyclic compound biosynthetic process (GO: 1901362) (Figure 10A, p≤0.01). It also identified the following molecular functions, including retinal dehydrogenase activity (GO: 0001758) and tRNA (cytosine) methyltransferase activity (GO: 0016427) (Figure 10B, p≤0.01); and the following cellular components, including mitochondrion (GO: 0005739), mitochondrial matrix (GO: 0005759), and peroxisomal matrix (GO:0005782) (Figure 10C, p≤0.01). KEGG analysis further elucidated three enriched pathways: glyoxylate and dicarboxylate metabolism, tyrosine metabolism, and tryptophan metabolism (Figure 10D, p≤0.01). Diseases/drugs enrichment analysis found COVID-19-related gene sets, including downregulated by MHV-A59 in murine liver; COVID19-Nsp13 protein host PPI from Krogan, downregulated by SARS-CoV-2 infection in Vero E6 from GSE153940; COVID19-Nsp13 protein host PPI from Krogan (Figure 10E, p≤0.01); and OMIM disease, including glycogen storage disease (Figure 10F, p≤0.0238) and leigh syndrome (Figure 10F, p≤0.0274). Thus, from all this data, we concluded that the furin is mostly enriched in the genes for metabolic and biosynthetic processes, retinal dehydrogenase activity, tRNA methyltransferase activity, and diseases mainly involving COVID-19, demonstrating furin roles for COVID-19 and cancer metabolism.
Furin mutations in multiple cancer types and possible furin's effect on ACE2, HSPA5 and TMPRSS2 expressions. A. TIMER2.0 shows a bar plot presenting the furin mutant frequency in the indicated TCGA cancer types. B. The heatmap shows the log2 fold expression changes for genes of ACE2, furin, HSPA5 and TMPRSS2 for each cancer type. C. Differential ACE2 expression level in LUAD (sample numbers, n=511). D. Differential ACE2 expression level in UCEC (sample numbers, n=528). E. Differential HSPA5 expression level in PAAD (sample numbers, n=170). F. Differential TMPRSS2 expression level in BRCA (sample numbers, n=1017).
Overall survival (OS) analysis of furin expression from the cancer patients in LGG (A), LUSC (B), PCPG (C), COAD (D) and KIRC (E). The GENT2 databases were applied to evaluate on the TCGA cohort data and plot Kaplan-Meier curves. F. Comparison among the survival contribution of furin, TMPRSS2, ACE2 genes in multiple cancer types, estimated using Mantel-Cox test.
The results from analysis of GO enrichment and KEGG pathways. The enriched information for biological processes (A), molecular functions (B), cellular components (C) in analysis of ontologies, KEGG (D) in analysis of pathways, COVID-19 related gene sets (E), and OMIM disease (F) in diseases/drugs, were acquired via the Enrichr database, on the basis of the furin-associated genes. GO, gene ontology; KEGG, kyoto encyclopedia of genes and genomes.
Cordycepin (CD) inhibits furin expression in different cancer cells. A. Protein levels in lung cancer cell line H1975. B. Protein levels in breast cancer cell line BT549. C. RNA levels in breast cancer cell line BT549.
Cordycepin (CD) is one of the most functional components from traditional Chinese medicine (TCM) which showed a broad spectrum of biological activities including anticancer, anti-inflammatory, antidepressant, anti-virial replication, hepato-protective, and neuro-protective properties, etc. [46, 47]. To explore the potential in developing CD as an anti- SARS-CoV-2 drug, a furin expression analysis was conducted and the results showed that CD reduced furin expression of protein and/or mRNA in a dosage dependent manner on cancer cell lines H1975 and BT549 (Figure 11). These results implied that CD may have a potential to develop an anti-SARS-CoV-2 drug via inhibiting furin expression.
Patients with cancer are likely more susceptible to viral infection than those without cancer, and consequently more likely to become severely ill or die from COVID-19 [48-51]. Systematic review analyses showed that COVID-19 patients have more severe malignant cancer onset than that of patients without cancer (33.33% vs 16.09%) [50, 52]. Given this discrepancy, it is important to consider the expression of viral entry receptors and proteins in cancer tissues since cancer pathology may affect COVID-19 illness and susceptibility. However, it is unclear whether the patients' susceptibility to SARS-CoV-2 entry and disease severity are correlated with furin expression in different types of cancers.
The SARS-CoV-2 virus, responsible for COVID-19, contains an additional furin cleavage site that should promote viral infectivity, syncytia formation, and cell-cell fusion [10, 11, 53, 54]. This site is not presented in SARS-CoV and other coronaviruses [11, 13, 14]. Among these convertases, furin could play a predominant role in SARS-CoV-2 invasion. Furthermore, highly expressed SARS-CoV-2 entry receptor or proteins may facilitate viral infection [35, 55-57]. The expression levels for furin were also reported to be correlated to neutrophilic inflammation and inflammasome activation [56, 58]. Since the pandemic has been continuing at a rapid pace, understanding furin expressions and developing inhibitors as a form of COVID-19 therapy are urgently needed [20, 22]. Assessing furin expression in normal and tumor tissues will help foretell cancer patients' susceptibility to SARS-CoV-2 infections and disease outcomes. Even with its potentially vital prognostic and therapeutic use, furin expression and regulation in different malignant tumor tissues, in correlation with patient survival and susceptibility to SARS-CoV-2 entry, were unclear. In this study, using informatics analyzed from online databases, we revealed that furin is highly conserved in different species and highly expressed in normal tissues of humans, especially in comparison to other proteases like TMPRSS2 [44, 58], and found in the salivary gland, placenta, liver, pancreas. It is also increased in some cancer tissues. furin mRNA levels were 49.56-fold higher than that of ACE2 in lung normal tissues and 4-fold higher in lung cancer tissues, implying that SARS-CoV-2 might not only use TMPRSS2 or ACE2 to attack tissues from normal and cancerous organ tissues, but also use furin, particularly in lung cancer tissues. By analyzing the furin isoform expression and structures in 33 types of tissues, we found that furin expressed seven isoforms in different levels; for isoform usage, FURIN-201 is the highest in all cancer types. Furthermore, isoform structure prediction found that FURIN-001, FURIN-201 and FURIN-202 have the domains P_proprotein, Peptidase_S8, and S8_pro-domain, demonstrating their functional roles in tumorigenesis and SARS-CoV-2 infection in different tumor tissues. The GO and KEGG analyses found that the furin is mostly enriched in genes for metabolic and biosynthetic processes, retinal dehydrogenase activity, tRNA methyltransferase activity, and diseases mainly including COVID-19, demonstrating furin's roles in COVID-19 pathogenesis and cancer metabolism. Moreover, high furin expression is inversely correlated with long OS in LGG, LUSC and PCPG. Altogether, increased furin expression in cancer patients, such as LGG, LUSC and PCPG, might play a critical role in the susceptibility to SARS-CoV-2 uptake and severity of COVID-19 clinical symptoms.
The promoter methylation statuses in both ESCA and THYM were decreased in comparison to those of normal tissues, demonstrating furin expression and promoter methylation are significantly inversely correlated in ESCA, but not in THYM. Spearman and Pearson correlations disclosed an inverse correlation between furin expression and its promoter's methylation in ESCA and THYM tissues. Thus, furin promoter methylation may be one of the regulatory mechanisms for furin overexpression in ESCA tumors and pathogenesis. Tumor tissues usually have gene mutations which lead to recurrent malignancies. We found that UCES has the highest mutant frequency (29/531). Effect on ACE2, furin, HSPA6 and TMPRSS2 expression due furin mutation revealed that furin might elevate ACE2 expression in LUAD and UCEC, might reduce ACE2 expression in COAD, might elevate HSPA5 expression in PAAD, and might elevate TMPRSS2 expression in BRCA. These results indicated that furin mutations correlated with expression of ACE2, HSPA6 and TMPRSS2 in some cancers, implying furin's possible regulatory roles for the expression regulation of viral entry proteins.
Moreover, CD showed a broad spectrum of biological activities including anti-virial replication, anticancer, anti-inflammatory, antidepressant, hepato-protective, and neuro-protective properties [46, 47]. The recently study revealed that CD has a strong binding affinity with S protein and Mpro proteins of SARS-CoV-2 by molecular docking, and anti-SARS-CoV-2 testing further showed CD's antiviral action in vitro . Our furin expression analysis showed that CD reduced its expression in a dosage dependent manner, further supporting the possibility to develop an anti-SARS-CoV-2 drug. This also implies CD's role for anti-SARS-CoV-2 by inhibiting furin expression.
In conclusion, furin is highly expressed in normal tissues and increased significantly in some tumor tumors, implying the possibly higher susceptibility to SARS-CoV-2 uptake and possibly high disease severity of COVID-19 clinical symptoms for cancer patients. Our studies highlight the value of combating cancers by targeting therapeutic strategies for furin in the COVID-19 pandemic. CD might have a potential to develop an anti-SARS-CoV-2 drug through furin inhibition.
This work was supported by the Special Training Program for Young Science and Technology Talents from Southwest Medical University (grant no. 00031726), in part by the Translational Medicine Foundation of Southwest Medical University (grant no's. 00031476 and 00031477), and the National Natural Science Foundation of China (grant no's. 81672887 and 82073263).
The study was approved by the Ethical Committee of Southwest Medical University. The informed consent was obtained from tumor patients.
D. L., J. C., J. H., L. Z., X. L., X. C., S. L., Ji. F. did IHC, western blotting, collected and analyzed the data. J. F., J. C, and H. C designed and supervised the project. S. F, J.F. wrote and edited the manuscript.
The authors have declared that no competing interest exists.
1. Louagie E, Taylor NA, Flamez D, Roebroek AJ, Bright NA, Meulemans S. et al. Role of furin in granular acidification in the endocrine pancreas: identification of the V-ATPase subunit Ac45 as a candidate substrate. Proc Natl Acad Sci U S A. 2008;105:12319-24
2. Hendy GN, Bennett HP, Gibbs BF, Lazure C, Day R, Seidah NG. Proparathyroid hormone is preferentially cleaved to parathyroid hormone by the prohormone convertase furin. A mass spectrometric study. J Biol Chem. 1995;270:9517-25
3. Dubois CM, Laprise MH, Blanchette F, Gentry LE, Leduc R. Processing of transforming growth factor beta 1 precursor by human furin convertase. J Biol Chem. 1995;270:10618-24
4. Pesu M, Watford WT, Wei L, Xu L, Fuss I, Strober W. et al. T-cell-expressed proprotein convertase furin is essential for maintenance of peripheral immune tolerance. Nature. 2008;455:246-50
5. Izaguirre G. The Proteolytic Regulation of Virus Cell Entry by Furin and Other Proprotein Convertases. Viruses. 2019;11:837
6. Braun E, Sauter D. Furin-mediated protein processing in infectious diseases and cancer. Clin Transl Immunology. 2019;8:e1073
7. Wang C, Horby PW, Hayden FG, Gao GF. A novel coronavirus outbreak of global health concern. Lancet. 2020;395:470-3
8. Zheng J. SARS-CoV-2: an Emerging Coronavirus that Causes a Global Threat. International journal of biological sciences. 2020;16:1678-85
9. Deng CX. The global battle against SARS-CoV-2 and COVID-19. International journal of biological sciences. 2020;16:1676-7
10. Shang J, Wan Y, Luo C, Ye G, Geng Q, Auerbach A. et al. Cell entry mechanisms of SARS-CoV-2. Proc Natl Acad Sci U S A. 2020;117:11727-34
11. Coutard B, Valle C, de Lamballerie X, Canard B, Seidah NG, Decroly E. The spike glycoprotein of the new coronavirus 2019-nCoV contains a furin-like cleavage site absent in CoV of the same clade. Antiviral Res. 2020;176:104742
12. Walls AC, Park YJ, Tortorici MA, Wall A, McGuire AT, Veesler D. Structure, Function, and Antigenicity of the SARS-CoV-2 Spike Glycoprotein. Cell. 2020;181:281-92 e6
13. Cantuti-Castelvetri L, Ojha R, Pedro LD, Djannatian M, Franz J, Kuivanen S. et al. Neuropilin-1 facilitates SARS-CoV-2 cell entry and infectivity. Science. 2020;370:856-60
14. Daly JL, Simonetti B, Klein K, Chen KE, Williamson MK, Anton-Plagaro C. et al. Neuropilin-1 is a host factor for SARS-CoV-2 infection. Science. 2020;370:861-5
15. Johnson BA, Xie X, Kalveram B, Lokugamage KG, Muruato A, Zou J. et al. Furin Cleavage Site Is Key to SARS-CoV-2 Pathogenesis. bioRxiv. 2020. 2020 08.26.268854
16. Johnson BA, Xie X, Bailey AL, Kalveram B, Lokugamage KG, Muruato A. et al. Loss of furin cleavage site attenuates SARS-CoV-2 pathogenesis. Nature. 2021;591:293-9
17. Zhang L, Mann M, Syed Z, Reynolds HM, Tian E, Samara NL. et al. Furin cleavage of the SARS-CoV-2 spike is modulated by O-glycosylation. bioRxiv. 2021. 2021 02.05.429982
18. Hoffmann M, Kleine-Weber H, Schroeder S, Kruger N, Herrler T, Erichsen S. et al. SARS-CoV-2 Cell Entry Depends on ACE2 and TMPRSS2 and Is Blocked by a Clinically Proven Protease Inhibitor. Cell. 2020;181:271-80 e8
19. Li W, Moore MJ, Vasilieva N, Sui J, Wong SK, Berne MA. et al. Angiotensin-converting enzyme 2 is a functional receptor for the SARS coronavirus. Nature. 2003;426:450-4
20. Couture F, Kwiatkowska A, Dory YL, Day R. Therapeutic uses of furin and its inhibitors: a patent review. Expert Opin Ther Pat. 2015;25:379-96
21. Gioia M, Ciaccio C, Calligari P, De Simone G, Sbardella D, Tundo G. et al. Role of proteolytic enzymes in the COVID-19 infection and promising therapeutic approaches. Biochem Pharmacol. 2020;182:114225
22. He J, Hu L, Huang X, Wang C, Zhang Z, Wang Y. et al. Potential of coronavirus 3C-like protease inhibitors for the development of new anti-SARS-CoV-2 drugs: Insights from structures of protease and inhibitors. Int J Antimicrob Agents. 2020;56:106055
23. Jaaks P, Bernasconi M. The proprotein convertase furin in tumour progression. Int J Cancer. 2017;141:654-63
24. He Z, Khatib AM, Creemers JWM. Loss of the proprotein convertase Furin in T cells represses mammary tumorigenesis in oncogene-driven triple negative breast cancer. Cancer Lett. 2020;484:40-9
25. He Z, Khatib AM, Creemers JWM. Loss of Proprotein Convertase Furin in Mammary Gland Impairs proIGF1R and proIR Processing and Suppresses Tumorigenesis in Triple Negative Breast Cancer. Cancers (Basel). 2020 12
26. He Z, Thorrez L, Siegfried G, Meulemans S, Evrard S, Tejpar S. et al. The proprotein convertase furin is a pro-oncogenic driver in KRAS and BRAF driven colorectal cancer. Oncogene. 2020;39:3571-87
27. Mohanty A, Agnihotri S, Mehta A, Rawal S. COVID-19 and cancer: Sailing through the tides. Pathol Res Pract. 2021;221:153417
28. McKee DL, Sternberg A, Stange U, Laufer S, Naujokat C. Candidate drugs against SARS-CoV-2 and COVID-19. Pharmacological research. 2020;157:104859
29. Raghav PK, Kalyanaraman K, Kumar D. Human cell receptors: potential drug targets to combat COVID-19. Amino Acids. 2021;53:813-42
30. Uhlen M, Zhang C, Lee S, Sjostedt E, Fagerberg L, Bidkhori G. et al. A pathology atlas of the human cancer transcriptome. Science. 2017 357
31. Thul PJ, Akesson L, Wiking M, Mahdessian D, Geladaki A, Ait Blal H. et al. A subcellular map of the human proteome. Science. 2017 356
32. Uhlen M, Fagerberg L, Hallstrom BM, Lindskog C, Oksvold P, Mardinoglu A. et al. Proteomics. Tissue-based map of the human proteome. Science. 2015;347:1260419
33. Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic acids research. 2017;45:W98-W102
34. Zhang L, Yang M, Gan L, He T, Xiao X, Stewart MD. et al. DLX4 upregulates TWIST and enhances tumor migration, invasion and metastasis. International journal of biological sciences. 2012;8:1178-87
35. Fu J, Zhou B, Zhang L, Balaji KS, Wei C, Liu X. et al. Expressions and significances of the angiotensin-converting enzyme 2 gene, the receptor of SARS-CoV-2 for COVID-19. Molecular biology reports. 2020;47:4383-92
36. Fu J, Wei C, He J, Zhang L, Zhou J, Balaji KS. et al. Evaluation and characterization of HSPA5 (GRP78) expression profiles in normal individuals and cancer patients with COVID-19. International journal of biological sciences. 2021;17:897-910
37. Ding W, Chen J, Feng G, Chen G, Wu J, Guo Y. et al. DNMIVD: DNA methylation interactive visualization database. Nucleic acids research. 2020;48:D856-D62
38. Chandrashekar DS, Bashel B, Balasubramanya SAH, Creighton CJ, Ponce-Rodriguez I, Chakravarthi B. et al. UALCAN: A Portal for Facilitating Tumor Subgroup Gene Expression and Survival Analyses. Neoplasia. 2017;19:649-58
39. Fu J, Liao L, Balaji KS, Wei C, Kim J, Peng J. Epigenetic modification and a role for the E3 ligase RNF40 in cancer development and metastasis. Oncogene. 2021;40:465-74
40. Shen S, Wei C, Fu J. RNA-Sequencing Reveals Heat Shock 70-kDa Protein 6 (HSPA6) as a Novel Thymoquinone-Upregulated Gene That Inhibits Growth, Migration, and Invasion of Triple-Negative Breast Cancer Cells. Front Oncol. 2021;11:667995
41. Kuleshov MV, Jones MR, Rouillard AD, Fernandez NF, Duan Q, Wang Z. et al. Enrichr: a comprehensive gene set enrichment analysis web server 2016 update. Nucleic acids research. 2016;44:W90-7
42. Bourgonje AR, Abdulle AE, Timens W, Hillebrands JL, Navis GJ, Gordijn SJ. et al. Angiotensin-converting enzyme 2 (ACE2), SARS-CoV-2 and the pathophysiology of coronavirus disease 2019 (COVID-19). J Pathol. 2020;251:228-48
43. Blume C, Jackson CL, Spalluto CM, Legebeke J, Nazlamova L, Conforti F. et al. A novel ACE2 isoform is expressed in human respiratory epithelia and is upregulated in response to interferons and RNA respiratory virus infection. Nature genetics. 2021;53:205-14
44. Cheng J, Zhou J, Fu S, Fu J, Zhou B, Chen H. et al. Prostate adenocarcinoma and COVID-19: The possible impacts of TMPRSS2 expressions in susceptibility to SARS-CoV-2. J Cell Mol Med. 2021;25:4157-65
45. Bestle D, Heindl MR, Limburg H, Van Lam van T, Pilgram O, Moulton H. et al. TMPRSS2 and furin are both essential for proteolytic activation of SARS-CoV-2 in human airway cells. Life Sci Alliance. 2020: 3(9):e202000786.
46. Panya A, Songprakhon P, Panwong S, Jantakee K, Kaewkod T, Tragoolpua Y. et al. Cordycepin Inhibits Virus Replication in Dengue Virus-Infected Vero Cells. Molecules. 2021 26
47. Wei C, Yao X, Jiang Z, Wang Y, Zhang D, Chen X. et al. Cordycepin Inhibits Drug-resistance Non-small Cell Lung Cancer Progression by Activating AMPK Signaling Pathway. Pharmacological research. 2019;144:79-89
48. Penn I, Starzl TE. Immunosuppression and cancer. Transplant Proc. 1973;5:943-7
49. Yang K, Sheng Y, Huang C, Jin Y, Xiong N, Jiang K. et al. Clinical characteristics, outcomes, and risk factors for mortality in patients with cancer and COVID-19 in Hubei, China: a multicentre, retrospective, cohort study. Lancet Oncol. 2020;21:904-13
50. Zarifkar P, Kamath A, Robinson C, Morgulchik N, Shah SFH, Cheng TKM. et al. Clinical Characteristics and Outcomes in Patients with COVID-19 and Cancer: a Systematic Review and Meta-analysis. Clin Oncol (R Coll Radiol). 2021;33:e180-e91
51. Yang F, Shi S, Zhu J, Shi J, Dai K, Chen X. Clinical characteristics and outcomes of cancer patients with COVID-19. J Med Virol. 2020;92:2067-73
52. Zhang H, Han H, He T, Labbe KE, Hernandez AV, Chen H. et al. Clinical Characteristics and Outcomes of COVID-19-Infected Cancer Patients: A Systematic Review and Meta-Analysis. J Natl Cancer Inst. 2021;113:371-80
53. Papa G, Mallery DL, Albecka A, Welch LG, Cattin-Ortola J, Luptak J. et al. Furin cleavage of SARS-CoV-2 Spike promotes but is not essential for infection and cell-cell fusion. PLoS Pathog. 2021;17:e1009246
54. Mykytyn AZ, Breugem TI, Riesebosch S, Schipper D, van den Doel PB, Rottier RJ. et al. SARS-CoV-2 entry into human airway organoids is serine protease-mediated and facilitated by the multibasic cleavage site. Elife. 2021;10:e64508
55. Choudhary S, Sreenivasulu K, Mitra P, Misra S, Sharma P. Role of Genetic Variants and Gene Expression in the Susceptibility and Severity of COVID-19. Annals of laboratory medicine. 2021;41:129-38
56. Bradding P, Richardson M, Hinks TSC, Howarth PH, Choy DF, Arron JR. et al. ACE2, TMPRSS2, and furin gene expression in the airways of people with asthma-implications for COVID-19. J Allergy Clin Immunol. 2020;146:208-11
57. Torre-Fuentes L, Matias-Guiu J, Hernandez-Lorenzo L, Montero-Escribano P, Pytel V, Porta-Etessam J. et al. ACE2, TMPRSS2, and Furin variants and SARS-CoV-2 infection in Madrid, Spain. J Med Virol. 2021;93:863-9
58. Kermani NZ, Song WJ, Badi Y, Versi A, Guo Y, Sun K. et al. Sputum ACE2, TMPRSS2 and FURIN gene expression in severe neutrophilic asthma. Respir Res. 2021;22:10
59. Verma AK, Aggarwal R. Repurposing potential of FDA-approved and investigational drugs for COVID-19 targeting SARS-CoV-2 spike and main protease and validation by machine learning algorithm. Chem Biol Drug Des. 2021;97:836-53
Corresponding authors: Dr. Junjiang Fu, Key Laboratory of Epigenetics and Oncology, the Research Center for Preclinical Medicine, Southwest Medical University, Luzhou, Sichuan 646000, P R China. Tel/Fax: +86-830-3160283; E-mail: fujunjiangedu.cn; Jingliang Cheng, E-mail: jingliangcedu.cn; Dr. Hanchun Chen, E-mail: chenhanchunedu.cn.