Identification of genes associated with methotrexate resistance in methotrexate-resistant osteosarcoma cell lines

Exploring circulating MiRNA signature for osteosarcoma detection: Bioinformatics-based analyzing and validation

Early detection followed by efficient treatment still remain a considerable challenge for osteosarcoma (OS), indicating the importance of emerging innovative diagnostic methods. Circulating miRNAs offer a promising and non-invasive approach to assess the OS molecular landscapes. This study utilized RNAseq data from OS plasma miRNA expression profiles (PRJEB30542) and PCR Array data (GSE65071) from GEO and ENA databases. In total, 43 miRNAs demonstrated significant differential expression in OS samples of training dataset. A diagnostic model, including hsa-miR-30a-5p, hsa-miR-556-3p, hsa-miR-200a-3p, and hsa-miR-582-5p was identified through multivariate logistic regression analysis and demonstrated significant efficacy in differentiating OS patients from healthy controls in the validation group (AUC: 0.917, sensitivity: 1, specificity: 0.85). The result of target gene prediction and functional enrichment analyses revealed significant associations with terms such as epithelial morphogenesis, P53 and Wnt signaling pathways, and neoplasm metastasis. Further bioinformatics-based evaluations showed that the down-regulation of these miRNAs significantly correlates with poor prognosis and lower survival rate in OS patients and propose their tumor suppressor function in pathogenesis of OS. Furthermore, the study developed a miRNA-mRNA subnetwork that connects these miRNAs to the P53 and Wnt signaling pathways, which are critical pathways with oncogenic effects on OS progression. This comprehensive approach not only presents a promising diagnostic model but also proposes potential molecular markers for OS early diagnosis, making prognosis, and targeted therapy. The identified miRNA-mRNA functional axis holds promise as a valuable resource for further research in understanding OS pathogenesis and establishing therapeutic modalities.

Clinical Significance and Potential Mechanisms of the RNA Methyltransferase KIAA1429 in Osteosarcoma

Background: KIAA1429, a member of the RNA methyltransferase complex, is involved in cancer progression; however, the clinical significance and underlying mechanism of KIAA1429 in osteosarcoma (OS) remains to be reported. Methods: We evaluated the clinical significance of KIAA1429 in OS by performing RT-qPCR, microarray, and RNA sequencing and using published data as a reference. Two KIAA1429-targeting siRNA constructs were transfected into SW1353 cells. CCK-8 assay, colony formation assays, flow cytometry and the xenograft mouse model were conducted to investigate the biological function of KIAA1429 in OS. Results: The mRNA expression of KIAA1429 was markedly upregulated in 250 OS samples as compared to that in 71 non-cancer samples (standardized mean difference = 0.67). Summary receiver operating characteristic curve analysis revealed that KIAA1429 exhibited reliable diagnostic capacity to differentiate OS samples from non-cancer samples (area under the curve = 0.83). Further, survival analysis indicated that KIAA1429 overexpression was associated with shorter overall survival time. Knocking down KIAA1429 reduced m6A methylation levels, inhibited proliferation, prevented the growth of tumors in vivo and accelerated apoptosis of OS cells. In total, 395 KIAA1429-related genes were identified among co-expressed genes and differentially expressed genes, which were enriched in the cell cycle pathway. Protein-protein interaction network analysis showed that CDK1, CCNA2, and CCNB1 were KIAA1429-related genes, serving as major network hubs in OS. Conclusions: Our findings indicate that KIAA1429 plays an oncogenic role in OS and potentially facilitates OS progression via a mechanism that involves regulating CDK1, CCNA2, and CCNB1.

AURORA KINASE A and related downstream molecules: A potential network for cancer therapy

Aurora-A kinase (AURKA) belongs to the serine/threonine kinase family specific to cell division. In normal cells, activation of the AURKA protein is essential for regulating chromosomal segregation and centrosome maturation. The physiological concentration of AURKA accumulation has utmost importance during cell division. AURKA starts accumulating during the S phase of the cell cycle, gets functionally activated during the G2/M phase, attaches to the microtubule, and gets degraded during mitotic exit. Overexpression of AURKA could lead to deregulated cell cycle division, which is intrinsic to numerous cancers. Moreover, dysregulated AURKA affects various downstream molecules that aid in cancer pathogenesis. AURKA phosphorylates its substrates, including oncoproteins, transcriptional factors, tumor suppressor proteins, or other kinases central to various oncogenic signaling pathways critical to cancer. Considering the central role of AURKA in cell proliferation and tumorigenesis, targeting AURKA can be a novel alternative to cancer management. Several AURKA inhibitors have shown promising responses against different cancers either as a single agent or combined with various therapies. This chapter briefly discusses the role of AURKA and its downstream molecules in cancer vis-à-vis the role of AURKA inhibitor in chemoprevention.

Mitochondrial-Related Transcriptome Feature Correlates with Prognosis, Vascular Invasion, Tumor Microenvironment, and Treatment Response in Hepatocellular Carcinoma

Background

Hepatocellular carcinoma (HCC) is the most common subtype of primary liver cancer, which was highly correlated with metabolic dysfunction. Nevertheless, the association between nuclear mitochondrial-related transcriptome and HCC remained unclear.

Materials and Methods

A total of 147 nuclear mitochondrial-related genes (NMRGs) were downloaded from the MITOMAP: A Human Mitochondrial Genome Database. The training dataset was downloaded from The Cancer Genome Atlas (TCGA), while validation datasets were retrieved from the International Cancer Genome Consortium (ICGC) and Gene Expression Omnibus (GEO). The univariate and multivariate, and least absolute shrinkage and selection operator (LASSO) Cox regression analyses were applied to construct a NMRG signature, and the value of area under receiver operating characteristic curve (AUC) was utilized to assess the signature and nomogram. Then, data from the Genomics of Drug Sensitivity in Cancer (GDSC) were used for the evaluation of chemotherapy response in HCC.

Results

Functional enrichment of differentially expressed genes (DEGs) between HCC and paired normal tissue samples demonstrated that mitochondrial dysfunction was significantly associated with HCC development. Survival analysis showed a total of 35 NMRGs were significantly correlated with overall survival (OS) of HCC, and the LASSO Cox regression analysis further identified a 25-NMRG signature and corresponding prognosis score based on their transcriptional profiling. HCC patients were divided into high- and low-risk groups according to the median prognosis score, and high-risk patients had significantly worse OS (median OS: 27.50 vs. 83.18 months, P < 0.0001). The AUC values for OS at 1, 3, and 5 years were 0.79, 0.77, and 0.77, respectively. The prognostic capacity of NMRG signature was verified in the GSE14520 dataset and ICGC-HCC cohort. Besides, the NMRG signature outperformed each NMRG and clinical features in prognosis prediction and could also differentiate whether patients presented with vascular invasions (VIs) or not. Subsequently, a prognostic nomogram (C-index: 0.753, 95% CI: 0.703~0.804) by the integration of age, tumor metastasis, and NMRG prognosis score was constructed with the AUC values for OS at 1, 3, and 5 years were 0.82, 0.81, and 0.82, respectively. Notably, significant enrichment of regulatory and follicular helper T cells in high-risk group indicated the potential treatment of immune checkpoint inhibitors for these patients. Interestingly, the NMRG signature could also identify the potential responders of sorafenib or transcatheter arterial chemoembolization (TACE) treatment. Additionally, HCC patients in high-risk group appeared to be more sensitive to cisplatin, vorinostat, and methotrexate, reversely, patients in low-risk group had significantly higher sensitivity to paclitaxel and bleomycin instead.

Conclusions

In summary, the development of NMRG signature provided a more comprehensive understanding of mitochondrial dysfunction in HCC, helped predict prognosis and tumor microenvironment, and provided potential targeted therapies for HCC patients with different NMRG prognosis scores.

Genomic and transcriptomic analyses reveal a tandem amplification unit of 11 genes and mutations in mismatch repair genes in methotrexate-resistant HT-29 cells

DHFR gene amplification is commonly present in methotrexate (MTX)-resistant colon cancer cells and acute lymphoblastic leukemia. In this study, we proposed an integrative framework to characterize the amplified region by using a combination of single-molecule real-time sequencing, next-generation optical mapping, and chromosome conformation capture (Hi-C). We identified an amplification unit spanning 11 genes, from the DHFR gene to the ATP6AP1L gene position, with high adjusted interaction frequencies on chromosome 5 (~2.2 Mbp) and a twenty-fold tandemly amplified region, and novel inversions at the start and end positions of the amplified region as well as frameshift insertions in most of the MSH and MLH genes were detected. These mutations might stimulate chromosomal breakage and cause the dysregulation of mismatch repair. Characterizing the tandem gene-amplified unit may be critical for identifying the mechanisms that trigger genomic rearrangements. These findings may provide new insight into the mechanisms underlying the amplification process and the evolution of drug resistance.

Sequencing a large region of DNA containing many surplus copies of genes linked to drug resistance in colon cancer cells may illuminate how these genomic rearrangements arise. Such regions of gene amplification are highly repetitive, making them impossible to sequence using ordinary methods, and little is known about how they are generated. Using advanced methods, Jeong-Sun Seo at Seoul National University Bundang Hospital in South Korea and co-workers sequenced a region of gene amplification in colon cancer cells. The amplified region was approximately 20 times the length of that in healthy cells and contained many copies of an eleven-gene segment, including a gene implicated in drug resistance. The region also contained mutations in chromosomal repair genes which would disrupt repair pathways. These results illuminate the genetic changes that lead to gene amplification and drug resistance in cancer cells.

High expression levels of centromere protein A plus upregulation of the phosphatidylinositol 3-kinase/Akt/mammalian target of rapamycin signaling pathway affect chemotherapy response and prognosis in patients with breast cancer

Centromere proteins (CENPs) are involved in mitosis, and CENP gene expression levels are associated with chemotherapy responses in patients with breast cancer. The present study aimed to examine the roles and underlying mechanisms of the effects of CENP genes on chemotherapy responses and breast cancer prognosis. Using data obtained from the Gene Expression Omnibus (GEO) database, correlation and Cox multivariate regression analyses were used to determine the CENP genes associated with chemotherapy responses and survival in patients with breast cancer. Weighted gene co-expression network and correlation analyses were used to determine the gene modules co-expressed with the identified genes and the differential expression of gene modules associated with the pathological complete response (PCR) and residual disease (RD) subgroups. CENPA, CENPE, CENPF, CENPI, CENPJ and CENPN were associated with a high nuclear grade and low estrogen and progesterone receptor expression levels. In addition, CENPA, CENPB, CENPC and CENPO were independent factors affecting the distant relapse-free survival (DRFS) rates in patients with breast cancer. Patients with high expression levels of CENPA or CENPO exhibited poor prognoses, whereas those with high expression levels of CENPB or CENPC presented with favorable prognoses. For validation between databases, the Cancer Genome Atlas (TCGA) database analysis also revealed that CENPA, CENPB and CENPO exerted similar effects on overall survival. However, according to the multivariate analyses, only CENPA was an independent risk factor associated with DRFS in GEO database. In addition, in the RD subgroup, patients with higher CENPA expression levels had a worse prognosis compared with those with lower CENPA expression levels. Among patients with high expression levels of CENPA, the PI3K/Akt/mTOR pathway was more likely to be activated in the RD compared with the PCR subgroup. The same trend was observed in TCGA data. These results suggested that high CENPA expression levels plus upregulation of the PI3K/Akt/mTOR signaling pathway may affect DRFS in patients with breast cancer.

Novel application of single-cell next-generation sequencing for determination of intratumoral heterogeneity of canine osteosarcoma cell lines

Osteosarcoma (OSA) is a highly aggressive and metastatic neoplasm of both the canine and human patient and is the leading form of osseous neoplasia in both species worldwide. To gain deeper insight into the heterogeneous and genetically chaotic nature of OSA, we applied single-cell transcriptome (scRNA-seq) analysis to 4 canine OSA cell lines. This novel application of scRNA-seq technology to the canine genome required uploading the CanFam3.1 reference genome into an analysis pipeline (10X Genomics Cell Ranger); this methodology has not been reported previously in the canine species, to our knowledge. The scRNA-seq outputs were validated by comparing them to cDNA expression from reverse-transcription PCR (RT-PCR) and Sanger sequencing bulk analysis of 4 canine OSA cell lines (COS31, DOUG, POS, and HMPOS) for 11 genes implicated in the pathogenesis of canine OSA. The scRNA-seq outputs revealed the significant heterogeneity of gene transcription expression patterns within the cell lines investigated (COS31 and DOUG). The scRNA-seq data showed 10 distinct clusters of similarly shared transcriptomic expression patterns in COS31; 12 clusters were identified in DOUG. In addition, cRNA-seq analysis provided data for integration into the Qiagen Ingenuity Pathway Analysis software for canonical pathway analysis. Of the 81 distinct pathways identified within the clusters, 33 had been implicated in the pathogenesis of OSA, of which 18 had not been reported previously in canine OSA.

CENPE, PRC1, TTK, and PLK4 May Play Crucial Roles in the Osteosarcoma Progression

Osteosarcoma (OS) is a cancerous tumor in a bone. We aimed to identify the critical genes involved in OS progression, and then try to elucidate the molecular mechanisms of this disease. The microarray data of GSE32395 was used for the present study. We analyzed differentially expressed genes (DEGs) in OS cells compared with control group by Student’s t-test. The significant enriched gene ontology (GO) and kyoto encyclopedia of genes and genomes (KEGG) pathways were analyzed for upregulated genes and downregulated genes, respectively. In addition, a protein-protein interaction (PPI) network was constructed. GO and KEGG enrichment analyses were conducted for genes in the PPI network. In total, 183 DEGs, including 100 upregulated DEGs and 83 downregulated DEGs were screened. The upregulated DEGs were significantly enriched in 2 KEGG pathways, such as “Glycosaminoglycan biosynthesis-chondroitin sulfate” and the downregulated DEGs were significantly enriched in 12 pathways, including “cell adhesion molecules,” “pentose phosphate pathway” and “allograft rejection.” GO enrichment analysis indicated that the upregulated DEGs were significantly involved in biological process, such as “multicellular organismal metabolic process” and “limb morphogenesis,” while the downregulated DEGs were significantly enriched in biological process, such as “Positive regulation of pathway-restricted SMAD protein phosphorylation.” The PPI network included 84 interactions and 51 nodes. The “glycosaminoglycan biosynthesis-chondroitin sulfate pathway,” “microtubule motor activityfunction,” and “regulation of mitosis process” were significantly enriched by genes in PPI network. In particular, CENPE, PRC1, TTK, and PLK4 had higher degrees in the PPI network. The interactions between TTK and PLK4 as well as CENPE and PRC1 may involve in the OS development. These 4 genes might be possible biomarkers for the treatment and diagnosis of OS.

Targeting Molecular Mechanisms Underlying Treatment Efficacy and Resistance in Osteosarcoma: A Review of Current and Future Strategies

Osteosarcoma is the most common primary malignant bone tumour in children and adolescents. Due to micrometastatic spread, radical surgery alone rarely results in cure. Introduction of combination chemotherapy in the 1970s, however, dramatically increased overall survival rates from 20% to approximately 70%. Unfortunately, large clinical trials aiming to intensify treatment in the past decades have failed to achieve higher cure rates. In this review, we revisit how the heterogenous nature of osteosarcoma as well as acquired and intrinsic resistance to chemotherapy can account for stagnation in therapy improvement. We summarise current osteosarcoma treatment strategies focusing on molecular determinants of treatment susceptibility and resistance. Understanding therapy susceptibility and resistance provides a basis for rational therapy betterment for both identifying patients that might be cured with less toxic interventions and targeting resistance mechanisms to sensitise resistant osteosarcoma to conventional therapies.

High mobility group box 2 modulates the progression of osteosarcoma and is related with poor prognosis

Background

Increased expression of high mobility group box 2 (HMGB2) has been reported to promote the progression of several malignancies and be related to poor outcome. However, few studies have explored the relationship between HMGB2 and osteosarcoma. In this study, we aimed to obtain a better understanding of HMGB2 and its function in osteosarcoma.

Methods

Utilizing osteosarcoma paraffin sections and osteosarcoma cell lines, we observed the clinico-pathological relationship of osteosarcoma with HMGB2 expression and investigated the functions of HMGB2 in vitro. The possible pathways and regulation networks in which HMGB2 is involved were further explored through analysis of miRNA, mRNA and lncRNA micro array data sets.

Results

Strong expression of HMGB2 was found to be related with Enneking staging (P=0.002), tumor size (P=0.006), metastasis (P<0.001), and survival (P=0.011) in osteosarcoma. Multivariate analysis revealed that HMGB2 might have independent prognostic value in osteosarcoma (P=0.022). Kaplan-Meier curves and the log-rank test showed that survival time was significantly reduced in OS patients with strong HMGB2 expression (P=0.0056). In vitro experiments showed that HMGB2 overexpression promoted cell proliferation and enhanced the migration and invasion ability of osteosarcoma cells. Gene Ontology (GO) term analysis of osteosarcoma cell lines revealed HMGB2 to have various functions and to be mainly enriched in regulation of cell proliferation, cell death, and DNA binding. A competing endogenous RNA (ceRNA) network of miR-139-5p and six candidate lncRNAs was also suggested as targeting HMGB2 in osteosarcoma.

Conclusions

Our findings suggest that HMGB2 might have various functions in promoting the progression of osteosarcoma and may serve as a new target for osteosarcoma research.

Small Nucleolar RNAs Determine Resistance to Doxorubicin in Human Osteosarcoma

Doxorubicin (Dox) is one of the most important first-line drugs used in osteosarcoma therapy. Multiple and not fully clarified mechanisms, however, determine resistance to Dox. With the aim of identifying new markers associated with Dox-resistance, we found a global up-regulation of small nucleolar RNAs (snoRNAs) in human Dox-resistant osteosarcoma cells. We investigated if and how snoRNAs are linked to resistance. After RT-PCR validation of snoRNAs up-regulated in osteosarcoma cells with different degrees of resistance to Dox, we overexpressed them in Dox-sensitive cells. We then evaluated Dox cytotoxicity and changes in genes relevant for osteosarcoma pathogenesis by PCR arrays. SNORD3A, SNORA13 and SNORA28 reduced Dox-cytotoxicity when over-expressed in Dox-sensitive cells. In these cells, GADD45A and MYC were up-regulated, TOP2A was down-regulated. The same profile was detected in cells with acquired resistance to Dox. GADD45A/MYC-silencing and TOP2A-over-expression counteracted the resistance to Dox induced by snoRNAs. We reported for the first time that snoRNAs induce resistance to Dox in human osteosarcoma, by modulating the expression of genes involved in DNA damaging sensing, DNA repair, ribosome biogenesis, and proliferation. Targeting snoRNAs or down-stream genes may open new treatment perspectives in chemoresistant osteosarcomas.

Exploration and validation of downregulated microRNA-199a-3p, downstream messenger RNA targets and transcriptional regulation in osteosarcoma

