Abstract Telomere maintenance mechanisms (TMMs) are important for cell survival and homeostasis. However, most related cancer research studies have used heterogenous bulk tumor tissue, which consists of various single cells, and the cell type properties cannot be precisely recognized. In particular, cells exhibiting non-defined TMM (NDTMM) indicate a poorer prognosis than those exhibiting alternative lengthening of telomere (ALT)-like mechanisms. In this study, we used bioinformatics to classify TMMs by cell type in gastric cancer (GC) in single cells and compared the biological processes of each TMM. We elucidated the pharmacological vulnerabilities of NDTMM type cells, which are associated with poor prognosis, based on molecular mechanisms. We analyzed differentially expressed genes in cells exhibiting different TMMs in two single-cell GC cohorts and the pathways enriched in single cells. NDTMM type cells showed high stemness, epithelial–mesenchymal transition, cancer hallmark activity, and metabolic reprogramming with mitochondrial abnormalities. Nuclear receptor subfamily 4 group A member 1 (NR4A1) activated parkin-dependent mitophagy in association with tumor necrosis factor-alpha (TNFA) to maintain cellular homeostasis without TMM. NR4A1 overexpression affected TNFA-induced GC cell apoptosis by inhibiting Jun N-terminal kinase/parkin-dependent mitophagy. Our findings also revealed that NR4A1 is involved in cell cycle mediation, inflammation, and apoptosis to maintain cell homeostasis, and is a novel potential therapeutic target in recalcitrant GC. Keywords: gastric cancer, non-defined telomere maintenance mechanism, NR4A1, parkin-dependent mitophagy 1. Background Cell immortality, one of the main characteristics of cancer cells, is mediated through telomere maintenance mechanisms (TMMs) [[26]1], which are activated dependently or independently of the telomerase (TEL) enzyme [[27]2]. In approximately 85% of human tumors, TEL ensures telomere maintenance, whereas the ALT mechanism only occurs in 10–15% [[28]3]. The mechanism underlying the maintenance of telomere has been widely investigated in pan-cancer studies [[29]4]. A previous study using bulk RNA-sequencing (RNA-seq) analysis of pan-cancer data of the Cancer Genome Atlas (TCGA) to classify TMMs into four types was the first to identify non-defined TMM (NDTMM) [[30]5]. NDTMM is associated with a poor prognosis in glioblastoma [[31]6]. Although TMM is essential for the immortalization of cancer cells, NDTMM has be shown to frequently occur [[32]7] where a TMM was not observed in the actual bulk tumor tissue. In particular, alternative lengthening of telomere (ALT) is associated with poor prognosis of gastric cancer (GC) and a stem-like molecular mechanism [[33]8]. However, most previous studies used bulk RNA-seq, which limits their relevance and applicability because the specific cell types were not considered [[34]5]. Each TMM differs in its immortalization process and progression into aggressive drug-resistant cancer cells. Although ALT-positive (ALT+) cells have been reported in fibroblasts, few studies have investigated cell types and TMM in other immune or stromal cells. Previous studies investigating the maintenance of telomere in ALT cells analyzed the proliferative potential and telomere dynamics of GM847 ALT cells (SV40 immortalized human skin fibroblasts) co-cultured with normal fibroblasts or TEL+ immortalized human cells [[35]7]. The results revealed that ALT phenotypic repressors were present in normal and some TEL+ immortal cells [[36]9]. Another study identified APBs in TERC+ keratinocytes and squamous cell carcinomas of mice, demonstrating that ALT and TEL coexist in the same cell (SCC) [[37]10]. In the current study, we used a bioinformatics approach to comprehensively classify TMM by cell type at the single-cell level, and comparatively analyzed the biological processes of each TMM type. 2. Methods 2.1. Data Preprocessing and Pathway Enrichment Analysis Two GC single-cell cohorts were used in this study [[38]11]. Data for eight original tumor samples from single cells collected from six patients with stomach cancer were downloaded from the Ji research group website ([39]https://dna-discovery.stanford.edu/research/datasets/ (accessed on 1 January 2022) [[40]11]. Seurat software v. 4.0.3 was used for the single-cell analysis [[41]12]. Quality filters were applied to 12,422 cells from the eight samples. We screened cells that met the following criteria for downstream analysis: cells with a (1) unique feature count of ≥200 and (2) mitochondrial count >5%. We identified eight different cell types. Single-cell libraries were created using the Chromium Single-Cell 3 Library & Gel Bead kit v2 (10× Genomics) according to the manufacturer’s instructions, and then sequenced using Illumina sequencers with a single-cell RNA-seq (scRNAseq) technique (i.e., NovaSeq, HiSeq, and NextSeq). We examined 84 metabolic pathways from the Kyoto Encyclopedia of Genes and Genomes (KEGG) database to determine their activity. Markov affinity-based graph imputation of cells (MAGIC) was used to impute missing values and restore the structure of the single-cell data [[42]13]. We used the single-sample gene set enrichment plugin in the R package gene set variation analysis (GSVA) [[43]14] to acquire information regarding the pathway activity of each cell. We calculated the GSVA score to analyze the TMM pathway enrichment, and attempted 10,000,000 runs to judge the accuracy and significance (false discovery rate [FDR] < 0.01). The TMM-related gene signature was obtained from a previous study [[44]15]. The second single-cell cohort (5927 cells) was downloaded from the Gene Express Omnibus (GEO) database ([45]GSE150290) from tumor samples from 23 patients [[46]16]. The method used to select the marker gene for cell type classification was described in a previous study [[47]11,[48]16]. MSigDB ([49]http://software.broadinstitute.org/gsea/msigdb, accessed on 2 February 2022) includes epithelial–mesenchymal transition (EMT) signatures as well as 50 cancer hallmark gene sets [[50]17]. We also used microarray data generated from bulk samples of 497 stomach cancer patients from the Yonsei cohort (registration number: [51]GSE84437, accessed 2 February 2022) from the GEO database [[52]18]. TCGA-stomach adenocarcinoma (STAD) was used for Gene Expression Profiling Interactive Analysis 2 (GEPIA2) [[53]19] to validate patient survival. 2.2. Analysis of DEGs and Single-Cell Meta-Analysis DESingle was used to examine the differentially expressed genes (DEGs) in cells based on their TMM type [[54]20]. Pathways enriched with METASCAPE [[55]21] for significantly affected genes selected using DEG, protein–protein interaction (PPI), molecular complex detection (MCODE), transcription factors (TF), and other software were analyzed. The R packages RaceID and StemID were used to identify single-cell stemness [[56]22]. Signatures related to stemness were obtained from previous studies [[57]23] and cell types were categorized from microarray Y497 bulk samples [[58]18] using the CIBERSORT software [[59]24]. 2.3. Drug Target Analysis In order to investigate gene regulatory interactions with genes overexpressed in the NDTMM type as input, we analyzed the DEGs [[60]20] to generate a list of genes with high or low expression in each TMM type and ConsensusPathDB (CPDB; [61]http://cpdb.molgen.mpg.de/ (accessed on 29 March 2022)) [[62]25]. The target gene and drug were discovered through a search of the Cancer Cell Line Encyclopedia (CCLE) [[63]26]. 3. Results 3.1. Single-Cell Analysis Revealed TMM Type in GC We analyzed single-cell data from two GC cohorts (12,422 and 5927 cells). By comparing cell type categorization and TMM type in single-cell data subjected to imputation, we determined the TMM types of single cells. TMM types were classified into four categories in order to determine the most prevalent type in each cell, and six pathways were assessed using the gene profile from our previous pan-cancer TMM study [[64]5]. In adenocarcinoma cells, telomerase (TEL) and ALT activity were observed simultaneously, and ALT activity was high in B, T, and NK cells, which are related to innate immunity. Telomerase activity was relatively higher in endothelial cells and granulocytes than it was in other cells ([65]Figure 1A). This trend was similar to that observed in the second GC single-cell cohort ([66]Figure 1B), and telomerase activity was higher in endothelial cells than it was in other cells (FDR < 0.01). In gland mucous cells (GMC) and fibroblasts, ALT activity was relatively high. According to analysis of the frequency of TMM types in different cells, 39% and 10% of adenocarcinoma cells showed ALT-like and 10% NDTMM, respectively, whereas very few cells exhibited telomerase activity ([67]Figure 1C). Figure 1. [68]Figure 1 [69]Open in a new tab Single−cell analysis revealed telomere maintenance mechanisms (TMMs) in gastric cancer (GC). Uniform manifold approximation and projection (UMAP) and heat map of TMM activity of eight single−cell types in (A) adenocarcinoma, B cells, endothelial cells, fibroblasts, granulocytes, macrophages, natural killer (NK) cells, and T cells (12,422 cells, cohort 1) and (B) enteroendocrine cells, fibroblasts, gland mucous cells (GMC), mesenchymal stem cells, pit mucous cell (PMC), tumor cells, and endothelial cells (5927 cells, cohort 2). (C) Bar graph showing frequency of TMM types in eight single−cell types (12,422 cells, left); different cohort eight single−cell types (5927 cells, right). Boxplot of epithelial–mesenchymal transition (EMT) signature of (D) adenocarcinoma cells, fibroblasts, (E) PMC, and GMC. (F) Boxplot of MKI67 expression in adenocarcinoma cells (cohort 1) and PMCs. This was similar to a study that reported approximately 30–40% and 30% ALT activity in a GC cohort and in tumor cells of another cohort, respectively [[70]5]. Specifically, in endothelial cells, granulocytes, and macrophages, the number of cells exhibiting ALT and telomerase activities was low and relatively high, respectively. These results suggest that the TMM of each cell type was different and the preferred type was used for maintenance. In addition, the various cell types also included those with a previously unknown TMM (NDTMM). We also analyzed the relationship of the four TMM and specific cell types to EMT. NDTMM showed the highest EMT activity in adenocarcinoma cells. In fibroblasts from the same cohort, telomerase exhibited the highest EMT activity. In another cohort, pit mucous cells (PMC) showed a similar tendency to that of adenocarcinoma cells, and although it was not significantly different from the ALT-like type, EMT was highest in the NDTMM type. In GMCs, ALT-like and NDTMM types exhibited similar trends ([71]Figure 1D,E). We considered TMM to be closely related to cell growth, maintenance, differentiation, and proliferation, and a clear difference existed in MKI67 expression. The number of ALT-like and NDTMM type cells was relatively low in adenocarcinoma; however, telomerase (TEL) and TEL+ALT-like type cells showed high proliferation. A similar trend was observed in PMC cells, and TEL+ALT-like type cells predominately showed high proliferation ([72]Figure 1F). Next, we analyzed the correlation between immune cells and TMM in a 497 Yonsei Hospital cohort (Y497) from the bulk transcriptome data, and the correlation was clearly separated (positively/negatively). Our analysis revealed that ALT chromatin-like type-related cells were T regulatory, naïve B, and memory B cells, and the remaining the immune cells were positively correlated with both TEL and ALT pathways ([73]Supplementary Figure S1A). We stratified patients by TMM type using bulk transcriptome data. In both the 497 Yonsei [[74]18] and TCGA-STAD cohorts, the four TMM types were not significantly distinguishable; however, the same pattern showed the best prognosis in the TEL only group, and the poorest prognosis was observed in the NDTMM type patients ([75]Supplementary Figure S1B). Therefore, these results indicate that the TMM of each cell type, which was not identified in the bulk sample, was detected at the single-cell level. Thus, TMM heterogeneity was observed by cell type. 3.2. Cancer Hallmark Activity and TMM at Single-Cell Level We divided adenocarcinoma cells into four TMM types and analyzed the differences in each cell type with respect to the characteristics observed in the bulk ([76]Figure 2A). Initially, 924 adenocarcinoma cells from the first cohort were classified into four TMM types, and 50 cancer hallmark pathways were analyzed. Contrary to our expectations, >60% of the cancer hallmark activity was confirmed in the NDTMM type cells. In the ALT-like type cells, both the KRAS and Wnt-catenin signaling pathways were highly enriched ([77]Figure 2B). Figure 2. [78]Figure 2 [79]Figure 2 [80]Open in a new tab Cancer hallmark and telomere maintenance mechanism (TMM) at the single-cell level. (A) UMAP of adenocarcinoma cells with four TMM types (ALT-like, TEL+ALT-like, non-defined telomere maintenance mechanism [NDTMM], TEL). (B) Heat map of cancer hallmark of adenocarcinoma cells categorized by four TMM types. (C) Boxplot of PML expression in the adenocarcinoma cells categorized by four TMM types. (D) tSNE plot of PMC stemness by entropy. (E) Bar graph of frequency of TMM types by stemness cluster. (F) Heat map of stemness related genes of four TMM types in adenocarcinoma. (G) Heat map of oncogene and tumor suppressor gene by four TMM types in adenocarcinoma. (H) Boxplot of autophagy and necrosis signature of the four TMM types in adenocarcinoma. We predicted the telomere length based on the expression of the promyelocytic leukemia (PML) gene, and determined that the higher the PML expression, the shorter the telomere length and the higher the ability to maintain telomeres. In our adenocarcinoma data, TEL+ALT-like was the most common TMM type and, interestingly, the prevalence of the NDTMM type was higher than that of the ALT-like type. In the ALT-like type cells, the telomeres were predicted to be long, whereas they were assumed to be shortened in the NDTMM type cells ([81]Figure 1C). We analyzed the stemness of cells in the PMC using StemID [[82]22] through entropy calculation using the stemness analysis. Cluster 9 showed the highest entropy, whereas clusters 1 and 7 showed low entropy (FDR < 0.05, [83]Figure 2D). The classification of the four TMM types for each cluster according to the degree of entropy revealed numerous TEL+ALT-like type cells in the high-entropy cluster, and ALT-like and NDTMM type cells were relatively common in clusters 1 and 7. In cluster 1, ALT-like and NDTMM type cells accounted for 50% and 20%, respectively ([84]Figure 2E). We performed a stemness marker analysis of adenocarcinomas in order to determine the stemness marker genes associated with TMM. The Yamanaka factor was mainly high in the TEL+ALT-like type, whereas the expression level of the normal stem cell marker genes (PROM1, ABCG2, and CD34) was high in the NDTMM type. Furthermore, the expression of the cancer stem cell marker gene was high in the NDTMM and TEL+ALT-like types, with similar ratios ([85]Figure 2F). Next, we analyzed the expression of the oncogenes and tumor suppressor genes. In the case of ALT-like type cells, most genes exhibited low activity. In the NDTMM, TEL, and TEL+ALT-like type cells, the opposite trend was observed ([86]Figure 2G). TMM is related to cell survival, and autophagy signature analysis confirmed that the TEL type was high but the NDTMM type was low, similarly to findings with the ALT-like type ([87]Figure 2H). However, with the necrosis signature, the NDTMM and ATL-like types showed the highest and lowest activity, respectively ([88]Figure 2H). These results confirmed that the NDTMM type cells had short telomeres just before death, possessed high EMT and stemness, and were highly related to cancer hallmarks. Conversely, the ALT-like type cells possessed elongated telomeres, and exhibited low necrosis and entropy. 3.3. Differentially Expressed Biological Pathways for TMM Type Currently, information about other biological pathways for each TMM type in GC is limited, and we previously reported ALT-like type TMM and epigenetic characteristics on the basis of bulk RNA-seq analysis in GC [[89]8]. In this study, differently enriched pathways for each TMM type at the adenocarcinoma single-cell level were analyzed. The ALT-like type was enriched for NABA ECM-affiliated pathways and regulation of lymphocyte activation, and the NDTMM type was enriched for Parkinson’s disease, neutrophil degranulation, and vascular endothelial growth factor-alpha (VEGFA) signaling pathways ([90]Figure 3A). Interestingly, for NDTMM, pathways such as multicellular organismal and cellular maintenance homeostasis were enriched. Figure 3. [91]Figure 3 [92]Figure 3 [93]Open in a new tab Biological pathways related to differentially expressed genes (DEGs) of different telomere maintenance mechanism (TMM) types in adenocarcinoma. Enriched biological network for (A) ALT-like (left) and NDTMM (right) and (B) NDTMM (left) and TEL+ALT-like (right) types of adenocarcinoma. (C) PPI network of nine clusters for NDTMM type analyzed using MCODE algorithm. (D) Volcano plot of DEGs between NDTMM and TEL+ALT-like types of adenocarcinoma. (E) PPI network of 12 clusters of TEL+ALT-like type cells analyzed using MCODE. (F) Bar graph of TF for ALT-like and NDTMM types. In contrast, for the TEL+ALT-like type, pathways such as RNA metabolism, Huntington’s disease, and translation were enriched ([94]Figure 3B). We analyzed the PPI of genes that were differently expressed between the NDTMM and ALT-like type using MCODE, and the genes were classified into nine MCODE clusters. Genes were mainly enriched in neutrophil degranulation, generation of precursor metabolites and energy, and ATP metabolism ([95]Figure 3C). It has been recognized that these genes create an environment for NDTMM type cells to generate a large amount of energy for survival. The maintenance of the gastrointestinal epithelium and epithelial structure of adenocarcinoma cells to maintain survival in a TMM-deficient environment is also of interest. Gene ontology (GO) term analysis of the nine MCODE clusters revealed that neutrophil degradation, multicellular organismal homeostasis, and tissue homeostasis were enriched. REG3A, H1F0, and KRT7 were highly expressed in the ALT-like type cells ([96]Supplementary Table S1). We compared genes that were significantly different between NDTMM and TEL+ALT groups using DEG analysis. In the NDTMM type, the expression of INPP1, TRIM15, FAM83H, FUOM, and ASL was high, and in the TEL+ALT-like type, the expression of NOP56, C8orf59, NUDC, and CD44 was high ([97]Figure 3D). Based on these results, we conducted an MCODE analysis of the PPI of genes that were highly significantly expressed in the TEL+ALT-like type, which were classified into 12 clusters. These clusters were enriched for RNA metabolism, Huntington’s disease, and human immunodeficiency virus (HIV) infection ([98]Figure 3E). We determined that the TFs in ALT-like and NDTMM types, general transcription factor IIE subunit 2 (GTF2E2), and PSMB5, were ALT-like, whereas proteasome 20S subunit beta 5 (PSMB5) was also found in NDTMM ([99]Figure 3F). Specifically, early growth response 2 (EGR2) [[100]27] was identified as a TEL inhibitory TF. 3.4. Landscape of Metabolic Reprogramming of TMM Types In ALT+ cells, telomere length is increased and PGC-1, a key regulator of mitochondrial biogenesis and function, is amplified or overexpressed in ALT+ tumors, which are very sensitive to PGC-1 or SOD2 knockdown. In order to improve anti-TEL cancer therapy, genetic modeling of TEL elimination exposes vulnerabilities that accidentally suppress mitochondrial homeostasis and oxidative defense systems [[101]28]. However, currently, few studies have reported metabolic reprogramming based on other TMM types. We systematically analyzed the metabolic reprogramming of adenocarcinoma and other cell types at the single-cell level using 84 metabolic pathways from KEGG [[102]29]. In the ALT-like type of adenocarcinoma cells, linoleic acid metabolism, nitrogen, glycosphingolipid biosynthesis, and glycosaminoglycan-related pathway activities were high. Tricyclic acid (TCA) and oxidative phosphorylation (OXPOHS) showed high activity in the TEL type cells; whereas in the NDTMM type cells, high activity was observed in >50% of the pathways. The metabolic pathway activity was not as high as expected in the TEL+ALT-like group, which also showed higher activity in riboflavin metabolism, glycosaminoglycan biosynthesis, and heparan sulfate than in other processes. The innate immune cells also demonstrated a pattern of metabolic reprogramming that was clearly distinguished according to TMM type ([103]Figure 4A). NDTMM and TEL type B cells showed high metabolic reprogramming and used various energy sources, whereas the ALT-like type cells showed high ether lipid metabolism activity ([104]Figure 4B). NDTMM type T cells showed high nitrogen, linoleic acid, histidine, and phenylalanine metabolism, and the TEL+ALT-like type demonstrated high OXPHOS and TCA activity ([105]Figure 4C). NK cells showed three types of TMM. Figure 4. [106]Figure 4 [107]Figure 4 [108]Open in a new tab Landscape of metabolic reprogramming of and telomere maintenance mechanism (TMM) types. (A) Heatmap of metabolic reprogramming landscape of cells with four TMM types in adenocarcinoma. Spider plots of top enriched metabolic pathways for (B) B, (C) T, and (D) NK cells. (E) Box plot of PPARGC1A, PPARGC1B (mitochondria bioenergetics master regulators) for adenocarcinoma cells. (F) Boxplot of mitochondria bioenergetics signature score for adenocarcinoma cells and PMC. Spider plots of top enriched metabolic pathways for (G) GMC, (H) intestinal metaplasia (IM), metaplasia stem cell (MSC), and (I) PMC. In the ALT-like type cells, nitrogen metabolism, steroid metabolism, and arginine activity were high, and in the NDTMM type, glycosaminoglycan biosynthesis and heparan sulfate, thiamine, and biotin metabolism were high ([109]Figure 4D). The master regulators of mitochondrial bioenergetics, PGC1A and PGC1B, showed the highest activity in NDTMM type adenocarcinoma cells ([110]Figure 4E). However, in the mitochondrial bioenergetics signature analysis, the TEL+ALT-like type and TEL type showed high activity, whereas the ALT type exhibited the lowest activity. This trend was also observed in the PMC of other GC single-cell cohorts, showing the highest and lowest activity in the TEL+ALT-like and NDTMM types, respectively ([111]Figure 4F). We analyzed the metabolic reprogramming of GMC, intestinal metaplasia (IM, MSC), and PMC in single-cell cohorts of different GCs. The ALT-like type GMC showed high nitrogen metabolism activity, similar to that of T and NK cells. Taurine and hypotaurine metabolism, as well as mucin type o-glycan biosynthesis, were high in both the ALT-like and NDTMM type cells ([112]Figure 4G). In the IM (MSE), TCA was high in the TEL+ALT-like group, and the ALT-like and NDTMM groups showed similar trends ([113]Figure 4H). In PMC, NDTMM and ALT-like types showed similar trends, and TEL+ALT-like type showed high activity in one-carbon metabolism, fatty acid biosynthesis, phenylalanine metabolism, and cysteine and methionine metabolism, unlike other types ([114]Figure 4I). These results demonstrate that the energy metabolism pathways changed according to TMM type at the single-cell level, and differed within the same cell type. Cells select TMM to maintain cellular homeostasis and, thus, different metabolic reprogramming methods are used as energy sources for survival. In summary, these results may provide insights into cell type-specific TMM inhibition and the development of metabolic inhibitory drugs. 3.5. Vulnerabilities of NDTMM Type for Cancer Therapy Our categorization of adenocarcinoma cells into four TMM types and examination of their features confirmed that NDTMM and ALT-like types were vital to the development of aggressive GC cells. We analyzed gene regulatory and drug–target interactions using the induced network modules of CPDB to identify a drug target for the NDTMM type [[115]25] ([116]Figure 5A). Figure 5. [117]Figure 5 [118]Open in a new tab Vulnerabilities of non-defined telomere maintenance mechanism NDTMM type for cancer therapy. (A) Gene regulatory interaction network for NDTMM type of adenocarcinoma constructed using ConsensusPathDB. (CPDB) (B) Drug–target interaction for NDTMM type of adenocarcinoma. (C) Kaplan–Meier plots show overall survival rates for high and low NR4A1 expression (p < 0.011). (D) Drug and target genes for gastric cell line molecular subtypes from the Cancer Cell Line Encyclopedia (CCLE). High drug sensitivity is indicated in green (negative correlation). When the intermediate node z-score threshold was set to 100, the drug targets and drugs at the upper level were as follows: chembl1449836, nuclear receptor subfamily 4 group A member 1 (NR4A1), MCL1/sodium benzoate, sodium phenylacetate, ASL, ASS1/GRASSYSTATIN A, PMID8410973C3, chembl1088572, CTSD, CASP4, PGC, and PTP4A1 ([119]Figure 5B). We investigated NR4A1 as the drug target because it plays an important role in the development of mitochondrial abnormalities in NDTMM type cells, as well as in cell maintenance in the absence of TMM. We also confirmed the association of NR4A1 with survival in the TCGA-STAD dataset ([120]Figure 5C). We also identified prospective therapeutics for NR4A1 (nilotinib, AZD6244, PF-2341066, and paclitaxel; in stomach cancer cell line, [121]Figure 5D). Tumor necrosis factor (TNF) causes mitochondrial dysfunction in GC, as demonstrated by Cyt-c leakage, mitochondrial membrane potential (MMP) collapse, and energy metabolism disturbance, which all activate cellular death processes. TNF therapy resulted in GC cell death, but also induced protective mitophagy. TNF boosted the activity of c-JNK, which raised parkin expression, subsequently triggering mitophagy to remove damaged mitochondria and prevent cell death. Mitophagy was reduced by overexpression of NR4A1, making GC cells more susceptible to TNF-induced cell death [[122]30]. 4. Discussion Telomere maintenance is essential for cancer cell survival and proliferation [[123]31]. However, our study demonstrates that cancer cells maintain their own survival and proliferation by changing and using the cancer microenvironment [[124]32] to their advantage, even where TMM is absent. We defined four types of TMM, but the current knowledge about the NDTMM type in GC is sparse. Although the ALT-like mechanism was associated with a poorer survival prognosis in GC, the NDTMM type was associated with a more deleterious condition than the ALT-like type was. Overall, the NDTMM type displayed strong cancer hallmark activity, as well as the highest EMT and stemness among the four TMM types. This feature of the NDTMM type led us to infer the length of telomeres based on PML expression. NDTMM type cells had telomeres that were shorter than those of the ALT-like type cells, but longer than those of the TEL+ALT-like type cells. NDTMM type cells were previously thought to be cancer cells immediately before death due to autophagy or high necrosis that could not maintain their telomeres. However, contrary to our hypotheses, the pathways that were enriched in this cell type were mostly metabolic and cellular homeostasis and maintenance pathways. Despite a lack of TMM, cancer cells can continuously attempt to avoid death. ALT-like cells were previously considered to possess a strong mitochondrial bioenergetic activity. However, our findings revealed that NDTMM type cells which were more active than ALT-like cells were present in more than half of the metabolic pathways. This tendency was observed in adenocarcinoma and other innate immune cells, and cell-specific metabolic reprogramming of the NDTMM type was exploited as an energy source. NR4A1 was determined to be a potential therapeutic target based on our observation of the regulatory interactions of the differential genes that were substantially expressed in the NDTMM type cells. In TMM-free cancer cells, NR4A1 prevents cell death via a JNK/parkin-dependent mitophagy. Overexpression of NR4A1 causes a mitochondrial energy imbalance by suppressing the expression of the mitochondrial respiratory complex. Our findings also demonstrated that TNF therapy stimulated parkin-dependent mitophagy, which, in excess, inhibits mitochondrial apoptosis, reducing the cytotoxicity of TNF. In contrast, NR4A1 overexpression reduced parkin-dependent mitophagy by inhibiting JNK. In addition, our study demonstrates the potential of NR4A1 as a novel drug target for killing NDTMM type cancer cells by inhibiting JNK/parkin-dependent mitophagy, making GC cells more susceptible to TNF-induced apoptosis. To the best of our knowledge, this is the first study to show that the NDTMM type is more detrimental in GC than the ALT-like type. Our findings will contribute to the development of precision medicine strategies for patients with refractory GC and provide unique insights into the development of anticancer therapies. 5. Conclusions This study demonstrates the survival, maintenance, and proliferation of GC cells in absence of a TMM. Cancer cells were classified into four types on based on TMM types, and the NDTMM type has been sparingly studied. We comparatively analyzed TMM activity, DEGs, PPI, and enriched GO terms among individual TMM type cells in GC, and demonstrated that the NDTMM type cells maintained their survival, proliferation, and homeostasis by altering their cell environments. Furthermore, we identified the NR4A1 gene as a potential therapeutic target for the development of efficient anticancer therapies. Supplementary Materials The following supporting information can be downloaded at: [125]https://www.mdpi.com/article/10.3390/cells11213342/s1, Figure S1: Correlation of immune cells and TMM pathway in Y497 cohort; Table S1: List of differential expressed genes (DEGs) between ALT-like type and NDTMM type. [126]Click here for additional data file.^ (872.8KB, zip) Author Contributions J.-Y.S.; methodology: J.-Y.S.; data analysis: J.-Y.S.; writing—original draft: J.-Y.S.; writing—review and editing: J.-Y.S.; supervision; J.-Y.S.; project administration: J.-Y.S.; formal analysis: J.-H.C.; funding acquisition: J.-Y.S. Both authors contributed to the interpretation of results. Both authors read and agreed to the published version of the manuscript. All authors have read and agreed to the published version of the manuscript. Institutional Review Board Statement Not applicable. Informed Consent Statement Not applicable. Data Availability Statement All data are available in the main text or the [127]supplementary materials. Conflicts of Interest The authors have no conflicts of interest to declare. Funding Statement This research was supported by a grant from KHIDI, funded by the Ministry of Health and Welfare, Republic of Korea (HI14C1324). Footnotes Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. References