High expression of VARS promotes the growth of multiple myeloma cells by causing imbalance in valine metabolism

OBJECTIVES

Multiple myeloma (MM) is a plasma cell malignancy occurring in middle and old age. MM is still an incurable disease due to its frequent recurrence and drug resistance. However, its pathogenesis is still unclear. Abnormal amino acid metabolism is one of the important characteristics of MM, and the important metabolic pathway of amino acids participates in protein synthesis as basic raw materials. Aminoacyl transfer ribonucleic acid synthetase (ARS) gene is a key regulatory gene in protein synthesis. This study aims to explore the molecular mechanism for ARS, a key factor of amino acid metabolism, in regulating amino acid metabolism in MM and affecting MM growth.

METHODS

The corresponding gene number was combined with the gene expression profile GSE5900 dataset and GSE2658 dataset in Gene Expression Omnibus (GEO) database to standardize the gene expression data of ARS. GSEA_4.2.0 software was used to analyze the difference of gene enrichment between healthy donors (HD) and MM patients in GEO database. GraphPad Prism 7 was used to draw heat maps and perform data analysis. Kaplan-Meier and Cox regression model were used to analyze the expression of ARS gene and the prognosis of MM patients, respectively. Bone marrow samples from 7 newly diagnosed MM patients were collected, CD138+ and CD138- cells were obtained by using CD138 antibody magnetic beads, and the expression of ARS in MM clinical samples was analyzed by real-time RT-PCR. Human B lymphocyte GM12878 cells and human MM cell lines ARP1, NCI-H929, OCI-MY5, U266, RPMI 8266, OPM-2, JJN-3, KMS11, MM1.s cells were selected as the study objects. The expression of ARS in MM cell lines was analyzed by real-time RT-PCR and Western blotting. Short hairpin RNA (shRNA) lentiviruses were used to construct gene knock-out plasmids (VARS-sh group). No-load plasmids (scramble group) and gene knock-out plasmids (VARS-sh group) were transfected into HEK 293T cells with for virus packaging, respectively. Stable expression cell lines were established by infecting ARP1 and OCI-MY5 cells, and the effects of knockout valyl-tRNA synthetase (VARS) gene on proliferation and apoptosis of MM cells were detected by cell counting and flow cytometry, respectively. GEO data were divided into a high expression group and a low expression group according to the expression of VARS. Bioinformatics analysis was performed to explore the downstream pathways affected by VARS. Gas chromatography time-of-flight mass spectrometry (GC-TOF/MS) and high performance liquid chromatography (HPLC) were used to detect the valine content in CD138+ cells and ARP1, OCI-MY5 cells and supernatant of knockdown VARS gene in bone marrow samples from patients, respectively.

RESULTS

Gene enrichment analysis showed that tRNA processing related genes were significantly enriched in MM compared with HD (P<0.0001). Further screening of tRNA processing-pathway related subsets revealed that cytoplasmic aminoacyl tRNA synthetase family genes were significantly enriched in MM (P<0.0001). The results of gene expression heat map showed that the ARS family genes except alanyl-tRNA synthetase (AARS), arginyl-tRNA synthetase (RARS), seryl-tRNA synthetase (SARS) in GEO data were highly expressed in MM (all P<0.01). With the development of monoclonal gammopathy of undetermined significance (MGUS) to MM, the gene expression level was increased gradually. Kaplan-Meier univariate analysis of survival results showed that there were significant differences in the prognosis of MM patients in methionyl-tRNA synthetase (MARS), asparaginyl-tRNA synthetase (NARS) and VARS between the high expression group and the low expression group (all P<0.05). Cox regression model multivariate analysis showed that the high expression of VARS was associated with abnormal overall survival time of MM (HR=1.83, 95% CI 1.10 to 3.06, P=0.021). The high expression of NARS (HR=0.90, 95% CI 0.34 to 2.38) and MARS (HR=1.59, 95% CI 0.73 to 3.50) had no effect on the overall survival time of MM patients (both P>0.05). Real-time RT-PCR and Western blotting showed that VARS, MARS and NARS were highly expressed in CD138+ MM cells and MM cell lines of clinical patients (all P<0.05). Cell counting and flow cytometry results showed that the proliferation of MM cells by knockout VARS was significantly inhibited (P<0.01), the proportion of apoptosis was significantly increased (P<0.05). Bioinformatics analysis showed that in addition to several pathways including the cell cycle regulated by VARS, the valine, leucine and isoleucine catabolic pathways were upregulated. Non-targeted metabolomics data showed reduced valine content in CD138+ tumor cells in MM patients compared to HD (P<0.05). HPLC results showed that compared with the scramble group, the intracellular and medium supernatant content of ARP1 cells and the medium supernatant of OCI-MY5 in the VARS-shRNA group was increased (all P<0.05).

CONCLUSIONS

MM patients with abnormal high expression of VARS have a poor prognosis. VARS promotes the malignant growth of MM cells by affecting the regulation of valine metabolism.

Thymidine phosphorylase facilitates retinoic acid inducible gene-I induced endothelial dysfunction

Activation of nucleic acid sensors in endothelial cells (ECs) has been shown to drive inflammation across pathologies including cancer, atherosclerosis and obesity. We previously showed that enhancing cytosolic DNA sensing by inhibiting three prime exonuclease 1 (TREX1) in ECs led to EC dysfunction and impaired angiogenesis. Here we show that activation of a cytosolic RNA sensor, Retinoic acid Induced Gene 1 (RIG-I) diminishes EC survival, angiogenesis and triggers tissue specific gene expression programs. We discovered a RIG-I dependent 7 gene signature that affects angiogenesis, inflammation and coagulation. Among these, we identified the thymidine phosphorylase TYMP as a key mediator of RIG-I induced EC dysfunction via its regulation of a subset of interferon stimulated genes. Our RIG-I induced gene signature was also conserved in the context of human diseases - in lung cancer vasculature and herpesvirus infection of lung endothelial cells. Pharmacological or genetic inhibition of TYMP rescues RIG-I induced EC death, migration arrest and restores sprouting angiogenesis. Interestingly, using RNAseq we identified a gene expression program that was RIG-I induced but TYMP dependent. Analysis of this dataset indicated that IRF1 and IRF8 dependent transcription is diminished in RIG-I activated cells when TYMP is inhibited. Functional RNAi screen of our TYMP dependent EC genes, we found that a group of 5 genes - Flot1, Ccl5, Vars2, Samd9l and Ube2l6 are critical for endothelial cell death mediated by RIG-I activation. Our observations identify mechanisms by which RIG-I drives EC dysfunction and define pathways that can be pharmacologically targeted to ameliorate RIG-I induced vascular inflammation.

[Amino-tRNA-synthetase activity of rabbit muscles in fasting and aging]

The aminoacyl-tRNA-synthetase activity of muscular soluble fraction was studied in normal rabbits, in those fasting for a long period and old rabbits. A rigion of linear dependence is found for the initial rate of the tRNA aminoacylation reaction on the concentration of the enzymic preparation and also on the substrate amount. Having valuated the initial rate of the aminoacylation reaction, the asparagil- and valyl-tRNA-synthetase activities of muscles were determined in normal rabbits, in those fasting for a long period and in old rabbits. The data obtained evidence for an increase of the studied activities in the muscles of the fasting animals and for their decrease in old animals.

Ultracentrifugation studies of yeast valyl-tRNA synthetase and of its interaction with tRNAVal

Yeast valyl-tRNA synthetase and its complexes with yeast tRNAVal were investigated by means of analytical ultracentrifugation. A molecular weight of 125 700 +/- 1500 and a sedimentation coefficient (SO 20, w) of 6.3 +/- 0.3 were found for the native enzyme. When the enzyme (3--60 muM) was mixed with its cognate tRNA, several types of complex were observed, depending on the relative amounts of the two macromolecules. In the presence of equimolecular amounts of tRNA and enzyme, a complex formed by the association of one of each molecule was observed with a sedimentation coefficient of about 7.3 S. However, for tRNA/enzyme stoichiometries lower than one, beside the 1 : 1 complex, a complex of higher molecular weight was observed, with a sedimentation coefficient of about 10.0 S which fits with the association of two valyl-tRNA synthetase molecules with one tRNA molecule. This 2 : 1 complex was predominant from tRNA/enzyme stoichiometries lower than 0.3. It dissociated into the 1 : 1 complex upon addition of monovalent salts or MgCl2, suggesting the electrostatic nature of the interaction in this association. All these association and dissociation phenomena were detected over a large range of pH (6.0--7.5) and in various buffers.

[Differences in the thermostability of aminoacyl-mRNA-synthetases from intact and denervated rabbit muscles]

The stability of alanyl-tRNA synthetase and valyl-tRNA synthetase from the intact and denervated for 30 days rabbit skeletal muscles towards the inactivating effect of heat (42 degrees for 10 min) has been studied. The activity of the enzymes was measured in the supernatant fraction obtained by centrifugation of muscle homogenates at 105, 000 g for 1 h. The stability of both alanyl- and valyl-tRNA synthetases from the denervated m. soleus and m. gastrocnemius was shown to be considerably decreased as compared with the enzymes from the intact muscles.

A neutron investigation of yeast valyl-tRNA synthetase interaction with tRNAs

A new way of studying RNA-protein complexes, using neutron small angle scattering in solution, is described and was applied in the case of the system, yeast valyl-tRNA synthetase, interacting with its cognate and non cognate yeast tRNAs. It was shown that, when limited amounts of tRNA (either cognate or non cognate) are added to valyl-tRNA synthetase, a complex consisting of two enzyme molecules and one tRNA molecule is first formed. It is subsequently dissociated to a one to one complex when more tRNA is present in the solution. The association curve shows a maximum for a molecular ratio, enzyme over tRNA, equal to 2.

The infrastructure of valyl-transfer Ribonucleic Acid synthetase from yeast

Valyl-tRNA synthetase from yeast is shown to contain a significant proportion of repeated sequence, similar to the enzyme from bacterial sources.