Role of EHD2 in migration and invasion of human breast cancer cells

EH domain-containing protein 2 (EHD2): Overview, biological function, and therapeutic potential

EH domain-containing protein 2 (EHD2) is a member of the EHD protein family and is mainly located in the plasma membrane, but can also be found in the cytoplasm and endosomes. EHD2 is also a nuclear-cytoplasmic shuttle protein. After entering the cell nuclear, EHD2 acts as a corepressor of transcription to inhibit gene transcription. EHD2 regulates a series of biological processes. As a key regulator of endocytic transport, EHD2 is involved in the formation and maintenance of endosomal tubules and vesicles, which are critical for the intracellular transport of proteins and other substances. The N-terminal of EHD2 is attached to the cell membrane, while its C-terminal binds to the actin-binding protein. After binding, EHD2 connects with the actin cytoskeleton, forming the curvature of the membrane and promoting cell endocytosis. EHD2 is also associated with membrane protein trafficking and receptor signaling, as well as in glucose metabolism and lipid metabolism. In this review, we highlight the recent advances in the function of EHD2 in various cellular processes and its potential implications in human diseases such as cancer and metabolic disease. We also discussed the prospects for the future of EHD2. EHD2 has a broad prospect as a therapeutic target for a variety of diseases. Further research is needed to explore its mechanism, which could pave the way for the development of targeted treatments.

Endocytosis in cancer and cancer therapy

Endocytosis is a complex process whereby cell surface proteins, lipids and fluid from the extracellular environment are packaged, sorted and internalized into cells. Endocytosis is also a mechanism of drug internalization into cells. There are multiple routes of endocytosis that determine the fate of molecules, from degradation in the lysosomes to recycling back to the plasma membrane. The overall rates of endocytosis and temporal regulation of molecules transiting through endocytic pathways are also intricately linked with signalling outcomes. This process relies on an array of factors, such as intrinsic amino acid motifs and post-translational modifications. Endocytosis is frequently disrupted in cancer. These disruptions lead to inappropriate retention of receptor tyrosine kinases on the tumour cell membrane, changes in the recycling of oncogenic molecules, defective signalling feedback loops and loss of cell polarity. In the past decade, endocytosis has emerged as a pivotal regulator of nutrient scavenging, response to and regulation of immune surveillance and tumour immune evasion, tumour metastasis and therapeutic drug delivery. This Review summarizes and integrates these advances into the understanding of endocytosis in cancer. The potential to regulate these pathways in the clinic to improve cancer therapy is also discussed.

mRNA transcriptome profiling of human hepatocellular carcinoma cells HepG2 treated with Catharanthus roseus-silver nanoparticles

BACKGROUND

The demand for the development of cancer nanomedicine has increased due to its great therapeutic value that can overcome the limitations of conventional cancer therapy. However, the presence of various bioactive compounds in crude plant extracts used for the synthesis of silver nanoparticles (AgNPs) makes its precise mechanisms of action unclear.

AIM

To assessed the mRNA transcriptome profiling of human HepG2 cells exposed to Catharanthus roseus G. Don (C. roseus)-AgNPs.

METHODS

The proliferative activity of hepatocellular carcinoma (HepG2) and normal human liver (THLE3) cells treated with C. roseusAgNPs were measured using MTT assay. The RNA samples were extracted and sequenced using BGIseq500 platform. This is followed by data filtering, mapping, gene expression analysis, differentially expression genes analysis, Gene Ontology analysis, and pathway analysis.

RESULTS

The mean IC₅₀ values of C. roseusAgNPs on HepG2 was 4.38 ± 1.59 μg/mL while on THLE3 cells was 800 ± 1.55 μg/mL. Transcriptome profiling revealed an alteration of 296 genes. C. roseusAgNPs induced the expression of stress-associated genes such as MT, HSP and HMOX-1. Cellular signalling pathways were potentially activated through MAPK, TNF and TGF pathways that are responsible for apoptosis and cell cycle arrest. The alteration of ARF6, EHD2, FGFR3, RhoA, EEA1, VPS28, VPS25, and TSG101 indicated the uptake of C. roseus-AgNPs via both clathrin-dependent and clathrin-independent endocytosis.

CONCLUSION

This study provides new insights into gene expression study of biosynthesised AgNPs on cancer cells. The cytotoxicity effect is mediated by the aberrant gene alteration, and more interestingly the unique selective antiproliferative properties indicate the C. roseusAgNPs as an ideal anticancer candidate.

A methylomics-associated nomogram predicts the overall survival risk of stage III to IV ovarian cancer

Accumulating studies demonstrated that DNA methylation may be potential prognostic hallmarks of various cancers. However, few studies have focused on the power of DNA methylation for prognostic prediction in patients with stage III to IV ovarian cancer (OC). Therefore, constructing a methylomics-related indicator to predict overall survival (OS) of stage III to IV OC was urgently required. A total of 520 OC patients with 485,577 DNA methylation sites from TCGA database were selected to develop a robust DNA methylation signature. The 520 patients were clustered into a training group (70%, n = 364 samples) and an internal validation group (30%, n = 156). The training group was used for digging a prognostic predictor based on univariate Cox proportional hazard analysis, least absolute shrinkage and selection operator (LASSO) as well as multivariate Cox regression analysis. The internal and external validation group (ICGC OV-AU project) were used for validating the predictive robustness of the predictor based on receiver operating characteristic (ROC) analysis and Kaplan-Meier survival analysis. We identified a 21-DNA methylation signature-based classifier for stage III-IV OC patients' OS. According to ROC analysis in the internal validation, external validation and entire TCGA set, we proved the high power of the 21-DNA methylation signature for predicting OS (area under the curve [AUC] at 1, 3, 5 years in internal validation set (0.782, 0.739, 0.777, respectively), external validation set (0.828, 0.760, 0.741, respectively), entire TCGA set (0.741, 0.748, 0.781, respectively). Besides, a nomogram was developed via methylation risk score as well as a few clinical variables, and the result showed a high ability of the predictive nomogram. In summary, we used integrated bioinformatics approaches to successfully identified a DNA methylation-associated nomogram, which can predict effectively the OS of patients with stage III to IV OC.

EHD2 overexpression promotes tumorigenesis and metastasis in triple-negative breast cancer by regulating store-operated calcium entry

With nearly all cancer deaths a result of metastasis, elucidating novel pro-metastatic cellular adaptations could provide new therapeutic targets. Here, we show that overexpression of the EPS15-Homology Domain-containing 2 (EHD2) protein in a large subset of breast cancers (BCs), especially the triple-negative (TNBC) and HER2+ subtypes, correlates with shorter patient survival. The mRNAs for EHD2 and Caveolin-1/2, structural components of caveolae, show co-overexpression across breast tumors, predicting shorter survival in basal-like BC. EHD2 shRNA knockdown and CRISPR-Cas9 knockout with mouse Ehd2 rescue, in TNBC cell line models demonstrate a major positive role of EHD2 in promoting tumorigenesis and metastasis. Mechanistically, we link these roles of EHD2 to store-operated calcium entry (SOCE), with EHD2-dependent stabilization of plasma membrane caveolae ensuring high cell surface expression of the SOCE-linked calcium channel Orai1. The novel EHD2-SOCE oncogenic axis represents a potential therapeutic target in EHD2- and CAV1/2-overexpressing BC.

Integrative analysis of transcriptional profile reveals LINC00052 as a suppressor of breast cancer cell migration

BACKGROUND

Long-non-coding RNAs, a class of transcripts with lengths > 200 nt, play key roles in tumour progression. Previous reports revealed that LINC00052 (long intergenic non-coding RNA 00052) was strongly downregulated during breast cancer multicellular spheroids formation and suggested a role in cell migration and oxidative metabolism.

OBJECTIVE

To examine the function of LINC00052 in MCF-7 breast cancer cells.

METHODS

Loss-of-function studies were performed to evaluate LINC00052 role on MCF-7 breast cancer cells. Microarray expression assays were performed to determine genes and cellular functions modified after LINC00052 knockdown. Next, the impact of LINC00052 depletion on MCF-7 cell respiration and migration was evaluated.

RESULTS

1,081 genes were differentially expressed upon LINC00052 inhibition. Gene set enrichment analysis, Gene Ontology and Key Pathway Advisor analysis showed that signalling networks related to cell migration and oxidative phosphorylation were enriched. However, whereas LINC00052 knockdown in MCF-7 cells revealed marginal difference in oxygen consumption rates when compared with control cells, LINC00052 inhibition enhanced cell migration in vitro and in vivo, as observed using a Zebrafish embryo xenotransplant model.

CONCLUSION

Our data show that LINC00052 modulates MCF-7 cell migration. Genome-wide microarray experiments suggest that cancer cell migration is affected by LINC00052 through cytoskeleton modulation and Notch/β-catenin/NF-κB signalling pathways.

Caveola-forming proteins and prostate cancer

Caveolae are specialised and dynamic plasma membrane subdomains, involved in many cellular functions including endocytosis, signal transduction, mechanosensing and lipid storage, trafficking, and metabolism. Two protein families are indispensable for caveola formation and function, namely caveolins and cavins. Mutations of genes encoding these caveolar proteins cause serious pathological conditions such as cardiomyopathies, skeletal muscle diseases, and lipodystrophies. Deregulation of caveola-forming protein expression is associated with many types of cancers including prostate cancer. The distinct function of secretion of the prostatic fluid, and the unique metabolic phenotype of prostate cells relying on lipid metabolism as a main bioenergetic pathway further suggest a significant role of caveolae and caveolar proteins in prostate malignancy. Accumulating in vitro, in vivo, and clinical evidence showed the association of caveolin-1 with prostate cancer grade, stage, metastasis, and drug resistance. In contrast, cavin-1 was found to exhibit tumour suppressive roles. Studies on prostate cancer were the first to show the distinct function of the caveolar proteins depending on their localisation within the caveolar compartment or as cytoplasmic or secreted proteins. In this review, we summarise the roles of caveola-forming proteins in prostate cancer and the potential of exploiting them as therapeutic targets or biological markers.

EHD2 Overexpression Suppresses the Proliferation, Migration, and Invasion in Human Colon Cancer

Background: To investigate how EHD2 influences the development of colon cancer.Methods: Immunohistochemistry of 90 colon cancer tissue specimens were determined the expression of EHD2. The lentivirus-EHD2-transfected colon cancer cells were conducted to evaluate the biological behaviors.Results: EHD2 was closely associated with clinic pathological parameters (p < 0.001). EHD2 upregulation was relative with a longer overall survival. The results of the univariate and multivariate analyses indicated that EHD2 could be an independent prognosis marker. EHD2 overexpression suppressed cell invasion and proliferation, but enhanced cell apoptosis and cell cycle arrest.Conclusions: EHD2 might represent a therapeutic target of colon cancer.IMPACT STATEMENTWhat is already known on this subject? Membrane trafficking is crucial for cell proliferation, differentiation and apoptosis, especially tumorigenesis and development. EHD2 proteins play an important role in the regulation of membrane trafficking in endocytosis. EHD2 has been suggested to participate in the occurrence of some malignancies.What are the new findings? EHD2 could be an independent prognosis marker in colon cancer. EHD2 overexpression suppressed cell invasion and proliferation, but enhanced cell apoptosis and cell cycle arrest in vitro. EHD2 overexpression markedly increased the expression of EMT marker E-cadherin in colon cancer.How might it impact on clinical practice in the foreseeable future? EHD2 overexpression may inhibit tumorigenesis in colon cancer through the modulation of E-cadherin, the critical marker of EMT which is closely related to invasion and distant metastasis of tumor cells.

A Genetic Analysis of Tumor Progression in Drosophila Identifies the Cohesin Complex as a Suppressor of Individual and Collective Cell Invasion

Metastasis is the leading cause of death for patients with cancer. Consequently it is imperative that we improve our understanding of the molecular mechanisms that underlie progression of tumor growth toward malignancy. Advances in genome characterization technologies have been very successful in identifying commonly mutated or misregulated genes in a variety of human cancers. However, the difficulty in evaluating whether these candidates drive tumor progression remains a major challenge. Using the genetic amenability of Drosophila melanogaster we generated tumors with specific genotypes in the living animal and carried out a detailed systematic loss-of-function analysis to identify conserved genes that enhance or suppress epithelial tumor progression. This enabled the discovery of functional cooperative regulators of invasion and the establishment of a network of conserved invasion suppressors. This includes constituents of the cohesin complex, whose loss of function either promotes individual or collective cell invasion, depending on the severity of effect on cohesin complex function.

•

Screen identifies genes that affect tumor behavior in a wide variety of ways

•

A functionally validated network of invasion-suppressor genes was generated

•

Loss of cohesin complex function can promote individual or collective cell invasion

•

The fly pupal notum is an excellent in vivo system to study tumor progression

EHD2 is a Predictive Biomarker of Chemotherapy Efficacy in Triple Negative Breast Carcinoma

EHD2 is a mechanotransducing ATPase localized in caveolae invaginations at the plasma membrane. EHD2 has recently been associated with several human cancers, however the significance of EHD2 transcript levels in cancer prognosis remains debated. Breast cancer is the most commonly occurring cancer in women and prognosis is variable depending on the subtypes. Triple negative breast cancer (TNBC) often has a poor therapeutic response. The aim of this study was to assess the prognostic significance of EHD2 transcripts and protein expression levels in breast carcinomas. We found that low EHD2 levels were associated with enhanced proliferation, migration and invasion of TNBC cells. EHD2 expression was significantly reduced in TNBC tissues and the loss of EHD2 led to higher expression of the pro-tumoral cytokine IL-8. In apparent contradiction with in vitro data, multivariate analysis of two independent cohorts of breast cancer patients revealed that low EHD2 was in fact associated with good prognosis in the highly proliferative TNBC subtype. Accordingly, TNBC low EHD2 expressers were found to benefit the most from chemotherapy when compared to all subtypes of breast cancers. Our study validates EHD2 expression level as an independent prognostic factor of metastasis-free survival and as a new predictive marker of chemotherapy efficacy in TNBC patients.

Effect of EH domain containing protein 2 on the biological behavior of clear cell renal cell carcinoma

To investigate the effects of EH domain containing protein 2 (EHD2) on clear cell renal cell carcinoma (ccRCC) and provide new insights for the clinical treatment of rental cancer. Forty patients (26 males and 14 females, 62.4 ± 5.7 years old) with ccRCC were selected from January 2015 to December 2016 to serve as research subjects in this study. The EHD2 protein expression in the tumor tissues and adjacent healthy tissues of ccRCC patients were detected by Western Blot assay. The cells of ccRCC cell lines RLC-310 and 786-O were divided into normal control group (control), no-load control group (pLV), EHD2 overexpression group (pLV-EHD2), and EHD2 interference group (pLV-siEHD2). The expression levels of EHD2 protein in each group of cells were detected by western blot. The cell proliferation was detected by Cell Counting Kit-8 (CCK-8) assay. Wound healing assay was performed to check the cell migration ability. Transwell invasion assay was used to detect the cell invasion ability. Cell apoptosis was detected by flow cytometry. The expression level of EHD2 was significantly increased in pLV-EHD2 group and decreased in pLV-siEHD2 group compared with control group and pLV-siEHD2 group, indicating the successfully established EHD2 overexpression cell line and EHD2 RNA interference cell line. EHD2 overexpression enhanced the proliferation, invasion, and migration but inhibited the apoptosis of ccRCC cells, while EHD2 interference showed opposite functions. EHD2 interference can inhibit the development of ccRCC by inhibiting the proliferation, invasion, and migration, and EHD2 can potentially serve as a molecular target for the clinical treatment of ccRCC.

EHD2 is a mechanotransducer connecting caveolae dynamics with gene transcription

Caveolae are dynamic mechanosensors. Torrino et al. show that EHD2 plays a crucial role in the adaptation to mechanical perturbations by maintaining the caveolae reservoir at the plasma membrane after changes in membrane tension and connecting caveolae mechanosensing at the plasma membrane with the regulation of gene transcription.

Caveolae are small invaginated pits that function as dynamic mechanosensors to buffer tension variations at the plasma membrane. Here we show that under mechanical stress, the EHD2 ATPase is rapidly released from caveolae, SUMOylated, and translocated to the nucleus, where it regulates the transcription of several genes including those coding for caveolae constituents. We also found that EHD2 is required to maintain the caveolae reservoir at the plasma membrane during the variations of membrane tension induced by mechanical stress. Metal-replica electron microscopy of breast cancer cells lacking EHD2 revealed a complete absence of caveolae and a lack of gene regulation under mechanical stress. Expressing EHD2 was sufficient to restore both functions in these cells. Our findings therefore define EHD2 as a central player in mechanotransduction connecting the disassembly of the caveolae reservoir with the regulation of gene transcription under mechanical stress.

A muscle-specific protein 'myoferlin' modulates IL-6/STAT3 signaling by chaperoning activated STAT3 to nucleus

Myoferlin, a member of ferlin family of proteins, was first discovered as a candidate gene for muscular dystrophy and cardiomyopathy. Recently, myoferlin was shown to be also expressed in endothelial and cancer cells where it was shown to modulate vascular endothelial growth factor (VEGFR)-2 and epidermal growth factor receptor (EGFR) signaling by enhancing their stability and recycling. Based on these reports, we hypothesized that myoferlin might be regulating IL-6 signaling by modulating IL-6R stabilization and recycling. However, in our immunoprecipitation (IP) experiments, we did not observe myoferlin binding with IL-6R. Instead, we made a novel discovery that in resting cells myoferlin was bound to EHD2 protein and when cells were treated with IL-6, myoferlin dissociated from EHD2 and binds to activated STAT3. Interestingly, myoferlin depletion did not affect STAT3 phosphorylation, but completely blocked STAT3 translocation to nucleus. In addition, inhibition of STAT3 phosphorylation by phosphorylation-defective STAT3 mutants or JAK inhibitor blocked STAT3 binding to myoferlin and nuclear translocation. Myoferlin knockdown significantly decreased IL-6-mediated tumor cell migration, tumorsphere formation and ALDH-positive cancer stem cell population, in vitro. Furthermore, myoferlin knockdown significantly decreased IL-6-meditated tumor growth and tumor metastasis. Based on these results, we have proposed a novel model for the role of myoferlin in chaperoning phosphorylated STAT3 to the nucleus.

Prognostic implication of histological features associated with EHD2 expression in papillary thyroid carcinoma

Papillary thyroid carcinoma (PTC) is a heterogeneous tumor with various histological and molecular subtypes. EHD2 is involved in endocytosis and endosomal recycling. This study aimed to investigate the prognostic significance of EHD2 expression in PTC and develop a new model for predicting persistent/recurrent disease after thyroidectomy. Pathologic slides of 512 consecutive patients with PTC ≥ 1 cm were retrospectively reviewed. BRAF mutation analysis and immunohistochemistry for EHD2 were performed. Clinical significance of EHD2 mRNA expression was analyzed in 388 PTC patients using The Cancer Genome Atlas dataset. The presence of dyscohesive cells and psammoma bodies were found have significant association with persistent/recurrent disease (p = 0.049 and p = 0.038, respectively). The best discrimination of disease-free survival was found by dividing patients into three prognostic groups based on the following two risk factors according to the size category: psammoma bodies ≥ 4 and dyscohesive cells (≥ 1% and ≥ 20% in PTCs of < 2.0 cm and ≥ 2.0 cm, respectively). In PTCs of ≥ 2.0 cm, patients with the two risk factors had a hazard ratio of 13.303 (p = 0.005) compared to those without risk factors. High expression level of EHD2 was associated with BRAF V600E (p < 0.001), presence of dyscohesive cells (p = 0.010), and absence of psammoma bodies (p = 0.001). Increased EHD2 mRNA expression level was associated with extrathyroidal extension (p < 0.001), pT3-4 (p < 0.001), lymph node metastasis (p < 0.001), higher risk of recurrence (p < 0.001), and BRAF V600E (p < 0.001). Our prognostic model is useful for predicting persistent/recurrent disease after surgery of PTC. EHD2 mRNA expression could be a novel prognostic marker for PTC patients.

Decreased Expression of EHD2 Promotes Tumor Metastasis and Indicates Poor Prognosis in Hepatocellular Carcinoma

BACKGROUND

Metastasis remains the most common cause of lethal outcomes in hepatocellular carcinoma (HCC) after curative resection. Understanding molecular mechanisms that regulate metastasis process is crucial for improving treatment of hepatocellular carcinoma.

AIMS

In this article, we examined whether Eps15 homology domain-containing 2 (EHD2) played a critical role in hepatocellular carcinoma metastasis and explored the possible mechanism.

METHODS

EHD2 and E-cadherin expression levels in hepatocellular carcinoma patients were examined using Western blotting and immunohistochemistry. The cell migration and invasion were evaluated by wound-healing assay and trans-well assay. Epithelial-mesenchymal transition was analyzed by immunofluorescence, and the vital markers were detected by Western blotting. The correlation of EHD2 and E-cadherin was confirmed by co-immunoprecipitation.

RESULTS

EHD2 expression, along with the epithelial marker E-cadherin, was markedly reduced in tumor tissues than in adjacent noncancerous tissues. Moreover, EHD2 was positively correlated with E-cadherin, histological grade, tumor metastasis, and microvascular invasion. Kaplan-Meier survival analysis showed that hepatocellular carcinoma patients with decreased EHD2 expression had shorter overall survival times than those with higher EHD2 expression. Knockdown of EHD2 induced an increase in cell invasion and changes characteristic of epithelial-mesenchymal transition, while overexpression of EHD2 inhibited these processes.

CONCLUSIONS

Molecular data indicated that EHD2 inhibited migration and invasion of hepatocellular carcinoma probably by interacting with E-cadherin and it might be an independent, significant risk factor for survival after curative resection.

Over-expression of EPS15 is a favorable prognostic factor in breast cancer

As a crucial player in terminating growth factor signaling, EPS15 plays important roles in many malignancies including breast cancer. To explore the potential association of EPS15 with the clinical outcome of breast cancer, we conducted gene expression survival analysis using six independent datasets, checked its expression quantitative loci and their associated genes, and explored the networking of these genes with EPS15. Our results show that over-expression of EPS15 is significantly associated with a favorable clinical outcome of breast cancer, especially in tumors harbouring a positive estrogen receptor status. 21 unique SNPs were found to be associated with EPS15 expression. Among the neighboring genes of these SNPs, five (MTUS1, DOCK5, MSRA, SLIT3 and SKAP1) are genetically connected with EPS15 and its physical partners. These genes including EPS15 also show significant concurrent expressions, and four exhibit distinct relevance regarding patient survival. High expressions of EPS15 and MSRA show a distinct combinatorial favorable survival, suggesting the clinical relevance of their co-activation. In summary, over-expression of EPS15 is a potential favorable prognostic marker in breast cancer, which can be used clinically alone or together with other genes such as MSRA to avail therapeutic decision-making.

Phosphatidylinositol glycan anchor biosynthesis, class X containing complex promotes cancer cell proliferation through suppression of EHD2 and ZIC1, putative tumor suppressors

We identified phosphatidylinositol glycan anchor biosynthesis, class X (PIGX), which plays a critical role in the biosynthetic pathway of glycosylphosphatidylinositol (GPI)-anchor motif, to be upregulated highly and frequently in breast cancer cells. Knockdown of PIGX as well as reticulocalbin 1 (RCN1) and reticulocalbin 2 (RCN2), which we found to interact with PIGX and was indicated to regulate calcium-dependent activities, significantly suppressed the growth of breast cancer cells. We also identified PIGX to be a core protein in an RCN1/PIGX/RCN2 complex. Microarray analysis revealed that the expression of two putative tumor suppressor genes, Zic family member 1 (ZIC1) and EH-domain containing 2 (EHD2), were upregulated commonly in cells in which PIGX, RCN1, or RCN2 was knocked down, suggesting that this RCN1/PIGX/RCN2 complex could negatively regulate the expression of these two genes and thereby contribute to human breast carcinogenesis. Our results imply that PIGX may be a good candidate molecule for development of novel anticancer drugs for breast cancer.

Integrative radiogenomic analysis for genomic signatures in glioblastomas presenting leptomeningeal dissemination

Despite therapeutic advances, the prognosis for glioblastoma (GBM) remains poor. In particular, leptomeningeal dissemination (LMD) has a dismal prognosis. The aim of this study was to identify tumor molecular phenotype, which has a great propensity to develop LMD. Between May 2004 and December 2012, a total of 145 GBM tumor samples were obtained from data registry. A total of 20 of the 145 patients with GBM were found to develop LMD. A specialized radiologist confirmed the diagnosis of LMD on magnetic resonance imaging. To clarify the genomic signatures in GBM with LMD, we performed integrative analysis of whole transcriptome sequencing and copy number alteration in the radiological features indicating LMD phenotypes in GBM. Eleven newly diagnosed patients with GBM with LMD had worse prognosis than those without LMD (median 5.55 vs. 12.94 months, P < 0.0001). Integrating analysis using gene expression based on the change of copy number revealed that SPOCK1, EHD2, SLC2A3, and ANXA11 were highly expressed with the gain of copy number, compared with the gene expression in the non-LMD group. In addition, it was demonstrated that NME2, TMEM100, and SIVA1 were downregulated with the loss of copy number. We also found that mesenchymal subtype accounted for 50% in LMD group, whereas mesenchymal subtype consisted of 29% in non-LMD group, even though there was no statistical significance (P = 0.06). Through this radiogenomic analysis, we suggested the possibility of finding candidate genes associated with LMD and highlighted the significance of integrating approach to clarify the molecular characteristics in LMD.