Next Article in Journal
Associations between Adverse Childhood Experiences (ACEs) and Prenatal Mental Health and Substance Use
Previous Article in Journal
Syndromic Surveillance of Health Effects Due to Summer Sun Overexposure: Construction of an Indicator Based on Drug Sales in Pharmacies—Preliminary Study within the PRISME Project
Previous Article in Special Issue
COVID-19 Classification on Chest X-ray Images Using Deep Learning Methods
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Soft Tissue Ewing Sarcoma Cell Drug Resistance Revisited: A Systems Biology Approach

by
Seyedehsadaf Asfa
1,2,
Halil Ibrahim Toy
1,2,
Reza Arshinchi Bonab
1,2,
George P. Chrousos
3,4,
Athanasia Pavlopoulou
1,2,* and
Styliani A. Geronikolou
3,4,*
1
Izmir Biomedicine and Genome Center (IBG), 35340 Izmir, Turkey
2
Izmir International Biomedicine and Genome Institute, Dokuz Eylül University, 35340 Izmir, Turkey
3
Clinical, Translational and Experimental Surgery Research Centre, Biomedical Research Foundation Academy of Athens, Soranou Ephessiou 4, 11527 Athens, Greece
4
University Research Institute of Maternal and Child Health and Precision Medicine and UNESCO Chair on Adolescent Health Care, National and Kapodistrian University of Athens, Aghia Sophia Children’s Hospital, Levadeias 8, 11527 Athens, Greece
*
Authors to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2023, 20(13), 6288; https://doi.org/10.3390/ijerph20136288
Submission received: 9 March 2023 / Revised: 8 May 2023 / Accepted: 30 June 2023 / Published: 3 July 2023

Abstract

:
Ewing sarcoma is a rare type of cancer that develops in the bones and soft tissues. Drug therapy represents an extensively used modality for the treatment of sarcomas. However, cancer cells tend to develop resistance to antineoplastic agents, thereby posing a major barrier in treatment effectiveness. Thus, there is a need to uncover the molecular mechanisms underlying chemoresistance in sarcomas and, hence, to enhance the anticancer treatment outcome. In this study, a differential gene expression analysis was conducted on high-throughput transcriptomic data of chemoresistant versus chemoresponsive Ewing sarcoma cells. By applying functional enrichment analysis and protein–protein interactions on the differentially expressed genes and their corresponding products, we uncovered genes with a hub role in drug resistance. Granted that non-coding RNA epigenetic regulators play a pivotal role in chemotherapy by targeting genes associated with drug response, we investigated the non-coding RNA molecules that potentially regulate the expression of the detected chemoresistance genes. Of particular importance, some chemoresistance-relevant genes were associated with the autonomic nervous system, suggesting the involvement of the latter in the drug response. The findings of this study could be taken into consideration in the clinical setting for the accurate assessment of drug response in sarcoma patients and the application of tailored therapeutic strategies.

1. Introduction

Sarcomas represent one-fifth of pediatric cancers and 1% of adult solid malignant tumors. They rise most frequently from connective tissues and only in one-tenth of cases from osseus tissues. Based on their tissue of origin, there are two main classes of sarcomas. The pathogenesis of sarcomas is still unclear, as amid established risk factors reside genetic and epigenetic causes. Literature evidence includes smoking, age, gestational age and weight, parental and maternal health status, occupational exposure to drugs, chemicals, etc. Unfortunately, the survival prognosis of sarcomas is rather poor [1].
At least seventy types of sarcomas have been described so far [2]. Ewing sarcoma is a rare type that develops in bones (skull, spine, pelvis, chest, legs, arms, feet) and soft tissues (head and neck, legs, retroperitoneum) equally frequently and within 13–20 years [3]. More specifically, the Ewing sarcoma family of tumors includes neoplasms localized in the chest wall; this type of peripheral primitive neuroectodermal cancer has been termed an Askin tumor. Other localized sarcomas include an extraosseous Ewing sarcoma (rising in tissues other than bone). Medical record of hernias strongly contributes to Ewing sarcoma [4,5,6].
The main symptomatology of sarcomas includes swelling and pain around the tumor; these symptoms may coincide with febrile episodes and unexplained bone breaks, while their treatment is designed after gene studies in tumor tissue. Treatment involves surgery, palliative care, chemotherapy options, stem cells, CAR-T cells, targeted therapy with a NEDD-8 inhibitor or a tyrosine kinase inhibitor, immunotherapy with immune checkpoint inhibitors, radiotherapy, etc. Chemotherapy is a widely used modality for the effective treatment of diverse types of cancers, including sarcomas. However, in many cases, cancer cells acquire resistance to chemotherapy, which poses a major problem as to the efficacy of the treatment [7,8], as described in detail by Nikolaou and colleagues (2018) [7,8].
Differential gene expression profiles are important indicators of cell resistance to anticancer treatments [9,10,11,12]. The differential expression of genes is often associated with epigenetic regulation by DNA methylation, histone modifications, or non-coding RNA species (ncRNAs), which affect gene expression without modifying the underlying DNA sequence [13,14,15]. There is compelling evidence that gene expression is regulated by ncRNA molecules at different levels [16,17,18]. MicroRNAs (miRNAs) are endogenous, small (~22 nucleotides in length), non-protein-coding RNAs that play a critical role in regulating the expression of target mRNAs at the post-transcriptional level by binding to complementary sequences on target mRNAs, the so-called miRNA response elements (MREs) [19]. MiRNAs reduce the stability of the target mRNAs and/or inhibit their translation, thereby downregulating the expression of the corresponding genes [20]. A single miRNA can potentially regulate the expression of multiple genes, and, conversely, a gene can be targeted by numerous miRNAs [21]. MiRNAs have been reported to play important regulatory roles in cancer drug resistance [22,23,24].
Lately, accumulating evidence points towards the key role of competing endogenous RNAs (ceRNAs) in miRNA-mediated gene regulation. CeRNAs are ncRNAs that share MREs with miRNAs, and, therefore, can sequester miRNAs (acting like ‘sponges’), preventing thus miRNAs from binding to their MREs and reducing their regulatory effect on target mRNAs [25,26,27].
Post-transcriptional miRNA-mediated crosstalk between long ncRNAs (>200 bp long) (lncRNAs) and mRNAs has also been reported in chemoresistance [28,29]. The lncRNAs can be capped, polyadenylated, and spliced, but they lack a functional open reading frame. These versatile molecules are implicated in a range of biological functions and cellular processes, including regulation of gene expression, cell–cell signaling, genomic instability, and RNA decay [17,18]. Notably, lncRNAs are largely involved in diverse cellular processes [17], including response to chemotherapies [30,31]. In this study, an integrated bioinformatics approach was applied to identify competing endogenous RNA networks—i.e., the miRNAs that regulate the expression of chemoresistance-related genes, and the lncRNAs that act as molecular sponges of those miRNAs.
The above evidence highlights the importance of identifying genetic and epigenetic biomarkers for chemoresistance in sarcomas. To this end, we have designed an analysis of ‘-omics’ data relevant to the response of A673 Ewing sarcoma cells to drug treatment, towards the identification of novel genetic biomarkers of drug response in sarcomas, as well as pivotal epigenetic regulators.

2. Materials and Methods

2.1. High-Throughput Gene Expression Data

