MECOM/PRDM3 and PRDM16 Serve as Prognostic-Related Biomarkers and Are Correlated With Immune Cell Infiltration in Lung Adenocarcinoma

Background: The MDS1 and EVI1 complex locus (MECOM, also called PRDM3) and PR domain containing 16 (PRDM16) are two highly related zinc finger transcription factors associated with many malignancies. However, the mechanisms of MECOM and PRDM16 in prognosis and tumor immune infiltration in lung adenocarcinoma (LUAD) remain uncertain.

Methods: The Cancer Genome Atlas (TCGA), Oncomine, UALCAN, GEPIA, and TIMER databases were searched to determine the relationship between the expression of MECOM and PRDM16, clinicopathological features, immune infiltration, and prognosis in LUAD. Coexpressed genes of the two genes were investigated by CBioPortal, and the potential mechanism of MECOM- and PRDM16-related genes was elucidated by GO and KEGG analyses. STRING database was utilized to further construct the protein-protein interaction network of the coexpressed genes, and the hub genes were identified by Cytoscape. Finally, qRT-PCR was performed to identify the mRNA levels of the target genes in LUAD.

Results: mRNA levels of MECOM and PRDM16 were downregulated in LUAD (p < 0.05), and the low expression of the two genes was associated with the age, gender, smoking duration, tissue subtype, poor stage, nodal metastasis status, TP53 mutation, and prognosis in LUAD (p < 0.05). MECOM and PRDM16 were also found to be correlated with the expression of a variety of immune cell subsets and their markers. KEGG analysis showed that both of them were mainly enriched in the cell cycle, cellular senescence, DNA replication, and p53 signaling pathway. Importantly, the mRNA levels of the two genes were also found to be decreased in the clinical samples of LUAD by qRT-PCR.

Conclusion: MECOM and PRDM16 may serve as potential prognostic biomarkers which govern immune cell recruitment to LUAD.

Introduction

Lung cancer, a malignant disease, poses a serious threat to human health, especially in East Asia (1). Non-small cell lung cancer (NSCLC) accounts for approximately 85% of all lung carcinoma (2), and adenocarcinoma is the most prevalent subtype (3). In recent years, significant advances have been made in targeted therapy and immunotherapy in lung cancer treatment, some of which have replaced traditional chemoradiotherapy as the first-line treatment (4). Despite the advances in lung cancer-related treatment technology, the long-time survival remains poor, with a 5-year survival rate of less than 20% (5). Therefore, it is of great significance to investigate more potential molecular targets for its diagnosis, treatment, and prognostic evaluation.

The PR/SET domain (PRDM) gene family contains 19 members, which is structurally defined by the conserved N-terminal PR domain that relates to the SET methyltransferase domain and multiple zinc fingers that mediate sequence-specific DNA binding and protein-protein interactions (6, 7). PRDM proteins are involved in a spectrum of crucial cellular processes such as cell fate determination and development, and abnormal function of some members may contribute to malignant transformation (6–8). MDS1 and EVI1 complex locus (MECOM) encodes PRDM3 (also known as MDS1-EVI1). PR domain containing 16 (PRDM16, also called MEL1) was identified as a closely related paralog of PRDM3, which shares 53% sequence identity with the N-terminus of PRDM3 (9). Given the striking similarity between MECOM/PRDM3 and PRDM16, the two are thought to have a cooperative effect on certain biological functions, the occurrence, and the development of the disease. PRDM3 and PRDM16 were found to maintain the survival of normal hematopoietic stem cells with additive roles (10), and they were also demonstrated to have a direct interaction with the NuRD chromatin remodeling complex (11). Of note, MECOM and PRDM16 are the most mutated genes in the PRDM family in multiple human cancers (12, 13), which has been verified to promote the occurrence and progression of acute myeloid leukemia (14, 15). MECOM was found to be upregulated in well-differentiated pancreatic ductal adenocarcinoma (PDAC), which contributes to better survival in PDAC patients (16). By contrast, PRDM16 was observed with extensive epigenetic gene silencing and inhibits tumor growth by suppressing HIF-targeted gene semaphorin 5B in renal cancer. (17). PRDM16 was also found to function as a suppressor in lung cancer metastasis and downregulated in LUAD (18). The potential mechanism of MECOM and PRDM16, especially MECOM, in terms of tumor progression, prognosis, and immune cell infiltration in LUAD, however, remains unclear.

Herein, we analyzed the expression of MECOM and PRDM16 in LUAD and investigate the correlation between the two genes and clinicopathological features, prognosis, immune cell infiltration, and markers of immune cells in LUAD patients on the basis of various databases. Finally, the expression levels of MECOM and PRDM16 were further verified by Quantitative reverse transcription polymerase chain reaction (qRT-PCR) in LUAD patients.

Materials and Methods

Analysis of MECOM and PRDM16 Expression in Different Cancers

Oncomine (19) and TIMER databases (20) were utilized to explore the mRNA levels of MECOM and PRDM16 in tumor and adjacent normal tissues. Filter criteria of Oncomine: p < 1E−04; multiple changes = 2. Data type was mRNA. Gene ranking is top 10%.

Data Source and Processing

The gene expression data in HTSeq-FPKM of 535 LUAD and 59 normal lung tissues and the clinical information of 522 LUAD patients were obtained from The Cancer Genome Atlas (TCGA) database (21). A total of 57 paired LUAD and adjacent normal lung tissues among the 594 samples were identified in TCGA. Samples with incomplete parameters or lack of prognostic follow-up clinical data were excluded. The details are shown in Table S1.

Expression and Survival Analysis of MECOM and PRDM16 in LUAD

The downloaded data of TCGA were analyzed by the R “limma” and “survival” package to explore the expression of MECOM and PRDM16 in LUAD, normal tissues, and different clinical features, as well as the correlation between MECOM, PRDM16, and the overall survival (OS) in LUAD. Subsequently, the association between MECOM, PRDM16, clinicopathologic features, and disease-free survival (DFS) were further analyzed by UALCAN (22) and GEPIA (23). CBioPortal database (24) is an open-platform integrated data from multiple databases including TCGA, which was searched to investigate the correlation between MECOM and PRDM16 expression and its methylation levels.

Visualization, Expression, and Correlation Analysis of Immune Cells in LUAD

Immune cell expression in each LUAD sample of TCGA was obtained by the “CIBERSORT” R package, and the results were visualized by the “Barplot” and “Corrplot” R package. The expression of different immune cells between patients with high and low expression of MECOM and PRDM16 in LUAD were investigated by the “limma” package and visualized by the “Vioplot” package.

Correlation Analysis Between MECOM and PRDM16 Expression, Immune Cell Infiltration, and Immune Cell Biomarkers in LUAD

The correlations between MECOM and PRDM16 and the infiltration of the immune cells in LUAD were analyzed by the Spearman correlation analysis of the R package. Gene_Corr module of the TIMER database was performed to investigate the relationship between the expression of the two genes and biomarkers of immune cells.

Survival Analysis of MECOM and PRDM16 Expression and Immune Cell Infiltration in LUAD

The effect of gene expression and immune cell infiltration on survival of LUAD patients was investigated by the Outcome Module of the TIMER database.

Screening of MECOM- and PRDM16-Coexpressed Genes

The coexpression genes of PRDM members MECOM and PRDM16 were analyzed using the CBioPortal database (coexpression criteria: |r| > 0.4 and p < 0.001), and heatmaps were generated using the “heatmap” package of R. The intersection genes were then extracted from the MECOM- and PRDM16-coexpressed genes by the “VennDiagram” package of R.

GO and KEGG Analysis

The biological functions and signaling pathways affected by the coexpressed genes of MECOM and PRDM16 were investigated by GO and KEGG analyses by the “clusterProfiler” package of R (25).

Protein-Protein Interaction Network Construction and Hub Genes Identification

STRING database was performed to construct a protein-protein interaction (PPI) network of the coexpressed genes of MECOM and PRDM16 with an interaction score of >0.7. Subsequently, the top 10 hub genes were identified using the CytoHubba plugin of Cytoscape that uses the degree method. Finally, the expression of hub genes and the relationship between them, OS, and DFS in LUAD were investigated by the GEPIA database.

Patients and Tissue Samples

A total of 10 patients with pathologically confirmed primary LUAD in the First Affiliated Hospital of Xi’an Jiaotong University (Shaanxi, China) from July 2020 to Mar 2021 were collected in this study. All patients did not receive radiotherapy and/or chemotherapy and are without another cancer history.

Quantitative Reverse Transcription Polymerase Chain Reaction

