Next Article in Journal
The Molecular Genetics of Dissociative Symptomatology: A Transdiagnostic Literature Review
Next Article in Special Issue
Identification and Functional Analysis of Individual-Specific Subpathways in Lung Adenocarcinoma
Previous Article in Journal
Whole-Genome Sequencing Reveals Lignin-Degrading Capacity of a Ligninolytic Bacterium (Bacillus cereus) from Buffalo (Bubalus bubalis) Rumen
Previous Article in Special Issue
In Silico Analysis of the L-2-Hydroxyglutarate Dehydrogenase Gene Mutations and Their Biological Impact on Disease Etiology
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Correlation of NTRK1 Downregulation with Low Levels of Tumor-Infiltrating Immune Cells and Poor Prognosis of Prostate Cancer Revealed by Gene Network Analysis

1
Department of Biology, Faculty of Sciences, University of Mohaghegh Ardabili, Ardabil 56199-11367, Iran
2
Department of Biology, Jahrom Branch, Islamic Azad University, Jahrom 74147-85318, Iran
3
Department of Biology, Faculty of Sciences, Shahid Chamran University of Ahvaz, Ahvaz 61357-83151, Iran
4
Centre for Interdisciplinary Research in Basic Sciences, Jamia Millia Islamia, New Delhi 110025, India
5
Department and Clinic of Internal Medicine, Angiology and Physical Medicine, Faculty of Medical Sciences in Zabrze, Medical University of Silesia, Batorego 15 St., 41-902 Bytom, Poland
*
Authors to whom correspondence should be addressed.
These authors have contributed equally to this work.
Genes 2022, 13(5), 840; https://doi.org/10.3390/genes13050840
Submission received: 20 March 2022 / Revised: 5 May 2022 / Accepted: 6 May 2022 / Published: 8 May 2022
(This article belongs to the Special Issue Bioinformatics Analysis for Diseases)

Abstract

:
Prostate cancer (PCa) is a life-threatening heterogeneous malignancy of the urinary tract. Due to the incidence of prostate cancer and the crucial need to elucidate its molecular mechanisms, we searched for possible prognosis impactful genes in PCa using bioinformatics analysis. A script in R language was used for the identification of Differentially Expressed Genes (DEGs) from the GSE69223 dataset. The gene ontology (GO) of the DEGs and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were performed. A protein–protein interaction (PPI) network was constructed using the STRING online database to identify hub genes. GEPIA and UALCAN databases were utilized for survival analysis and expression validation, and 990 DEGs (316 upregulated and 674 downregulated) were identified. The GO analysis was enriched mainly in the “collagen-containing extracellular matrix”, and the KEGG pathway analysis was enriched mainly in “focal adhesion”. The downregulation of neurotrophic receptor tyrosine kinase 1 (NTRK1) was associated with a poor prognosis of PCa and had a significant positive correlation with infiltrating levels of immune cells. We acquired a collection of pathways related to primary PCa, and our findings invite the further exploration of NTRK1 as a biomarker for early diagnosis and prognosis, and as a future potential molecular therapeutic target for PCa.

1. Introduction

The most common malignancy diagnosed in men worldwide is prostate cancer (PCa) [1]. In the most frequent cancers (2018), this malignancy takes the fourth place among other types of cancers [2]. Mortality rates caused by this heterogeneous disease have been static, but it does not alter the increased chances of prostate cancer by a man in the later ages of his life or a man with a family history of prostate cancer [2,3,4]. Statistics have also shown that many black men suffer from this situation compared to white or Asian men [4]. One of the known possible factors influencing the risks of prostate cancer is dietary habits [5]. Prostate cancer is caused by a number of variables, including smoking, obesity, race/ethnicity, food, age, chemical and radiation exposure, sexually transmitted illnesses, and so on [6]. However, the fundamental shift at the molecular level is the confirmed diagnosis of PCa. The glandular development and expression of luminal differentiation markers androgen receptor (AR) and prostate-specific antigen (PSA) identify most prostatic malignancies that are adenocarcinomas [7]. A blood test and biopsy-based PSA efficiently diagnose prostate cancer in the early stages [8]. This is due to mutations and other changes in the AR or signaling pathways that lead to its increased expression [9], with recent research revealing that the knocking out of the AR gene reduces PCa cell invasion and migration [10]. Although PSA is the most widely used biomarker for prostate cancer, it has low specificity and substantial limitations. Therefore, screening new early diagnostic and prognostic markers for the pathogenesis and prognosis of prostate cancer is essential [11]. The mutations in BRCA1 and BRCA2 can give rise to prostate cancer [12,13,14]. It is not expected that a single gene directs the pathogenesis of the disease, and usually, alterations in the gene expression profile form a pathogenetic network built up from interactions between multiple genes.
The genes with equivalent effects in a pathogenetic network are sited in the same functional portion defined as a module, and they cooperate to fulfill their biological function [15,16,17]. Developments in bioinformatics technologies, such as microarrays, transcriptome sequencing, and proteomics, have provided potential advancements in cancer biomarker research and have been surveyed in several studies on different types of cancer [11,18]. Several studies have identified genes that have a significant role in the onset and development of PCa, such as forkhead box A1 (FoxA1) [19], kallikrein-related peptidase 3 (KLK3) [20], insulin-like growth factor 2 (IGF2) [21], and phosphatase and tensin homolog (PTEN) [22]. The essential genes discovered by the previous research, on the other hand, are quite diverse from one another and have nothing in common, which might be explained by the fact that PCa is a heterogeneous illness in general. We used a single microarray data set of publicly available human primary prostate tissue and screened the differentially expressed genes (DEGs) in the first step. Gene annotation and pathway analysis were performed using gene ontology (GO) enrichment and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. A protein–protein interaction (PPI) network was then constructed to identify the hub genes (HUBs). This study aimed to find an effective prognostic gene in primary prostate cancer among HUBs by bioinformatics analysis and validated findings using the cancer genome atlas (TCGA) data. We used the UALCAN and cBioPortal databases to verify the expression levels and mutational conditions. We also looked into immune cell infiltration utilizing the tumor immune estimation resource (TIMER) database.

2. Materials and Methods

2.1. Microarray Data Extraction

The GSE69223 CEL files were obtained from the National Center for Biotechnology
Information (NCBI) gene expression omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo/, accessed on 20 March 2021) [23]. The data set comprised 30 samples; 15 were primary prostate cancer tissue, and 15 were adjacent normal prostate tissue. The tumor staging for primary prostate cancer was pT2 or pT3. Sequencing was performed using the GPL570 platform (Affymetrix Human Genome U133 Plus 2.0 Array) platform with 54,675 probes.

2.2. Data Preprocessing and Screening of DEGs