Osteosarcoma (OS) is a primary bone tumor with a high incidence and mortality in children and adolescents. Emerging evidence shows that microRNAs (miRNAs) participate in biological tumor mechanisms by targeting downstream messenger RNAs (mRNAs). This article aimed to investigate the potential regulatory targets of microRNA-199a-3p (miR-199a-3p) in OS and to contribute to the understanding of miR-199a-3p-related OS regulatory mechanisms. MicroRNA-related Gene Expression Omnibus (GEO) chips, ArrayExpress chips and literature data were used to determine the expression of miR-199a-3p in OS and pooled to explore its potential clinical value. To investigate the target genes of miR-199a-3p further, we integrated the results from the following three-part gene study: Twelve online prediction tools were used to predict the target genes of miR-199a-3p; the GEO GSE89370 chip transfected with miRSelect pEP-miR-199a-3p was used to analyze the downregulated differentially expressed genes (DEGs) in OS cells; and highly expressed DEGs were derived from an in-house microarray generated from three pairs of clinical OS and normal tissue samples acquired through our department. Then, we analyzed the target genes using the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases and the protein-protein interaction (PPI) network to further identify the primary target genes. In addition, we constructed transcription factor (TF)-miRNA-joint gene feed-forward regulatory loops (FFLs) with Circuits DB using miR-199a-3p as the core. A comprehensive meta-analysis of a hub of miR-199a-3p targeted genes was performed to integrate expression level, summary ROC (sROC) curves and survival analysis results from the GEO data for verification and exploration. Finally, the expression levels of the hub genes were verified in OS tissues and U2OS cells by immunohistochemistry (IHC) and immunocytochemistry (ICC). Data on miR-199a-3p expression were obtained from three data sets (GSE65071, GSE69524, and PMID 21666078), which showed low miR-199a-3p expression levels in OS tissues. The combined data indicated the same tendency, with the SMD of the random effect model, as shown in forest plots, being -2.8 (95% CI: -4.49, -1.11). In addition, we determined that miR-199a-3p may serve as a molecular marker useful for distinguishing OS tissues from normal tissues with high sensitivity and specificity, with the measured outcomes being 0.94 (95% CI: 0.80, 0.99) and 0.96 (95% CI: 0.78, 1.00), respectively. In addition, 391 genes were considered targets of miR-199a-3p in OS, and the enrichment analysis indicated that these targets were mainly enriched in proteoglycans in cancer and in spliceosomes. Four genes, CDKI, CCNB1, AURKA and NEK2, were regarded as hub targets based on the PPI data. Subsequently, TF-miRNA-joint genes FFLs were constructed in Circuits DB and included 17 TFs and 82 joint targets. These joint targets were mainly enriched in spliceosomes. UBE2D1 and RBM25 were regarded as hub joint targets based on the enrichment analysis. All selected target genes were further verified to ensure that they were upregulated in OS and to determine their prognostic significance. At the experimental verification level, the CDK1 protein was confirmed to be positively expressed in the cytoplasm of OS tissues and the U2OS cell line. Our study verified that miR-199a-3p was obviously downregulated in OS. CDK1, CCNB1, NEK2, AURKA, UBE2D1 and RBM25 were identified as potential target genes of miR-199a-3p in OS.

SKA1 induces de novo MTX-resistance in osteosarcoma through inhibiting FPGS transcription

De novo methotrexate (MTX)-resistance, whose underlying mechanism remains largely unknown, usually leads to very poor prognosis in patients with osteosarcoma (OS). In this study, we established the de novo MTX-resistant OS cell line SF-86 and identified the candidate gene spindle and kinetochore associated complex subunit 1 (SKA1) as potentially related to de novo MTX-resistance. Analysis of a cohort of 95 OS patients demonstrated that SKA1 overexpression significantly correlated with de novo MTX-resistance and poor 5-year survival. Mechanistically, SKA1 overexpression lead to a downregulation of folylpoly-γ-glutamate synthetase (FPGS), a key enzyme that converts MTX into its active form, MTX-PG. We further demonstrated that SKA1 interacts with DNA-directed RNA polymerase II subunit RPB3. ChIP analysis revealed that RPB3 binds the promoter region of the FPGS gene and triggers FPGS transcription upon MTX treatment in SW1353, a MTX-sensitive OS cell line lacking endogenous SKA1 expression. On the contrary, this process is blocked in SF-86 cells due to the formation of an inhibitory SKA1-RPB3 complex. Furthermore, downregulation of SKA1 levels restores MTX sensitivity in SF-86. Collectively, our study has established the de novo MTX-resistant cell line SF-86 and identified SKA1 as a novel regulator of FPGS, playing a key role in the development of de novo MTX-resistance in OS.

Effects of Sepantronium Bromide (YM-155) on the Whole Transcriptome of MDA-MB-231 Cells: Highlight on Impaired ATR/ATM Fanconi Anemia DNA Damage Response

Sepantronium bromide (YM-155) is believed to elicit apoptosis and mitotic arrest in tumor cells by reducing (BIRC5, survivin) mRNA. In this study, we monitored changes in survivin mRNA and protein after treating MDA-MB-231 cells with YM-155 concurrent with evaluation of whole transcriptomic (WT) mRNA and long intergenic non-coding RNA at 2 time points: 8 h sub-lethal (83 ng/mL) and 20 h at the LC50 (14.6 ng/mL). The data show a tight association between cell death and the precipitating loss of survivin protein and mRNA (-2.67 fold-change (FC), p<0.001) at 20 h, questioning if the decline in survivin is attributed to cell death or drug impact. The meager loss of survivin mRNA was overshadowed by enormous differential change to the WT in both magnitude and significance for over 2000 differentially up/down-regulated transcripts: (+22 FC to -12 FC, p<0.001). The data show YM-155 to up-regulate transcripts in control of circadian rhythm (NOCT, PER, BHLHe40, NFIL3), tumor suppression (SIK1, FOSB), histone methylation (KDM6B) and negative feedback of NF-kappa B signaling (TNFAIP3). Down-regulated transcripts by YM-155 include glucuronidase (GUSBP3), numerous micro-RNAs, DNA damage repair elements (CENPI, POLQ, RAD54B) and the most affected system was the ataxia-telangiectasia mutated (ATM)/Fanconi anemia E3 monoubiquitin ligase core complexes (FANC transcripts - A/B/E/F/G/M), FANC2, FANCI, BRCA1, BRCA2, RAD51, PALB2 gene and ATR (ATM- and Rad3-Related) pathway. In conclusion, these findings suggest that a primary target of YM-155 is the loss of replicative DNA repair systems.

Exploring the key genes and pathways of side population cells in human osteosarcoma using gene expression array analysis

BACKGROUND

Human osteosarcoma (OS) is one of the most common primary bone sarcoma, because of early metastasis and few treatment strategies. It has been reported that the tumorigenicity and self-renewal capacity of side population (SP) cells play roles in human OS via regulating of target genes. This study aims to complement the differentially expressed genes (DEGs) that regulated between the SP cells and the non-SP cells from primary human OS and identify their functions and molecular pathways associated with OS.

METHODS

The gene expression profile GSE63390 was downloaded, and bioinformatics analysis was made.

RESULTS

One hundred forty-one DEGs totally were identified. Among them, 72 DEGs (51.06%) were overexpressed, and the remaining 69 DEGs (48.94%) were underexpressed. Gene ontology (GO) and pathway enrichment analysis of target genes were performed. We furthermore identified some relevant core genes using gene-gene interaction network analysis such as EIF4E, FAU, HSPD1, IL-6, and KISS1, which may have a relationship with the development process of OS. We also discovered that EIF4E/mTOR signaling pathway could be a potential research target for therapy and tumorigenesis of OS.

CONCLUSION

This analysis provides a comprehensive understanding of the roles of DEGs coming from SP cells in the development of OS. However, these predictions need further experimental validation in future studies.

Plasma metabolic profiling on postoperative colorectal cancer patients with different traditional Chinese medicine syndromes

OBJECTIVES

This study aims to investigate the metabolic profiles of postoperative colorectal cancer (PCRC) patients with different traditional Chinese medicine (TCM) syndromes and to discuss the metabolic mechanism under PCRC progression and TCM syndrome classification.

METHODS

Fifty healthy controls (HC) and 70 PCRC patients, including 10 Dampness and heat syndrome (DHS), 33 Spleen deficiency syndrome (SDS), 19 Liver and kidney Yin deficiency syndrome (LKYDS) and 8 with non-TCM syndrome (NS) were enrolled. Plasma metabolic profiles were detected by Gas chromatography-mass spectrometry (GC-MS) and analyzed by principal component analysis (PCA) and partial least squares-discriminate analysis (PLS-DA). Furthermore, pathway enrichment was analyzed based on KEGG and DAVID databases and metabolic network was constructed via metaboanalyst and cytoscape.

RESULTS

The top-3 metabolites with higher abundance in PCRC compared with HC were terephthalic acid (165.417-fold), ornithine (24.484-fold) and aminomalonic acid (21.346-fold). And the cholesterol (0.588-fold) level was decreased in PCRC. l-Alanine, 1, 2-ethanediamine, urea, glycerol, glycine, aminomalonic acid, creatinine and palmitic acid were specifically altered in the DHS, while d-tryptophan was exclusively changed in SDS, and l-proline, 1, 2, 3-propanetricarboxylic acid, d-galactose and 2-indolecarboxylic acids in LKYDS.

CONCLUSIONS

The plasma metabolic profiles were perturbed in PCRC patients. Increased levels of terephthalic acid might indicate high risk of relapse and elevated ornithine may contribute to the post-operational recovery or may raise the susceptibility to PCRC recurrence. The metabolic profiles of DHS, SDS, LKYDS and NS were almost separately clustered, indicating the possibility of explaining TCM syndromes classification using metabolomics. Furthermore, creatinine and aminomalonic acid alternation might correlate with the formation of DHS, while d-tryptophan may associate with SDS and d-galactose and 1, 2, 3-propanetricarboxylic acid may relate to LKYDS. As numbers of patients in each TCM syndrome are small, further study is needed to verify those results.

Identification of differential co-expressed gene networks in early rheumatoid arthritis achieving sustained drug-free remission after treatment with a tocilizumab-based or methotrexate-based strategy

BACKGROUND

Methotrexate is endorsed to be used as first-line treatment in rheumatoid arthritis (RA). However, a large proportion of patients need additional treatment with a biological disease-modifying anti-rheumatic drug (DMARD) to adequately suppress their disease activity. A better understanding of genotypes could help to distinguish between patients with different pathogenic mechanisms. The aim of this study was therefore to identify networks of genes within DMARD-naive early RA patients associated with achieving sustained drug-free remission (sDFR) after initiating tocilizumab plus methotrexate, tocilizumab, or methotrexate therapy.

METHODS

Samples were used from 60 patients from the U-Act-Early study who received tocilizumab plus methotrexate, tocilizumab, or methotrexate therapy, and who achieved sDFR (≥3 months in drug-free remission until the end of the study, n = 37) after therapy was tapered and subsequently stopped, or who were not able to discontinue the therapy as controls (n = 23). Whole blood samples were collected and ribonucleic acid (RNA) was isolated from positive cluster of differentiation 4 (CD4⁺) and CD14⁺ cells and analysed using high-throughput sequencing. Weighted gene co-expression network analyses were performed to identify clusters (i.e. modules) of differently expressed genes associated with achieving sDFR and which were subsequently used for pathway analyses.

RESULTS

Network analyses within CD4⁺ cells identified two significant modules in the tocilizumab plus methotrexate arm and four modules in the tocilizumab and methotrexate arms, respectively (p ≤ 0.039). Important pathways in the module best correlating with achieving sDFR were in the tocilizumab plus methotrexate arm related to processes involved with transcription and translation; in the tocilizumab arm, pathways were related to migration of white blood cells and G-protein coupled receptors, and in the methotrexate arm pathways were involved with the response to a bacterial or biotic (i.e. biological material)-related stimulus. No relevant networks could be identified in the sequenced CD14⁺ cells.

CONCLUSIONS

Within networks of co-expressed genes, several pathways were found related to achieving sDFR after initiating therapy with tocilizumab, methotrexate, or the combination. Between the three strategy arms, we identified different networks of predisposing genes which indicates that specific gene expression profiles, depending on the treatment strategy chosen, are associated with a higher chance of achieving sDFR.

TRIAL REGISTRATION

Clinicaltrials.gov, NCT01034137 . Registered on 16 December 2009.