Total RNA was extracted from lung tissues by RNA-easyIsolation Reagent (Vazyme #R701). SuperScript IV Reverse Transcriptase (cat. no. 18090010, Thermo Fisher, Waltham, MA, USA) was performed for reverse transcription, and SYBR^® Premix Ex Taq™ II (Tli RNaseH Plus) (cat. no. RR820Q, Takara, Kusatsu, Japan) was utilized for the qRT-PCR in accordance with the instructions. The primer sequences are shown in Table S2. The thermocycling program of qRT-PCR is shown below: initial denaturation for 30 s at 95°C, 40 cycles of denaturation for 5 s at 95°C, and 40 cycles of amplification for 30 s at 60°C. The mRNA expression of MECOM and PRDM16 in cells was represented as the 2^−ΔΔCt, and β-actin was used as the internal reference.

Statistical Analysis

Perl, R(4.0.1), and SPSS 18.0 were used for data processing and statistical analysis. Wilcoxon test and Mann-Whitney U test were performed to evaluate the expression of MECOM and PRDM16 and their association with the clinicopathological features of LUAD patients. Spearman correlation analysis was conducted to assess the correlations between MECOM and PRDM16 expression and the immune infiltration in LUAD. p < 0.05 was statistically significant.

Results

The Expression of MECOM and PRDM16 Were Downregulated in LUAD

As shown in Figure 1A, the results of Oncomine revealed that both MECOM and PRDM16 expression were downregulated in lung, breast, gastric, sarcoma, and kidney tumors compared with normal tissues (all p < 0.05). In addition, downregulation of MECOM and PRDM16 expression was also found in LUAD and lung squamous cell carcinomas (LUSC) in the TIMER database (Figures 1B, C, all p < 0.05). Further analysis of Oncomine elucidated that MECOM and PRDM16 were downexpressed in LUAD patients with a fold change of −2.081 and −3.069, respectively (Table S3). In addition, the expression of MECOM and PRDM16 in LUAD tissues were found to be significantly lower in the normal tissues by evaluating the 594 samples and 57 paired LUAD and adjacent noncancerous tissues in the TCGA database (Figures 2A, B, D, E, p < 0.05), and the low expression of the two genes were supposed to be related to the hypermethylation level of its DNA by CBioPortal analysis (Figures 2C, F).

FIGURE 1

Figure 1 Expression of MECOM and PRDM16 in various tumors. (A) Expression of MECOM and PRDM16 in various tumors in the Oncomine database. Note: red and blue represent upregulation and downregulation of the target gene, respectively. (B, C) Expression of MECOM (B) and PRDM16 (C) in various tumors in the TIMER database. ^*p < 0.05; ^**p < 0.01; ^***p < 0.001. MECOM, MDS1 and EVI1 complex locus, also called PRDM3; PRDM16, PR domain containing 16, also called MEL1.

FIGURE 2

Figure 2 Expression of MECOM and PRDM16 in LUAD. (A, D) mRNA levels of MECOM (A) and PRDM16 (D) in LUAD tissues in TCGA. (B, E) MECOM (B) and PRDM16 (E) mRNA expression levels in paired tumor and adjacent normal tissues of 57 patients with LUAD in TCGA. (C, F) Correlation scatter plot of MECOM (C) and PRDM16 (F) mRNA expression and its methylation levels in LUAD in the cBioPortal database. MECOM, MDS1 and EVI1 complex locus, also called PRDM3; PRDM16, PR domain containing 16, also called MEL1; LUAD, lung adenocarcinoma.

MECOM and PRDM16 Were Correlated With Clinicopathological Features and Prognosis in LUAD Patients

The results showed that low expression of MECOM and PRDM16 were correlated with poor pathological stage, high lymph node status, and poor OS in LUAD (Figures 3A, B and 4A, C). Whereas, the expression of MECOM and PRDM16 had no effect on the DFS of LUAD (Figures 4B, D). In addition, MECOM and PRDM16 were also found to be associated with age, gender, race, stage, smoking status, tissue subtypes, TP53 mutation, and lymph node metastasis status in the UALCAN database (Tables S4 and S5, all p < 0.05).

FIGURE 3

Figure 3 Expression levels of MECOM (A) and PRDM16 (B) in different clinicopathological characteristics of LUAD patients. MECOM, MDS1 and EVI1 complex locus, also called PRDM3; PRDM16, PR domain containing 16, also called MEL1; LUAD, lung adenocarcinoma.

FIGURE 4

Figure 4 Correlation between the expression of MECOM and PRDM16 and prognosis in LUAD. (A, B) Correlation between MECOM, OS (A), and DFS (B) in LUAD. (C, D) Correlation between PRDM16, OS (C), and DFS (D) in LUAD. MECOM, MDS1 and EVI1 complex locus, also called PRDM3; PRDM16, PR domain containing 16, also called MEL1; LUAD, lung adenocarcinoma.

Expression of MECOM and PRDM16 Correlated With Immune Cell Infiltration in LUAD

We visualized the expression ratio of the immune cells in each LUAD sample from TCGA database and analyzed the correlation between different immune cells (Figures S1A, B). Some immune cells were differentially expressed between LUAD patients with high and low expression of MECOM and PRDM16 (Figures 5A, B). In addition, a significant positive relationship was found between MECOM and naive B cells, resting mast cells, activated NK cells, activated dendritic cells, plasma cells, and resting memory CD4 T cells (Figure 6A; Figure S2A; p < 0.05). Macrophages M0 and M1, activated mast cells, resting NK cells, activated memory CD4 T cells, and Tregs were proved to be negatively correlated with MECOM (Figure 6A). As shown in Figure 6B and Figure S2B, there was a positive relationship between PRDM16 and monocytes, resting mast cells, resting and activated dendritic cells, activated NK cells, resting memory CD4 T cells, and Tregs. Whereas, macrophages M0 and M1, resting NK cells, activated memory CD4 T cells, eosinophils, and CD8 T cells were negatively correlated with PRDM16 (p < 0.05).

FIGURE 5

Figure 5 Expression of immune cells in LUAD patients with high and low expression of MECOM (A) and PRDM16 (B). MECOM, MDS1 and EVI1 complex locus, also called PRDM3; PRDM16, PR domain containing 16, also called MEL1; LUAD, lung adenocarcinoma.

FIGURE 6

Figure 6 Correlation between MECOM (A), PRDM16 (B) expression, and the level of immune cell infiltration in LUAD. MECOM, MDS1 and EVI1 complex locus, also called PRDM3; PRDM16, PR domain containing 16, also called MEL1; LUAD, lung adenocarcinoma.

MECOM- and PRDM16-Related Immune Cells Correlated With Prognosis in LUAD

Kaplan-Meier plots were constructed by the TIMER database to investigate the relationship between the expression of MECOM and PRDM16 and their related immune cells in the prognosis of LUAD. Figure 7A reveals that patients with high MECOM expression and high B-cell, mast cell, and myeloid dendritic cell levels were correlated with better prognosis (B cells: HR = 0.548, p = 0.015; mast cells: HR = 0.551, p = 0.015; myeloid dendritic cells: HR = 0.552, p = 0.020). However, patients with low levels of MECOM and high CD4⁺ T-cell levels seem to have a worse prognosis (HR = 1.720, p = 0.009). Among patient with high levels of PRDM16 and Tregs, mast cells and myeloid dendritic cells correlated with favorable outcomes (Figure 7B, Tregs: HR = 0.568, p = 0.047; mast cells: HR = 0.549, p = 0.026; myeloid dendritic cells: HR = 0.404, p = 0.001), and low expression of PRDM16 and high CD4⁺ T-cell levels was found to be associated with poor outcomes in LUAD (HR = 1.730, p =0.005).

FIGURE 7

Figure 7 Correlation between MECOM- (A) and PRDM16-related (B) immune cell expression and prognosis in LUAD. MECOM, MDS1 and EVI1 complex locus, also called PRDM3; PRDM16, PR domain containing 16, also called MEL1; LUAD, lung adenocarcinoma.

The Correlation Between MECOM, PRDM16, and Immune Cell Marker Expression

TIMER database was searched to assess the correlation between the expression of MECOM and PRDM16 and the levels of particular immune cell subset markers. As shown in Table 1, the results adjusted based on tumor purity revealed a significant correlation between MECOM and markers of monocyte (CD86), M1 macrophage (NOS2), M2 macrophage (COX2) neutrophils (CD66b, CCR7), NK (KIR2DL4), DC (HLA-DPB1, CD1C), Th1 (STAT1, IFNG), Th2 (STAT6, STAT5A), Tfh (BCL6), Th17 (STAT3), Treg (STAT5B), and T-cell exhaustion (PD-1, LAG3, HAVCR2, GZMB) in LUAD. PRDM16 was also demonstrated to be closely related to the markers of some immune cells including CD8⁺ T, monocyte, M2 macrophage, neutrophils, DC, NK, Th1, Th2, Tfh, Th17, Treg, and T-cell exhaustion in LUAD (Table 1).

TABLE 1

Table 1 Correlation analysis between MECOM, PRDM16, and gene markers of immune cells in TIMER.

Identification of MECOM and PRDM16 Coexpressed Genes

In the CBioPortal database, 103 and 156 genes were positively and negatively correlated with MECOM, respectively. In addition, 215 positively and 276 negatively correlated genes with PRDM16 were also identified. The heatmap (Figures 8A, B) showed the top 15 positively and negatively correlated genes with MECOM and PRDM16. A total of 137 intersections of the two coexpressed genes were further explored by the Venn diagram (Figure 8C).

FIGURE 8

Figure 8 The top 15 coexpression genes of MECOM (A) and PRDM16 (B) in LUAD in the cBioPortal database. (C) Intersection coexpression genes of MECOM and PRDM16. Note: Genes with |r| > 0.4 and p < 0.001 were considered coexpression genes. MECOM, MDS1 and EVI1 complex locus, also called PRDM3; PRDM16, PR domain containing 16, also called MEL1; LUAD, lung adenocarcinoma.

The Pathways Regulated by MECOM and PRDM16 in LUAD Were Revealed by GO and KEGG Analyses

GO analysis elucidated that MECOM- and PRDM16-coexpressed genes were mainly involved in mitosis, DNA replication, and organelle division (Figures 9A, C, E). Coexpressed genes were also found enriched in regulating cell cycle, cellular senescence, DNA replication, and p53 signaling pathway by KEGG analysis (Figures 9B, D, F).

FIGURE 9

Figure 9 GO and KEGG analyses of MECOM (A, B) and PRDM16 (C, D) coexpression genes and intersection of coexpression genes (E, F). MECOM, MDS1 and EVI1 complex locus, also called PRDM3; PRDM16, PR domain containing 16, also called MEL1; LUAD, lung adenocarcinoma.

MECOM- and PRDM16-Related Hub Genes in PPI Network

A total of 10 hub genes (CDK1, CDC20, BUB1, CCNA2, CCNB2, AURKB, CCNB1, KIF2C, CDCA8, TOP2A) identified by PPI network analysis (Figures 10A, B; Table 2) were observed to be elevated in LUAD and related to poor OS by the GEPIA database (Figures 11A–J and 12A–J). Additionally, CDK1, AURKB, BUB1, CCNB1, and CCNB2 were also demonstrated to correlate with the DFS of LUAD patients (Figure S3).

FIGURE 10

Figure 10 PPI network constructed by MECOM and PRDM16 intersection coexpression genes. (A) PPI network. (B) Ten hub genes in the PPI network. The top 10 hub genes are shown in different colors except blue. MECOM, MDS1 and EVI1 complex locus, also called PRDM3; PRDM16, PR domain containing 16, also called MEL1; LUAD, lung adenocarcinoma.

TABLE 2

Table 2 The top 10 hub genes in the protein-protein interaction (PPI) network ranked by degree method.

FIGURE 11

Figure 11 Expression of the top 10 hub genes in LUAD in the GEPIA database. The expressions of CDK1 (A), CDC20 (B), BUB1 (C), CCNA2 (D), CCNB2 (E), AURKB (F), CCNB1 (G), KIF2C (H), CDCA8 (I), and TOP2A (J) in LUAD. LUAD, lung adenocarcinoma. *p < 0.05.

FIGURE 12

Figure 12 The correlation between the top 10 hub genes and the OS of LUAD patients in the GEPIA database. CDK1 (A), CDC20 (B), BUB1 (C), CCNA2 (D), CCNB2 (E), AURKB (F), CCNB1 (G), KIF2C (H), CDCA8 (I), and TOP2A (J). LUAD, lung adenocarcinoma.

Verification of MECOM and PRDM16 Expression in LUAD Tissue Samples

The mRNA expression of MECOM and PRDM16 was evaluated in LUAD patients’ tissue by qRT-PCR. As shown in Figures 13A–D, both MECOM and PRDM16 were found to be downregulated in primary tumorous tissues of LUAD patients compared with the adjacent noncancerous tissues (p < 0.05), which is consistent with the results of the bioinformatic analysis.

FIGURE 13

Figure 13 The mRNA expression of MECOM (A, B) and PRDM16 (C, D) in clinical samples of LUAD. MECOM, MDS1 and EVI1 complex locus, also called PRDM3; PRDM16, PR domain containing 16, also called MEL1; LUAD, lung adenocarcinoma. **p < 0.01; ***p < 0.001.

Discussion

The PRDM family was illustrated to participate in a series of important cellular processes, and the abnormal function of some members may cause malignant transformation. In fact, we systematically investigate the expression of all members of the PRDM family and the relationship between them and prognosis in LUAD at the early stage of the study based on the above background. We found that only two members, MECOM and PRDM16, were differentially expressed and correlated with the prognosis in LUAD. Many kinds of literature further confirm that MECOM and PRDM16 are two highly correlated transcription factors and play an accumulative role in some biological functions and the pathogenesis of the disease (10, 11). MECOM encodes PRDM3, which is involved in the process of hematopoiesis, apoptosis, development, differentiation, and proliferation (26–29) and considered an oncoprotein in the hematopoietic system and has been found to promote leukemogenesis by overexpression of the AML1 gene through the translocation of AML1 (30). Whereas, MECOM has been found to function as both oncogene and tumor suppressor gene (TSG) in some solid tumors (13), which was observed to be overexpressed in colon and breast cancer specimens and correlated with worse outcomes (31, 32), by contrast, the favorable prognosis of ovarian cancer patients were identified to be associated with amplification of MECOM gene (33). PRDM16 is highly homologous to MECOM, and they have been found to be together involved in the maintenance of hematopoietic stem cell function and heterochromatin integrity (10, 34). Additionally, amplification of PRDM16 was found in leukemia (35), osteosarcoma (36), and gastric cancer (37), and reduced expression of PRM16 was observed in uterine leiomyosarcoma (38) and NSCLC (18, 39), which elucidated PRDM16 and also functions as an oncogene or TSG in several cancers as MECOM. Together, these data suggest that alterations of MECOM and PRDM16 are possibly involved in tumorigenesis by either upregulation or downregulation. Although research demonstrated that PRDM16 can inhibit metastasis of LUAD (18), the relationship between PRDM16, prognosis, and immune infiltration in LUAD were not further evaluated. In addition, the role of MECOM in LUAD has not been studied.

In the present study, we found that MECOM and PRDM16 were downregulated and related to age, gender, race, smoking status, pathological stage, N stage, TP53 mutation status, and poor OS in LUAD. Previous studies also found that MECOM and PRDM16 were identified to be involved in the development and drug resistance in leukemia by coupling with p53 (40, 41). In addition, the coexpressed genes of the MECOM and PRDM16 were further found to be enriched in regulating cell cycle, cellular senescence, DNA replication, and p53 signaling pathway; these are the classic pathways involved in tumorigenesis.

The tumor immune microenvironment (TIME), including the composition, distribution, and expression of various biomarkers of tumor immune cells, plays a crucial role in the progression and treatment response of LUAD (42). Hence, immunotherapy has become a new direction of tumor therapy. Herein, the correlation between the MECOM and PRDM16 expression and the infiltration of the immune cells were further investigated to evaluate whether the two genes influence the prognosis by participating in the immune response to LUAD. We found that MECOM and PRDM16 were closely associated with multiple immune cell infiltration such as macrophages, mast cells, T cells, B cells, and NK cells. The expression of MECOM, PRDM16, and immune cells together affect the prognosis of patients with LUAD. Similar to our findings, PRDM3 decreases pancreatic tumorigenesis by regulating immune cell activation and infiltration (16). Whereas, the function of PRDM16 in the immune response of cancers was not systematically investigated, which needs to be further studied.

We also identified 10 hub genes (CDK1, CDC20, BUB1, CCNA2, CCNB2, AURKB, CCNB1, KIF2C, CDCA8, TOP2A) that were overexpressed in LUAD and associated with poor survival, among which, CDK1, AURKB, BUB1, CCNB2, and CCNB were also correlated with the DFS of LUAD patients. Notably, we further validated the expression of MECOM and PRDM16 in clinical specimens, which is consistent with the results of the bioinformatics analysis.

The advantage of this study is that the expression of the two highly related PRDM members, MECOM and PRDM16, and the correlation between them, prognosis, and immune cell infiltration in the LUAD were explored for the first time. Importantly, the results were further verified in the clinical samples. Whereas, there are still some limitations: (i) the sample size of the clinical sample is insufficient and (ii) the molecular biology experiment was not conducted to investigate the function of the two genes in LUAD. Considering the importance of the two genes in the prognosis and immune cell infiltration of LUAD, it is worthwhile to further study the mechanism of the two genes in regulating the occurrence and development of LUAD in the future.

Conclusion

These findings together support the notion that MECOM and PRDM16 are potential immune-related and prognostic biomarkers for LUAD, and the mechanism of MECOM and PRDM16 in LUAD deserves in-depth investigation by molecular biology experiment.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

The studies involving human participants were reviewed and approved by the Ethics Committee of the First Affiliated Hospital of Xi’an Jiaotong University. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

Conception and design: PS and MC. Administrative support: PS and MC. Provision of study materials or patients: ML, HR, YZ, NL, MF, KW, and TY. Experiment operations, collection, and assembly of data: ML, HR, YZ, NL, MF, KW, and TY. Data analysis and interpretation: ML, HR, YZ, NL, MF, KW, TY, MC, and PS. Manuscript writing: all authors. Final approval of manuscript: all authors. All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.

Funding

The research was funded by the following projects: the Key research and development projects of Shaanxi Province (No. 2018KW-039; No. 2021ZDLSF02-05; No. 2019-KW-034); the National Natural Science Foundation of China (No. 81802291; No. 82103467); Fundamental Research Funds for the Central Universities in Xi’an Jiaotong University (No. 1191320149; No. xzy012021064); and China Postdoctoral Science Foundation Grant (No. 2019M653666).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2022.772686/full#supplementary-material

Supplementary Figure S1 | Proportion of immune cells (A) and correlation between immune cells (B) in LUAD samples in TCGA database. LUAD, Lung adenocarcinoma.

Supplementary Figure S2 | Correlation between MECOM (A), PRDM16 (B) expression and the level of immune cells infiltration in LUAD. MECOM, MDS1 and EVI1 complex locus, also called PRDM3; PRDM16, PR domain containing 16, also called MEL1; LUAD, Lung adenocarcinoma.

Supplementary Figure S3 | The correlations of the top 10 Hub genes with DFS of LUAD patients in the GEPIA database. CDK1 (A), CDC20 (B), BUB1 (C), CCNA2 (D), CCNB2 (E), AURKB (F), CCNB1 (G), KIF2C (H), CDCA8 (I) and TOP2A (J). DFS, Disease-free survival.

References

1. Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global Cancer Statistics 2018: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J Clin (2018) 68(6):394–424. doi: 10.3322/caac.21492

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Duma N, Santana-Davila R, Molina JR. Non-Small Cell Lung Cancer: Epidemiology, Screening, Diagnosis, and Treatment. Mayo Clin Proc (2019) 94(8):1623–40. doi: 10.1016/j.mayocp.2019.01.013

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Chen Z, Fillmore CM, Hammerman PS, Kim CF, Wong KK. Non-Small-Cell Lung Cancers: A Heterogeneous Set of Diseases. Nat Rev Cancer (2014) 14(8):535–46. doi: 10.1038/nrc3775

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Anagnostou VK, Brahmer JR. Cancer Immunotherapy: A Future Paradigm Shift in the Treatment of Non-Small Cell Lung Cancer. Clin Cancer Res (2015) 21(5):976–84. doi: 10.1158/1078-0432.Ccr-14-1187

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Lu T, Yang X, Huang Y, Zhao M, Li M, Ma K, et al. Trends in the Incidence, Treatment, and Survival of Patients With Lung Cancer in the Last Four Decades. Cancer Manag Res (2019) 11:943–53. doi: 10.2147/cmar.S187317

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Fog CK, Galli GG, Lund AH. PRDM Proteins: Important Players in Differentiation and Disease. Bioessays (2012) 34(1):50–60. doi: 10.1002/bies.201100107

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Hohenauer T, Moore AW. The Prdm Family: Expanding Roles in Stem Cells and Development. Development (2012) 139(13):2267–82. doi: 10.1242/dev.070110

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Casamassimi A, Rienzo M, Di Zazzo E, Sorrentino A, Fiore D, Proto MC, et al. Multifaceted Role of PRDM Proteins in Human Cancer. Int J Mol Sci (2020) 21(7):2648. doi: 10.3390/ijms21072648

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Yoshida M, Nosaka K, Yasunaga J, Nishikata I, Morishita K, Matsuoka M. Aberrant Expression of the MEL1S Gene Identified in Association With Hypomethylation in Adult T-Cell Leukemia Cells. Blood (2004) 103(7):2753–60. doi: 10.1182/blood-2003-07-2482

PubMed Abstract | CrossRef Full Text | Google Scholar

10. McGlynn KA, Sun R, Vonica A, Rudzinskas S, Zhang Y, Perkins AS. Prdm3 and Prdm16 Cooperatively Maintain Hematopoiesis and Clonogenic Potential. Exp Hematol (2020) 85:20–32.e23. doi: 10.1016/j.exphem.2020.04.010

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Ivanochko D, Halabelian L, Henderson E, Savitsky P, Jain H, Marcon E, et al. Direct Interaction Between the PRDM3 and PRDM16 Tumor Suppressors and the NuRD Chromatin Remodeling Complex. Nucleic Acids Res (2019) 47(3):1225–38. doi: 10.1093/nar/gky1192

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Mzoughi S, Tan YX, Low D, Guccione E. The Role of PRDMs in Cancer: One Family, Two Sides. Curr Opin Genet Dev (2016) 36:83–91. doi: 10.1016/j.gde.2016.03.009

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Sorrentino A, Federico A, Rienzo M, Gazzerro P, Bifulco M, Ciccodicola A, et al. PR/SET Domain Family and Cancer: Novel Insights From the Cancer Genome Atlas. Int J Mol Sci (2018) 19(10):3250. doi: 10.3390/ijms19103250

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Zhang Y, Owens K, Hatem L, Glass CH, Karuppaiah K, Camargo F, et al. Essential Role of PR-Domain Protein MDS1-EVI1 in MLL-AF9 Leukemia. Blood (2013) 122(16):2888–92. doi: 10.1182/blood-2012-08-453662

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Yamato G, Yamaguchi H, Handa H, Shiba N, Kawamura M, Wakita S, et al. Clinical Features and Prognostic Impact of PRDM16 Expression in Adult Acute Myeloid Leukemia. Genes Chromosomes Cancer (2017) 56(11):800–9. doi: 10.1002/gcc.22483

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Ye J, Huang A, Wang H, Zhang AMY, Huang X, Lan Q, et al. PRDM3 Attenuates Pancreatitis and Pancreatic Tumorigenesis by Regulating Inflammatory Response. Cell Death Dis (2020) 11(3):187. doi: 10.1038/s41419-020-2371-x

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Kundu A, Nam H, Shelar S, Chandrashekar DS, Brinkley G, Karki S, et al. PRDM16 Suppresses HIF-Targeted Gene Expression in Kidney Cancer. J Exp Med (2020) 217(6):e20191005. doi: 10.1084/jem.20191005

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Fei LR, Huang WJ, Wang Y, Lei L, Li ZH, Zheng YW, et al. PRDM16 Functions as a Suppressor of Lung Adenocarcinoma Metastasis. J Exp Clin Cancer Res (2019) 38(1):35. doi: 10.1186/s13046-019-1042-1

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Rhodes DR, Kalyana-Sundaram S, Mahavisno V, Varambally R, Yu J, Briggs BB, et al. Oncomine 3.0: Genes, Pathways, and Networks in a Collection of 18,000 Cancer Gene Expression Profiles. Neoplasia (2007) 9(2):166–80. doi: 10.1593/neo.07112

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Li T, Fan J, Wang B, Traugh N, Chen Q, Liu JS, et al. TIMER: A Web Server for Comprehensive Analysis of Tumor-Infiltrating Immune Cells. Cancer Res (2017) 77(21):e108–10. doi: 10.1158/0008-5472.Can-17-0307

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Hutter C, Zenklusen JC. The Cancer Genome Atlas: Creating Lasting Value Beyond Its Data. Cell (2018) 173(2):283–5. doi: 10.1016/j.cell.2018.03.042

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Chandrashekar DS, Bashel B, Balasubramanya SAH, Creighton CJ, Ponce-Rodriguez I, Chakravarthi B, et al. UALCAN: A Portal for Facilitating Tumor Subgroup Gene Expression and Survival Analyses. Neoplasia (2017) 19(8):649–58. doi: 10.1016/j.neo.2017.05.002

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z. GEPIA: A Web Server for Cancer and Normal Gene Expression Profiling and Interactive Analyses. Nucleic Acids Res (2017) 45(W1):W98–102. doi: 10.1093/nar/gkx247

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Cerami E, Gao J, Dogrusoz U, Gross BE, Sumer SO, Aksoy BA, et al. The Cbio Cancer Genomics Portal: An Open Platform for Exploring Multidimensional Cancer Genomics Data. Cancer Discov (2012) 2(5):401–4. doi: 10.1158/2159-8290.Cd-12-0095

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Zhou Y, Zhou B, Pache L, Chang M, Khodabakhshi AH, Tanaseichuk O, et al. Metascape Provides a Biologist-Oriented Resource for the Analysis of Systems-Level Datasets. Nat Commun (2019) 10(1):1523. doi: 10.1038/s41467-019-09234-6

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Hoyt PR, Bartholomew C, Davis AJ, Yutzey K, Gamer LW, Potter SS, et al. The Evi1 Proto-Oncogene Is Required at Midgestation for Neural, Heart, and Paraxial Mesenchyme Development. Mech Dev (1997) 65(1-2):55–70. doi: 10.1016/s0925-4773(97)00057-9

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Kurokawa M, Mitani K, Yamagata T, Takahashi T, Izutsu K, Ogawa S, et al. The Evi-1 Oncoprotein Inhibits C-Jun N-Terminal Kinase and Prevents Stress-Induced Cell Death. EMBO J (2000) 19(12):2958–68. doi: 10.1093/emboj/19.12.2958

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Wieser R. The Oncogene and Developmental Regulator EVI1: Expression, Biochemical Properties, and Biological Functions. Gene (2007) 396(2):346–57. doi: 10.1016/j.gene.2007.04.012

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Karakaya K, Herbst F, Ball C, Glimm H, Krämer A, Löffler H. Overexpression of EVI1 Interferes With Cytokinesis and Leads to Accumulation of Cells With Supernumerary Centrosomes in G0/1 Phase. Cell Cycle (2012) 11(18):3492–503. doi: 10.4161/cc.21801

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Mitani K. Molecular Mechanisms of Leukemogenesis by AML1/EVI-1. Oncogene (2004) 23(24):4263–9. doi: 10.1038/sj.onc.1207777

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Choi EJ, Kim MS, Song SY, Yoo NJ, Lee SH. Intratumoral Heterogeneity of Frameshift Mutations in MECOM Gene Is Frequent in Colorectal Cancers With High Microsatellite Instability. Pathol Oncol Res (2017) 23(1):145–9. doi: 10.1007/s12253-016-0112-3

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Wang H, Schaefer T, Konantz M, Braun M, Varga Z, Paczulla AM, et al. Prominent Oncogenic Roles of EVI1 in Breast Carcinoma. Cancer Res (2017) 77(8):2148–60. doi: 10.1158/0008-5472.Can-16-0593

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Nanjundan M, Nakayama Y, Cheng KW, Lahad J, Liu J, Lu K, et al. Amplification of MDS1/EVI1 and EVI1, Located in the 3q26.2 Amplicon, Is Associated With Favorable Patient Prognosis in Ovarian Cancer. Cancer Res (2007) 67(7):3074–84. doi: 10.1158/0008-5472.Can-06-2366

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Pinheiro I, Margueron R, Shukeir N, Eisold M, Fritzsch C, Richter FM, et al. Prdm3 and Prdm16 Are H3K9me1 Methyltransferases Required for Mammalian Heterochromatin Integrity. Cell (2012) 150(5):948–60. doi: 10.1016/j.cell.2012.06.048

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Mochizuki N, Shimizu S, Nagasawa T, Tanaka H, Taniwaki M, Yokota J, et al. A Novel Gene, MEL1, Mapped to 1p36.3 is Highly Homologous to the MDS1/EVI1 Gene and Is Transcriptionally Activated in T(1;3)(P36;Q21)-Positive Leukemia Cells. Blood (2000) 96(9):3209–14.

PubMed Abstract | Google Scholar

36. Man TK, Lu XY, Jaeweon K, Perlaky L, Harris CP, Shah S, et al. Genome-Wide Array Comparative Genomic Hybridization Analysis Reveals Distinct Amplifications in Osteosarcoma. BMC Cancer (2004) 4:45. doi: 10.1186/1471-2407-4-45

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Takahata M, Inoue Y, Tsuda H, Imoto I, Koinuma D, Hayashi M, et al. SKI and MEL1 Cooperate to Inhibit Transforming Growth Factor-Beta Signal in Gastric Cancer Cells. J Biol Chem (2009) 284(5):3334–44. doi: 10.1074/jbc.M808989200

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Cuppens T, Moisse M, Depreeuw J, Annibali D, Colas E, Gil-Moreno A, et al. Integrated Genome Analysis of Uterine Leiomyosarcoma to Identify Novel Driver Genes and Targetable Pathways. Int J Cancer (2018) 142(6):1230–43. doi: 10.1002/ijc.31129

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Tan SX, Hu RC, Xia Q, Tan YL, Liu JJ, Gan GX, et al. The Methylation Profiles of PRDM Promoters in Non-Small Cell Lung Cancer. Onco Targets Ther (2018) 11:2991–3002. doi: 10.2147/ott.S156775

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Shing DC, Trubia M, Marchesi F, Radaelli E, Belloni E, Tapinassi C, et al. Overexpression of Sprdm16 Coupled With Loss of P53 Induces Myeloid Leukemias in Mice. J Clin Invest (2007) 117(12):3696–707. doi: 10.1172/jci32390

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Cai SF, Chu SH, Goldberg AD, Parvin S, Koche RP, Glass JL, et al. Leukemia Cell of Origin Influences Apoptotic Priming and Sensitivity to LSD1 Inhibition. Cancer Discovery (2020) 10(10):1500–13. doi: 10.1158/2159-8290.Cd-19-1469

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Chi A, He X, Hou L, Nguyen NP, Zhu G, Cameron RB, et al. Classification of Non-Small Cell Lung Cancer's Tumor Immune Micro-Environment and Strategies to Augment Its Response to Immune Checkpoint Blockade. Cancers (Basel) (2021) 13(12):2924. doi: 10.3390/cancers13122924

PubMed Abstract | CrossRef Full Text | Google Scholar