The unprocessed CEL files were background-corrected, normalized, and converted to an expression set by the “affy” package (https://www.r-project.org/, accessed on 20 March 2021, Version 4.0.5) using the MAS 5.0 expression measure (mas5) function and then log 2   transformation was applied to the expression values. The expression data were scaled using the scale function in R, followed by applying the principal component analysis (PCA) to remove the outlier samples. The probe IDs were translated into their HUGO Gene Nomenclature Committee (HGNC) gene symbols corresponding to the official sequencing platform. The maximum value across the probes was used to compute the expression value for a particular gene symbol represented by multiple probe IDs [24], and probes that did not have gene information were excluded. Finally, the Escape Excel plug-in [25] was used to prevent gene name mangling using Microsoft Excel (2019). Empirical Bayes statistics (eBayes, p-value < 0.05) was applied through the “limma” package [26] in R to uncover the DEGs between normal and cancer tissue, according to the following criteria: log 2 FC > 1.5 and adjusted p-value (FDR) < 0.05. The adjusted p-value was calculated using the Benjamini–Hochberg (BH) method by limma package.

2.3. Pathway and Functional Enrichment Analysis

The “clusterProfiler” package in R [27] was applied to annotate and visualize the functional profiles for the DEGs. GO term enrichment [28] and KEGG pathway analyses [29] were performed using this package. The cell component (CC), biological process (BP), and molecular function (MF) terms associated with the DEGs were characterized by the GO enrichment analysis. The KEGG pathway analysis uncovered biological pathways correlated with the DEGs. The threshold was set at an adjusted p-value < 0.05 (obtained using the BH procedure).

2.4. PPI Network Construction and HUBs Selection

The species Homo sapiens was chosen, and the PPI network was constructed using the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING, https://string-db.org/, accessed on 12 May 2021) v11.0 b database [30] to investigate the interactions between the DEGs. Only experimentally validated interactions with the minimum required interaction score 0.4 (default) were selected and retained. Cytoscape (v 3.9.0) [31] was used to analyze the network parameters for further HUBs identification and sub-network visualization. The eigenvector was computed using the “CentiScaPe 2.2” plug-in, while other topological parameters were computed using the Cytoscape network analyzer tool.

2.5. Overall Survival (OS) Analysis of HUBs

In this case, we used the Gene Expression Profiling and Interactive Analysis (GEPIA, http://gepia2.cancer-pku.cn/, accessed on 2 August 2021) v2 database [32], which is an RNA-Seq web server based on the UCSC Xena project, calculated by a standard pipeline. We used the GEPIA to see if the expressions of the hub genes were relevant to the survival of PCa patients in the TCGA cohort. Patients with a high level of expression (>median expression value) and patients with a low level of expression (<median expression value) were defined. The Kaplan–Meier (KM) method was used to evaluate overall survival [33] using a log-rank test (statistically significant: p-value < 0.05).

2.6. Validation of Prognostic HUBs Using cBioPortal and UALCAN Databases

The cBioPortal for Cancer Genomics (https://www.cbioportal.org/, accessed on 23 December 2021) [34] and UALCAN (http://ualcan.path.uab.edu/, accessed on 20 March 2022) [35] was tested for expression validation and investigation of genomic alteration frequencies, including mutations and CNA (amplifications and homozygous deletions) in our prognostic PCa HUBs, respectively. The prostate adenocarcinoma (PRAD) dataset of the TCGA was selected in both web-based tools for analysis.

2.7. Tumor Infiltration Analysis

The tumor immune estimation resource (TIMER) web-based tool (http://timer.cistrome.org/, accessed on 23 December 2021) [36] was queried to explore the correlation between tumor-infiltrating immune cells in PRAD patients and the expression levels of the prognostic PCa HUBs. The Spearman’s test was used, and the p-value < 0.05 was regarded as the statistically significant threshold.

3. Results

3.1. Data Preprocessing and Identification of DEGs

Because outlying microarray samples can dramatically bias further analysis, preprocessing procedures to identify and eliminate such samples from each dataset prior to network construction were critical [37]. Normalization, background correction, and perfect match (PM)/mismatch (MM) correction were applied using the “mas” algorithm of the “affy” package. Additionally, log2 (expression values + 1) was generated. Sample outlier detection (using PCA) was subsequent to the normalization and logarithmic transformation of all the probe sets. Box plots of 30 samples were plotted before and after normalization, which can be seen in Figure 1A,B, respectively. After normalization, each box’s median gene expression values became approximately centric, showing an appropriate normalization and a significant data distribution. Accordingly, no sample with considerable deference in Interquartile Range (IQR) was considered an outlier. In addition, Figure 1C,D depicts the PCA plot for the samples in our dataset before and after the omission of outliers, respectively. Between two groups (tumor and normal), four samples were farther away from their clusters, because of which there was an inherent overlap between the two clusters (tumor and normal), which were considered outliers. On removal of these four samples, we saw that the two clusters were distinct and independent from each other. In conclusion, the outlier samples were: GSM1695588_K400_PN (K400_PN), GSM1695600_K815_PN (K815_PN), GSM1695606_L083_PN (L083_PN), and GSM1695593_K643_PC (K643_PC).
The probe IDs were translated into HGNC gene symbols corresponding to the hgu133plus2.db platform using the annotate R package. Utilizing limma, 990 DEGs were identified in PCa (threshold: log 2 FC > 1.5 and adjusted p-value < 0.05), including 316 upregulated and 674 downregulated genes. Table 1 represents the top ten DEGs that were upregulated and downregulated based on the order of fold changes. Table 2 shows the known PCa biomarkers available within our list of DEGs. The volcano plot was used to visualize the DEGs dispersions (Figure 2A). The expression pattern of the DEGs between samples was revealed by hierarchical clustering analysis, indicating that tumor tissues have substantially different gene expression patterns compared to adjacent normal tissues (Figure 2B). Supplementary Tables S1 and S2 provided detailed information on the upregulated and downregulated DEGs.

3.2. DEGs Enrichment Analysis

As reported by the KEGG pathway analysis obtained from the clusterProfiler package, the DEGs were linked to pathways, including focal adhesion, dilated cardiomyopathy, protein digestion and absorption, hypertrophic cardiomyopathy, calcium signaling pathway, arrhythmogenic right ventricular cardiomyopathy, ECM-receptor interaction, PI3K-Akt signaling pathway, proteoglycans in cancer, and arginine and proline metabolism (Table 3). The top 20 enriched KEGG pathways are depicted in Figure 3A.
The results of the GO analysis, also performed using the clusterProfiler package, revealed mainly the enrichment of DEG in the collagen-containing extracellular matrix (belonging to the cellular component), the extracellular matrix (a member of the biological process), and the organization of the extracellular structure (part of the biological process) (Table 4). The 10 enriched GO terms for the cellular component, the biological process, and the molecular function are shown in Figure 3B. Supplementary Tables S3 and S4 provide comprehensive information on the GO and KEGG enrichment findings.

3.3. PPI Network Construction and HUBs Selection

Based on the STRING interaction score > 0.4, our network consisted of 886 nodes and 3579 edges built and visualized by Gephi (Figure 4). Various topological/centrality distributions of the PPI network can be seen in Figure 5. The top 50 genes ranked on the basis of four topological algorithms (i.e., Degree, Betweenness, Closeness, and Eigenvector) were carried out to detect the potential key genes in the PPI network. As shown by the Venn plot in Figure 6, 22 genes were overlapped within four topological algorithms and were considered HUBs in our study. The epidermal growth factor (EGF), transforming growth factor β 3 (TGFB3), brain-derived neurotrophic factor (BDNF), caveolin 1 (CAV1), cadherin 2 (CDH2), collagen type I α 2 chain (COL1A2), decorin (DCN), fibrillin 1 (FBN1), fibroblast growth factor 2 (FGF2), C-X-C motif chemokine ligand 12 (CXCL12), insulin-like growth factor 1 (IGF1), laminin subunit β 1 (LAMB1), matrix metallopeptidase 14 (MMP14), 5′-nucleotidase ecto (NT5E), neurotrophic receptor tyrosine kinase 1 (NTRK1), peroxisome proliferator activated receptor γ (PPARG), epithelial cell adhesion molecule (EPCAM), snail family transcriptional repressor 2 (SNAI2), bone morphogenetic protein 4 (BMP4), neural cell adhesion molecule 1 (NCAM1), secreted protein acidic and cysteine rich (SPARC), and vinculin (VCL) were our hub genes. EGF and EPCAM were upregulated, and others were downregulated. Among the DEGs, EGF had the highest node degree (99).

3.4. OS Analysis of HUBs

We used GEPIA to explore the prognostic worthiness of HUBs and performed an OS analysis. As shown in Figure 7, only NTRK1 (downregulated) had a statistically significant (log-rank p-value < 0.05) impact on the PCa patients’ OS. As a result, the relatively low expression of NTRK1 was significantly related to a poor prognosis of PCa, while the remaining 21 hub genes were not reported to be correlated to PCa prognosis (Supplementary Figure S1).

3.5. Validation of Prognostic HUBs Using cBioPortal and UALCAN Databases

The cBioPortal was used to investigate the specific genetic alterations of NTRK1 (prognostically significant) and 21 other hub genes across the TCGA-PRAD messenger RNA (mRNA) cohort comprising 501 patient samples. All these 22 HUBs were altered in 34% (i.e., 172 cases) of the patient samples, most of which were deep deletions (19.24%: altered in 96 cases) followed by amplifications (6.61%: altered in 33 cases), missense mutations (4.61%: altered in 23 cases), and multiple alterations (4.01%: altered in 20 cases). OncoPrint results for all these HUBs, as shown in Figure S2, revealed genetic alterations, including amplification, mRNA upregulation, mRNA downregulation, and various others in 34% of patient samples. The top five highly mutated HUBs were NT5E, SNAI2, FBN1, FGF2, and NTRK1. As only NTRK1 was prognostically significant, unlike other HUBs, we further analyzed only this gene and eliminated the rest. Figure 8A shows the overall alteration frequency barplot of NTRK1 in the PRAD dataset with a deep deletion frequency of 1.8% (9 samples), a missense mutation frequency of 1.2% (6 samples), and an amplification frequency of 0.4% (2 samples). As shown in Figure 8B, the Lollipop plot displayed the frequency and location of all possible mutations in the Pfam protein domains where NTRK1 had a somatic mutation frequency of 1.2% (i.e., eight missense mutations + 1 inframe mutation).
The NTRK1 expression in the TCGA-PRAD cohort was validated using the UALCAN database, based on the sample types, the nodal metastasis, the molecular signature, and the TP53 mutation status. As shown in Figure 8C, the expression level of NTRK1 in the primary tumor and the normal cells is significantly different (p-value = 1.73 × 10 3 ). In Figure 8D,F, a significant correlation of NTRK1 expression with nodal metastasis [normal vs. N0 (p-value = 1.74 × 10 3 ), normal vs. N1 (p-value = 1.62 × 10 3 )] and TP53 mutation status [normal vs. TP53-Mutant (p-value = 6.01 × 10 4 ), normal vs.TP53-nonmutant (p-value = 2.00 × 10 3 )] is depicted. Figure 8E represents significant differential expression based on molecular signatures [normal NTRK1 vs. ERG fusion (p-value = 6.62 × 10 4 ), normal vs. ETV1 fusion (p-value = 4.22 × 10 4 ), normal vs. FLI1 fusion (p-value = 2.12 × 10 3 ), and normal vs. SPOP mutation (p-value = 1.06 × 10 3 )].

3.6. Tumor Infiltration Analysis

The TIMER database was used to evaluate the significant correlation of NTRK1 with tumor-infiltrating immune cells across TCGA-PRAD cohort. The NTRK1 expression had a significant positive correlation with the infiltrating levels of T cell CD8+ ( r = 0.127 , p-value = 9.71 × 10 3 ), T cell CD4+ ( r = 0.275 , p-value = 1.17 × 10 8 ), dendritic cell ( r = 0.347 , p-value = 3.34 × 10 13 ), macrophage ( r = 0.123 , p-value = 1.18 × 10 2 ), neutrophil ( r = 0.264 , p-value = 4.47 × 10 8 ), and T cell NK ( r = 0.103 , p-value = 3.49 × 10 2 ), as shown in Figure 9, respectively.

4. Discussion

Finding new potential biomarkers to achieve the early detection of PCa in the primary stages and identify new molecular drug targets is necessary and plays a considerable role in helping patients most likely to benefit from treatment. Our study concentrated on a single cohort profile dataset through a microarray analysis. Therefore, we tried to achieve the highest possible quality in terms of gene expression levels by removing the outlier samples and using an efficient normalization method. We identified 990 DEGs between the primary PCa and the adjacent normal tissues, including 316 upregulated and 674 downregulated genes from the GEO dataset-GSE69223. It is apparent that the number of downregulated genes is notably higher than the upregulated genes. We performed a functional enrichment analysis using the clusterProfiler package to better understand the interactions among the DEGs. The DEGs were involved mainly in the organization of the collagen-containing extracellular structure, the organization of the extracellular matrix, the structural constituent of the extracellular matrix, the extracellular matrix, and the sarcolemma, according to the current enrichment analysis of the GO term (Table 4). The enrichment analysis of the KEGG pathway showed that the DEGs were involved in several cancer-related pathways: the calcium signaling pathway, the PI3K-Akt signaling pathway, the Wnt signaling pathway, the cGMP-PKG signaling pathway, and the focal adhesion, which are important in tumor growth and carcinogenesis [38,39,40,41,42,43]. In addition, the DEGs were mainly involved in the digestion of proteins and the absorption and cardiomyopathy pathways (Table 3). Our analysis represents NTRK1 as a new protein involved in PCa prognosis. Moreover, as we have shown in Supplementary Tables S1 and S2 (summarized in Table 2), we obtained some known PCa biomarkers such as α-methylacyl-CoA racemase (AMACR) [44,45], forkhead box A1 (FOXA1) [46,47,48], two members of the kallikrein related peptidase (KLK) gene family KLK2 [49,50] and KLK4 [51,52,53], prostate-cancer-associated 3 (PCA3) [54,55,56], distal-less homeobox 1 (DLX1), and a member of the homeobox (HOX) family HOXC6 [57], which play specific roles in PCa development.
AMACR plays a key role in the peroxisomal β-oxidation of dietary branched fatty acids. Previous studies have shown that AMACR is upregulated in prostate cancer. However, the mechanism underlying the correlation between AMACR and prostate cancer has not been clarified yet. In a meta-analysis of 22 studies on 4385 participants from various geographic regions, the results show the association between PCa risk and AMACR. In this study, AMACR expression is significantly associated with an increased diagnosis of PCa [58]. FOXA1 helps to shape AR signaling through direct interactions with the AR and drives the growth and survival of normal prostate and prostate cancer cells. FOXA1 also possesses an AR-independent role in regulating epithelial-to-mesenchymal transition [47]. Previous in vitro studies have shown that FOXA1 increases pro-angiogenic factors, including EGF, endothelin-1, and endoglin in prostate cancer cells and promotes endothelial cell proliferation, migration, and tube formation. Moreover, in vivo studies and a clinical samples investigation have shown that FOXA1 facilitates prostate cancer angiogenesis [59].
KLKs are involved in the regulation of tumor growth, neoplastic progression, tumor angiogenesis, and metastasis. Tailor et al. exhibited significant gene upregulation of KLK2 and KLK4 in PRAD. KLK2 can increase ECM degradation due to its proteolytic effects on fibronectin, laminin, gelatin, fibrinogen, and collagenases, leading to metastasis. In addition, KLK4 has been shown to increase the activation of plasmin via the activation of the urokinase plasminogen activator, which helps with the angiogenesis, invasion, and metastasis of the tumor [60].
PCA3 is significantly elevated in patients with prostate cancer, and several available studies show the utility of PCA3, as a urinary biomarker, for the diagnosis of early prostate cancer with reasonable specificity and sensitivity [61,62]. Liang et al. have shown that DLX1 is upregulated in the prostate clinical samples and exerts its oncogenic roles on PCa by activating β-catenin/TCF signaling and promoting the growth and migration of PCa cells [63]. In earlier studies, it has been reported that HOXC6 is involved in PCa development. Recently, Zhou et al. have shown that upregulated HOXC6 could participate in the progression of PCa and function as an independent prognostic marker for cancer [64].
Subsequently, after bioinformatical analysis and the construction of the PPI network, we identified potential HUBs in primary PCa (Figure 2C). One of the key genes with the highest degree (99) in the network was the EGF, which was upregulated in this study, and only NTRK1 (downregulated) was significantly (p-value < 0.05) associated with PCa prognosis (Figure 7). A UALCAN-based analysis in PRAD supported that the downregulation of NTRK1 was significantly associated with tumorigenesis (Figure 8C), and this expression decreases when the status of the nodal metastasis increases (Figure 8D). The TP53 mutation status does not affect the NTRK1 downregulation in tumor cells, according to (Figure 8F). As represented in (Figure 8A,B), there was an intermediate tumor mutational burden (TMB) in NTRK1. NTRK1 had a somatic mutation frequency of 1.2%, and most of the mutations are accumulated in the tyrosine kinase domain (Pkinase_Tyr).
NTRK1 or TrkA is a nerve growth factor (NGF) receptor that is part of the tyrosine kinase receptor family. Protein kinases play a critical role in PCa tumor proliferation, development, and metastasis [65], and several malignancies have been shown to cause changes in NTRK1 expression or mutations in this gene [66]. NTRK1 is actively involved in developing, protecting, and maintaining neurons [67,68,69,70]. This study shows the downregulation of NTRK1 in PCa cells, and thus its tumor-suppressive role is expected. However, the tumor-suppressive behavior of NTRK1 in PCa is controversial as many studies have delineated the pro-tumorigenic role of NTRK1 in PCa. The alternative splicing of NTRK1 could result in numerous different protein isoforms, and three of them (TrkAI, TrkAII, and TrkAIII) have been described in humans [66]. The human NTRK1, which is 25 kb in length, is located on chromosome 1q21-q22 and is made up of 17 exons. The full-length isoform is TrkAII, is mostly found in neuronal tissues, and if we remove exon 9 there will be TrkAI in most non-neuronal tissues. Exons 6, 7, and 9 of TrkAIII are missing, making it unable to bind to NGF; therefore, TrkAIII is autophosphorylated, does not bind to NGF, and antagonizes NGF/TrkAI signaling [66,71].
As we understood, the oncogenic and tumor-suppressive nature of NTRK1 depends on several things, such as the tumor environment and NTRK1 splicing patterns. Studies have shown that the upregulation of regular NTRK1 isoforms (TrkAI/II) in normoxia condition will have a good prognosis in neuroblastoma (NB), so it plays an antioncogenic role in NB. Moreover, under hypoxic conditions, the upregulation of TrkAIII occurs, which provides tumor progression and metastasis promotion; it plays an oncogenic role in NB and indicates relevance to the NB-regulated tumor-promoting switch by generating an angiogenic, stress-resistant, and tumorigenic NB phenotype via IP/Akt signaling [71,72,73,74,75]. At the molecular level, the activation of NTRK1 confers pro-differentiation programs by binding the specific ligand, binding the NGF, inhibiting angiogenesis, increasing immunogenicity, inducing the differentiation and growth arrest, and mediating apoptosis. On the contrary, downregulating NTRK1 results in proliferation and angiogenesis and thus tumor growth and aggressiveness [73]. Future studies on the splicing patterns of NTRK1 in PCa are needed, and we think that the downregulation of NTRK1 can still be meaningful in the poor prognosis of PCa patients.
According to (Figure 9), it can be inferred that there is a positive correlation between the suppression of the NTRK1 gene and the decrease in immune cell infiltration, such as T cell CD8+. CD8+ T cells and monocytes have been known to suppress PCa. The presence of iNKT cells has been shown to delay prostate cancer progression; however, M1 macrophages and neutrophils infiltration are associated with a poor prognosis [76]. The higher number of T cells, especially CD8+ T cells, memory cells, and CD4+ Th1 cells, can produce a better prognosis in some cancers [77]. Moreover, CD4+ T cells have been linked to human cancers, and they are thought to play a role in PCa growth and promotion [76,78]. As mentioned before, some immune cells such as CD8+ T, monocytes, and NKT cells have a cancer suppression pattern; thus, low levels of these cells’ infiltration could lead to a poor prognosis. However, other immune cells (such as M1 macrophages and neutrophils) have a cancer amplifier pattern. However, these are in line with our findings, which show that intermediate TMB in NTRK1 causes low levels of tumor infiltration of immune cells, notably CD8+ T cells, and the role of CD8+ T cells is more established over other immune cells in the killing of cancerous cells and is used in cancer immunotherapy [79]. The stage of the disease where the host immune response may decrease with increased tumor development has been proposed as a primary factor of immune cell infiltration [80]. Pajtler et al. have provided some evidence that NTRK1 leads to increased immunogenicity in neuroblastoma, which may contribute to a less malignant phenotype and/or the spontaneous regression of neuroblastoma cells [81]. It can be suspected that the suppression of the NTRK1 gene might ignite prostate cancer with the decrease in the level of infiltration of CD8+ T cells.
Finally, in the NTRK fusion-positive tumors, genomic co-alterations and DNA rearrangements were often found, notably in genes implicated in cell-cycle-associated regulators, PI3K signaling, the MAPK pathway, and the tyrosine kinase families [70]. Fusions have been regularly found in rare malignancies, including mammary analog secretory carcinoma, congenital infantile fibrosarcoma, secretory breast carcinoma, and congenital mesoblastic nephroma, as well as in several pediatric case malignancies. NTRK gene fusions can be found in a small percentage of frequent adult cancers, including head and neck cancer, non-small-cell lung cancer, colorectal cancer, salivary gland cancer, thyroid cancer, bladder cancer, brain tumors, and soft tissue sarcomas [66,82,83]. Figure 8E shows a decrease in the expression of the NTRK1- FLI1/ETV1/ERG gene fusions compared to normal tissue. The downregulation of NTRK1 has been detected in several cancers, such as PCa and breast cancer, which lead to tumor progression [84].
Note that (1) this research was entirely based on public databases. As a result, more research is required to confirm the validity of the results; (2) although our objective was to identify a trustworthy candidate associated with the prognosis of primary prostate cancer, there is always the probability of having missed some details.

5. Conclusions

Consequently, we recognized the key genes and pathways correlated with the pathogenesis and the prognosis of primary PCa by a bioinformatics analysis. The downregulation of NTRK1 was linked with the poor prognosis of PCa and may be used as a prognostic marker of primary PCa. Nevertheless, NTRK1 expression was shown to significantly correlate with immune cell infiltration levels, such as the CD8+ T cells currently used in cancer immunotherapy. Eventually, further molecular biological studies and clinical research are imperative to confirm these results and their specific functional roles in PCa.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/genes13050840/s1, Figure S1 an OS analysis of PCa HUBs based on the GEPIA database determined by the KM curve. The expression of other 21 HUBs was not shown to significantly impact the prognosis of PCa in using KM estimates (their Log-rank p-value ≥ 0.05); Figure S2 the OncoPrint summarizes the genomic alterations in 22 PCa HUBs across the TCGA-PRAD cohort, comprising 501 patient samples. The bottom row represents the frequency of the genomic alterations in these HUBs, with green, orange, grey, red, blue, and golden bars signifying missense, splice, truncating, amplification, deep deletion, and inframe mutations, respectively. The first, second, third, fourth, and fifth rows depict the clinical annotation bars such as samples per patient, profiled in mutations, and putative copy-number alterations from GISTIC, the mutation spectrum, and the mutation count; Table S1 a detailed information of the upregulated probes; Table S2 detailed information on the downregulated probes; Table S3 detailed information on the KEGG enrichment results; Table S4 detailed information on the GO enrichment results.

Author Contributions

Conceptualization, A.B.; methodology, A.B. and A.H.; validation, S.Z., A.S. and R.D.; formal analysis, A.B.; investigation, A.H. and N.S.; data curation, A.B. and A.H.; writing—original draft preparation, A.B., A.H., N.S. and P.S.; review and editing, S.Z., A.S. and R.D.; visualization, A.B., A.H., S.Z., A.S. and R.D.; supervision, S.Z., A.S. and R.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

This study created or analyzed no new data. Data sharing does not apply to this article.

Acknowledgments

The authors gratefully thank Parvaneh Nikpour for her facilitating guidance regarding the pathogenesis and prognosis concepts. We would like to thank Peyman Choopanian for his valuable advice about outlier sample detection. The technical assistance in R programming provided by Fayyaz Bahari is greatly appreciated. Thanks to Ehsan Bastami for his efforts and to the scientists who helped us with tumor immunology concepts.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jemal, A.; Center, M.M.; DeSantis, C.; Ward, E.M. Global patterns of cancer incidence and mortality rates and trends. Cancer Epidemiol. Prev. Biomark. 2010, 19, 1893–1907. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Mattiuzzi, C.; Lippi, G. Current Cancer Epidemiology. J. Epidemiol. Glob. Health 2019, 9, 217–222. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Frame, F.M.; Maitland, N.J. Epigenetic Control of Gene Expression in the Normal and Malignant Human Prostate: A Rapid Response Which Promotes Therapeutic Resistance. Int. J. Mol. Sci. 2019, 20, 2437. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Merriel, S.W.D.; Funston, G.; Hamilton, W. Prostate Cancer in Primary Care. Adv. Ther. 2018, 35, 1285–1294. [Google Scholar] [CrossRef] [Green Version]
  5. Chen, F.Z.; Zhao, X.K. Prostate cancer: Current treatment and prevention strategies. Iran. Red Crescent Med. J. 2013, 15, 279–284. [Google Scholar] [CrossRef] [Green Version]
  6. Shen, M.M.; Abate-Shen, C. Molecular genetics of prostate cancer: New prospects for old challenges. Genes Dev. 2010, 24, 1967–2000. [Google Scholar] [CrossRef] [Green Version]
  7. Tai, S.; Sun, Y.; Squires, J.M.; Zhang, H.; Oh, W.K.; Liang, C.-Z.; Huang, J. PC3 is a cell line characteristic of prostatic small cell carcinoma. Prostate 2011, 71, 1668–1679. [Google Scholar] [CrossRef] [Green Version]
  8. Kohaar, I.; Petrovics, G.; Srivastava, S. A Rich Array of Prostate Cancer Molecular Biomarkers: Opportunities and Challenges. Int. J. Mol. Sci. 2019, 20, 1813. [Google Scholar] [CrossRef] [Green Version]
  9. Watson, P.A.; Arora, V.K.; Sawyers, C.L. Emerging mechanisms of resistance to androgen receptor inhibitors in prostate cancer. Nat. Rev. Cancer 2015, 15, 701–711. [Google Scholar] [CrossRef] [Green Version]
  10. Lin, C.-Y.; Jan, Y.-J.; Kuo, L.-K.; Wang, B.-J.; Huo, C.; Jiang, S.S.; Chen, S.-C.; Kuo, Y.-Y.; Chang, C.-R.; Chuu, C.-P. Elevation of androgen receptor promotes prostate cancer metastasis by induction of epithelial-mesenchymal transition and reduction of KAT5. Cancer Sci. 2018, 109, 3564–3574. [Google Scholar] [CrossRef]
  11. Liu, S.; Wang, W.; Zhao, Y.; Liang, K.; Huang, Y. Identification of Potential Key Genes for Pathogenesis and Prognosis in Prostate Cancer by Integrated Analysis of Gene Expression Profiles and the Cancer Genome Atlas. Front. Oncol. 2020, 10, 809. [Google Scholar] [CrossRef]
  12. Cavanagh, H.; Rogers, K.M.A. The role of BRCA1 and BRCA2 mutations in prostate, pancreatic and stomach cancers. Hered. Cancer Clin. Pract. 2015, 13, 16. [Google Scholar] [CrossRef] [Green Version]
  13. Stone, L. The IMPACT of BRCA2 in prostate cancer. Nat. Rev. Urol. 2019, 16, 639. [Google Scholar] [CrossRef]
  14. Nyberg, T.; Frost, D.; Barrowdale, D.; Evans, D.G.; Bancroft, E.; Adlard, J.; Ahmed, M.; Barwell, J.; Brady, A.F.; Brewer, C.; et al. Prostate Cancer Risks for Male BRCA1 and BRCA2 Mutation Carriers: A Prospective Cohort Study. Eur. Urol. 2020, 77, 24–35. [Google Scholar] [CrossRef] [Green Version]
  15. Min, F.; Gao, F.; Liu, Z. Screening and further analyzing differentially expressed genes in acute idiopathic pulmonary fibrosis with DNA microarray. Eur. Rev. Med. Pharmacol. Sci. 2013, 17, 2784–2790. [Google Scholar]
  16. Tan, S.H.; Petrovics, G.; Srivastava, S. Prostate Cancer Genomics: Recent Advances and the Prevailing Underrepresentation from Racial and Ethnic Minorities. Int. J. Mol. Sci. 2018, 19, 1255. [Google Scholar] [CrossRef] [Green Version]
  17. Yang, G.; Chen, S.; Ma, A.; Lu, J.; Wang, T. Identification of the difference in the pathogenesis in heart failure arising from different etiologies using a microarray dataset. Clinics 2017, 72, 600–608. [Google Scholar] [CrossRef]
  18. Deng, S.P.; Zhu, L.; Huang, D.S. Predicting Hub Genes Associated with Cervical Cancer through Gene Co-Expression Networks. IEEE/ACM Trans. Comput. Biol. Bioinform. 2016, 13, 27–35. [Google Scholar] [CrossRef]
  19. Katoh, M.; Igarashi, M.; Fukuda, H.; Nakagama, H.; Katoh, M. Cancer genetics and genomics of human FOX family genes. Cancer Lett. 2013, 328, 198–206. [Google Scholar] [CrossRef]
  20. Hong, S.K. Kallikreins as Biomarkers for Prostate Cancer. BioMed Res. Int. 2014, 2014, 526341. [Google Scholar] [CrossRef] [Green Version]
  21. Livingstone, C. IGF2 and cancer. Endocr. Relat. Cancer 2013, 20, R321–R339. [Google Scholar] [CrossRef] [Green Version]
  22. Khan, M.M.; Mohsen, M.T.; Malik, M.Z.; Bagabir, S.A.; Alkhanani, M.F.; Haque, S.; Serajuddin, M.; Bharadwaj, M. Identification of Potential Key Genes in Prostate Cancer with Gene Expression, Pivotal Pathways and Regulatory Networks Analysis Using Integrated Bioinformatics Methods. Genes 2022, 13, 655. [Google Scholar] [CrossRef]
  23. Barrett, T.; Wilhite, S.E.; Ledoux, P.; Evangelista, C.; Kim, I.F.; Tomashevsky, M.; Marshall, K.A.; Phillippy, K.H.; Sherman, P.M.; Holko, M.; et al. NCBI GEO: Archive for functional genomics data sets—update. Nucleic Acids Res. 2013, 41, D991–D995. [Google Scholar] [CrossRef] [Green Version]
  24. Li, L.; Zhu, Z.; Zhao, Y.; Zhang, Q.; Wu, X.; Miao, B.; Cao, J.; Fei, S. FN1, SPARC, and SERPINE1 are highly expressed and significantly related to a poor prognosis of gastric adenocarcinoma revealed by microarray and bioinformatics. Sci. Rep. 2019, 9, 7827. [Google Scholar] [CrossRef] [Green Version]
  25. Welsh, E.A.; Stewart, P.A.; Kuenzi, B.M.; Eschrich, J.A. Escape Excel: A tool for preventing gene symbol and accession conversion errors. PLoS ONE 2017, 12, e0185207. [Google Scholar] [CrossRef] [Green Version]
  26. Ritchie, M.E.; Phipson, B.; Wu, D.; Hu, Y.; Law, C.W.; Shi, W.; Smyth, G.K. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015, 43, e47. [Google Scholar] [CrossRef]
  27. Yu, G.; Wang, L.-G.; Han, Y.; He, Q.-Y. clusterProfiler: An R package for comparing biological themes among gene clusters. Omics: A J. Integr. Biol. 2012, 16, 284–287. [Google Scholar] [CrossRef]
  28. Chen, L.; Zhang, Y.-H.; Wang, S.; Zhang, Y.; Huang, T.; Cai, Y.-D. Prediction and analysis of essential genes using the enrichments of gene ontology and KEGG pathways. PLoS ONE 2017, 12, e0184129. [Google Scholar] [CrossRef] [Green Version]
  29. Kanehisa, M.; Goto, S. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 2000, 28, 27–30. [Google Scholar] [CrossRef]
  30. Szklarczyk, D.; Gable, A.L.; Nastou, K.C.; Lyon, D.; Kirsch, R.; Pyysalo, S.; Doncheva, N.T.; Legeay, M.; Fang, T.; Bork, P.; et al. The STRING database in 2021: Customizable protein–protein networks, and functional characterization of user-uploaded gene/measurement sets. Nucleic Acids Res. 2021, 49, D605–D612. [Google Scholar] [CrossRef]
  31. Shannon, P.; Markiel, A.; Ozier, O.; Baliga, N.S.; Wang, J.T.; Ramage, D.; Amin, N.; Schwikowski, B.; Ideker, T. Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res 2003, 13, 2498–2504. [Google Scholar] [CrossRef] [PubMed]
  32. Tang, Z.; Kang, B.; Li, C.; Chen, T.; Zhang, Z. GEPIA2: An enhanced web server for large-scale expression profiling and interactive analysis. Nucleic Acids Res. 2019, 47, W556–W560. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Guyot, P.; Ades, A.E.; Ouwens, M.J.N.M.; Welton, N.J. Enhanced secondary analysis of survival data: Reconstructing the data from published Kaplan-Meier survival curves. BMC Med. Res. Methodol. 2012, 12, 9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Cerami, E.; Gao, J.; Dogrusoz, U.; Gross, B.E.; Sumer, S.O.; Aksoy, B.A.; Jacobsen, A.; Byrne, C.J.; Heuer, M.L.; Larsson, E.; et al. The cBio Cancer Genomics Portal: An Open Platform for Exploring Multidimensional Cancer Genomics Data. Cancer Discov. 2012, 2, 401. [Google Scholar] [CrossRef] [Green Version]
  35. Chandrashekar, D.S.; Bashel, B.; Balasubramanya, S.A.H.; Creighton, C.J.; Ponce-Rodriguez, I.; Chakravarthi, B.V.S.K.; Varambally, S. UALCAN: A Portal for Facilitating Tumor Subgroup Gene Expression and Survival Analyses. Neoplasia 2017, 19, 649–658. [Google Scholar] [CrossRef]
  36. Li, T.; Fu, J.; Zeng, Z.; Cohen, D.; Li, J.; Chen, Q.; Li, B.; Liu, X.S. TIMER2.0 for analysis of tumor-infiltrating immune cells. Nucleic Acids Res. 2020, 48, W509–W514. [Google Scholar] [CrossRef]
  37. Oldham, M.; Horvath, S.; Konopka, G.; Iwamoto, K.; Langfelder, P.; Kato, T.; Geschwind, D. Identification and Removal of Outlier Samples (Illumina) Supplement for: “Functional Organization of the Transcriptome in Human Brain”. Dim (Dat1) 2021, 1, 197. [Google Scholar]
  38. Bong, A.H.L.; Monteith, G.R. Calcium signaling and the therapeutic targeting of cancer cells. Biochim. Et Biophys. Acta (BBA) Mol. Cell Res. 2018, 1865, 1786–1794. [Google Scholar] [CrossRef]
  39. Fajardo, A.M.; Piazza, G.A.; Tinsley, H.N. The Role of Cyclic Nucleotide Signaling Pathways in Cancer: Targets for Prevention and Treatment. Cancers 2014, 6, 436–458. [Google Scholar] [CrossRef] [Green Version]
  40. Jiang, N.; Dai, Q.; Su, X.; Fu, J.; Feng, X.; Peng, J. Role of PI3K/AKT pathway in cancer: The framework of malignant behavior. Mol. Biol. Rep. 2020, 47, 4587–4629. [Google Scholar] [CrossRef]
  41. Jung, Y.-S.; Park, J.-I. Wnt signaling in cancer: Therapeutic targeting of Wnt signaling beyond β-catenin and the destruction complex. Exp. Mol. Med. 2020, 52, 183–191. [Google Scholar] [CrossRef] [Green Version]
  42. Maziveyi, M.; Alahari, S.K. Cell matrix adhesions in cancer: The proteins that form the glue. Oncotarget 2017, 8, 48471–48487. [Google Scholar] [CrossRef] [Green Version]
  43. Nallanthighal, S.; Heiserman, J.P.; Cheon, D.-J. The Role of the Extracellular Matrix in Cancer Stemness. Front. Cell Dev. Biol. 2019, 7, 86. [Google Scholar] [CrossRef]
  44. Abou El-Kasem, F.; Abulkheir, I.; Sidhom, N.; Ismail, A.; Habashy, H.; Elsayed, H. 275. Role of immunohistochemical expression of AMACR as a prognostic and predictive biologic marker in advanced prostatic carcinoma. Eur. J. Surg. Oncol. 2016, 42, S139. [Google Scholar] [CrossRef]
  45. Box, A.; Alshalalfa, M.; Hegazy, S.A.; Donnelly, B.; Bismar, T.A. High alpha-methylacyl-CoA racemase (AMACR) is associated with ERG expression and with adverse clinical outcome in patients with localized prostate cancer. Tumor Biol. 2016, 37, 12287–12299. [Google Scholar] [CrossRef]
  46. Gerhardt, J.; Montani, M.; Wild, P.; Beer, M.; Huber, F.; Hermanns, T.; Müntener, M.; Kristiansen, G. FOXA1 Promotes Tumor Progression in Prostate Cancer and Represents a Novel Hallmark of Castration-Resistant Prostate Cancer. Am. J. Pathol. 2012, 180, 848–861. [Google Scholar] [CrossRef]
  47. Teng, M.; Zhou, S.; Cai, C.; Lupien, M.; He, H.H. Pioneer of prostate cancer: Past, present and the future of FOXA1. Protein Cell 2021, 12, 29–38. [Google Scholar] [CrossRef]
  48. Labbé, D.P.; Brown, M. Transcriptional Regulation in Prostate Cancer. Cold Spring Harb. Perspect. Med. 2018, 8. [Google Scholar] [CrossRef]
  49. Chao, J.; Chen, L.-M.; Chai, K.X. Chapter 608—Human Kallikrein-related Peptidase 2. In Handbook of Proteolytic Enzymes, 3rd ed.; Rawlings, N.D., Salvesen, G., Eds.; Academic Press: Cambridge, MA, USA, 2013; pp. 2762–2765. [Google Scholar]
  50. Williams, S.A.; Xu, Y.; De Marzo, A.M.; Isaacs, J.T.; Denmeade, S.R. Prostate-specific antigen (PSA) is activated by KLK2 in prostate cancer ex vivo models and in prostate-targeted PSA/KLK2 double transgenic mice. Prostate 2010, 70, 788–796. [Google Scholar] [CrossRef] [Green Version]
  51. Fizazi, K. The role of Src in prostate cancer. Ann. Oncol. 2007, 18, 1765–1773. [Google Scholar] [CrossRef]
  52. Seiz, L.; Kotzsch, M.; Grebenchtchikov, N.I.; Geurts-Moespot, A.J.; Fuessel, S.; Goettig, P.; Gkazepis, A.; Wirth, M.P.; Schmitt, M.; Lossnitzer, A.; et al. Polyclonal antibodies against kallikrein-related peptidase 4 (KLK4): Immunohistochemical assessment of KLK4 expression in healthy tissues and prostate cancer. Biol. Chem. 2010, 391, 391–401. [Google Scholar] [CrossRef]
  53. Zhou, H.-J.; Yan, J.; Luo, W.; Ayala, G.; Lin, S.-H.; Erdem, H.; Ittmann, M.; Tsai, S.Y.; Tsai, M.-J. SRC-3 Is Required for Prostate Cancer Cell Proliferation and Survival. Cancer Res. 2005, 65, 7976–7983. [Google Scholar] [CrossRef] [Green Version]
  54. Groskopf, J.; Aubin, S.M.; Deras, I.L.; Blase, A.; Bodrug, S.; Clark, C.; Brentano, S.; Mathis, J.; Pham, J.; Meyer, T.; et al. APTIMA PCA3 Molecular Urine Test: Development of a Method to Aid in the Diagnosis of Prostate Cancer. Clin. Chem. 2006, 52, 1089–1095. [Google Scholar] [CrossRef] [Green Version]
  55. Hessels, D.; van Gils, M.P.M.Q.; van Hooij, O.; Jannink, S.A.; Witjes, J.A.; Verhaegh, G.W.; Schalken, J.A. Predictive value of PCA3 in urinary sediments in determining clinico-pathological characteristics of prostate cancer. Prostate 2010, 70, 10–16. [Google Scholar] [CrossRef]
  56. Roobol, M.J.; Schröder, F.H.; van Leeuwen, P.; Wolters, T.; van den Bergh, R.C.N.; van Leenders, G.J.L.H.; Hessels, D. Performance of the Prostate Cancer Antigen 3 (PCA3) Gene and Prostate-Specific Antigen in Prescreened Men: Exploring the Value of PCA3 for a First-line Diagnostic Test. Eur. Urol. 2010, 58, 475–481. [Google Scholar] [CrossRef]
  57. Fujita, K.; Nonomura, N. Urinary biomarkers of prostate cancer. Int. J. Urol. 2018, 25, 770–779. [Google Scholar] [CrossRef] [Green Version]
  58. Jiang, N.; Zhu, S.; Chen, J.; Niu, Y.; Zhou, L. A-Methylacyl-CoA Racemase (AMACR) and Prostate-Cancer Risk: A Meta-Analysis of 4,385 Participants. PLoS ONE 2013, 8, e74386. [Google Scholar] [CrossRef]
  59. Su, Y.; Zhang, Y.; Zhao, J.; Zhou, W.; Wang, W.; Han, B.; Wang, X. FOXA1 promotes prostate cancer angiogenesis by inducing multiple pro-angiogenic factors expression. J. Cancer Res. Clin. Oncol. 2021, 147, 3225–3243. [Google Scholar] [CrossRef]
  60. Tailor, P.D.; Kodeboyina, S.K.; Bai, S.; Patel, N.; Sharma, S.; Ratnani, A.; Copland, J.A.; She, J.-X.; Sharma, A. Diagnostic and prognostic biomarker potential of kallikrein family genes in different cancer types. Oncotarget 2018, 9, 17876–17888. [Google Scholar] [CrossRef] [Green Version]
  61. Gunelli, R.; Fragalà, E.; Fiori, M. PCA3 in Prostate Cancer. Methods Mol. Biol. (Clifton N.J.) 2021, 2292, 105–113. [Google Scholar] [CrossRef]
  62. Jiang, Z.; Zhao, Y.; Tian, Y. Comparison of diagnostic efficacy by two urine PCA3 scores in prostate cancer patients undergoing repeat biopsies. Minerva Urol. E Nefrol. = Ital. J. Urol. Nephrol. 2019, 71, 373–380. [Google Scholar] [CrossRef] [PubMed]
  63. Liang, M.; Sun, Y.; Yang, H.-L.; Zhang, B.; Wen, J.; Shi, B.-K. DLX1, a binding protein of beta-catenin, promoted the growth and migration of prostate cancer cells. Exp. Cell Res. 2018, 363, 26–32. [Google Scholar] [CrossRef] [PubMed]
  64. Zhou, J.; Yang, X.; Song, P.; Wang, H.; Wang, X. HOXC6 in the prognosis of prostate cancer. Artif. Cells Nanomed. Biotechnol. 2019, 47, 2715–2720. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Bagheri, S.; Rahban, M.; Bostanian, F.; Esmaeilzadeh, F.; Bagherabadi, A.; Zolghadri, S.; Stanek, A. Targeting Protein Kinases and Epigenetic Control as Combinatorial Therapy Options for Advanced Prostate Cancer Treatment. Pharmaceutics 2022, 14, 515. [Google Scholar] [CrossRef]
  66. Luberg, K.; Park, R.; Aleksejeva, E.; Timmusk, T. Novel transcripts reveal a complex structure of the human TRKA gene and imply the presence of multiple protein isoforms. BMC Neurosci. 2015, 16, 78. [Google Scholar] [CrossRef] [Green Version]
  67. Amatu, A.; Sartore-Bianchi, A.; Siena, S. NTRK gene fusions as novel targets of cancer therapy across multiple tumour types. ESMO Open 2016, 1. [Google Scholar] [CrossRef] [Green Version]
  68. Khotskaya, Y.B.; Holla, V.R.; Farago, A.F.; Mills Shaw, K.R.; Meric-Bernstam, F.; Hong, D.S. Targeting TRK family proteins in cancer. Pharmacol. Ther. 2017, 173, 58–66. [Google Scholar] [CrossRef]
  69. Nakagawara, A. Trk receptor tyrosine kinases: A bridge between cancer and neural development. Cancer Lett. 2001, 169, 107–114. [Google Scholar] [CrossRef]
  70. Okamura, R.; Boichard, A.; Kato, S.; Sicklick, J.K.; Bazhenova, L.; Kurzrock, R. Analysis of NTRK Alterations in Pan-Cancer Adult and Pediatric Malignancies: Implications for NTRK-Targeted Therapeutics. JCO Precis. Oncol. 2018, 2, 1–20. [Google Scholar] [CrossRef]
  71. Tacconelli, A.; Farina, A.R.; Cappabianca, L.; DeSantis, G.; Tessitore, A.; Vetuschi, A.; Sferra, R.; Rucci, N.; Argenti, B.; Screpanti, I.; et al. TrkA alternative splicing: A regulated tumor-promoting switch in human neuroblastoma. Cancer Cell 2004, 6, 347–360. [Google Scholar] [CrossRef] [Green Version]
  72. Tacconelli, A.; Farina, A.R.; Cappabianca, L.; Gulino, A.; Mackay, A.R. Alternative TrkAIII splicing: A potential regulated tumor-promoting switch and therapeutic target in neuroblastoma. Future Oncol. 2005, 1, 689–698. [Google Scholar] [CrossRef]
  73. Funke, L.; Bracht, T.; Oeck, S.; Schork, K.; Stepath, M.; Dreesmann, S.; Eisenacher, M.; Sitek, B.; Schramm, A. NTRK1/TrkA Signaling in Neuroblastoma Cells Induces Nuclear Reorganization and Intra-Nuclear Aggregation of Lamin A/C. Cancers 2021, 13, 5293. [Google Scholar] [CrossRef]
  74. Eggert, A.; Ikegaki, N.; Liu, X.-g.; Brodeur, G.M. Prognostic and Biological Role of Neurotrophin- Receptor TrkA and TrkB in Neuroblastoma. Klin Padiatr 2000, 212, 200–205. [Google Scholar] [CrossRef]
  75. Iraci, N.; Diolaiti, D.; Papa, A.; Porro, A.; Valli, E.; Gherardi, S.; Herold, S.; Eilers, M.; Bernardoni, R.; Valle, G.D.; et al. A SP1/MIZ1/MYCN Repression Complex Recruits HDAC1 at the TRKA and p75NTR Promoters and Affects Neuroblastoma Malignancy by Inhibiting the Cell Response to NGF. Cancer Res. 2011, 71, 404–412. [Google Scholar] [CrossRef] [Green Version]
  76. Wu, Z.; Chen, H.; Luo, W.; Zhang, H.; Li, G.; Zeng, F.; Deng, F. The Landscape of Immune Cells Infiltrating in Prostate Cancer. Front. Oncol. 2020, 10, 517637. [Google Scholar] [CrossRef]
  77. Abbas, A.K.; Lichtman, A.H.; Pillai, S. Cellular and Molecular Immunology; Saunders/Elsevier: Philadelphia, PA, USA, 2021. [Google Scholar]
  78. Miller, A.M.; Lundberg, K.; Özenci, V.; Banham, A.H.; Hellström, M.; Egevad, L.; Pisa, P. CD4+CD25high T Cells Are Enriched in the Tumor and Peripheral Blood of Prostate Cancer Patients. J. Immunol. 2006, 177, 7398. [Google Scholar] [CrossRef] [Green Version]
  79. Raskov, H.; Orhan, A.; Christensen, J.P.; Gögenur, I. Cytotoxic CD8+ T cells in cancer and cancer immunotherapy. Br. J. Cancer 2021, 124, 359–367. [Google Scholar] [CrossRef]
  80. Zidlik, V.; Brychtova, S.; Uvirova, M.; Ziak, D.; Dvorackova, J. The Changes of Angiogenesis and Immune Cell Infiltration in the Intra- and Peri-Tumoral Melanoma Microenvironment. Int. J. Mol. Sci. 2015, 16, 7876–7889. [Google Scholar] [CrossRef] [Green Version]
  81. Pajtler, K.W.; Rebmann, V.; Lindemann, M.; Schulte, J.; Schulte, S.; Stauder, M.; Leuschner, I.; Schmid, K.; Koehl, U.; Schramm, A.; et al. Expression of NTRK1/TrkA affects immunogenicity of neuroblastoma cells. Int. J. Cancer 2013, 133, 908–919. [Google Scholar] [CrossRef]
  82. Frattini, M.; Ferrario, C.; Bressan, P.; Balestra, D.; De Cecco, L.; Mondellini, P.; Bongarzone, I.; Collini, P.; Gariboldi, M.; Pilotti, S.; et al. Alternative mutations of BRAF, RET and NTRK1 are associated with similar but distinct gene expression patterns in papillary thyroid cancer. Oncogene 2004, 23, 7436–7440. [Google Scholar] [CrossRef] [Green Version]
  83. Gatalica, Z.; Xiu, J.; Swensen, J.; Vranic, S. Molecular characterization of cancers with NTRK gene fusions. Mod. Pathol. 2019, 32, 147–153. [Google Scholar] [CrossRef]
  84. Sreenivasan, S.; Thirumalai, K.; Krishnakumar, S. Expression profile of genes regulated by curcumin in Y79 retinoblastoma cells. Nutr. Cancer 2012, 64, 607–616. [Google Scholar] [CrossRef]
Figure 1. An overview of the preprocessing plots. The box plots show the expression value of the total of 30 samples before and after normalization. The ordinate reflects the gene expression values, while the abscissa displays the distinct samples. (A) The median gene expression values within each box (shown by the black line) were not equal before normalization. (B) The median of the expression values were nearly on the same line after normalizing, indicating satisfactory normalization performance. (C) The PCA plot for 30 samples. (D) The PCA plot for the 26 samples selected in this study (K400_PN, K815_PN, L083_PN, and K643_PC excluded). PCA, principal component analysis.
Figure 1. An overview of the preprocessing plots. The box plots show the expression value of the total of 30 samples before and after normalization. The ordinate reflects the gene expression values, while the abscissa displays the distinct samples. (A) The median gene expression values within each box (shown by the black line) were not equal before normalization. (B) The median of the expression values were nearly on the same line after normalizing, indicating satisfactory normalization performance. (C) The PCA plot for 30 samples. (D) The PCA plot for the 26 samples selected in this study (K400_PN, K815_PN, L083_PN, and K643_PC excluded). PCA, principal component analysis.
Genes 13 00840 g001
Figure 2. (A) The volcano plot highlighting the DEGs. A plot of log2FC vs. −log10 (adjusted-p-value) for the DEGs in PCa shows that the downregulated DEGs highlighted in blue are more than the upregulated DEGs shown in red. The X-axis represents the Log2FC, whereas the Y-axis displays the adjusted p-value (−log10 scale). (B) The hierarchical clustering analysis of the 990 DEGs. Each column indicates a sample, and each row demonstrates the gene expression level. The color scale ranges from red to green to represent high to low expression. (C) The expression heatmap of 22 HUBs. FC, fold change; DEGs, differentially expressed genes; and HUBs, hub genes.
Figure 2. (A) The volcano plot highlighting the DEGs. A plot of log2FC vs. −log10 (adjusted-p-value) for the DEGs in PCa shows that the downregulated DEGs highlighted in blue are more than the upregulated DEGs shown in red. The X-axis represents the Log2FC, whereas the Y-axis displays the adjusted p-value (−log10 scale). (B) The hierarchical clustering analysis of the 990 DEGs. Each column indicates a sample, and each row demonstrates the gene expression level. The color scale ranges from red to green to represent high to low expression. (C) The expression heatmap of 22 HUBs. FC, fold change; DEGs, differentially expressed genes; and HUBs, hub genes.
Genes 13 00840 g002
Figure 3. The scatter plots of the KEGG pathway and the GO enrichment. (A) The top 20 KEGG pathways are shown in a scatter plot. (B) The scatter plot of the top 10 separate GO terms. The Rich factor is the proportion of DEGs numbers among all gene numbers (in the database), indicated in the pathway term that is linked to it. KEGG, Kyoto Encyclopedia of Genes and Genomes; GO, gene ontology; DEGs, differentially expressed genes; BP, biological process; MF, molecular function; and CC, cellular component.
Figure 3. The scatter plots of the KEGG pathway and the GO enrichment. (A) The top 20 KEGG pathways are shown in a scatter plot. (B) The scatter plot of the top 10 separate GO terms. The Rich factor is the proportion of DEGs numbers among all gene numbers (in the database), indicated in the pathway term that is linked to it. KEGG, Kyoto Encyclopedia of Genes and Genomes; GO, gene ontology; DEGs, differentially expressed genes; BP, biological process; MF, molecular function; and CC, cellular component.
Genes 13 00840 g003
Figure 4. The PPI network of the DEGs. The circle size and color are set for degree and betweenness scores, respectively. The larger the circles, the greater the degree score. In addition, richly colored circles have a higher betweenness score. PPI network, protein–protein interaction network; DEGs, differentially expressed genes.
Figure 4. The PPI network of the DEGs. The circle size and color are set for degree and betweenness scores, respectively. The larger the circles, the greater the degree score. In addition, richly colored circles have a higher betweenness score. PPI network, protein–protein interaction network; DEGs, differentially expressed genes.
Genes 13 00840 g004
Figure 5. The topological property/centrality distribution plots showing the node degree distribution, the betweenness centrality, the closeness centrality, the clustering coefficient, the neighborhood connectivity, the average shortest path length distribution, the topological coefficient, and the eccentricity of PPI network. PPI network, protein–protein interaction network.
Figure 5. The topological property/centrality distribution plots showing the node degree distribution, the betweenness centrality, the closeness centrality, the clustering coefficient, the neighborhood connectivity, the average shortest path length distribution, the topological coefficient, and the eccentricity of PPI network. PPI network, protein–protein interaction network.
Genes 13 00840 g005
Figure 6. The Venn plot shows 22 overlapping HUBs between the top 50 genes ranked based on four topological algorithms: degree, eigenvector, closeness, and betweenness. The green, yellow, blue, and red areas represent the top 50 genes ranked based on closeness, betweenness, degree, and eigenvector. HUBs, hub genes.
Figure 6. The Venn plot shows 22 overlapping HUBs between the top 50 genes ranked based on four topological algorithms: degree, eigenvector, closeness, and betweenness. The green, yellow, blue, and red areas represent the top 50 genes ranked based on closeness, betweenness, degree, and eigenvector. HUBs, hub genes.
Genes 13 00840 g006
Figure 7. The KM curve was used to estimate OS in PCa patients according to the GEPIA database. The expression of NTRK1 was shown to have a significant impact (log-rank p-value < 0.05) on the PCa prognosis using KM estimates (log-rank test). The solid line shows the survival curve, and the dotted line represents the 95% CI. Patients with expression levels above the median are shown with red lines, while those with levels below the median are marked with blue lines. GEPIA, gene expression profiling interactive analysis. OS, overall survival; KM, Kaplan–Meier; PCa, prostate cancer; and CI, confidence interval.
Figure 7. The KM curve was used to estimate OS in PCa patients according to the GEPIA database. The expression of NTRK1 was shown to have a significant impact (log-rank p-value < 0.05) on the PCa prognosis using KM estimates (log-rank test). The solid line shows the survival curve, and the dotted line represents the 95% CI. Patients with expression levels above the median are shown with red lines, while those with levels below the median are marked with blue lines. GEPIA, gene expression profiling interactive analysis. OS, overall survival; KM, Kaplan–Meier; PCa, prostate cancer; and CI, confidence interval.
Genes 13 00840 g007
Figure 8. The validation of NTRK1 using the UALCAN and cBioPortal databases. (A) The barplot shows the alteration frequency of NTRK1 (3% of 499 PRAD cases) across the TCGA-PRAD dataset (TCGA, Firehose Legacy). The blue bar depicts 1.8% deep deletions, the green bar depicts 1.2% missense mutations, and the red bar depicts 0.4% amplifications. (B) The lollipop plot showing nine nonsynonymous mutations in NTRK1 protein domains. The grey horizontal bar represents the whole length of the NTRK1 protein, and the number of amino acids is displayed below the grey bar. The protein domains are shown with the red, blue, and green colored solid boxes. The locations and the frequencies of the mutations were denoted by the solid vertical lines and lollipop-like dots at their ends, respectively. The green and brown lollipops represent eight missense mutations and one inframe mutation. The NTRK1 expression levels in the TCGA-PRAD cohort are shown by box-and-whisker plots based on (C) the sample types, and (D) the nodal metastasis statuses (N0 means no regional lymph node metastasis; N1 means metastases in 1 to 3 axillary lymph nodes). (E) The molecular signatures, and (F) the TP53 mutation status via the UALCAN database. ** p-value < 0.01, and *** p-value < 0.001 vs. normal.
Figure 8. The validation of NTRK1 using the UALCAN and cBioPortal databases. (A) The barplot shows the alteration frequency of NTRK1 (3% of 499 PRAD cases) across the TCGA-PRAD dataset (TCGA, Firehose Legacy). The blue bar depicts 1.8% deep deletions, the green bar depicts 1.2% missense mutations, and the red bar depicts 0.4% amplifications. (B) The lollipop plot showing nine nonsynonymous mutations in NTRK1 protein domains. The grey horizontal bar represents the whole length of the NTRK1 protein, and the number of amino acids is displayed below the grey bar. The protein domains are shown with the red, blue, and green colored solid boxes. The locations and the frequencies of the mutations were denoted by the solid vertical lines and lollipop-like dots at their ends, respectively. The green and brown lollipops represent eight missense mutations and one inframe mutation. The NTRK1 expression levels in the TCGA-PRAD cohort are shown by box-and-whisker plots based on (C) the sample types, and (D) the nodal metastasis statuses (N0 means no regional lymph node metastasis; N1 means metastases in 1 to 3 axillary lymph nodes). (E) The molecular signatures, and (F) the TP53 mutation status via the UALCAN database. ** p-value < 0.01, and *** p-value < 0.001 vs. normal.
Genes 13 00840 g008
Figure 9. The scatter plots exhibiting significant positive correlations of NTRK1 with the infiltrating levels of T cell CD4+, T cell CD8+, neutrophil, dendritic cell, macrophage, and T cell NK across TCGA-PRAD cohort. In addition, the Spearman’s correlation value and the estimated statistical significance are shown as the legends for each scatter plot. PRAD, prostate adenocarcinoma.
Figure 9. The scatter plots exhibiting significant positive correlations of NTRK1 with the infiltrating levels of T cell CD4+, T cell CD8+, neutrophil, dendritic cell, macrophage, and T cell NK across TCGA-PRAD cohort. In addition, the Spearman’s correlation value and the estimated statistical significance are shown as the legends for each scatter plot. PRAD, prostate adenocarcinoma.
Genes 13 00840 g009
Table 1. The top 10 DEGs that were upregulated and downregulated between tumor and normal tissues. DEGs, differentially expressed genes; FC, fold change.
Table 1. The top 10 DEGs that were upregulated and downregulated between tumor and normal tissues. DEGs, differentially expressed genes; FC, fold change.
Genes SymbolProbe IDLog2FCAdjusted p-ValueState
OR51E2236121_at4.244412 4.76 × 10 4 Upregulated
DLX1242138_at3.898046 5.66 × 10 5 Upregulated
B3GAT1219521_at3.883768 4.70 × 10 7 Upregulated
LINC00992239319_at3.556998 2.40 × 10 5 Upregulated
DNASE1210165_at3.435794 4.77 × 10 7 Upregulated
FFAR2221345_at3.332304 1.67 × 10 6 Upregulated
FOLH1B211303_x_at3.277545 9.01 × 10 5 Upregulated
CLDN3203954_x_at3.057903 2.79 × 10 6 Upregulated
HPN204934_s_at3.054792 3.34 × 10 8 Upregulated
TRPM4219360_s_at3.014302 3.79 × 10 8 Upregulated
CXCL13205242_at−6.36684 3.98 × 10 9 Downregulated
SMTNL2229730_at−4.85304 1.42 × 10 9 Downregulated
SMR3B207441_at−4.66651 1.09 × 10 4 Downregulated
FBXL21P1555412_at−4.07877 1.59 × 10 6 Downregulated
NELL2203413_at−3.99192 6.67 × 10 13 Downregulated
SERPINA5209443_at−3.56837 2.17 × 10 7 Downregulated
BMP5205431_s_at−3.52423 6.47 × 10 10 Downregulated
RBP4219140_s_at−3.51651 5.50 × 10 6 Downregulated
FOXF2206377_at−3.49783 6.61 × 10 10 Downregulated
CA3204865_at−3.33564 2.15 × 10 8 Downregulated
Table 2. The previously known PCa biomarkers among our DEGs.
Table 2. The previously known PCa biomarkers among our DEGs.
Genes SymbolProbe IDLog2FCAdjusted p-ValueState
AMACR217111_at2.272572347 1.59 × 10 2 Upregulated
FOXA1204667_at1.988720507 2.92 × 10 6 Upregulated
KLK2210339_s_at1.577097564 1.51 × 10 4 Upregulated
KLK4224062_x_at1.84883585 1.64 × 10 3 Upregulated
PCA3232575_at2.979943566 6.63 × 10 3 Upregulated
DLX1242138_at3.89804597 5.66 × 10 5 Upregulated
HOXC6206858_s_at2.908075942 1.42 × 10 3 Upregulated
Table 3. The top 10 DEGs enrichment analysis of the KEGG pathway. KEGG, Kyoto Encyclopedia of Genes and Genomes. The Rich factor is the proportion of selected gene numbers (DEGs) compared to all gene numbers involved in each pathway term. The degree of pathway enrichment increases as the Rich factor increases.
Table 3. The top 10 DEGs enrichment analysis of the KEGG pathway. KEGG, Kyoto Encyclopedia of Genes and Genomes. The Rich factor is the proportion of selected gene numbers (DEGs) compared to all gene numbers involved in each pathway term. The degree of pathway enrichment increases as the Rich factor increases.
KEGG IDDescriptionCategoryGene CountRich FactorBH-p-Value
has04510Focal adhesionKEGG pathway2914.42% 2.06 × 10 5
has05414Dilated cardiomyopathyKEGG pathway1818.75% 6.87 × 10 5
has04974Protein digestion and absorptionKEGG pathway1817.47% 1.35 × 10 4
has05410Hypertrophic cardiomyopathyKEGG pathway1617.77% 3.10 × 10 4
has04020Calcium signaling pathwayKEGG pathway2711.25% 1.48 × 10 3
has05412Arrhythmogenic right ventricular cardiomyopathyKEGG pathway1316.88% 2.56 × 10 3
has04512ECM–receptor interactionKEGG pathway1415.90% 2.56 × 10 3
has04151PI3K-Akt signaling pathwayKEGG pathway339.32% 5.47 × 10 3
has05205Proteoglycans in cancerKEGG pathway2210.73% 9.18 × 10 3
has00330Arginine and proline metabolismKEGG pathway917.64% 1.69 × 10 2
Table 4. Top 10 GO enrichment analyses of the DEGs. GO, gene ontology.
Table 4. Top 10 GO enrichment analyses of the DEGs. GO, gene ontology.
GO TermDescriptionCategoryGene CountRich FactorBH-p-Value
GO:0062023collagen-containing extracellular matrixCC7518.47% 8.56 × 10 24
GO:0030198extracellular matrix organizationBP6216.84% 3.03 × 10 16
GO:0043062extracellular structure organizationBP6216.80% 3.03 × 10 16
GO:0005201extracellular matrix structural constituentMF4024.53% 9.28 × 10 16
GO:0042383sarcolemmaCC3022.05% 4.84 × 10 11
GO:0005539glycosaminoglycan bindingMF3917.03% 4.73 × 10 10
GO:0060485mesenchyme developmentBP4315.41% 1.91 × 10 9
GO:0048762mesenchymal cell differentiationBP3716.81% 3.84 × 10 9
GO:0048565digestive tract developmentBP2820.89% 6.63 × 10 9 .
GO:0042692muscle cell differentiationBP5012.98% 1.02 × 10 8
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bagherabadi, A.; Hooshmand, A.; Shekari, N.; Singh, P.; Zolghadri, S.; Stanek, A.; Dohare, R. Correlation of NTRK1 Downregulation with Low Levels of Tumor-Infiltrating Immune Cells and Poor Prognosis of Prostate Cancer Revealed by Gene Network Analysis. Genes 2022, 13, 840. https://doi.org/10.3390/genes13050840

AMA Style

Bagherabadi A, Hooshmand A, Shekari N, Singh P, Zolghadri S, Stanek A, Dohare R. Correlation of NTRK1 Downregulation with Low Levels of Tumor-Infiltrating Immune Cells and Poor Prognosis of Prostate Cancer Revealed by Gene Network Analysis. Genes. 2022; 13(5):840. https://doi.org/10.3390/genes13050840

Chicago/Turabian Style

Bagherabadi, Arash, Amirreza Hooshmand, Nooshin Shekari, Prithvi Singh, Samaneh Zolghadri, Agata Stanek, and Ravins Dohare. 2022. "Correlation of NTRK1 Downregulation with Low Levels of Tumor-Infiltrating Immune Cells and Poor Prognosis of Prostate Cancer Revealed by Gene Network Analysis" Genes 13, no. 5: 840. https://doi.org/10.3390/genes13050840

APA Style

Bagherabadi, A., Hooshmand, A., Shekari, N., Singh, P., Zolghadri, S., Stanek, A., & Dohare, R. (2022). Correlation of NTRK1 Downregulation with Low Levels of Tumor-Infiltrating Immune Cells and Poor Prognosis of Prostate Cancer Revealed by Gene Network Analysis. Genes, 13(5), 840. https://doi.org/10.3390/genes13050840

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop