Crosstalk between microRNA30a/b/c/d/e-5p and the canonical Wnt pathway: implications for multiple myeloma therapy

BCL9 is a Risk Factor of Neck Lymph Nodes Metastasis and Correlated with Immune Cell Infiltration in Papillary Thyroid Carcinoma

Purpose

B-cell lymphoma 9 (BCL9), a key transcription co-activator of the Wnt pathway, contributed to tumor progression and metastasis in various tumors, whereas, the role of BCL9 in papillary thyroid cancer (PTC) has not been investigated.

Methods

We acquired PTC gene expression data from The Cancer Genome Atlas (TCGA) database. Fifty-nine PTC tissues were applied to validate the clinical significance of BCL9. Cell experiments were applied to investigate the role of BCL9. Bioinformatics analysis was employed to investigate the biological functions of BCL9.

Results

We found that BCL9 was higher expressed (P < 0.05) and an independent risk factor for lymph node metastasis (OR = 3.770, P = 0.025), as well as associated with poorer progression-free survival (PFS) (P = 0.049) in PTC. BCL9 knockdown inhibited proliferation and invasion of PTC cells. BCL9 was positively associated with the key genes of Wnt/β-catenin and MAPK pathway by co-expression analysis. GO, KEGG and GSEA analysis showed BCL9 might participated in PPAR, cAMP, and focal adhesion pathway. CIBERSORT analysis found BCL9 was negatively associated with CD8+ T cells and NK cell infiltration and positively with PD-L1 expression.

Conclusion

Therefore, BCL9 was associated with lymph node metastasis and shorter PFS of PTC, due to promotion of PTC cell proliferation and invasion, activation of Wnt/β-catenin and MAPK pathway, inhibition of CD8+ T and NK cell infiltration, and promotion of PD-L1 expression.

Chromosome 1 instability in multiple myeloma: Aberrant gene expression, pathogenesis, and potential therapeutic target

Multiple myeloma (MM), the terminally differentiated B cells malignancy, is widely considered to be incurable since many patients have either developed drug resistance or experienced an eventual relapse. To develop precise and efficient therapeutic strategies, we must understand the pathogenesis of MM. Thus, unveiling the driver events of MM and its further clonal evolution will help us understand this complicated disease. Chromosome 1 instabilities are the most common genomic alterations that participate in MM pathogenesis, and these aberrations of chromosome 1 mainly include copy number variations and structural changes. The chromosome 1q gains/amplifications and 1p deletions are the most frequent structural changes of chromosomes in MM. In this review, we intend to focus on the genes that are affected by chromosome 1 instability: some tumor suppressors were lost or down regulated in 1p deletions, and others that contributed to tumorigenesis were upregulated in 1q gains/amplifications. We have summarized their biological function as well as their roles in the MM pathogenesis, hoping to uncover potential novel therapeutical targets and promote the development of future therapeutic approaches.

A novel β-catenin/BCL9 complex inhibitor blocks oncogenic Wnt signaling and disrupts cholesterol homeostasis in colorectal cancer

Dysregulated Wnt/β-catenin signaling is implicated in the pathogenesis of many human cancers, including colorectal cancer (CRC), making it an attractive clinical target. With the aim of inhibiting oncogenic Wnt activity, we developed a high-throughput screening AlphaScreen assay to identify selective small-molecule inhibitors of the interaction between β-catenin and its coactivator BCL9. We identified a compound that consistently bound to β-catenin and specifically inhibited in vivo native β-catenin/BCL9 complex formation in CRC cell lines. This compound inhibited Wnt activity, down-regulated expression of the Wnt/β-catenin signature in gene expression studies, disrupted cholesterol homeostasis, and significantly reduced the proliferation of CRC cell lines and tumor growth in a xenograft mouse model of CRC. This study has therefore identified a specific small-molecule inhibitor of oncogenic Wnt signaling, which may have value as a probe for functional studies and has important implications for the development of novel therapies in patients with CRC.

Identification of candidate targets and mechanisms involved in miRNA regulation in multiple myeloma

Background

Multiple myeloma (MM) is a complex disease affected by many factors. The recognition of miRNA networks is helpful for specific detection and personalised treatment.

Methods