The publicly available repository NCBI GEO (Gene Expression Omnibus) DataSets (https://www.ncbi.nlm.nih.gov/gds/; accessed on 23 October 2022) [32,33] was thoroughly searched for gene expression datasets related to sarcomas and drug treatment using the keywords: “drug” AND “sarcoma” AND (“human” or “homo sapiens”). The criteria for selecting datasets were: (i) gene expression data from treated and untreated human sarcoma tissues or cell lines, (ii) more than 5000 genes were included in the transcriptomic dataset. In this way, one eligible RNA-Seq dataset was obtained. The GEO series GSE118871 [34] (Table S1) includes the genome-wide gene expression of Ewing sarcoma A673 cells treated with SP-2509. This dataset contains cell lines both responsive and resistant to SP-2509. The drug-resistant cell lines (herein referred to as “resistant”) were established by exposing the corresponding parental A673 cells (untreated) to increasing concentrations of SP-2509 over a 7 month period (in increments of 100 nM) or 48 h (2 uM). Both the resistant cells and those responsive to drug treatment (hereinafter called “responsive”) were assessed by Pishas et al. [34] based on experimental cell viability assays; the resistant cells demonstrated a significantly increased viability as compared to parental A673 cells treated with SP-2509, whereas the responsive cells displayed a reduced viability following drug treatment (Table S1). The Illumina HiSeq 2000 (Homo sapiens) GPL11154 platform was employed.

2.2. RNA-Seq Data Processing and Identification of Differentially Expressed Genes

The FASTQ files that contain raw 2 × 50-bp paired-end RNA-Seq reads were downloaded from the respective Sequence Read Archive (SRA) using the SRA Tool Kit v.2.9.0 [35] with the fastq-dump –gzip –skip-technical –readids –dumpbase –clip –split-3 command. Raw RNA-Seq reads were mapped to the human reference genome GRCh38 (Ensembl version 104) with the usage of the splice junction aligner HISAT2 v.2.1.0 [36] with the “hisat2 -p -dta -x {input.index} -U {input.fq} -S {out.sam}” parameters. The output SAM files were converted to the compressed binary BAM file with the usage of SAMtools v.1.14 [37] with the “samtools sort -@ 10 -o {output.bam} {input.sam}” commands. Transcriptome normalization, reconstruction, and quantification was conducted by employing the StringTie version 1.3.5 [38], with the “stringtie -e -B -p -G {input.gtf} -A {output.tab} -o {output.gtf} -l {input.label}{input.bam}” parameters. The assembled transcripts and their estimated abundances were included in the output GTF file.
To identify DEGs between resistant versus responsive, and responsive versus parental Ewing sarcoma cells, the EdgeR package version 3.32.0 [39] of the R statistical computation environment v.3.6.1 (https://www.r-project.org; accessed on 28 October 2022) was used. The negative binomial (NB) distribution was used to model the RNA-Seq read counts per gene per sample in EdgeR. Then, the estimating dispersion was calculated with the estimateDisp function. Differential expression analysis between the two RNA-Seq groups was performed using the exactTest function of the EdgeR package v3.32.0. For determining statistically significant DEGs, the cutoff for the absolute log2-fold change (FC) was set at two (|log2FC≥2|), and the Benjamini and Hochberg (BH)-adjusted p-value [40] was ≤0.05.
The official HUGO Gene Nomenclature Committee (HGNC) [41] gene symbols and gene names were used.

2.3. Gene Set Enrichment Analysis

Gene set enrichment analysis (GSEA), a method to identify biological terms that are enriched in a large gene set, was conducted to functionally annotate the drug resistance-associated genes detected in this study. To this end, the ‘resistant’ DEGs were provided as input to WebGestalt (WEB-based GEne SeT AnaLysis Toolkit) [42,43] to identify statistically significant over-represented Gene Ontology (GO) terms; the non-redundant Biological Process subontology of GO was selected, the threshold for the BH-corrected p-value [40] was set at 10−3; the hypergeometric distribution was applied.

2.4. Functional Association Network

The associations among the ‘resistant’ genes/proteins under study were investigated and visualized with the usage of the STRING (Search Tool for Retrieval of Interacting Genes/Proteins) v11.5 [44], a database of both experimentally derived or predicted, direct or indirect, association data among genes/proteins extracted from diverse resources. A relatively high confidence score (≥0.6) for displaying interactions was chosen. The associations were further investigated, analyzed, and visualized through the open-source platform Cytoscape (v.3.8.2) (https://cytoscape.org/; accessed on 11 December 2022) [45].

2.5. Epigenetic Regulators of Chemoresistant Genes

The ncRNAs—i.e., miRNAs and lncRNAs—likely regulating the ‘resistant’ DEGs under study were investigated by employing state-of-the-art software tools.
Both the experimentally supported and predicted miRNAs targeting the DEGs were collected with the usage of: (i) microT_CDS (http://diana.imis.athena-innovation.gr/DianaTools/index.php?r=microT_CDS/index; accessed on 7 January 2022) has implemented the target prediction algorithm DIANA-microT-CDS [46]; (ii) TargetScan (http://www.targetscan.org/vert_72/; accessed on 7 January 2022) searches for the presence of conserved 8mer, 7mer, and 6mer sites that match the seed region of each miRNA [47]; (iii) miRDB (http://mirdb.org/; accessed on 8 January 2022) predicts functional miRNA gene targets by machine learning modeling of target gene expression data [48,49]; (iv) PicTar (https://pictar.mdc-berlin.de/; accessed on 8 January 2022) includes published, experimentally validated miRNA targets [50]. The miRNA/mRNA associations from the four databases were integrated and the duplicates were removed.
The DIANA-LncBase v.3 (https://diana.e-ce.uth.gr/lncbasev3; accessed on 24 January 2022) [51], a comprehensive database of experimentally supported miRNA targets on non-coding transcripts, was employed to identify those lncRNAs that potentially interact with the miRNAs detected in the previous step. To increase the accuracy of our analysis, only the miRNA–lncRNA interactions with high confidence detected in cancer/malignant cell types were considered.
An in-house Python script was used to retrieve miRNA–mRNA and miRNA–lncRNA interactions from the respective databases.

3. Results

The overall procedure followed in this study is outlined in Figure 1.

3.1. Gene Expression Profiles of Ewing Sarcoma Drug Resistance

The differentially expressed genes (DEGs) detected between the chemoresistant and chemoresponsive, as well as the chemoresponsive versus the untreated parental Ewing sarcoma cells [34], were 1601 and 1196, respectively (Table S1). Among the genes found upregulated in the drug resistant cells were those encoding “drug resistance-associated membrane proteins” or “DRAMPs”, which affected the transport of drugs across membranes. The ABCA2/8/9 genes belong to the broad ATP-binding cassette (ABC) transporter superfamily, the members of which pump drug molecules out of the cell, thereby decreasing the net accumulation of drugs within cancer cells. ABCA2 and ABCA9 were down-regulated in the responsive cells. The expression of a great number of genes coding for another class of DRAMPs, the solute carrier (SLC) transporters, which interfere with the translocation of drug molecules across membranes relying on physicochemical processes [52,53], was also dysregulated (Table S1). A key player in angiogenic signaling in cancer, the vascular epithelial growth factor, VEGFB, was overexpressed in resistant cells (Table S1). Treatment with specific VEGF inhibitors results in transient tumor vascular normalization and increased cancer cell response to chemotherapeutic drugs [54]. An increased expression of ALDH1A3 and ALDH1B1—encoding aldehyde dehydrogenases (ALDH), detoxifying enzymes that catalyze the oxidation of intracellular aldehydes—was detected in the chemoresistant cell lines; ALDHs have been proposed as biomarkers of chemoresistance [55]. ALDH1B1 was also under-expressed in the chemoresponsive cells.
Of note, non-common DEGs (with the same direction, either up- or down-regulated) were detected when the resistance versus responsiveness and responsiveness versus parental (untreated) cells were compared (Table S1). These findings indicate that our results were biologically meaningful, as different molecular factors and mechanisms were implicated in these two phenomena.

3.2. Functional Annotation Analysis

A total of 641 ‘chemoresistant’ genes products were found to be implicated in known biological processes (Figure 2 and Table S2). The protein products of 369 (out of 641) genes formed a highly interconnected network (Figure 3 and Table S1), suggesting a functional association among these genes leading to their drug resistance effect. One of the over-represented biological processes was ‘drug transport’, which included 27 genes, 22 of which were up-regulated in the drug resistant cells (Figure 2 and Table S1). Given that an increased drug transport activity plays a critical role in drug resistance [7,56] and the bioentities that participate in networks are of a higher biological significance [57], we selected these eight over-expressed genes (DRD1, DRD2, GHRL, KCNA2, SLC7A10, SLC25A13, STRA6, SYT2), the protein products of which participate in the constructed network (Figure 3). Of note, the genes DRD2, KCNA2, SLC7A10, STRA6, and SYT2 were under-expressed in the chemoresponsive versus untreated parental cells (Table S1).

3.3. Epigenetic Regulatory Network of Drug Resistance-Relevant Genes

To gain a glimpse into the post-transcriptional regulatory mechanism(s) of the eight SP-2509-resistant genes, a ceRNA network was constructed, consisting of the miRNAs that regulate the expression of the drug resistance genes, and the lncRNAs, capable of acting as molecular sponges of these miRNAs. For instance, several studies suggest that the downregulation of cancer-relevant miRNAs through sponging alleviates the suppression of downstream mRNAs, affecting different aspects of carcinogenesis [58,59,60].
The miRNAs that potentially target the eight genes in Figure 3 were predicted using different methods. The computational methods for miRNA/mRNA target prediction usually depend on sequence-based features, thermodynamic stability, evolutionary information, or probabilistic models, etc. [61,62,63]. In our study, we used four tools based on different algorithms, so as to extract the maximum possible information. To enhance the prediction accuracy, only those miRNA–target gene interactions predicted by more than three tools were included in the study. Moreover, to obtain more robust results, miRNAs targeting more than two genes were selected for analysis. Collectively, 30 miRNAs were found to target more than 2 genes, and, conversely, 5 genes were targeted by more than 2 miRNAs (DRD1/2, SLC25A13, STRA6, SYT2), suggesting possible co-regulation at the post-transcriptional level. Of these miRNAs, based on our stringent selection criteria, 8 (hsa-miR-149-5p, hsa-miR-29a/b/c-3p, hsa-miR-330-5p, hsa-miR-501-3p, hsa-miR-760 and hsa-miR-766-3p) interacted with the lncRNAs in DIANA-LncBase (Table S2). We retained only the lncRNAs that were upregulated in the drug resistant cells—that is, co-upregulated with the target genes (Tables S1 and S2). Subsequently, a competing endogenous RNA network of the selected six lncRNAs (EDRF1-DT, HAGLR, LINC00997, LOXL1-AS1, SRRM2-AS1, TMPO-AS1), their sponged miRNAs, and the genes targeted by the miRNAs was constructed; however, the mirnas that interact with the DRD2 were not sponged by any lncRNAs (Figure 4).
The three members of the mir-29 family, miR-29a/b/c, were the ones mostly targeted by the lncRNAs (Figure 4). In particular, miR-29s have been reported to act primarily as tumor suppressors in numerous cancers through the upregulation of tumor suppressor genes or/and downregulation of oncogenes; they can, thus, elicit apoptosis and inhibit invasion and proliferation of cancer cells [64,65]. In a similar manner, they suppressed the activity of the SP-2509 resistance genes, resulting in the Ewing sarcoma cells responding to drugs. Conversely, the lncRNAs identified in this study could act as drivers of chemoresistance in sarcomas through their sponging role regarding miRNAs (Figure 4). In agreement with this, accumulating evidence suggests that LOXL1-AS1 [66,67,68,69], HAGLR [70,71], LINC00997 [72], and TMPO-AS1 [73,74,75] contribute to carcinogenesis by acting as molecular sponges of miRNAs.

4. Discussion

In this work, we investigated the molecular determinants of the Ewing sarcoma cell’s line A673 response to diverse anti-neoplastic agents through an in silico approach. To this end, a relevant, publicly available gene expression dataset was used, in which Pishas and colleagues [34] treated Ewing sarcoma cells with increasing concentrations of SP-2509 to investigate the mechanisms underlying the resistance to SP-2509, a small molecular reversible inhibitor of LSD1 [76], in sarcoma. The lysine-specific demethylase 1 (LSD1), also known as KDM1A, demethylates histone 3 lysine 4 (H3K4), thereby acting either as a transcriptional co-repressor or as a transcriptional co-activator by catalyzing the demethylation of histone 3 lysine 9 (H3K9) [77,78]. LSD1 has been demonstrated to contribute to carcinogenesis via chromatin modification, as it promotes or represses the transcription of oncogenes or tumor suppressor genes, respectively [79,80,81,82]. The majority (76%) of the genes found to be differentially expressed between the SP-2509 resistant and responsive Ewing sarcoma cells were up-regulated (Table S1), suggesting that the overexpressed genes contributed the most to drug resistance. Of those, eight genes were implicated in drug transport and their corresponding protein products were components of a highly connected network (Figure 2 and Figure 3). These genes/proteins included the signaling receptor and transporter of retinol STRA6, the synaptotagmin 2, the potassium voltage-gated channel KCNA2, as well as the receptor ligands ghrelin and obestatin prepropeptide GHRL. The dopamine receptors DRD1 and DRD2, a class of G protein-coupled receptors, activated adenylyl cyclase to convert ATP into cyclic AMP (cAMP), a second messenger, and were implicated in the dopaminergic regulation of essential neurophysiological processes [83,84]. Two solute carrier (SLC) transporters, SLC25A13 and SLC7A10, were also upregulated. SLC transporters facilitated the influx of drug molecules into cells [85]. However, it has also been suggested that several SLC transporters can mediate bi-directional transport (i.e., both influx and efflux) [86].

Autonomic Nervous System and Drug Resistance

The autonomic nervous system (ANS) regulates many bodily functions, participating in a major fashion in the maintenance of homeostasis and the body’s response to stress, and may, thus, influence carcinogenesis. On the other hand, members of the dopamine receptor (DR) family are upregulated in diverse cancers. A positive correlation was observed between a reduced cancer risk and neurodegenerative disorders, such as Parkinson’s disease or schizophrenia, where DR-targeting drugs are administered. Moreover, DR antagonists have displayed anticancer efficacy [87,88]. Amid the nodes involved in our network, SYT2 has been suggested to play a regulatory role in synaptic vesicle trafficking and in promoting metastasis in ovarian cancer [89]. In addition, polymorphisms in the gene encoding the neuropeptide GHRL have been associated with non-Hodgkin lymphoma [90]. Three more of the identified major hubs in our constructed interactome: (1) the ion transporter KCNA2, normally expressed in the central nervous system (CNS) [91]; (2) the STRA6, the expression of which was reported in the differentiating nervous system [92]; and (3), the SLC7A10, also called the astrocytic transporter (Asc-1), which has been proposed as a primary driver of the D-serine uptake in the CNS [93]. Taken together, the above findings suggested a link between drug resistance and nervous system function.
Of note, in a study by Gaynes et al., the CNS niche was shown to enhance chemoresistance in leukemia [94]. Additionally, denervation enhanced the effectiveness of chemotherapy in gastric cancer [95]. Finally, Logotheti and colleagues had suggested that genes implicated in neuronal function are activated in cancer cells [96]. As little is known on the potential effects of the ANS on cancer chemoresistance, we compared the up-regulated drug resistance-associated molecules identified herein (Table S1) with murine knockout genes associated with abnormalities in neuronal development, as reported by Logotheti et al. [96]. Of the 119 molecules found in common (Table 1), 40 were oncogenes, including 2 Ewing sarcoma specific genes—namely, the NK2 Homeobox 2 (NKX2) and FEV transcription factor (FEV) [97,98,99]. It has been reported that the NKX2 expression as a biomarker has a high sensitivity (100%), but a moderate specificity in cytologic specimens [97]. Yet, the Ewing sarcoma specificity increased when combined with CD99 [100]. The NKX2 is selectively expressed in the brain, thyroid gland, parathyroid glands, lungs, skin, and enteric ganglia, from prenatal development to adulthood, and has key functions at the interface of the brain, the endocrine, and the immune systems. It is highly expressed in mature limbic circuits related to context-dependent goal-directed patterns of behavior, social interaction and reproduction, as well as in fear responses [101].
The Nkx2-1 is indirectly involved in autonomic regulation as a major central target for angiotensinogen produced in the subfornical organ and directed to the hypothalamic paraventricular nucleus, where it is converted to angiotensin II [101]. Angiotensin II in the paraventricular nucleus contributes to the autonomic output and sympathetic nervous system excitation [102] and the so termed “enteric branch” (third autonomic branch beyond the sympathetic and parasympathetic ones) of the ANS, as well [103,104,105,106]. The peculiar behavior of the latter branch was initially described by Sternini (1997) [107]. Its association to obesity was first observed in animal models, followed by a meta-analysis in human bariatric surgery data [108], and then described in a preliminary molecular network involved in the “third type of obesity” [109]. Yet, concerning cancer, although this protein has been proposed as a cytostatic factor in certain cancer types, this has not been confirmed on a large scale [110].
FEV gene rearrangements have been observed in a fraction of Ewing sarcoma patients (3.5%) being linked to axial extraskeletal locations (mainly soft tissue), “older age at diagnosis and aggressive clinical behavior”, but retaining the classic Ewing sarcoma image [99].
Most of the rest are molecules associated with sensory (deafness, eye disorders) and motor impairments, cardiovascular diseases (i.e., valvulopathies, Brugada syndrome, hypertension), metabolic disorders, neurodegenerative disorders (Parkinson’s, Alzheimer, myoskeletal disabilities), encepalopathies and various syndromes, as shown in Table 1. In this table, some molecules are of unidentified function, meriting further research.
Furthermore, five genes involved in drug transport (DRD1/2, SLC7A10, STRA6, and SYT2) and the SLC25A13 paralog, SLC25A1, (Figure 3), are related to neuronal development (Table 1, italicized), further supporting the emerging role of the ANS in drug resistance in chemotherapy.

5. Conclusions

By employing an interdisciplinary in silico methodology, we identified genes involved in drug resistance in sarcoma cells. A non-negligible fraction of these genes was also associated with the nervous system function and/or development. Of note, the study of the chemoresistance modification effect of the autonomic nervous system in cancer is an ongoing research field with interesting and unpredictable results. Furthermore, this study suggested that future investigations should be directed towards deciphering the crosstalk between the ceRNA and epigenetic regulation of interconnected protein-coding genes. LncRNAs may be considered potential targets for drug resistance based on their capability to attract and inhibit downstream miRNAs, influencing their ability to suppress chemoresistance—associated mRNAs. Both the protein–protein interaction and ceRNA networks comprise of an orchestrated multi-molecular mechanism, the dissection of which would enhance our understanding of the molecular determinants of cancer cell chemoresistance. Therefore, by disrupting the lncRNAs–miRNA interplay and by repressing the corresponding mRNAs, as well as by manipulating interactions between proteins encoded by chemoresistance-relevant genes, the drug resistance of cancer cells could be attenuated. This information can be utilized in the clinical setting for the design of combinatorial therapies, targeting both proteins and epigenetic regulatory factors, such as miRNAs and lncRNAs.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijerph20136288/s1, Table S1: Gene expression data analyzed in this study. Samples from the GSE118871 transcriptomic dataset; differentially expressed genes (resistant versus responsive, responsive versus parental); Table S2: Functional analysis of the Ewing chemoresistant genes.

Author Contributions

Conceptualization, A.P. and S.A.G.; methodology, S.A., H.I.T., R.A.B. and A.P.; software, S.A. and H.I.T.; validation, S.A., H.I.T., A.P. and S.A.G.; formal analysis, S.A., H.I.T., R.A.B., A.P. and S.A.G.; investigation, S.A., H.I.T., A.P. and S.A.G.; resources, S.A., H.I.T., R.A.B. and A.P.; data curation, S.A., H.I.T. and A.P.; visualization, S.A., G.P.C. and R.A.B.; supervision, A.P. and S.A.G.; project administration, A.P. and S.A.G; writing—original draft preparation, all authors; writing—review and editing, all authors. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data and analysis methodologies are contained in the manuscript. Any additional data requests can be addressed to the corresponding authors.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yang, J.; Ren, Z.; Du, X.; Hao, M.; Zhou, W. The role of mesenchymal stem/progenitor cells in sarcoma: Update and dispute. Stem Cell Investig. 2014, 1, 18. [Google Scholar] [CrossRef]
  2. Burningham, Z.; Hashibe, M.; Spector, L.; Schiffman, J.D. The Epidemiology of Sarcoma. Clin. Sarcoma Res. 2012, 2, 14. [Google Scholar] [CrossRef] [Green Version]
  3. Khan, S.; Abid, Z.; Haider, G.; Bukhari, N.; Zehra, D.; Hashmi, M.; Abid, M.; Ibrahim, U. Incidence of Ewing’s Sarcoma in Different Age Groups, Their Associated Features, and Its Correlation with Primary Care Interval. Cureus 2021, 13, e13986. [Google Scholar] [CrossRef] [PubMed]
  4. WHO. Undifferentiated Small Round Cell Sarcoma of Bone and Soft Tissue: Ewing Sarcoma; International Agency for Research on Cancer: Lyon, France, 2020; Volume 3.
  5. Valery, P.C.; Holly, E.A.; Sleigh, A.C.; Williams, G.; Kreiger, N.; Bain, C. Hernias and Ewing’s sarcoma family of tumours: A pooled analysis and meta-analysis. Lancet Oncol. 2005, 6, 485–490. [Google Scholar] [CrossRef] [PubMed]
  6. Cope, J.U.; Tsokos, M.; Helman, L.J.; Gridley, G.; Tucker, M.A. Inguinal hernia in patients with Ewing sarcoma: A clue to etiology. Med. Pediatr. Oncol. 2000, 34, 195–199. [Google Scholar] [CrossRef]
  7. Alfarouk, K.O.; Stock, C.M.; Taylor, S.; Walsh, M.; Muddathir, A.K.; Verduzco, D.; Bashir, A.H.; Mohammed, O.Y.; Elhassan, G.O.; Harguindey, S.; et al. Resistance to cancer chemotherapy: Failure in drug response from ADME to P-gp. Cancer Cell Int. 2015, 15, 71. [Google Scholar] [CrossRef] [Green Version]
  8. Nikolaou, M.; Pavlopoulou, A.; Georgakilas, A.G.; Kyrodimos, E. The challenge of drug resistance in cancer treatment: A current overview. Clin. Exp. Metastasis 2018, 35, 309–318. [Google Scholar] [CrossRef]
  9. Gyorffy, B.; Surowiak, P.; Kiesslich, O.; Denkert, C.; Schafer, R.; Dietel, M.; Lage, H. Gene expression profiling of 30 cancer cell lines predicts resistance towards 11 anticancer drugs at clinically achieved concentrations. Int. J. Cancer 2006, 118, 1699–1712. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Kang, H.C.; Kim, I.J.; Park, J.H.; Shin, Y.; Ku, J.L.; Jung, M.S.; Yoo, B.C.; Kim, H.K.; Park, J.G. Identification of genes with differential expression in acquired drug-resistant gastric cancer cells using high-density oligonucleotide microarrays. Clin. Cancer Res. Off. J. Am. Assoc. Cancer Res. 2004, 10, 272–284. [Google Scholar] [CrossRef] [Green Version]
  11. Kim, S.C.; Shin, Y.K.; Kim, Y.A.; Jang, S.G.; Ku, J.L. Identification of genes inducing resistance to ionizing radiation in human rectal cancer cell lines: Re-sensitization of radio-resistant rectal cancer cells through down regulating NDRG1. BMC Cancer 2018, 18, 594. [Google Scholar] [CrossRef]
  12. Serafim, R.B.; da Silva, P.; Cardoso, C.; Di Cristofaro, L.F.M.; Netto, R.P.; de Almeida, R.; Navegante, G.; Storti, C.B.; de Sousa, J.F.; de Souza, F.C.; et al. Expression Profiling of Glioblastoma Cell Lines Reveals Novel Extracellular Matrix-Receptor Genes Correlated with the Responsiveness of Glioma Patients to Ionizing Radiation. Front. Oncol. 2021, 11, 668090. [Google Scholar] [CrossRef]
  13. Gibney, E.R.; Nolan, C.M. Epigenetics and gene expression. Heredity 2010, 105, 4–13. [Google Scholar] [CrossRef] [Green Version]
  14. Golbabapour, S.; Abdulla, M.A.; Hajrezaei, M. A concise review on epigenetic regulation: Insight into molecular mechanisms. Int. J. Mol. Sci. 2011, 12, 8661–8694. [Google Scholar] [CrossRef]
  15. Holoch, D.; Moazed, D. RNA-mediated epigenetic regulation of gene expression. Nat. Rev. Genet. 2015, 16, 71–84. [Google Scholar] [CrossRef] [Green Version]
  16. Kaikkonen, M.U.; Lam, M.T.; Glass, C.K. Non-coding RNAs as regulators of gene expression and epigenetics. Cardiovasc. Res. 2011, 90, 430–440. [Google Scholar] [CrossRef] [Green Version]
  17. Statello, L.; Guo, C.J.; Chen, L.L.; Huarte, M. Gene regulation by long non-coding RNAs and its biological functions. Nat. Rev. Mol. Cell Biol. 2021, 22, 96–118. [Google Scholar] [CrossRef] [PubMed]
  18. Zhang, X.; Wang, W.; Zhu, W.; Dong, J.; Cheng, Y.; Yin, Z.; Shen, F. Mechanisms and Functions of Long Non-Coding RNAs at Multiple Regulatory Levels. Int. J. Mol. Sci. 2019, 20, 5573. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. O’Brien, J.; Hayder, H.; Zayed, Y.; Peng, C. Overview of MicroRNA Biogenesis, Mechanisms of Actions, and Circulation. Front. Endocrinol. 2018, 9, 402. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Ambros, V. The functions of animal microRNAs. Nature 2004, 431, 350–355. [Google Scholar] [CrossRef] [PubMed]
  21. Lewis, B.P.; Burge, C.B.; Bartel, D.P. Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell 2005, 120, 15–20. [Google Scholar] [CrossRef] [Green Version]
  22. Si, W.; Shen, J.; Zheng, H.; Fan, W. The role and mechanisms of action of microRNAs in cancer drug resistance. Clin. Epigenet. 2019, 11, 25. [Google Scholar] [CrossRef] [PubMed]
  23. Allegra, A.; Ettari, R.; Innao, V.; Bitto, A. Potential Role of microRNAs in inducing Drug Resistance in Patients with Multiple Myeloma. Cells 2021, 10, 448. [Google Scholar] [CrossRef] [PubMed]
  24. Hu, X.Y.; Song, Z.; Yang, Z.W.; Li, J.J.; Liu, J.; Wang, H.S. Cancer drug resistance related microRNAs: Recent advances in detection methods. Analyst 2022, 147, 2615–2632. [Google Scholar] [CrossRef] [PubMed]
  25. Beylerli, O.; Gareev, I.; Sufianov, A.; Ilyasova, T.; Guang, Y. Long noncoding RNAs as promising biomarkers in cancer. Non-Coding RNA Res. 2022, 7, 66–70. [Google Scholar] [CrossRef] [PubMed]
  26. Chen, L.L. The expanding regulatory mechanisms and cellular functions of circular RNAs. Nat. Rev. Mol. Cell Biol. 2020, 21, 475–490. [Google Scholar] [CrossRef]
  27. Diener, C.; Keller, A.; Meese, E. Emerging concepts of miRNA therapeutics: From cells to clinic. Trends Genet. TIG 2022, 38, 613–626. [Google Scholar] [CrossRef]
  28. Kong, X.; Hu, S.; Yuan, Y.; Du, Y.; Zhu, Z.; Song, Z.; Lu, S.; Zhao, C.; Yan, D. Analysis of lncRNA, miRNA and mRNA-associated ceRNA networks and identification of potential drug targets for drug-resistant non-small cell lung cancer. J. Cancer 2020, 11, 3357–3368. [Google Scholar] [CrossRef] [Green Version]
  29. Liu, H.; Wang, S.; Zhou, S.; Meng, Q.; Ma, X.; Song, X.; Wang, L.; Jiang, W. Drug Resistance-Related Competing Interactions of lncRNA and mRNA across 19 Cancer Types. Mol. Ther. Nucleic Acids 2019, 16, 442–451. [Google Scholar] [CrossRef] [Green Version]
  30. Qu, Y.; Tan, H.Y.; Chan, Y.T.; Jiang, H.; Wang, N.; Wang, D. The functional role of long noncoding RNA in resistance to anticancer treatment. Ther. Adv. Med. Oncol. 2020, 12, 1758835920927850. [Google Scholar] [CrossRef]
  31. Liu, K.; Gao, L.; Ma, X.; Huang, J.J.; Chen, J.; Zeng, L.; Ashby, C.R., Jr.; Zou, C.; Chen, Z.S. Long non-coding RNAs regulate drug resistance in cancer. Mol. Cancer 2020, 19, 54. [Google Scholar] [CrossRef] [Green Version]
  32. Barrett, T.; Wilhite, S.E.; Ledoux, P.; Evangelista, C.; Kim, I.F.; Tomashevsky, M.; Marshall, K.A.; Phillippy, K.H.; Sherman, P.M.; Holko, M.; et al. NCBI GEO: Archive for functional genomics data sets--update. Nucleic Acids Res. 2013, 41, D991–D995. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Clough, E.; Barrett, T. The Gene Expression Omnibus Database. Methods Mol. Biol. 2016, 1418, 93–110. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Pishas, K.I.; Lessnick, S.L. Ewing sarcoma resistance to SP-2509 is not mediated through KDM1A/LSD1 mutation. Oncotarget 2018, 9, 36413–36429. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Alnasir, J.; Shanahan, H.P. Investigation into the annotation of protocol sequencing steps in the sequence read archive. GigaScience 2015, 4, 23. [Google Scholar] [CrossRef] [Green Version]
  36. Kim, D.; Langmead, B.; Salzberg, S.L. HISAT: A fast spliced aligner with low memory requirements. Nat. Methods 2015, 12, 357–360. [Google Scholar] [CrossRef] [Green Version]
  37. Li, H.; Handsaker, B.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.; Durbin, R.; Genome Project Data Processing, S. The Sequence Alignment/Map format and SAMtools. Bioinformatics 2009, 25, 2078–2079. [Google Scholar] [CrossRef] [Green Version]
  38. Pertea, M.; Pertea, G.M.; Antonescu, C.M.; Chang, T.C.; Mendell, J.T.; Salzberg, S.L. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat. Biotechnol. 2015, 33, 290–295. [Google Scholar] [CrossRef] [Green Version]
  39. Robinson, M.D.; McCarthy, D.J.; Smyth, G.K. edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 2010, 26, 139–140. [Google Scholar] [CrossRef] [Green Version]
  40. Benjamini, Y.; Hochberg, Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J. R. Stat. Soc. Ser. B 1995, 57, 289–300. [Google Scholar] [CrossRef]
  41. Tweedie, S.; Braschi, B.; Gray, K.; Jones, T.E.M.; Seal, R.L.; Yates, B.; Bruford, E.A. Genenames.org: The HGNC and VGNC resources in 2021. Nucleic Acids Res. 2021, 49, D939–D946. [Google Scholar] [CrossRef]
  42. Kirov, S.; Ji, R.; Wang, J.; Zhang, B. Functional annotation of differentially regulated gene set using WebGestalt: A gene set predictive of response to ipilimumab in tumor biopsies. Methods Mol. Biol. 2014, 1101, 31–42. [Google Scholar] [CrossRef] [PubMed]
  43. Liao, Y.; Wang, J.; Jaehnig, E.J.; Shi, Z.; Zhang, B. WebGestalt 2019: Gene set analysis toolkit with revamped UIs and APIs. Nucleic Acids Res. 2019, 47, W199–W205. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Szklarczyk, D.; Gable, A.L.; Nastou, K.C.; Lyon, D.; Kirsch, R.; Pyysalo, S.; Doncheva, N.T.; Legeay, M.; Fang, T.; Bork, P.; et al. The STRING database in 2021: Customizable protein-protein networks, and functional characterization of user-uploaded gene/measurement sets. Nucleic Acids Res. 2021, 49, D605–D612. [Google Scholar] [CrossRef] [PubMed]
  45. Shannon, P.; Markiel, A.; Ozier, O.; Baliga, N.S.; Wang, J.T.; Ramage, D.; Amin, N.; Schwikowski, B.; Ideker, T. Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13, 2498–2504. [Google Scholar] [CrossRef]
  46. Paraskevopoulou, M.D.; Georgakilas, G.; Kostoulas, N.; Vlachos, I.S.; Vergoulis, T.; Reczko, M.; Filippidis, C.; Dalamagas, T.; Hatzigeorgiou, A.G. DIANA-microT web server v5.0: Service integration into miRNA functional analysis workflows. Nucleic Acids Res. 2013, 41, W169–W173. [Google Scholar] [CrossRef] [Green Version]
  47. Agarwal, V.; Bell, G.W.; Nam, J.W.; Bartel, D.P. Predicting effective microRNA target sites in mammalian mRNAs. eLife 2015, 4, e05005. [Google Scholar] [CrossRef]
  48. Chen, Y.; Wang, X. miRDB: An online database for prediction of functional microRNA targets. Nucleic Acids Res. 2020, 48, D127–D131. [Google Scholar] [CrossRef] [Green Version]
  49. Liu, W.; Wang, X. Prediction of functional microRNA targets by integrative modeling of microRNA binding and target expression data. Genome Biol. 2019, 20, 18. [Google Scholar] [CrossRef]
  50. Krek, A.; Grun, D.; Poy, M.N.; Wolf, R.; Rosenberg, L.; Epstein, E.J.; MacMenamin, P.; da Piedade, I.; Gunsalus, K.C.; Stoffel, M.; et al. Combinatorial microRNA target predictions. Nat. Genet. 2005, 37, 495–500. [Google Scholar] [CrossRef]
  51. Karagkouni, D.; Paraskevopoulou, M.D.; Tastsoglou, S.; Skoufos, G.; Karavangeli, A.; Pierros, V.; Zacharopoulou, E.; Hatzigeorgiou, A.G. DIANA-LncBase v3: Indexing experimentally supported miRNA targets on non-coding transcripts. Nucleic Acids Res. 2020, 48, D101–D110. [Google Scholar] [CrossRef]
  52. Sherlach, K.S.; Roepe, P.D. Drug resistance associated membrane proteins. Front. Physiol. 2014, 5, 108. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Anderson, J.T.; Huang, K.M.; Lustberg, M.B.; Sparreboom, A.; Hu, S. Solute Carrier Transportome in Chemotherapy-Induced Adverse Drug Reactions. Rev. Physiol. Biochem. Pharmacol. 2022, 183, 177–215. [Google Scholar] [CrossRef] [PubMed]
  54. Elebiyo, T.C.; Rotimi, D.; Evbuomwan, I.O.; Maimako, R.F.; Iyobhebhe, M.; Ojo, O.A.; Oluba, O.M.; Adeyemi, O.S. Reassessing vascular endothelial growth factor (VEGF) in anti-angiogenic cancer therapy. Cancer Treat. Res. Commun. 2022, 32, 100620. [Google Scholar] [CrossRef]
  55. Pavlopoulou, A.; Oktay, Y.; Vougas, K.; Louka, M.; Vorgias, C.E.; Georgakilas, A.G. Determinants of resistance to chemotherapy and ionizing radiation in breast cancer stem cells. Cancer Lett. 2016, 380, 485–493. [Google Scholar] [CrossRef] [PubMed]
  56. Joyce, H.; McCann, A.; Clynes, M.; Larkin, A. Influence of multidrug resistance and drug transport proteins on chemotherapy drug metabolism. Expert Opin. Drug Metab. Toxicol. 2015, 11, 795–809. [Google Scholar] [CrossRef] [PubMed]
  57. Barabasi, A.L.; Gulbahce, N.; Loscalzo, J. Network medicine: A network-based approach to human disease. Nat. Rev. Genet. 2011, 12, 56–68. [Google Scholar] [CrossRef] [Green Version]
  58. Sun, J.; Jiang, Z.; Li, Y.; Wang, K.; Chen, X.; Liu, G. Downregulation of miR-21 inhibits the malignant phenotype of pancreatic cancer cells by targeting VHL. OncoTargets Ther. 2019, 12, 7215–7226. [Google Scholar] [CrossRef] [Green Version]
  59. Zhang, X.; Wang, S.; Wang, H.; Cao, J.; Huang, X.; Chen, Z.; Xu, P.; Sun, G.; Xu, J.; Lv, J.; et al. Circular RNA circNRIP1 acts as a microRNA-149-5p sponge to promote gastric cancer progression via the AKT1/mTOR pathway. Mol. Cancer 2019, 18, 20. [Google Scholar] [CrossRef] [Green Version]
  60. Zhou, L.; Jiang, F.; Chen, X.; Liu, Z.; Ouyang, Y.; Zhao, W.; Yu, D. Downregulation of miR-221/222 by a microRNA sponge promotes apoptosis in oral squamous cell carcinoma cells through upregulation of PTEN. Oncol. Lett. 2016, 12, 4419–4426. [Google Scholar] [CrossRef] [Green Version]
  61. Ekimler, S.; Sahin, K. Computational Methods for MicroRNA Target Prediction. Genes 2014, 5, 671–683. [Google Scholar] [CrossRef] [Green Version]
  62. Hamzeiy, H.; Allmer, J.; Yousef, M. Computational methods for microRNA target prediction. Methods Mol. Biol. 2014, 1107, 207–221. [Google Scholar] [CrossRef] [PubMed]
  63. Quillet, A.; Saad, C.; Ferry, G.; Anouar, Y.; Vergne, N.; Lecroq, T.; Dubessy, C. Improving Bioinformatics Prediction of microRNA Targets by Ranks Aggregation. Front. Genet. 2019, 10, 1330. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Jiang, H.; Zhang, G.; Wu, J.H.; Jiang, C.P. Diverse roles of miR-29 in cancer (review). Oncol. Rep. 2014, 31, 1509–1516. [Google Scholar] [CrossRef] [Green Version]
  65. Muniyappa, M.K.; Dowling, P.; Henry, M.; Meleady, P.; Doolan, P.; Gammell, P.; Clynes, M.; Barron, N. MiRNA-29a regulates the expression of numerous proteins and reduces the invasiveness and proliferation of human carcinoma cell lines. Eur. J. Cancer 2009, 45, 3104–3118. [Google Scholar] [CrossRef]
  66. Li, W.; Zhang, B.; Jia, Y.; Shi, H.; Wang, H.; Guo, Q.; Li, H. LncRNA LOXL1-AS1 regulates the tumorigenesis and development of lung adenocarcinoma through sponging miR-423-5p and targeting MYBL2. Cancer Med. 2020, 9, 689–699. [Google Scholar] [CrossRef] [Green Version]
  67. Wu, C.; Zhang, J. Long non-conding RNA LOXL1-AS1 sponges miR-589-5p to up-regulate CBX5 expression in renal cell carcinoma. Biosci. Rep. 2020, 40, BSR20200212. [Google Scholar] [CrossRef] [PubMed]
  68. Xie, N.; Fei, X.; Liu, S.; Liao, J.; Li, Y. LncRNA LOXL1-AS1 promotes invasion and proliferation of non-small-cell lung cancer through targeting miR-324-3p. Am. J. Transl. Res. 2019, 11, 6403–6412. [Google Scholar] [PubMed]
  69. Yu, W.; Dai, Y. lncRNA LOXL1-AS1 promotes liver cancer cell proliferation and migration by regulating the miR-377-3p/NFIB axis. Oncol. Lett. 2021, 22, 624. [Google Scholar] [CrossRef] [PubMed]
  70. Hu, J.; Huang, L.; Ding, Q.; Lv, J.; Chen, Z. Long noncoding RNA HAGLR sponges miR-338-3p to promote 5-Fu resistance in gastric cancer through targeting the LDHA-glycolysis pathway. Cell Biol. Int. 2021, 46, 173–184. [Google Scholar] [CrossRef]
  71. Yang, C.; Shen, S.; Zheng, X.; Ye, K.; Sun, Y.; Lu, Y.; Ge, H. Long noncoding RNA HAGLR acts as a microRNA-143-5p sponge to regulate epithelial-mesenchymal transition and metastatic potential in esophageal cancer by regulating LAMP3. FASEB J. Off. Publ. Fed. Am. Soc. Exp. Biol. 2019, 33, 10490–10504. [Google Scholar] [CrossRef]
  72. Zou, J.; Wu, K.; Lin, C.; Jie, Z.G. LINC00319 acts as a microRNA-335-5p sponge to accelerate tumor growth and metastasis in gastric cancer by upregulating ADCY3. Am. J. Physiol. Gastrointest. Liver Physiol. 2020, 318, G10–G22. [Google Scholar] [CrossRef] [PubMed]
  73. Chang, H.; Yao, Y. lncRNA TMPO antisense RNA 1 promotes the malignancy of cholangiocarcinoma cells by regulating let-7g-5p/high-mobility group A1 axis. Bioengineered 2022, 13, 2889–2901. [Google Scholar] [CrossRef] [PubMed]
  74. Liu, G.; Yang, H.; Cao, L.; Han, K.; Li, G. LncRNA TMPO-AS1 Promotes Proliferation and Invasion by Sponging miR-383-5p in Glioma Cells. Cancer Manag. Res. 2020, 12, 12001–12009. [Google Scholar] [CrossRef] [PubMed]
  75. Wang, Y.; Ma, J.; Li, R.; Gao, X.; Wang, H.; Jiang, G. LncRNA TMPO-AS1 serves as a sponge for miR-4731-5p modulating breast cancer progression through FOXM1. Am. J. Transl. Res. 2021, 13, 11094–11106. [Google Scholar]
  76. Welch, D.; Kahen, E.; Fridley, B.; Brohl, A.S.; Cubitt, C.L.; Reed, D.R. Small molecule inhibition of lysine-specific demethylase 1 (LSD1) and histone deacetylase (HDAC) alone and in combination in Ewing sarcoma cell lines. PLoS ONE 2019, 14, e0222228. [Google Scholar] [CrossRef]
  77. Kozub, M.M.; Carr, R.M.; Lomberk, G.L.; Fernandez-Zapico, M.E. LSD1, a double-edged sword, confers dynamic chromatin regulation but commonly promotes aberrant cell growth. F1000Research 2017, 6, 2016. [Google Scholar] [CrossRef] [Green Version]
  78. Maiques-Diaz, A.; Somervaille, T.C. LSD1: Biologic roles and therapeutic targeting. Epigenomics 2016, 8, 1103–1116. [Google Scholar] [CrossRef] [Green Version]
  79. Huang, Z.; Li, S.; Song, W.; Li, X.; Li, Q.; Zhang, Z.; Han, Y.; Zhang, X.; Miao, S.; Du, R.; et al. Lysine-specific demethylase 1 (LSD1/KDM1A) contributes to colorectal tumorigenesis via activation of the Wnt/beta-catenin pathway by down-regulating Dickkopf-1 (DKK1). PLoS ONE 2013, 8, e70077. [Google Scholar] [CrossRef]
  80. Kashyap, V.; Ahmad, S.; Nilsson, E.M.; Helczynski, L.; Kenna, S.; Persson, J.L.; Gudas, L.J.; Mongan, N.P. The lysine specific demethylase-1 (LSD1/KDM1A) regulates VEGF-A expression in prostate cancer. Mol. Oncol. 2013, 7, 555–566. [Google Scholar] [CrossRef]
  81. Lim, S.; Janzer, A.; Becker, A.; Zimmer, A.; Schule, R.; Buettner, R.; Kirfel, J. Lysine-specific demethylase 1 (LSD1) is highly expressed in ER-negative breast cancers and a biomarker predicting aggressive biology. Carcinogenesis 2010, 31, 512–520. [Google Scholar] [CrossRef]
  82. Magliulo, D.; Bernardi, R.; Messina, S. Lysine-Specific Demethylase 1A as a Promising Target in Acute Myeloid Leukemia. Front. Oncol. 2018, 8, 255. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Rondou, P.; Haegeman, G.; Van Craenenbroeck, K. The dopamine D4 receptor: Biochemical and signalling properties. Cell. Mol. Life Sci. CMLS 2010, 67, 1971–1986. [Google Scholar] [CrossRef]
  84. Undieh, A.S. Pharmacology of signaling induced by dopamine D(1)-like receptor activation. Pharmacol. Ther. 2010, 128, 37–60. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  85. Lin, L.; Yee, S.W.; Kim, R.B.; Giacomini, K.M. SLC transporters as therapeutic targets: Emerging opportunities. Nat. Rev. Drug Discov. 2015, 14, 543–560. [Google Scholar] [CrossRef] [Green Version]
  86. Li, Q.; Shu, Y. Role of solute carriers in response to anticancer drugs. Mol. Cell. Ther. 2014, 2, 15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  87. Rosas-Cruz, A.; Salinas-Jazmin, N.; Velazquez, M.A.V. Dopamine Receptors in Cancer: Are They Valid Therapeutic Targets? Technol. Cancer Res. Treat. 2021, 20, 15330338211027913. [Google Scholar] [CrossRef]
  88. Weissenrieder, J.S.; Neighbors, J.D.; Mailman, R.B.; Hohl, R.J. Cancer and the Dopamine D2 Receptor: A Pharmacological Perspective. J. Pharmacol. Exp. Ther. 2019, 370, 111–126. [Google Scholar] [CrossRef]
  89. Sung, H.Y.; Han, J.; Ju, W.; Ahn, J.H. Synaptotagmin-like protein 2 gene promotes the metastatic potential in ovarian cancer. Oncol. Rep. 2016, 36, 535–541. [Google Scholar] [CrossRef] [Green Version]
  90. Skibola, D.R.; Smith, M.T.; Bracci, P.M.; Hubbard, A.E.; Agana, L.; Chi, S.; Holly, E.A. Polymorphisms in ghrelin and neuropeptide Y genes are associated with non-Hodgkin lymphoma. Cancer Epidemiol. Biomark. Prev. 2005, 14, 1251–1256. [Google Scholar] [CrossRef] [Green Version]
  91. Syrbe, S.; Hedrich, U.B.S.; Riesch, E.; Djemie, T.; Muller, S.; Moller, R.S.; Maher, B.; Hernandez-Hernandez, L.; Synofzik, M.; Caglayan, H.S.; et al. De novo loss- or gain-of-function mutations in KCNA2 cause epileptic encephalopathy. Nat. Genet. 2015, 47, 393–399. [Google Scholar] [CrossRef]
  92. Bouillet, P.; Sapin, V.; Chazaud, C.; Messaddeq, N.; Decimo, D.; Dolle, P.; Chambon, P. Developmental expression pattern of Stra6, a retinoic acid-responsive gene encoding a new type of membrane protein. Mech. Dev. 1997, 63, 173–186. [Google Scholar] [CrossRef]
  93. Rutter, A.R.; Fradley, R.L.; Garrett, E.M.; Chapman, K.L.; Lawrence, J.M.; Rosahl, T.W.; Patel, S. Evidence from gene knockout studies implicates Asc-1 as the primary transporter mediating d-serine reuptake in the mouse CNS. Eur. J. Neurosci. 2007, 25, 1757–1766. [Google Scholar] [CrossRef] [PubMed]
  94. Gaynes, J.S.; Jonart, L.M.; Zamora, E.A.; Naumann, J.A.; Gossai, N.P.; Gordon, P.M. The central nervous system microenvironment influences the leukemia transcriptome and enhances leukemia chemo-resistance. Haematologica 2017, 102, e136–e139. [Google Scholar] [CrossRef] [Green Version]
  95. Zhao, C.M.; Hayakawa, Y.; Kodama, Y.; Muthupalani, S.; Westphalen, C.B.; Andersen, G.T.; Flatberg, A.; Johannessen, H.; Friedman, R.A.; Renz, B.W.; et al. Denervation suppresses gastric tumorigenesis. Sci. Transl. Med. 2014, 6, 250ra115. [Google Scholar] [CrossRef] [Green Version]
  96. Logotheti, S.; Marquardt, S.; Richter, C.; Sophie Hain, R.; Murr, N.; Takan, I.; Pavlopoulou, A.; Pützer, B.M. Neural Networks Recapitulation by Cancer Cells Promotes Disease Progression: A Novel Role of p73 Isoforms in Cancer-Neuronal Crosstalk. Cancers 2020, 12, 3789. [Google Scholar] [CrossRef] [PubMed]
  97. Russell-Goldman, E.; Hornick, J.L.; Qian, X.; Jo, V.Y. NKX2.2 immunohistochemistry in the distinction of Ewing sarcoma from cytomorphologic mimics: Diagnostic utility and pitfalls. Cancer Cytopathol. 2018, 126, 942–949. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  98. Peter, M.; Couturier, J.; Pacquement, H.; Michon, J.; Thomas, G.; Magdelenat, H.; Delattre, O. A new member of the ETS family fused to EWS in Ewing tumors. Oncogene 1997, 14, 1159–1164. [Google Scholar] [CrossRef] [Green Version]
  99. Tsuda, Y.; Dickson, B.C.; Swanson, D.; Sung, Y.S.; Zhang, L.; Meyers, P.; Healey, J.H.; Antonescu, C.R. Ewing sarcoma with FEV gene rearrangements is a rare subset with predilection for extraskeletal locations and aggressive behavior. Genes Chromosomes Cancer 2020, 59, 286–294. [Google Scholar] [CrossRef]
  100. Shibuya, R.; Matsuyama, A.; Nakamoto, M.; Shiba, E.; Kasai, T.; Hisaoka, M. The combination of CD99 and NKX2.2, a transcriptional target of EWSR1-FLI1, is highly specific for the diagnosis of Ewing sarcoma. Virchows Arch. 2014, 465, 599–605. [Google Scholar] [CrossRef]
  101. Malt, E.A.; Juhasz, K.; Malt, U.F.; Naumann, T. A Role for the Transcription Factor Nk2 Homeobox 1 in Schizophrenia: Convergent Evidence from Animal and Human Studies. Front. Behav. Neurosci. 2016, 10, 59. [Google Scholar] [CrossRef] [Green Version]
  102. Ferguson, A.V.; Bains, J.S. Actions of angiotensin in the subfornical organ and area postrema: Implications for long term control of autonomic output. Clin. Exp. Pharmacol. Physiol. 1997, 24, 96–101. [Google Scholar] [CrossRef]
  103. Leon, T.Y.; Ngan, E.S.; Poon, H.C.; So, M.T.; Lui, V.C.; Tam, P.K.; Garcia-Barcelo, M.M. Transcriptional regulation of RET by Nkx2-1, Phox2b, Sox10, and Pax3. J. Pediatr. Surg. 2009, 44, 1904–1912. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  104. Camilleri, M.; Wieben, E.; Eckert, D.; Carlson, P.; Hurley O’Dwyer, R.; Gibbons, D.; Acosta, A.; Klee, E.W. Familial chronic megacolon presenting in childhood or adulthood: Seeking the presumed gene association. Neurogastroenterol. Motil. 2019, 31, e13550. [Google Scholar] [CrossRef] [PubMed]
  105. Danková, M.; Tóth, Š.; Holodová, M.; Fagová, Z.; Čurgali, K.; Mechírová, E.; Maretta, M.; Nemcová, R.; Gancarčíková, S.; Polák, Š. Immunohistochemical visualisation of the enteric nervous system architecture in the germ-free piglets. J. Mol. Histol. 2022, 53, 773–780. [Google Scholar] [CrossRef]
  106. Costes, L.M.; Boeckxstaens, G.E.; de Jonge, W.J.; Cailotto, C. Neural networks in intestinal immunoregulation. Organogenesis 2013, 9, 216–223. [Google Scholar] [CrossRef] [PubMed]
  107. Sternini, C. Organization of the peripheral nervous system: Autonomic and sensory ganglia. J. Investig. Dermatol. Symp. Proc. 1997, 2, 1–7. [Google Scholar] [CrossRef]
  108. Geronikolou, S.A.; Albanopoulos, K.; Chrousos, G.; Cokkinos, D. Evaluating the Homeostasis Assessment Model Insulin Resistance and the Cardiac Autonomic System in Bariatric Surgery Patients: A Meta-Analysis. Adv. Exp. Med. Biol. 2017, 988, 249–259. [Google Scholar] [CrossRef]
  109. Geronikolou, S.A.; Pavlopoulou, A.; Cokkinos, D.; Chrousos, G. Interactome of Obesity: Obesidome: Genetic Obesity, Stress Induced Obesity, Pathogenic Obesity Interaction. Adv. Exp. Med. Biol. 2017, 987, 233–241. [Google Scholar] [CrossRef]
  110. Molteni, A.; Ward, W.F.; Ts’ao, C.H.; Taylor, J.; Small, W., Jr.; Brizio-Molteni, L.; Veno, P.A. Cytostatic properties of some angiotensin I converting enzyme inhibitors and of angiotensin II type I receptor antagonists. Curr. Pharm. Des. 2003, 9, 751–761. [Google Scholar] [CrossRef]
Figure 1. Graphical illustration of the overall methodology of this study. Genes differentially expressed between the chemoresistant and chemoresponsive, as well as the chemoresponsive versus parental Ewing sarcoma cells were detected. Functional annotation analysis of the genes upregulated in the chemoresistant cells and their protein products was performed. Subsequently, a ceRNA was generated from the drug resistance genes, their corresponding regulating miRNAs, and the lncRNAs that act as sponges of these miRNAs.
Figure 1. Graphical illustration of the overall methodology of this study. Genes differentially expressed between the chemoresistant and chemoresponsive, as well as the chemoresponsive versus parental Ewing sarcoma cells were detected. Functional annotation analysis of the genes upregulated in the chemoresistant cells and their protein products was performed. Subsequently, a ceRNA was generated from the drug resistance genes, their corresponding regulating miRNAs, and the lncRNAs that act as sponges of these miRNAs.
Ijerph 20 06288 g001
Figure 2. Bar plot depicting the over-represented biological processes (GO terms) in the chemoresistant genes. The x-axis corresponds to the number of genes associated with the GO terms; the length of the bars is proportional to the number of genes participating in the given biological process.
Figure 2. Bar plot depicting the over-represented biological processes (GO terms) in the chemoresistant genes. The x-axis corresponds to the number of genes associated with the GO terms; the length of the bars is proportional to the number of genes participating in the given biological process.
Ijerph 20 06288 g002
Figure 3. The STRING output interaction network of the resistance genes/proteins. The nodes represent genes/gene products, and the connecting lines denote functional associations. Pink and cyan fill colors indicate upregulation and downregulation of the corresponding genes, respectively.
Figure 3. The STRING output interaction network of the resistance genes/proteins. The nodes represent genes/gene products, and the connecting lines denote functional associations. Pink and cyan fill colors indicate upregulation and downregulation of the corresponding genes, respectively.
Ijerph 20 06288 g003
Figure 4. Competing endogenous RNA network in Ewing sarcoma resistance. The lncRNAs are represented by polygons, the miRNAs sponged by lncRNAs are denoted by rectangles, and the miRNA target genes are indicated by circles.
Figure 4. Competing endogenous RNA network in Ewing sarcoma resistance. The lncRNAs are represented by polygons, the miRNAs sponged by lncRNAs are denoted by rectangles, and the miRNA target genes are indicated by circles.
Ijerph 20 06288 g004
Table 1. Common genes between the SP-2509 drug resistance associated genes and murine genes affecting nervous system development. Associated disease syndromes for each gene may be crosschecked in the GeneCards/disorders platform. Genes associated to neurological disorders are marked in blue, genes associated to tumorogenesis or malignancies are marked in red, whilst genes involved in neural and/or cancer entities are marked in purple.
Table 1. Common genes between the SP-2509 drug resistance associated genes and murine genes affecting nervous system development. Associated disease syndromes for each gene may be crosschecked in the GeneCards/disorders platform. Genes associated to neurological disorders are marked in blue, genes associated to tumorogenesis or malignancies are marked in red, whilst genes involved in neural and/or cancer entities are marked in purple.
Gene SymbolGene NameAssociated Disease/Syndrome
AADATalpha-aminoadipate aminotransferaseDetrusor Sphincter Dyssynergia and Huntington Disease
ACVRL1activin A receptor like type 1Telangiectasia, Hereditary Hemorrhagic, Type 2; Idiopathic Pulmonary Arterial Hypertension
ARandrogen receptorSpinal And Bulbar Muscular Atrophy
ARAP3Ankyrin Repeat And Plekstrin Homology Domains-Containing Protein 3Neurofibromatosis 1
ARHGAP44Rho GTPase Activating Protein 44Hyperalphalipoproteinemia 2
ASS1Argininosuccinate Synthase 1Citrullinemia
ATP11AATPase Phospholipid Transporting 11ACOVID-19; leukodystrophy
BMP7Bone Morphogenetic Protein 7renal fibrosis, spondylolistheis,
CACNB3Calcium Voltage-Gated Channel Auxiliary Subunit Beta 3distal hereditary motor neuropathy 2
CC2D1ACoiled-Coil And C2 Domain Containing 1AIntellectual Developmental Disorder, Autosomal Recessive 3; Autosomal Recessive Non-Syndromic Intellectual Disability; cerebral palsy
CCKARCholecystokinin A Receptorpanic disorder, functionless pituitary adenoma
CCNFcyclin Famyotrophic lateral sclerosis
CHATCholine O-Acetyltransferasemyasthenic syndrome, central sleep apnea, respiratory failure
CNTNAP1ontactin Associated Protein 1lethal congenital contracture syndrome, neuropathy
Collagen Type II Alpha 1 Chainepiphysial dysplasia
CRB2Crumbs Cell Polarity Complex Component 2Genetic Steroid-Resistant Nephrotic Syndrome; Ventriculomegaly With Cystic Kidney Disease
CST3Cystatin Ckidney disease, Cerebral Amyloid Angiopathy, Cst3-Relate
CTSFCathepsin FNeuronal Ceroid Lipofuscinosis; Spinocerebellar Ataxia
CXCR4C-X-C Motif Chemokine Receptor 4dementia; whim syndrome
DAB1DAB Adaptor Protein 1epilepsy; spinocerebral ataxia
DLX6Distal-Less Homeobox 6Isolated Split Hand-Split Foot Malformation
DMRTA2DMRT Like Family A2Uncertain
DNM1Dynamin 1Developmental And Epileptic Encephalopathy 31
DRD1Dopamine Receptor D1cerebral meningioma; Early-Onset Schizophrenia; Attention Deficit-Hyperactivity Disorder; heroin dependence
DRD2Dopamine Receptor D2Cocaine Dependence; drug deeondence; Antisocial Personality Disorder
DZIP1DAZ Interacting Zinc Finger Protein 1Orthostatic Intolerance; Mitral Valve Prolapse 3; spermatogenic failure
EEF1A2Eukaryotic Translation Elongation Factor 1 Alpha 2Developmental And Epileptic Encephalopathy
EEF2KEukaryotic Elongation Factor 2 Kinasecolorectal adenocarcinoma
EFNB3ephrin beta3Craniofrontonasal Syndrome
ENTPD1Ectonucleoside Triphosphate Diphosphohydrolase 1spastic paraplegia; Radiation Proctitis; beta-thalassemia
FASFas Cell Surface Death ReceptorAutoimmune Lymphoproliferative Syndrome
FASNFatty Acid SynthaseNAFLD, prostate Ca
FBXL21PF-Box And Leucine Rich Repeat Protein 21, Pseudogeneuncertain
FEVFEV Transcription FactorEwing sarcoma; anxiety
FLT4Fms Related Receptor Tyrosine Kinase 4Congenital Heart Defects, Multiple Types, 7; Hereditary Lymphedema I; Hemangioma
FN1fibronectin 1soft tissue chondroma; Spondylometaphyseal Dysplasia; Spondyloepimetaphyseal Dysplasia
FRZBFrizzled Related Proteinosteogenic sarcoma; retinal degeneration
GAL3ST1Galactose-3-O-Sulfotransferase 1Renal Cell Carcinoma, Nonpapillary
GAS1Growth Arrest Specific 1Septopreoptic Holoprosencephaly
GAS7Growth Arrest Specific 7Open-Angle Glaucoma; deafness; glaucoma
GBAGlucosylceramidase BetaGaucher Disease, Type I,2,3,3c
GDF11Growth Differentiation Factor 11 Vertebral Hypersegmentation And Orofacial Anomalies; aging
GLI2GLI Family Zinc Finger 2Holoprosencephaly 9; Combined Pituitary Hormone Deficiencies
HAND2Heart And Neural Crest Derivatives Expressed 2Cardiomyopathy, Dilated, 1a/h; Dilated Cardiomyopathy
HHEXHematopoietically Expressed HomeoboxHeart Defects, Congenital, And Other Congenital Anomalies; Diabetes Mellitus, Neonatal, With Congenital Hypothyroidism
HPCAHippocalcinDystonia 2, Torsion, Autosomal Recessive
HSPG2Heparan Sulfate Proteoglycan 2Dyssegmental Dysplasia, Silverman-Handmaker Type Tardive Dyskinesia
IL33Interleukin 33Chronic Asthma
INKA1Inka Box Actin Regulator 1Visual Cortex Disease
IQSEC2IQ Motif And Sec7 Domain ArfGEF 2Intellectual Developmental Disorder, X-Linked 1
ISL2ISL LIM Homeobox 2Amme complex
KIRREL3Kirre Like Nephrin Family Adhesion Molecule 3Autosomal Dominant Non-Syndromic Intellectual Disability; Familial Nephrotic Syndrome; Granulomatous Disease
KNDC1Kinase Non-Catalytic C-Lobe Domain Containing 1Brugada syndrome
L1CAML1 Cell Adhesion MoleculeCorpus Callosum, Partial Agenesis Of, X-Linked; hydrocephaly
LAMA5Laminin Subunit Alpha 5Vitreous Detachment; Presynaptic Congenital Myasthenic Syndromes; Lama5-Related Multisystemic Syndrome
LHX6LIM Homeobox 6Waardenburg Syndrome, Type 2c; tooth ageenesis
LIMK1LIM Domain Kinase 1Supravalvular Aortic Stenosis
LMNB1Lamin B1Hutchinson-Gilford Progeria Syndrome; Leukodystrophy, Demyelinating, Adult-Onset
LYNX1Ly6/Neurotoxin 1Malde Meleda; Ovarian Serous Cystadenocarcinoma; amblyopia
MAPTMicrotubule Associated Protein TauFrontotemporal Dementia; Supranuclear Palsy,; parkinsson disease; parkinsson-dementia syndrome
MARVELD1MARVEL Domain Containing 1Facial Nerve Disease; Facial Paralysis; Heart Fibrosarcoma; Cranial Nerve Disease
MAVSMitochondrial Antiviral Signaling Proteinhepatitis; Newcastle diseases; Oropouche Fever
MMP9Matrix Metallopeptidase 9Metaphyseal Anadysplasia; Central Nervous System Tuberculosis; Brain Glioblastoma Multiforme; internal hemoroids
NAT8LN-Acetyltransferase 8 LikeN-Acetylaspartate Deficiency; Canavan Disease; microcephaly; Miliaria Rubra
NCKIPSDNCK Interacting Protein With SH3 DomainWiskott-Aldrich Syndrome; thrombocytopenia; leucemia, spheroid syndrome
NDRG2NDRG Family Member 2Charcot-Marie-Tooth Disease, Glioblastoma, Meningioma
NESnestinMedulloepithelioma; Central Neurocytoma; Optic Nerve Glioma; Periventricular Leukomalacia
NFIXNuclear Factor I XMalan Syndrome; Marshall-Smith Syndrome; megalocornea
NKX2-2NK2 Homeobox 2Maturity-Onset Diabetes Of The Young; Oligodendroglioma; Ewing Sarcoma
NODALNodal Growth Differentiation FactorWolff–Parkinson–White Syndrome; heart disease
NOS3Nitric Oxide Synthase 3stroke; alzheimer; hypertention, preeclampsia
NOTCH3Notch Receptor 3Cerebral Arteriopathy, Autosomal Dominant, With Subcortical Infarcts And Leukoencephalopathy
NPY1RNeuropeptide Y Receptor Y1Body Mass Index Quantitative Trait Locus 11
NTN1netrin 1Mirror Movements 4/1, Superior Semicircular Canal Dehiscence, Colorectal Ca
OTX1Orthodenticle Homeobox 1Agnathia-Otocephaly Complex, epilepsy, Bardet-Biedl Syndrome 15; dyslexia
PDGFRBPlatelet Derived Growth Factor Receptor BetaKosaki Overgrowth Syndrome, premature aging syndrome,
PIP5K1CPhosphatidylinositol-4-Phosphate 5-Kinase Type 1 Gamma Neurogenic Bladder; alcohol use disorder; Cerebellar Ataxia, Cayman Type
PITX1Paired Like Homeodomain 1Clubfoot, Congenital, With Or Without Deficiency Of Long Bones And/Or Mirror-Image Polydactyly, Clubfoot
PMP22Peripheral Myelin Protein 22Charcot-Marie-Tooth Disease And Deafness, Neuropathy, Hereditary, With Liability To Pressure Palsies, Guillain-Barre Syndrome
PRPHperipherinFrontotemporal Dementia And/Or Amyotrophic Lateral Sclerosis 3; Amyotrophic Lateral Sclerosis 19
PRRT2Proline Rich Transmembrane Protein 2Episodic Kinesigenic Dyskinesia 1, Prrt2-Associated Paroxysmal Movement Disorders, Convulsions, Familial Infantile, With Paroxysmal Choreoathetosis
RASGRF1Ras Protein Specific Guanine Nucleotide Releasing Factor 1myopia, Degenerative Myopia, Transient Neonatal Diabetes Mellitus
RETRet Proto-OncogenePheochromocytoma, Multiple Endocrine Neoplasia, Thyroid Carcinoma
RNF165Ring Finger Protein 165uncertaın
RTN4RReticulon 4 Receptorspinal cord injury, schizophrenia, psychotic disorder, leucemia
S1PR2Sphingosine-1-Phosphate Receptor 2deafness, pulmonary edema
SARM1Sterile Alpha And TIR Motif Containing 1Wallerian Degeneration, Retinoschisis 1
SCTRSecretin ReceptorGastrinoma, Jejunal Somatostatinoma, pancreas disease, Primary Biliary Cholangiti
SEMA3Fsemaphorin 3Fneuroma, megacolon, lung Ca, tongue carcinoma
SEMA6BSemaphorin 6Bepilepsy
SH3TC2SH3 Domain And Tetratricopeptide Repeats 2Genetic Motor Neuron Disease, Charcot-Marie-Tooth
SLC18A3Solute Carrier Family 18 Member 3Myasthenic Syndrome, Fetal Akinesia Deformation Sequence 1
SLC19A3Solute Carrier Family 19 Member 3Thiamine Metabolism Dysfunction Syndrome 2, Infantile Spasms-Psychomotor Retardation-Progressive Brain Atrophy-Basal Ganglia Disease Syndrome
SLC25A1Solute Carrier Family 25 Member 1Myasthenic Syndrome,, Combined D-2- And L-2-Hydroxyglutaric Aciduria
SLC7A10Solute Carrier Family 7 Member 10Hypotonia-Cystinuria Syndrome
SMARCC2SWI/SNF Related, Matrix Associated, Actin Dependent Regulator Of Chromatin Subfamily C Member 2Neurilemmomatosis
SRCIN1SRC Kinase Signaling Inhibitor 1Kohlschutter-Tonz Syndrome
SRRM3Serine/Arginine Repetitive Matrix 3deafness, burst Ca
STRA6Signaling Receptor And Transporter Of Retinol STRA6Microphthalmia, Syndromic 9
SYT2Synaptotagmin 2Myasthenic Syndrome, Congenital, 7a, Presynaptic, And Distal Motor Neuropathy; Lambert-Eaton Myasthenic Syndrome
TCF7L1Transcription Factor 7 Like 1Brust Ca, colorectal Ca, heaptocellular Ca, Arrhythmogenic Right Ventricular Cardiomyopathy
TEAD2TEA Domain Transcription Factor 2Multiple Acyl-Coa Dehydrogenase Deficiency; spastic paraplegia
THRAThyroid Hormone Receptor AlphaHypothyroidism, Resistance To Thyroid Hormone Due To A Mutation In Thyroid Hormone Receptor Alpha, Hyperthyroxinemia, bone giant cell tumor
TLE5TLE Family Member 5, Transcriptional ModulatorArthrogryposis, Distal
TNFRSF1ATNF Receptor Superfamily Member 1ACharge syndrome, periodic fever
TNFRSF1BTNF Receptor Superfamily Member 1BPsoriatic Arthritis, Rheumatoid Arthritis, Juvenile Rheumatoid Arthritis
TP73Tumor Protein P73respiratory failure, Oligodendroglioma, neuroblastoma
TRIM71Tripartite Motif Containing 71all types hydrocephalus, Autosomal Recessive Limb-Girdle Muscular Dystrophy Type 2h
TRPV4Transient Receptor Potential Cation Channel Subfamily V Member 4metatropic dysplasia, Hereditary Motor And Sensory Neuropathy
TSKUTsukushi, Small Leucine Rich Proteoglycanmiliaria
TUBB4ATubulin Beta 4A Class IVacerebral palsy, torsion dystonia, nervous system disease
UCP2Uncoupling Protein 2Hyperinsulinism Due To Ucp2 Deficiency
VAV3Vav Guanine Nucleotide Exchange Factor 3ovarian Ca, meningioma, glaucoma
VAX1Ventral Anterior Homeobox 1microphalmia
VEGFBVascular Endothelial Growth Factor BMacular Retinal Edema; Macular Degeneration, Age-Related, 1
VWA1Von Willebrand Factor A Domain Containing 1Neuropathy, Hereditary Motor, With Myopathic Features; Neuromuscular Disease
WDR62WD Repeat Domain 62Congenital Nervous System Abnormality
ZIC1Zic Family Member 1Craniosynostosis 6
ZMIZ1Zinc Finger MIZ-Type Containing 1Neurodevelopmental Disorder With Dysmorphic Facies And Distal Skeletal Anomalies; Syndromic Intellectual Disability; Brain Stem Infarction
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Asfa, S.; Toy, H.I.; Arshinchi Bonab, R.; Chrousos, G.P.; Pavlopoulou, A.; Geronikolou, S.A. Soft Tissue Ewing Sarcoma Cell Drug Resistance Revisited: A Systems Biology Approach. Int. J. Environ. Res. Public Health 2023, 20, 6288. https://doi.org/10.3390/ijerph20136288

AMA Style

Asfa S, Toy HI, Arshinchi Bonab R, Chrousos GP, Pavlopoulou A, Geronikolou SA. Soft Tissue Ewing Sarcoma Cell Drug Resistance Revisited: A Systems Biology Approach. International Journal of Environmental Research and Public Health. 2023; 20(13):6288. https://doi.org/10.3390/ijerph20136288

Chicago/Turabian Style

Asfa, Seyedehsadaf, Halil Ibrahim Toy, Reza Arshinchi Bonab, George P. Chrousos, Athanasia Pavlopoulou, and Styliani A. Geronikolou. 2023. "Soft Tissue Ewing Sarcoma Cell Drug Resistance Revisited: A Systems Biology Approach" International Journal of Environmental Research and Public Health 20, no. 13: 6288. https://doi.org/10.3390/ijerph20136288

APA Style

Asfa, S., Toy, H. I., Arshinchi Bonab, R., Chrousos, G. P., Pavlopoulou, A., & Geronikolou, S. A. (2023). Soft Tissue Ewing Sarcoma Cell Drug Resistance Revisited: A Systems Biology Approach. International Journal of Environmental Research and Public Health, 20(13), 6288. https://doi.org/10.3390/ijerph20136288

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop