Abstract

Background. Myelin and lymphocyte, T cell differentiation protein 2 (MAL2) is highly expressed in various cancers and associated with the development and prognosis of cancer. However, the relationship between MAL2 and breast cancer requires further investigation. This study aimed to explore the prognostic significance of MAL2 in breast cancer. Methods. MAL2 expression was initially assessed using the Oncomine database and The Cancer Genome Atlas (TCGA) database and verified by quantitative real-time polymerase chain reaction (RT-qPCR). The chi-square test or Fisher’s exact test was used to explore the association between clinical characteristics and MAL2 expression. The prognostic value of MAL2 in breast cancer was assessed by the Kaplan–Meier method and Cox regression analysis. Gene set enrichment analysis (GSEA) was performed to identify the biological pathways correlated with MAL2 expression in breast cancer. Besides, a single-sample GSEA (ssGSEA) was used to assess the relationship between the level of immune infiltration and MAL2 in breast cancer. Results. Both bioinformatics and RT-qPCR results showed that MAL2 was expressed at high levels in breast cancer tissues compared with the adjacent tissues. The chi-square test or Fisher’s exact test indicated that MAL2 expression was related to stage, M classification, and vital status. Kaplan–Meier curves implicated that high MAL2 expression was significantly associated with the poor prognosis. Cox regression models showed that high MAL2 expression could be an independent risk factor for breast cancer. GSEA showed that 14 signaling pathways were enriched in the high-MAL2-expression group. Besides, the MAL2 expression level negatively correlated with infiltrating levels of eosinophils and plasmacytoid dendritic cells in breast cancer. Conclusion. Overexpression of MAL2 correlates with poor prognosis and lower immune infiltrating levels of eosinophils and plasmacytoid dendritic cells in breast cancer and may become a biomarker for breast cancer prognosis.

1. Introduction

Among women, breast cancer is the most commonly diagnosed cancer and the main cause of cancer death [1]. The incidence of breast cancer has increased every year [2]. Nearly 3.8 million women in the United States have been diagnosed with breast cancer, including 268,600 new cases in 2019 [3].

Despite the advances in surgery, chemotherapy, radiotherapy, and endocrine therapy for breast cancer, more than 20% of patients still develop metastatic disease with a poor prognosis [4]. Since the molecular mechanisms of breast cancer remain unclear, the identification of novel prognostic biomarkers for breast cancer is necessary, which could contribute to the early detection and treatment of breast cancer [5]. MAL2, as an essential member of the MAL proteolipid family, is a four-pass membrane protein consisting of 176 amino acid residues. By encoding a multi-transmembrane protein, MAL2 mainly participates in endocytosis under physiological conditions and mediates the transport of substances between cells [6]. MAL2 is located in chromosome 8q23, an area where the copy number often increases in various types of cancer [7]. Previous studies have confirmed the increased expression of MAL2 in ovarian cancer [6], pancreatic cancer [8], thyroid cancer [9], and colorectal cancer [10]. Moreover, MAL2 expression was associated with pancreatic cancer and colorectal cancer overall survival [8], which suggested that MAL2 might be an important molecule involved in the progression and prognosis of tumors. Besides, MAL2 was recently identified as a breast cancer immunology target. Reduction of MAL2 in breast tumor cells can enhance CD8+ T cell-mediated cytotoxicity and inhibit the growth of breast tumors [11]. Although some previous studies suggested that MAL2 was overexpressed in breast cancer [7, 12, 13], the prognostic value of the MAL2 expression and its correlation with clinical features and immunotherapy in breast cancer requires further investigation, and whether MAL2 could be a specific marker for breast cancer still needs to be elucidated. In the present study, we applied a comprehensive strategy to uncover the importance of MAL2 in breast cancer.

2. Materials and Methods

2.1. Oncomine Database

The Oncomine database [14] (http://www.oncomine.org), a gene chip-based database and data mining platform, served to analyze the expression of MAL2 in various types of cancers. The filter conditions set were as follows: gene, MAL2; cancer type, breast cancer; differential analysis, cancer vs normal analysis; and data type, mRNA. Besides, we selected value = 1E–4, twofold change, and top 10% gene rank as the threshold [15]. All statistical methods and results were obtained from Oncomine.

2.2. The Cancer Genome Atlas

The Cancer Genome Atlas (TCGA) is a landmark cancer genomics program that comprised over 20,000 primary cancer data and matched normal samples spanning 33 cancer types. Raw counts of RNA-sequencing data and corresponding clinical information were obtained from The Cancer Genome Atlas (TCGA) data set (https://www.cancer.gov/tcga).

2.3. The Verification of MAL2 by qRT-PCR

Tumor tissues and paired adjacent tissues were collected from patients diagnosed with breast cancer at the First Affiliated Hospital of Guangzhou University of Chinese Medicine between 2019 and 2020. None of these patients received therapy before surgery. Specimens from surgery were stored at −80°C. Study approval was obtained from the Research Ethics Committee of the First Affiliated Hospital of Guangzhou University of Chinese Medicine, and patients signed informed consent forms before operation.

Total RNA was extracted from 32 pairs of breast cancer tissues and adjacent nontumor tissues frozen in liquid nitrogen using TRIzol reagent (Invitrogen, United States). The isolated RNA was measured at 260/280 nm using a NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, USA). Evo M-MLV RT Premix for qPCR (AG11706, China) was used to reverse-transcribe RNA into cDNA. According to the manufacturer’s instructions, cDNAs were amplified using the SYBR Green Premix Pro Taq HS qPCR Kit (AG11701, China). β-actin was used as an internal control for mRNA expression. The primer sequences were as follows: MAL2 Forward: 5′-ACGTAGCAGCCTCAATTTTTGC-3′ and Reverse: 5′-CATCTTCGTAAAGCCAGACCC-3′; β-actin, forward: 5′-TGGCACCCAGCACAATGAA-3′ and reverse: 5′-CTAAGTCATAGTCCGCCTAGAAGCA-3′. Each sample was carried out three times, and data were calculated using the 2−ΔΔCt (2 to the power of minus Delta Delta CT) method.

2.4. Statistical Analyses

Data analysis was performed using SPSS software 26.0 (IBM Corporation, Armonk, NY, USA) and R version 3.6.3 [16]. The ggplot2 package [17] was used to draw boxplots of clinical features according to MAL2 expression variation. The chi-square test or Fisher’s exact test was used to explore the association between clinical characteristics and MAL2 expression. The Wilcoxon signed rank sum test and Kruskal–Wallis test were utilized to measure the differential expression of MAL2 in the subgroup, including age; tumor (T), node (N), and metastasis (M) classification; estrogen receptor (ER) status; progesterone receptor (PR) status; human epidermal growth factor receptor-2 (HER2) status; and vital status. The impact of MAL2 expression levels on the overall survival (OS) of patients with breast cancer was analyzed by Kaplan–Meier curves using an R package named survminer [16].

Since we aimed to analyze the relationship of MAL2 with clinical features and its significance on the survival of patients with breast cancer, the median expression value of MAL2 was applied as a cut-off value for further analysis according to the previous studies [1821]. Therefore, the median value of MAL2 expression was utilized as a cut-off value to divide patients with complete clinical data from the TCGA database into the high-MAL2-expression group and low-MAL2-expression group.

Then, the univariate and multivariate Cox analyses were performed to determine the related variables.

2.5. Gene Set Enrichment Analysis

Gene set enrichment analysis (GSEA) [22] was conducted to determine the biological processes activated in the low-MAL2-expression and high-MAL2-expression groups. GSEA software (version 4.0) was downloaded from their official website. Tumor tissue samples were divided into high-expression groups and low-expression groups according to the median value of MAL2 expression level. The gene set of h.all.v7.4.symbols.gmt and c2.cp.kegg.v7.4.symbols.gmt, containing a large number of tumors signaling pathway gene sets, were obtained from the MSigDB database. Depending on the default parameters of GSEA software, 1,000 times a random combination was used for enrichment analysis. Sorted by false discovery rates (FDRs), the gene set with FDR <0.25 and normal value <0.05 was the significantly enriched gene set [22].

2.6. Immune Infiltration Analysis of MAL2

The transcripts per million (TPM) normalized RNA-seq data of breast cancer were downloaded from UCSC Xena (https://xena.ucsc.edu/) [23]. Twenty-eight immune gene sets, including different classic immune cell types, were obtained from the study by Charoentong et al. [24]. Next, the ssGSEA was performed by R package GSVA [25] to assess the immune infiltration level of each sample of breast cancer and calculate the responding immune infiltration score based on the TPM data and gene sets [25]. Subsequently, Spearman’s correlation analysis between the MAL2 expression data and the immune infiltration score of 28 immune gene sets was performed. The threshold was set as value < 0.05 and |r| ≥ 0.3 to filter the immune cell types that have significantly infiltrating level correlation with the MAL2 expression.

3. Results

3.1. Overexpression of MAL2 in Breast Cancer Based on Oncomine

We used Oncomine to analyze the difference between MAL2 tumor and normal tissues in different cancer types. Altogether, 214 different types of research results were collected in the Oncomine database. Among them, 25 results showed statistically different expressions of MAL2, including 13 with increased expression and 12 with decreased expression. Results revealed that the MAL2 expression level increased in breast, ovarian, prostate, and pancreatic cancers. On the contrary, MAL2 expression decreased in some cancer types, such as sarcoma, brain and CNS, esophageal, kidney, and lymphoma cancer (Figure 1(a)) [14]. To assess further the MAL2 expression in breast cancer, we analyzed MAL2 across six public expressions, containing 10 analyses from the Oncomine database. As a result, MAL2 was found to significantly upregulate in breast carcinoma tissues. Details are shown in Table 1. Based on the results of the meta-analysis of six data sets, including 10 analyses using the Oncomine database, the results indicated MAL2 overexpression in tumor tissues (median rank = 476.5, ) (Figure 1(b)).

3.2. MAL2 Was Upregulated in Breast Cancer Tissues Compared with Adjacent Tissues Based on TCGA

The MAL2 expression data of 1098 breast cancer samples and 113 adjacent controls were retrieved from the TCGA database. We measured the differences in MAL2 expression in tumor and adjacent tissues using the independent-sample Wilcoxon rank sum test and the paired Wilcoxon rank sum test. Results showed that MAL2 expression was significantly higher in different tissues (tumor tissues vs. adjacent tissues, , Figure 2(a); paired tumor vs. adjacent tissues, , Figure 2(b)).

3.3. Verification of MAL2 Upregulation in Breast Cancer by RT-qPCR

To certify the difference in MAL2 expression in the TCGA data set, we performed RT-qPCR to detect the expression of MAL2 in 32 pairs of breast cancer tissues (including 32 tumor tissues and 32 adjacent tissues). The result showed that MAL2 mRNA was significantly increased in breast cancer tumors compared with adjacent tissues (, Figure 2(c)), which was consistent with the Oncomine data sets analyses.

3.4. Association of MAL2 Expression with Clinical Features in Breast Cancer

We downloaded mRNA-seq data and clinical information of breast cancer patients from TCGA. There were 1211 cases of mRNA-seq data (1098 tumor samples; 113 adjacent tissues) derived from 1091 breast cancer patients in the TCGA database of 1090 equipped with complete clinical data. 1090 patients were finally enrolled in the study. The detailed clinical characteristics—age, gender, T, N, M classification, stage, ER, PR, HER2 status, vital status, and MAL2 expression—are presented in Table 2. In Figure 3, boxplots presented the MAL2 was expressed with the significant difference in the subgroup by M classification (), ER status (), HER2 status (), and vital status ().

3.5. Relationship between MAL2 Expression and Clinical Features in Breast Cancer

To determine correlations of MAL2 expression with clinical factors, we conducted a chi-square test. The median value of MAL2 expression was used to divide patients into a MAL2 high-expression group and a low-expression group. We observed that MAL2 expression was significantly associated with stage (), M classification (), and vital status () (Table 3).

3.6. MAL2 Overexpression Independently Predicted Poor Overall Survival

Kaplan–Meier curves with the log-rank test were applied to explore the prognostic value of MAL2 expression and the overall breast cancer survival rate. As shown in Figure 4, breast cancer patients with high MAL2 expression were associated with worse overall survival (). In addition, high MAL2 expression was associated with poor overall survival in old patients (); clinical stage I/II and III/IV ( and ); subgroup analysis of T classification (T1/T2 and T3/T4) ( and ); patients with lymphatic invasion (); patients with nondistant metastasis (); and patients with positive ER, PR, and HER2 status (, , and , respectively).

The univariate Cox proportional hazards model showed that MAL2 expression, stage, T classification, N classification, and M classification of breast cancer patients were significantly correlated with the prognosis of patients. Moreover, the multivariate Cox proportional hazards model indicated that high MAL2 expression (hazard ratio = 1.792, ), stage (hazard ratio = 1.473, ), and M classification (hazard ratio = 3.093, ) were independent risk factors for overall survival. The median overall survival of the high-expression group was 9.5 years, while that of the low-expression group was 18 years. Both the log-rank test and Cox proportional hazards model showed that the expression of MAL2 was significantly correlated with the prognosis of breast cancer. These results are described in Table 4 and Figure 5.

In conclusion, both in the Kaplan–Meier model and the Cox proportional hazard regression model, the results indicated that MAL2 expression in breast cancer was significantly correlated with the prognosis of patients.

3.7. Identification of MAL2-Related Signaling Pathways by GSEA

Data sets from GSEA showed significant differences (|NES| > 1, FDR < 0.25, NOM ) in MSigDB Collection. The details are described in Figure 6 and Table 5.

The significant pathways by GSEA included MYC targets V1, mTORC1 signaling pathway, insulin signaling pathway, E2F targets, UA response, G2M checkpoint, oocyte meiosis, mitotic spindle, peroxisome, spliceosome, cell cycle, and ubiquitin-mediated proteolysis enriched differentially in MAL2 high-expression phenotype.

3.8. High-Expressed MAL2 Correlates with Reduced Immune Infiltration in Breast Cancer

By Spearman’s correlation test, the MAL2 expression level was found to be negatively correlated with infiltrating levels of eosinophils (r = −0.38, ) and plasmacytoid dendritic cells (r = −0.33, ) (Figures 7(a)7(c)). Additionally, the abundance of infiltration of both eosinophils and plasmacytoid dendritic cells was significantly lower in tumor tissues than in adjacent tissues (Figure 7(d)).

4. Discussion

Breast cancer is the dominating cause of cancer-related mortality among females worldwide. With the rapid development of genomics and molecular biology, identifying the crucial biomarkers for the diagnosis and treatment of breast cancer has become an important research tendency. Recently, numerous research studies showed that MAL2 could be a promising biomarker in various solid tumors. Jennifer et al. showed that MAL2 was increased in ovarian carcinoma, and the overexpression of D52, the binding partner of MAL2, was linked to low overall survival in breast cancer [13]. Chen et al. found that MAL2 was utilized to distinguish pancreatic ductal adenocarcinoma from chronic pancreatitis [26]. Also, the high expression of MAL2 contributed to the short survival time and high distant metastasis rate of postoperative pancreatic cancer patients [8]. Besides, it was reported that MAL2 expression increased in colon cancer tissues and lymph nodes with metastasis with high accuracy and specificity in the diagnosis of colon cancer [10].

As for breast cancer, previous studies showed that MAL2 could induce proliferation and invasion of breast cancer cell lines by adjusting the epithelial-mesenchymal transition [12]. MAL2 was also proved to promote immune evasion by suppressing tumor antigen presentation in breast cancer [11]. However, a further study about the correlation between MAL2 and breast cancer was needed. Here, we conducted this study to consider further the significance of MAL2 expression in the prognosis of breast cancer.

Our study proved the value of MAL2 in breast cancer. MAL2 was upregulated significantly in breast cancer based on the results of the Oncomine and TCGA database analyses. Furthermore, the RT-qPCR outcome verified the high expression of MAL2 in breast cancer, which coincides with the results of the bioinformatics assay and previous studies in other tumors. Moreover, MAL2 expression was correlated with the prognosis of breast cancer. MAL2 could be a potential biomarker for breast cancer.

In clinical factors’ analysis, MAL2 expression was found to be associated with vital status, ER, PR status, and M classification. The MAL2 level of dead patients was higher than that of alive patients, implying that patients with high MAL2 expression were more aggressive. Additionally, patients without metastasis had a higher MAL2 expression level than patients with metastasis, and high MAL2 expression was associated with ER- and HER2-positive patients.

Both the chi-square test and multivariate analysis showed that MAL2 expression is associated with stage and M classification. In our study, breast cancer patients with high expression had shorter overall survival. In the multivariate model, stage, M classification, and MAL2 expression were significantly related to OS (all ), especially high MAL2 expression, advanced breast cancer patients, and those with distant metastasis. This finding may serve as a basis for the proper selection of specific and personalized treatment for breast cancer. Moreover, we found that MAL2 was an independent prognostic factor and might become a biomarker for breast cancer.

To identify the biological function of MAL2 in breast cancer, we used GSEA analysis to predict the pathway associated with MAL2. MYC target V1, mTORC1 signaling pathway, and E2F targets were closely related to the progression of breast cancer.

The high-MAL2-expression phenotype was associated with the activated MYC targets v1 and E2F targets v1 gene sets based on the enrichment score. MYC, a prominent gene in MYC targets v1 gene set, is common in aggressive tumors and contributes to cancer development [27, 28]. Also, MYC gene is overexpressed in triple-negative breast cancer and targeting the gene provides a new treatment [29]. Schulze’s study showed that high MYC Targets v1 enrichment scores were associated with high mutation load, increased infiltration of pro-and anticancerous immune cells, tumor aggressiveness, and poor prognosis of ER-positive cancer [30].

Besides, the E2F transcription factors are downstream effectors of the retinoblastoma protein (pRB) pathway, which is essential to regulate numerous genes essential for DNA replication and cell cycle progression [31]. A previous study showed that E2F transcription factors played key roles in mediating tumor development and metastasis for knockout E2F, leading to the decreased tumor angiogenesis and metastatic capacity of breast cancer. E2Fs could control the expression of genes critical to angiogenesis, remodeling the extracellular matrix, tumor cell survival, and tumor cell interactions with vascular endothelial cells that boost breast cancer metastasis to the lungs [32]. Furthermore, the E2F enrichment score is a marker of breast cancer aggressiveness and predicts the responsiveness of ER-positive/HER2-negative patients to neoadjuvant chemotherapy [33]. Similar research also found that ER-dependent E2F transcription enhanced endocrine resistance in breast cancer [34]. Taken together, speculated high MAL2 expression might be associated with proliferation, metastasis, and prognosis of breast cancer by regulating the genes in these two gene sets.

We finally identified the relationship between the infiltrating level of immune cells and MAL2 using ssGSEA. Interestingly, MAL2 expression was found to be significantly associated with the negatively infiltrating level of eosinophils and plasmacytoid dendritic cells.

A prior study showed eosinophil infiltration was considered a favorable prognosis in breast cancer [35]. In addition, a previous study indicated that low baseline eosinophil count was related to a higher recurrence rate in 419 patients diagnosed with breast cancer [36]. Similar results could be found in recent studies. Low blood eosinophilic relative count was associated with a worse prognosis in 930 breast cancer patients [37]. Moreover, there was a positive correlation between eosinophilic relative counts and both pathological complete remission and survival rate in triple-negative and hormone receptor-negative/HER2-positive breast cancer patients [38].

As for plasmacytoid dendritic cells, it was reported that the proportion of plasmacytoid dendritic cells in triple-negative breast cancer was higher than other subtypes of breast cancer [39]. Besides, high plasmacytoid dendritic cells in triple-negative breast cancer were bound up with a favorable immune response and predicted better survival in 2968 breast cancer patients [39]. Interestingly, in the present study, we found the overexpression of MAL2 was significantly associated with the low infiltrating level of eosinophilic relative counts and plasmacytoid dendritic cells. Interestingly, in the present study, we found that the high MAL2 expression level had a significantly negative correlation with the infiltrating level of eosinophils and plasmacytoid dendritic cells. Moreover, high MAL2 expression was associated with poor OS in patients with breast cancer. Since high eosinophils and plasmacytoid dendritic cells predicted better breast cancer survival in previous research, we speculate that high‐expressed MAL2 might impact the eosinophils and plasmacytoid dendritic cells, which triggered a disadvantageous immune response, leading to poor prognosis in breast cancer.

To the best of our knowledge, this is the first research further analyzing the relationship between MAL2 expression and clinical features and immune infiltration. MAL2 expression may be an independent predictor of a poor disease survival prognosis in breast cancer patients. However, there are a few shortcomings in this study. The sample size included in the experimental validation part was small. And the collected breast cancer tissues were fresh samples, which could not be followed up according to the prognosis of the patients. Besides, further studies are needed to explore the mechanism of MAL2 in breast cancer.

5. Conclusion

In conclusion, our study showed MAL2 expression increased markedly in breast cancer patients and was related to overall survival. MAL2 might be a novel prognostic biomarker of breast cancer. Moreover, high MAL2 expression correlates with reduced immune infiltration of eosinophils and plasmacytoid dendritic cells in breast cancer. However, further studies are warranted to verify the value of MAL2 in breast cancer prognosis evaluation.

Data Availability

Publicly available datasets were analyzed in this study. The data can be found at TCGA repository, https://www.cancer.gov/tcga, and Oncomine repository, http://www.oncomine.org. The data used to support the findings are available from the corresponding author upon request.

Conflicts of Interest

Yue Zhong and Zhenjie Zhuang are co-first authors.

Authors’ Contributions

Yue Zhong and Zhenjie Zhuang contributed equally to this work and are co-first authors.

Acknowledgments

The authors also thank the Lingnan Medical Research Center of Guangzhou University of Chinese Medicine for providing associated facilities.