mRNA expression profiles were obtained from GSE39754 and GSE87830, and differentially expressed mRNAs (DEmRs) between MM and controls were identified. The intersection of the two sets of DEmRs in GSE39754 and GSE87830 was identified as common mRNAs, and enrichment analysis was subsequently performed. Moreover, we analysed differentially expressed miRNAs (DEmiRs) between MM and controls in GSE87830. A regulatory network of target mRNAs related to the overall survival of MM patients was then constructed.

Results

In this study, a total of 356 common mRNAs were identified that were significantly enriched in neutrophil-mediated immunity, Th17 cell differentiation and PI3K-Akt signalling pathways. Moreover, we identified 103 DEmiRs and predicted 91 differentially expressed mRNAs as target mRNAs. Cox regression analysis was used to screen 14 target mRNAs that significantly affected the survival of MM patients. In the constructed integrated regulatory network, HIF1A and THBS1 were found to participate in Th17 cell differentiation and PI3K-Akt signalling pathways.

Conclusion

These findings improve the understanding of the regulatory mechanisms of MM. Genes that are part of integrated regulatory networks may represent candidate targets for MM treatment.

The role of Wnt/β-catenin signaling pathway in the pathogenesis and treatment of multiple myeloma (review)

Multiple myeloma (MM) is a refractory hematological malignancy characterized by aberrant accumulation of plasma cells. Patients with MM are susceptible to becoming resistant to chemotherapy, eventually leading to relapse. Progression of MM is largely dependent on the bone marrow microenvironment. Stromal cells in the bone marrow microenvironment secrete Wnt ligands to activate Wnt signaling in MM, which is mediated through the transcription regulator β-catenin. In addition, Wnt/β-catenin pathway encourages osteoblast differentiation and bone formation, dysregulation of which is responsible for proliferation and drug resistance of MM cells. As a result, direct inhibition or silencing of β-catenin or associated genes in the Wnt/β-catenin pathway has been proposed to be an effective therapeutic anti-MM strategy. However, the underlying regulatory mechanism of the Wnt/β-catenin pathway in MM remains to be fully elucidated. Herein, we summarized research advances on the specific genes and molecular biology process of Wnt/β-catenin pathway involved in tumorigenesis of MM, as well as the interaction with bone marrow microenvironment. Additionally, comprehensive summaries of drugs or small molecule inhibitors acting on Wnt/β-catenin pathway and targeting MM were introduced. This review intends to provide an overview of theoretical supports for novel Wnt/β-catenin pathway based treatment strategies in MM.

Current Status of MicroRNAs that Target the Wnt Signaling Pathway in Regulation of Osteogenesis and Bone Metabolism: A Review

The directional differentiation of bone mesenchymal stem cells (BMSCs) is regulated by a variety of transcription factors and intracellular signaling pathways. In the past, it was thought that the directional differentiation of BMSCs was related to transforming growth factors, such as bone morphogenetic protein (BMP) and MAPK pathway. However, in recent years, some scholars have pointed out that the Wnt signaling pathway, which is a necessary complex network of protein interactions for biological growth and development, takes a significant role in this process and plays a major part in regulating the development of osteoblasts by exerting signal transduction into cells. Also, they have proved the Wnt protein therapeutic truly have positive effects on the viability and osteogenic capacity of bone graft. Recent studies have shown that microRNAs (miRNAs) play an important regulatory role in this process. MiRNAs such as miRNA-218, miRNA-335, miRNA-29, microRNA-30 and other miRNAs exert negative or positive effects on some crucial molecules in the Wnt/β-catenin pathway, which in turn affect bone metabolism and osteopathy. Thus, miRNAs have been suggested as therapeutic targets for some metabolic bone diseases. This article aims to provide an update on the current status of microRNAs that target the Wnt signaling pathway in the regulation of osteogenesis and bone metabolism and includes a discussion of future areas of research, which can be a theoretical basis for bone metabolism-related diseases.

Identification of MicroRNAs With In Vivo Efficacy in Multiple Myeloma-related Xenograft Models

BACKGROUND/AIM

Multiple myeloma is a B-cell neoplasm, which can spread within the marrow of the bones forming many small tumors. In advanced disease, multiple myeloma can spread to the blood as plasma cell leukemia. In some cases, a localized tumor known as plasmacytoma is found within a single bone. Despite the approval of several agents such as melphalan, corticosteroids, proteasome inhibitors, thalidomide-based immuno-modulatory agents, histone deacetylase inhibitors, a nuclear export inhibitor and monoclonal antibodies daratuzumab and elatuzumab, the disease presently remains uncurable.

MATERIALS AND METHODS

In order to define new targets and treatment modalities we searched the literature for microRNAs, which increase or inhibit in vivo efficacy in multiple-myeloma-related xenograft models.

RESULTS AND CONCLUSION

We identified six up-regulated and twelve down-regulated miRs, which deserve further preclinical validation.

Targeting of CD38 by the Tumor Suppressor miR-26a Serves as a Novel Potential Therapeutic Agent in Multiple Myeloma

Multiple myeloma is an incurable refractory hematologic malignancy arising from plasma cells in the bone marrow. Here we investigated miR-26a function in multiple myeloma and tested single-wall carbon nanotube delivery of miR-26a in vitro and in vivo. miR-26a was downregulated in patients with multiple myeloma cells compared with plasma cells from healthy donors. miR-26a overexpression inhibited proliferation and migration and induced apoptosis in multiple myeloma cell lines. To identify the targets of miR-26a, RPMI8226-V-miR-26-GFP and RPMI8226-V-GFP cells were cultured using stable isotope labeling by amino acids in cell culture (SILAC) medium, followed by mass spectrometry analysis. In multiple myeloma cells overexpressing miR-26a, CD38 protein was downregulated and subsequently confirmed to be a direct target of miR-26a. Depletion of CD38 in multiple myeloma cells duplicated the multiple myeloma inhibition observed with exogenous expression of miR-26a, whereas restoration of CD38 overcame the inhibition of miR-26a in multiple myeloma cells. In a human multiple myeloma xenograft mouse model, overexpression of miR-26a inhibited CD38 expression, provoked cell apoptosis, and inhibited cell proliferation. Daratumumab is the first CD38 antibody drug for monotherapy and combination therapy for patients with multiple myeloma, but eventually resistance develops. In multiple myeloma cells, CD38 remained at low level during daratumumab treatment, but a high-quality response is sustained. In daratumumab-resistant multiple myeloma cells, CD38 expression was completely restored but failed to correlate with daratumumab-induced cell death. Therefore, a therapeutic strategy to confer selection pressure to maintain low CD38 expression in multiple myeloma cells may have clinical benefit. SIGNIFICANCE: These results highlight the tumor suppressor function of miR-26a via its targeting of CD38 and suggest the therapeutic potential of miR-26a in patients with multiple myeloma.

BCL9 provides multi-cellular communication properties in colorectal cancer by interacting with paraspeckle proteins

Colorectal cancer (CRC) is the third most commonly diagnosed cancer, which despite recent advances in treatment, remains incurable due to molecular heterogeneity of tumor cells. The B-cell lymphoma 9 (BCL9) oncogene functions as a transcriptional co-activator of the Wnt/β-catenin pathway, which plays critical roles in CRC pathogenesis. Here we have identified a β-catenin-independent function of BCL9 in a poor-prognosis subtype of CRC tumors characterized by expression of stromal and neural associated genes. In response to spontaneous calcium transients or cellular stress, BCL9 is recruited adjacent to the interchromosomal regions, where it stabilizes the mRNA of calcium signaling and neural associated genes by interacting with paraspeckle proteins. BCL9 subsequently promotes tumor progression and remodeling of the tumor microenvironment (TME) by sustaining the calcium transients and neurotransmitter-dependent communication among CRC cells. These data provide additional insights into the role of BCL9 in tumor pathogenesis and point towards additional avenues for therapeutic intervention.

MicroRNA-26a inhibits multiple myeloma cell growth by suppressing cyclin-dependent kinase 6 expression

MicroRNA-26a (miR-26a) has been reported to be involved in the tumorigenesis of several tumors, but its biological function and molecular mechanism in multiple myeloma (MM) are still unknown. In this study, we found that overexpression of miR-26a obviously inhibited MM cell growth, and delayed tumor growth in xenografts. Further studies showed that overexpression of miR-26a induced cell cycle arrest at G0/G1 phase in MM cells. MiR-26a mimic down-regulated the expression levels of CDK6 and E2F1, but up-regulated p53 and p21 expression. In contrast, overexpression of CDK6 decreased the effect of miR-26a mimic on MM cell survival. Moreover, miR-26a targeted CDK6 mRNA and thus suppressed CDK6 protein expression. Overexpression of miR-26a also enhanced the cytotoxic action of doxorubicin against MM. These results demonstrated that miR-26a was involved in the development of MM through regulating CDK6 signaling pathway, and indicated that miR-26a could be as a novel target for anti-tumor therapy in clinic as a single strategy or in combination with other anti-tumor drugs in MM.

Serum miR-30d as a novel biomarker for multiple myeloma and its antitumor role in U266 cells through the targeting of the MTDH/PI3K/Akt signaling pathway

Multiple myeloma (MM) is a hematological tumor and is characterized by the infiltration of malignant clonal plasma cells (PCs) in bone marrow. MicroRNAs (miRNAs or miRs) have been reported to play an important role in the genesis and progression of MM. However, little is known about the clinical diagnostic value and biological functions of miR-30d in MM. In this study, to investigate the role of miR-30d in MM, we used reverse transcription-quantitative polymerase chain reaction quantitative (RT-qPCR) to detect the relative expression level of miR-30d in the serum of 81 patients with primary MM and 78 healthy donors (HDs). The biological functions of miR-30d were then assessed by CCK-8 assay, flow cytometric analysis of apoptosis and western blot (WB) analysis in U266 cells. Moreover, the confirmation of the target gene of miR-30d was conducted by luciferase reporter assay. Our results indicated that miR-30d expression was significantly downregulated in the serum of patients with primary MM compared with that of the HDs and that it was significantly associated with several clinical indicators of MM. Further cell functional analyses using the U266 cells revealed that miR-30d functions as a tumor suppressor gene in MM by inhibiting cell viability and promoting cell apoptosis. Moreover, miR-30d was confirmed to directly bind to the 3'UTR of its target gene, metadherin (MTDH) and inhibit the activation of the downstream PI3K/Akt signaling pathway. On the whole, the findings of this study indicate that the serum expression level of miR-30d is of great significance to the diagnosis and treatment monitoring of patients with MM. Moreover, miR-30d carries out its antitumor role in U266 cells through the inhibition of the activation of the PI3K/Akt signaling pathway by negatively regulating MTDH, which reveals its potential for use as a therapeutic target for MM.

Circulating miR-30a-5p as a prognostic biomarker of left ventricular dysfunction after acute myocardial infarction

Left ventricular (LV) dysfunction after acute myocardial infarction (AMI) is associated with an increased risk of heart failure (HF) development. Diverse microRNAs (miRNAs) have been shown to appear in the bloodstream following various cardiovascular events. The aim of this study was to identify prognostic miRNAs associated with LV dysfunction following AMI. Patients were divided into subgroups comprising patients who developed or not LV dysfunction within six months of the infarction. miRNA profiles were determined in plasma and serum samples of the patients on the first day of AMI. Levels of 14 plasma miRNAs and 16 serum miRNAs were significantly different in samples from AMI patients who later developed LV dysfunction compared to those who did not. Two miRNAs were up-regulated in both types of material. Validation in an independent group of patients, using droplet digital PCR (ddPCR) confirmed that miR-30a-5p was significantly elevated on admission in those patients who developed LV dysfunction and HF symptoms six months after AMI. A bioinformatics analysis indicated that miR-30a-5p may regulate genes involved in cardiovascular pathogenesis. This study demonstrates, for the first time, a prognostic value of circulating miR-30a-5p and its association with LV dysfunction and symptoms of HF after AMI.

The potential function of microRNAs as biomarkers and therapeutic targets in multiple myeloma

Multiple myeloma (MM), accounting for ~1% of all types of human cancer and 13% of all hematological malignancies, is characterized by the malignant proliferation of monoclonal plasma cells (PCs) in the bone marrow. MM leads to end stage organ impairment, including bone lesions, renal dysfunction, hypercalcemia and anemia. So far, the specific pathogenesis of MM remains unclear and no early-stage sensitive biomarker of MM has been well characterized. Furthermore, treating MM is difficult, as the majority of patients eventually relapse or become refractory following treatment using presently available methods. To date, a number of studies have demonstrated that microRNAs (miRNAs) may serve crucial functions in the progression of numerous cancers, including MM. During the tumorigenesis and pathogenesis of MM, there are multiple carcinogenic events that involve the pernicious transformation from normal to malignant PCs. miRNAs, as oncogenes or tumor suppressors, regulate MM progression-related signaling pathways. In the present review, the up-to-date preliminary basic studies and associated clinical works on the underlying mechanisms of aberrant miRNA profiling in MM have been summarized, including an evaluation of its value as a potential biomarker and a novel therapeutic strategy for MM.

Molecular targets of luteolin in cancer

Many food-derived phytochemical compounds and their derivatives represent a cornucopia of new anticancer compounds. Despite extensive study of luteolin, the literature has no information on the exact mechanisms or molecular targets through which it deters cancer progression. This review discusses existing data on luteolin's anticancer activities and then offers possible explanations for and molecular targets of its cancer-preventive action. Luteolin prevents tumor development largely by inactivating several signals and transcription pathways essential for cancer cells. This review also offers insights into the molecular mechanisms and targets through which luteolin either prevents cancer or mediates cancer cell death.

Inhibition of miR-21 restores RANKL/OPG ratio in multiple myeloma-derived bone marrow stromal cells and impairs the resorbing activity of mature osteoclasts

miR-21 is an oncogenic microRNA (miRNA) with an emerging role as therapeutic target in human malignancies, including multiple myeloma (MM). Here we investigated whether miR-21 is involved in MM-related bone disease (BD). We found that miR-21 expression is dramatically enhanced, while osteoprotegerin (OPG) is strongly reduced, in bone marrow stromal cells (BMSCs) adherent to MM cells. On this basis, we validated the 3′UTR of OPG mRNA as miR-21 target. Constitutive miR-21 inhibition in lentiviral-transduced BMSCs adherent to MM cells restored OPG expression and secretion. Interestingly, miR-21 inhibition reduced RANKL production by BMSCs. Overexpression of protein inhibitor of activated STAT3 (PIAS3), which is a direct and validated target of miR-21, antagonized STAT3-mediated RANKL gene activation. Finally, we demonstrate that constitutive expression of miR-21 inhibitors in BMSCs restores RANKL/OPG balance and dramatically impairs the resorbing activity of mature osteoclasts. Taken together, our data provide proof-of-concept that miR-21 overexpression within MM-microenviroment plays a crucial role in bone resorption/apposition balance, supporting the design of innovative miR-21 inhibition-based strategies for MM-related BD.

Meta-Analysis of Differential Connectivity in Gene Co-Expression Networks in Multiple Sclerosis

Differential gene expression analyses to investigate multiple sclerosis (MS) molecular pathogenesis cannot detect genes harboring genetic and/or epigenetic modifications that change the gene functions without affecting their expression. Differential co-expression network approaches may capture changes in functional interactions resulting from these alterations. We re-analyzed 595 mRNA arrays from publicly available datasets by studying changes in gene co-expression networks in MS and in response to interferon (IFN)-β treatment. Interestingly, MS networks show a reduced connectivity relative to the healthy condition, and the treatment activates the transcription of genes and increases their connectivity in MS patients. Importantly, the analysis of changes in gene connectivity in MS patients provides new evidence of association for genes already implicated in MS by single-nucleotide polymorphism studies and that do not show differential expression. This is the case of amiloride-sensitive cation channel 1 neuronal (ACCN1) that shows a reduced number of interacting partners in MS networks, and it is known for its role in synaptic transmission and central nervous system (CNS) development. Furthermore, our study confirms a deregulation of the vitamin D system: among the transcription factors that potentially regulate the deregulated genes, we find TCF3 and SP1 that are both involved in vitamin D3-induced p27Kip1 expression. Unveiling differential network properties allows us to gain systems-level insights into disease mechanisms and may suggest putative targets for the treatment.

Role of microRNA in regulation of myeloma-related angiogenesis and survival

Multiple myeloma (MM) is a malignant disease caused by clonal proliferation of plasma cells that result in monoclonal gammopathy and severe end organ damage. Despite the uniform clinical signs, the disease is very diverse in terms of the nature and sequence of the underlying molecular events. Multiple cellular processes are involved in helping the malignant cells to remain viable and maintain proliferative properties in the hypoxic microenvironment of the bone marrow. Specifically, the process of angiogenesis, triggered by the interactions between the malignant MM cells and the stroma cells around them, was found to be critical for MM progression. In this review we highlight the current understanding about the epigenetic regulation of the proliferation and apoptosis of MM cells and its dependency on angiogenesis in the bone marrow that is carried out by different microRNAs.

miR-29s: a family of epi-miRNAs with therapeutic implications in hematologic malignancies

A wealth of studies has highlighted the biological complexity of hematologic malignancies and the role of dysregulated signal transduction pathways. Along with the crucial role of genetic abnormalities, epigenetic aberrations are nowadays emerging as relevant players in cancer development, and significant research efforts are currently focusing on mechanisms by which histone post-translational modifications, DNA methylation and noncoding RNAs contribute to the pathobiology of cancer. As a consequence, these studies have provided the rationale for the development of epigenetic drugs, such as histone deacetylase inhibitors and demethylating compounds, some of which are currently in advanced phase of pre-clinical investigation or in clinical trials. In addition, a more recent body of evidence indicates that microRNAs (miRNAs) might target effectors of the epigenetic machinery, which are aberrantly expressed or active in cancers, thus reverting those epigenetic abnormalities driving tumor initiation and progression. This review will focus on the broad epigenetic activity triggered by members of the miR-29 family, which underlines the potential of miR-29s as candidate epi-therapeutics for the treatment of hematologic malignancies.

MicroRNA Regulation of Epithelial to Mesenchymal Transition

Epithelial to mesenchymal transition (EMT) is a central regulatory program that is similar in many aspects to several steps of embryonic morphogenesis. In addition to its physiological role in tissue repair and wound healing, EMT contributes to chemo resistance, metastatic dissemination and fibrosis, amongst others. Classically, the morphological change from epithelial to mesenchymal phenotype is characterized by the appearance or loss of a group of proteins which have come to be recognized as markers of the EMT process. As with all proteins, these molecules are controlled at the transcriptional and translational level by transcription factors and microRNAs, respectively. A group of developmental transcription factors form the backbone of the EMT cascade and a large body of evidence shows that microRNAs are heavily involved in the successful coordination of mesenchymal transformation and vice versa, either by suppressing the expression of different groups of transcription factors, or otherwise acting as their functional mediators in orchestrating EMT. This article dissects the contribution of microRNAs to EMT and analyzes the molecular basis for their roles in this cellular process. Here, we emphasize their interaction with core transcription factors like the zinc finger enhancer (E)-box binding homeobox (ZEB), Snail and Twist families as well as some pluripotency transcription factors.

Diverse microRNAs with convergent functions regulate tumorigenesis

MicroRNAs (miRNAs) regulate several biological processes, including tumorigenesis. In order to comprehend the roles of miRNAs in cancer, various screens were performed to investigate the changes in the expression levels of miRNAs that occur in different types of cancer. The present review focuses on the results of five recent screens, whereby a number of overlapping miRNAs were identified to be downregulated or differentially regulated, whereas no miRNAs were observed to be frequently upregulated. Furthermore, the majority of the miRNAs that were common to >1 screen were involved in signaling networks, including wingless-related integration site, receptor tyrosine kinase and transforming growth factor-β, or in cell cycle checkpoint control. The present review will discuss the aforementioned miRNAs implicated in cell cycle checkpoint control and signaling networks.

Integrative Analysis with Monte Carlo Cross-Validation Reveals miRNAs Regulating Pathways Cross-Talk in Aggressive Breast Cancer

In this work an integrated approach was used to identify functional miRNAs regulating gene pathway cross-talk in breast cancer (BC). We first integrated gene expression profiles and biological pathway information to explore the underlying associations between genes differently expressed among normal and BC samples and pathways enriched from these genes. For each pair of pathways, a score was derived from the distribution of gene expression levels by quantifying their pathway cross-talk. Random forest classification allowed the identification of pairs of pathways with high cross-talk. We assessed miRNAs regulating the identified gene pathways by a mutual information analysis. A Fisher test was applied to demonstrate their significance in the regulated pathways. Our results suggest interesting networks of pathways that could be key regulatory of target genes in BC, including stem cell pluripotency, coagulation, and hypoxia pathways and miRNAs that control these networks could be potential biomarkers for diagnostic, prognostic, and therapeutic development in BC. This work shows that standard methods of predicting normal and tumor classes such as differentially expressed miRNAs or transcription factors could lose intrinsic features; instead our approach revealed the responsible molecules of the disease.